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'); 

將產生：