科技矢量图相关问题研究程序及图像(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;得到图像原图像还原拟合后数据图像还原拟合前后图像对比多项式拟合参数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+