Question

matlab二值图像 拟合直线

各位大侠，我这里有张二值图，希望将每行黑色像素的中点取出来，然后进行拟合，希望大神指点程序

Answer 1

clear all;close all;
I0=imread('line.jpg');
I=rgb2gray(I0);
I=im2bw(I);
subplot(131);imshow(I);
[m,n]=size(I);
Imid=zeros(m,n);
x = zeros(1,m);
y = zeros(1,n);
for i=1:m  %扫描行
    I1=I(i,:);  %提取一行
    dI1 = diff(I1);  %微分,找到边界点
    edge = find(dI1~=0);
    if length(edge)==2  
        y(i) = fix((edge(1)+edge(2))/2);
        x(i) = i;
        Imid(x(i),y(i)) = 1;
    end
end
%显示提取的中点
subplot(132);imshow(Imid)

%直线拟合
%去掉异常点
x(find(y==0))=[];
y(find(y==0))=[];
p = polyfit(x,y,1);

x1 = 1:m;
y1 = polyval(p,x1);
%画出直线
L=I0;
for i=1:m
     L(x1(i),round(y1(i)),1)=255;
end
subplot(133);imshow(L);title(['y=' num2str(p(1)) '+' num2str(p(2)) 'x']);

追问

十分感谢！x(find(y==0))=[];这句报错请问是什么原因？

追答

报什么错?

追问

Index of element to remove exceeds matrix dimensions.

Error in Untitled (line 51)
x(find(y==0))=[];

追答

去掉这两行看看
x(find(y==0))=[];
y(find(y==0))=[];

追问

还是有问题，给你私信了，麻烦看一下~这边追问太多了，我的财富值也没几个，剩下都给你了~