nihecaishangxian.m

用于直线拟合

clear all
clc
e=imread('129.jpg');
a=rgb2gray(e);
bw=edge(a,'canny',0.08); 
[x1,x2]=size(bw);
x=1;
for i=1:x1
    for j=1:x2
          if bw(i,j)==0
              continue
          else
             d(x)=i;
             c(x)=j;
             x=x+1;
           
          end
    end
end
 tu1=plot(d,c,'g.','markersize',2)
 grid on
 [xx,yy]=ginput(2);
 gtext('\fontsize{15}\fontname{隶书}关键点')
 gtext('\fontsize{15}\fontname{隶书}关键点')
 figure
 imshow (tu1);
 x1=1;
 xx1=round(xx(1))
 xx2=round(xx(2))
 yy1=round(yy(1))
 yy2=round(yy(2))

%for i1=xx1:xx2
    if yy1>yy2
           %t=yy2;yy2=yy1;yy1=t;
          %end
        for i1=xx1:1:xx2 
          for j1=yy1:-1:yy2
                 if bw(i1,j1)==0
                 continue
                 else
                 d1(x1)=i1;
                 c1(x1)=j1;
                 x1=x1+1;
                 break
                 end
          end
       end
    
    
    else
        for i1=xx2:-1:xx1
           for j1=yy2:-1:yy1
                 if bw(i1,j1)==0
                 continue
                 else
                 d1(x1)=i1;
                 c1(x1)=j1;
                 x1=x1+1;
                 break
                 end
           end
        end 
    end
huigui(d1,c1)