Matlab - 从边缘创建x，y曲线

Question

我想用Matlab解释我的问题。我的目标是在从照片坦克获得自由液面而要做到这一点，我已经使用这个算法：

A = 'C:\foto\img3.bmp'; 
B = imread(A, 'bmp'); 
figure(1), imshow(B); 
C = rgb2gray(B); 
level = graythresh(C); 
bw = im2bw(C,level); 
bw = bwareaopen(bw, 50); 
figure, imshow(bw); 
BW1 = edge(bw,'canny'); 
figure(2), imshow(BW1); 
imwrite(BW1, 'C:\foto\im1_edge.bmp', 'bmp') 

现在我有表面边缘为白色，与黑色的背景。 接下来我也只能检测白色像素的位置：

I= imread('C:\foto\img3_edge.bmp'); 
imshow(I); 
[r c] =size(I); 
for j=1:c 
    for i=1:r 
    if(I(i,j)==1)     
    [i j] 
    end 
    end 
end 

在这一点上，我怎么可以报告（含宏，自动可能）每对夫妇坐标的笛卡尔（x，y）平面？我的目标是从边缘获得重建的类型“y = f（x）”的函数。 我和另一个边缘尝试，并用油漆修改，删除所有无用像素的图像，例子是这样的一个：

http://i59.tinypic.com/20jh1g2.jpg

与代码：

I = im2bw(I); 

它返回我的错误“警告：输入图像已经是二进制。”接下来使用的代码：

[r c] = find(I), output = [r c]; 
plot(r,c,'.') 

我得到这一个：

http://i59.tinypic.com/2iqmqvl.jpg

此外，当我尝试在cftool插入R作为XDATA和c如YDATA，我得到了同样的问题，并且当我使用“Interpolant”装配，它会返回一个错误.....为什么？

请帮助我，谢谢大家！ :-) PS：对不起，太长的帖子:)

Answer 1

我不知道你想适用于你的数据是什么类型的适合。以下代码将适合数据上的三阶多项式。

I = imread('20jh1g2.jpg'); 
I = im2bw(I); 
imshow(I); 

[r, c] = find(I); 
figure; 
plot(c,r,'.'); 
hold on; 
f = fit(c, r, 'poly3'); 
plot((min(c):max(c)),f(min(c):max(c)), 'red', 'LineWidth', 3); 

将产生：

旋转可以由轴是如何定义进行说明。在您的图像中，Y轴从顶部的0开始到底部的374。您可以使用以下代码将您的拟合结果转换回二进制映像;

x = (min(c):max(c))'; 
y = round(f(x)); 
I = zeros(size(I)); 
I(y +((x-1)*size(I,1))) = 1; 
figure 
imshow(I); 

将产生：

拟合的结果，f存储在cfit对象。您可以通过为其提供x的值来评估此功能，如上所示。函数系数可以通过在命令窗口中打印cfit对象的字段来找到;

f = 

    Linear model Poly3: 
    f(x) = p1*x^3 + p2*x^2 + p3*x + p4 
    Coefficients (with 95% confidence bounds): 
     p1 = -6.252e-06 (-6.542e-06, -5.963e-06) 
     p2 = 0.001753 (0.001588, 0.001918) 
     p3 =  -0.3667 (-0.3932, -0.3401) 
     p4 =  290.4 (289.3, 291.6) 

要翻转参考框架内的函数并计算质心，您可以使用它;

I = imread('20jh1g2.jpg'); 
I = im2bw(I); 

[r, c] = find(I); 
r = -r + size(I,1); 
f = polyfit(c, r, 3); 
plot((min(c):max(c)),polyval(f,(min(c):max(c))), 'red', 'LineWidth', 3); 
hold on; 

xf = [f 0]; 
fx2 = sym2poly(poly2sym(f)^2); 
centroid = 1/polyval(polyint(f),size(I,2)) * [polyval(polyint(xf),size(I,2)) 1/2 * polyval(polyint(fx2),size(I,2))]; 
plot(centroid(1),centroid(2),'X'); 

将产生：