[ÇóÖú] ÔËÐÐʱºò³öÏÖÁËÒÔÏÂÇé¿ö£¬ÎÒÏëµÃµ½Ö±ÏßµÄбÂʺͽؾà´óÉñ°ïæÐÞ¸ÄÏÂ

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

2Â¥: Originally posted by ÔÂÖ»À¶ at 2017-04-13 13:44:27
¶Ô׎Øͼ£¬ÈâÑÛʶ±ð¡¢ÊÖ¶¯ÇÃÆäÖеĴúÂ룬ÊÇÒ»¼þÈÃÈ˺ܲ»Êæ·þµÄÊÂÇé¡£»»Î»Ë¼¿¼Ò»Ï£¬Äã»á·¢ÏÖÕâÊÇÒ»ÖֺܵÍЧµÄÇóÖú·½Ê½¡£
Äã°Ñ´úÂëÌù³öÀ´²»ºÃÂð£¿