科技矢量圖相關(guān)問題研究程序及圖像(1)

1、主程序：clear;yuchuli;hangxiantiqu;hangshujuzhengli;liexiantiqu;lieshujuzhengli;shanhang;shanlie;quzhi;nihe;子程序yuchulia=imread('C:UsersliangDesktopmatlab2010b623.jpg');figure(1);imshow(a);b=rgb2gray(a);figure(2);imshow(b);m n=size(b);ma=max(max(b);mi=min(min(b);limen=(ma+mi)/2;c=(b>(limen-50)

2、;figure(3);imshow(c);得到：原圖像灰度圖像二值化圖像子程序hangxiantiqui=1;k=1;j=620;for i=1:m if c(i,j)=0 u1(k)=i; k=k+1; endend子程序hangshujuzhenglij=1;x=k-2;v1(1)=u1(1);for i=1:x if u1(i+1)-v1(j)>25 j=j+1; v1(j)=u1(i+1); endendsum1=0;for i=1:j-2 sum1=sum1+v1(i+1)-v1(i); endbiy=sum1/(j-2);biliy=round(biy)求得縱坐標比例bili

3、y = 33子程序liexiantiquj=1;k=1;i=307;for j=1:n if c(i,j)=0 u2(k)=j; k=k+1; endend子程序lieshujuzhenglij=1;x=k-2;v2(1)=u2(1);for i=1:x if u2(i+1)-v2(j)>30 j=j+1; v2(j)=u2(i+1); endendsum2=0;for i=1:j-2 sum2=sum2+v2(i+1)-v2(i); endbix=sum2/(j-2);bilix=round(bix)得到橫坐標比例bilix =50子程序shanhangfor i=1:17 for x

4、=(v1(i)-3):(v1(i)+2); for j=1:n; c(x,j)=1; end endendimshow(c);得到去橫線圖像子程序shanliefor j=1:15 for x=(v2(j)-2):(v2(j)+3); for i=1:m; c(i,x)=1; end endendimshow(c);得到圖像去橫線豎線圖像子程序quzhid=c;x,y=find(d);i=1;t=size(x);for k1=2:t if x(k1-1)-x(k1)>=0; x1(i)=x(k1-1); y1(k1-1)=y(k1-1); i=i+1; end endy1(find(y1

5、=0)=;x2=m-x1;子程序nihey1=y1/bilix;x2=5*x2/biliy;p=polyfit(y1,x2,10);p,s=polyfit(y1,x2,10);figure(1);plot(y1,x2);grid on;x3=0:0.1:n/bilix;y3=polyval(p,x3);figure(2);plot(x3,y3);grid on;figure(3);plot(x3,y3,y1,x2);grid on;得到圖像原圖像還原擬合后數(shù)據(jù)圖像還原擬合前后圖像對比多項式擬合參數(shù)p = Columns 1 through 10 0.0000 -0.0000 0.0006 -0.0100 0.0895 -0.4594 1.1684 -0.3772 -6.8635 34.3664 Column 11-2.8896得到公式：y=p1xn+p2xn-1+pnx+