数学建模b题一等奖优秀论文

d=d+1; j=j-1; end end

A(i,2)=b; A(i,3)=c+d; A(i,4)=i+791; end

xlswrite('D:\\question1.xls',A,'A210:D418'); 附录六：

%%%%%%%%%%%%%注意！！！运行时将附件解压至D盘%%%%%%%%%%%%%% %%%%%%%%%%%%%该程序用来对图片的二值矩阵进行分类 A=zeros(209,5); for i=1:209

str='D:\\附件\\附件3\\';

I=imread([str,num2str(i),'.jpg']);%依次读取每一幅图像 i1=rgb2gray(I); %i1灰度图像

i2=im2bw(i1); % i2是二值图像，不需要 a=0;b=0;c=0;d=0; j=1;

if any(i2(1,:)==0)

while (any(i2(j,:)==0)) a=a+1; c=c+1; j=j+1; end else

while (all(i2(j,:)==1)) a=a-1; c=c+1; j=j+1; end end

A(i,1)=a;

if any(i2(j,:)==0)

while (any(i2(j,:)==0)) b=b+1; d=d+1; j=j+1; end else

while (all(i2(j,:)==1)) b=b-1; d=d+1;

j=j+1; end end

A(i,2)=b; A(i,3)=c+d; A(i,5)=i; if A(i,3)<30 A(i,4)=0;

elseif A(i,3)>30&&A(i,3)<=34 A(i,4)=1;

elseif A(i,3)>=40&&A(i,3)<=43 A(i,4)=2;

elseif A(i,3)>=44&&A(i,3)<=48 if A(i,1)>0 A(i,4)=3; else

A(i,4)=4; end

elseif A(i,3)>=50&&A(i,3)<=53 if A(i,1)>0 A(i,4)=5; else

A(i,4)=6; end

elseif A(i,3)>55&&A(i,3)<=60 if A(i,1)>0 A(i,4)=7; else

A(i,4)=8; end

elseif A(i,3)>61&&A(i,3)<=66 if A(i,1)>0 A(i,4)=9; else

A(i,4)=10; end

elseif A(i,3)>67&&A(i,3)<=70 A(i,4)=11;

elseif A(i,3)>76&&A(i,3)<=78 A(i,4)=12; elseif A(i,3)>80 A(i,4)=13; end end

xlswrite('D:\\photo3.xls',A); %将矩阵导入excel表里 附录七：

%以下程序的运行，请注意文件存放的位置！！！ %此程序运行需要时间较长，请耐心等待！！！ %请使用经格式转化器转化后的JPG格式图片 A=zeros(209,209);

a=A(:,1);%共209个纸条 l=A(:,1);%共209个纸条 for j=1:209

str='D:\\附件\\附件3\\';

I=imread([str,num2str(j),'.jpg']);%依次读取每一幅图像 i1=rgb2gray(I); %i1灰度图像

i2=im2bw(i1); % i2是二值图像，一个180*72的矩阵 a=i2(:,1); %取纸片左边缘 for i=1:209

I=imread([str,num2str(i),'.jpg']); %依次读取每一幅图像 i1=rgb2gray(I); %i1灰度图像 i2=im2bw(i1); % i2是二值图像 mi=i2(:,72); %取纸片左边缘 ni=a-mi;

A(j,i)=sqrt(dot(ni,ni)); %用以判断差距 end end

xlswrite('D:\\photo31.xls',A); %将矩阵元素导入excel表格 for i=1:209 for j=1:209

if(A(i,j)==0) l(i)=l(i)+1; end end end

for i=1:209

if(l(i)>=11) i end end

附录八：

%以下程序的运行，请注意文件存放的位置！！！ %请使用经格式转化器转化后的JPG格式图片。 %%%%%%选取以008作为开始第9组的行排序%%%%%%% A=zeros(19,19);

a=[1,8,33,46,54,57,69,71,94,127,138,139,154,159,167,175,176,197,209;]%共十九个纸条 for j=1:19

str='D:\\附件\\附件3\\'; s=a(j)

I=imread([str,num2str(a(j)),'.jpg']);%依次读取每一幅图像 i1=rgb2gray(I); %i1灰度图像 i2=im2bw(i1); % i2是二值图像 t=i2(:,72); %取纸片右边缘 for i=1:19

I=imread([str,num2str(a(i)),'.jpg']); %依次读取每一幅图像 i1=rgb2gray(I); %i1灰度图像 i2=im2bw(i1); % i2是二值图像 mi=i2(:,1); %取纸片左边缘 ni=t-mi;

A(j,i)=sqrt(dot(ni,ni)); end end

%%%%%%%%判断相邻图片并自动连接%%%%%%%%%%

r=[0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0];%记录和纸片右相邻的纸片 g=[1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1];%记录和纸片右相邻的纸片 for i=1:19 k=A(i,1); r(i)=1; for j=1:19

if(k>A(i,j)&&j~=2) k=A(i,j); r(i)=a(j); g(i)=j

end end end first=2;

str='D:\\附件\\附件3\\';

z= imread([str,num2str(a(first)),'.jpg']); for m=1:18 x=z;

y=imread([str,num2str(r(first)),'.jpg']);%记录下一张纸条 z = [x,y]; first=g(first); end

imshow(z) %请打开图片显示窗口查看拼接结果 附件九：

图15 附件三的拼接图片

附录十：

图16 附件4的拼接图片

附录十一：

%以下程序的运行，请注意文件存放的位置！！！ %请使用经格式转化器转化后的JPG格式图片。 A=zeros(19,19);

a=[791 1121 1260 1410 1560 1510 1841 1751 1110 670 1090 191 300 1901 821 1651 210 480 1371 900 111 370 771 1790 450 261 1930 1251 230 1211 1450 800 150 600 611 1480 1530 60 1871 1540 851 431 310 390 1220 990 951 621 1381 460 1390 571 1321 1881 871 2011 1441 2001 121 1620 1701 1951 1741 2071 1570 350 1821 1991 880 1331 940 731 1760 980 401 840 891 1080 1501 1810 381 1920 661 1161 1671 21 1521 1711 420 711 1401 30 1631 2041 910 1150 1851 1801 1171 2080 590 1590 1980 1551 291 130 181 1031 651 2090 1430 580 250 140 1470 1721 320 2020 510 1911 931 201 171 1781 541 2030 221 1310 1640 1941 741 1600 360

1661 1960 1290 1580 1690 470 680 641 761 1680 1181 91 691 1890 1280 410 1831 1230 1730 41 81 861 1491 780 50 700 330 751 1271 1770 1860 11 811 280 1361 1420 2051 1060 241 1340 490 521 960 1611 1200 341 721 530 630 1301 1191 1020 161 2060 831 1460 101 1000 440 971 1100 1240 70 1050 1350 1140 271 501 921 1071 1011 561 1040 1130 1971 551 ];%共十九个纸条

str='D:\\附件\\附件5\\';

I=imread([str,num2str(a(1,1)),'.jpg']);%依次读取每一幅图像 for j=2:19

str='C:\\Users\\Administrator\\Desktop\\附件5\\';

J=imread([str,num2str(a(1,j)),'.jpg']);%依次读取每一幅图像 I=[I,J]; end for i=2:11

str='D:\\附件\\附件5\\';

K=imread([str,num2str(a(i,1)),'.jpg']);%依次读取每一幅图像 for j=2:19

str='C:\\Users\\Administrator\\Desktop\\附件5\\';

J=imread([str,num2str(a(i,j)),'.jpg']);%依次读取每一幅图像 K=[K,J]; end

I=[I;K]; end

imshow(I) %请打开图片显示窗口查看拼接结果