[求助] 运行时候出现了以下情况，我想得到直线的斜率和截距大神帮忙修改下

2楼: Originally posted by 月只蓝 at 2017-04-13 13:44:27
对着截图，肉眼识别、手动敲其中的代码，是一件让人很不舒服的事情。换位思考一下，你会发现这是一种很低效的求助方式。
你把代码贴出来不好吗？

2楼: Originally posted by 月只蓝 at 2017-04-13 13:44:27
对着截图，肉眼识别、手动敲其中的代码，是一件让人很不舒服的事情。换位思考一下，你会发现这是一种很低效的求助方式。
你把代码贴出来不好吗？

clear ;
close all;
I0=imread('mid.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']);
k=num2str(p(2));
b=num2str(p(1));