2013B题碎纸片拼接

Photo1_1:

clear;clc;

path='E:\B\附件1\';

files=dir('E:\B\附件1\*.bmp');

% objdir='E:\B\附件1\';

% bgfile=[objdir,['00' int2str(i)],'.bmp'];

% lengthf=length(files);

% for i=1:lengthf

%     img=imread(strcat('E:\B\附件1\',files(i).name));

%     imshow(img);

% end

% bg=imread([path,files(1).name]);

% imshow(bg);

% imageinfo([path,files(1).name])

xsmat=ones(19,19);

for i=1:19

    imi=imread([path,files(i).name]);

    [mrow,mcol]=size(imi);

    for j=1:19

        imj=imread([path,files(j).name]);

        timi=im2double(imi);

        timj=im2double(imj);

        xssum=sum(abs(timi(:,mcol)-timj(:,1)));%计算图像i和图像j的相似度。

        xsmat(i,j)=xssum;

    end

end

first=0;last=0;

for i=1:19

    img=imread([path,files(i).name]);

   [ro,co]=size(img);

   tot=0;

   for j=1:ro

       if(img(j,1)==255)

           tot=tot+1;

       end

   end

   if tot==ro

       first=i;

       break;

   end

end

que=[];

que=[que first];

flag=1;

pos=1;

tf=first;

while(flag==1)

    t=255*255*254;

    tp=0;

    for i=1:19

        if i~=tf

            if t>xsmat(tf,i)

                t=xsmat(tf,i);

                tp=i;

            end

        end

    end

    tf=tp;

    pos=pos+1;

    que=[que tp];

    if pos==19

        flag=0;

    end

end

tim=[];

for i=1:19

    im=imread([path files(que(i)).name]);

    tim=[tim im];

end

imshow(tim);

　　Photo1_2:

clear;clc;

path='E:\B\附件2\';

files=dir('E:\B\附件2\*.bmp');

xsmat=ones(19,19);

for i=1:19

    imi=imread([path,files(i).name]);

    [mrow,mcol]=size(imi);

    for j=1:19

        imj=imread([path,files(j).name]);

        timi=im2double(imi);

        timj=im2double(imj);

        xssum=sum((timi(:,mcol)-timj(:,1)).^2);%计算图像i和图像j的相似度。

        xsmat(i,j)=xssum;

    end

end

first=0;last=0;

for i=1:19

    img=imread([path,files(i).name]);

   [ro,co]=size(img);

   tot=0;

   for j=1:ro

       if(img(j,1)==255)

           tot=tot+1;

       end

   end

   if tot==ro

       first=i;

       break;

   end

end

que=[]

que=[que first];

flag=1;

pos=1;

tf=first;

while(flag==1)

    t=255*255*254;

    tp=0;

    for i=1:19

        if i~=tf

            if t>xsmat(tf,i)

                t=xsmat(tf,i);

                tp=i;

            end

        end

    end

    tf=tp;

    pos=pos+1;

    que=[que tp];

    if pos==19

        flag=0;

    end

end

tim=[];

for i=1:19

    im=imread([path files(que(i)).name]);

    tim=[tim im];

end

imshow(tim);

photo2_1_a:

clear;clc;

path='E:\B\附件3\';

files=dir('E:\B\附件3\*.bmp');

fl=length(files);

eig=[];%特征向量矩阵

for i=1:fl

    imi=imread([path,files(i).name]);

    [ro,co]=size(imi);

    thresh=graythresh(imi);%图像二值化;

    imj=im2bw(imi,thresh);

    t=[];

    for j=1:ro

        if min(imj(j,:))==0

            t=[  t 0];

        else

            t=[t 1];

        end

    end

    eig=[eig t'];

end

P=eig;

net=newc(minmax(P),11,0.1);

net.trainParam.epochs=200;

net=init(net);

% net=newsom(P,[11]);

net=train(net,P);

Y=sim(net,P);

y=vec2ind(Y);

% y=kmeans(P',11);

% 测试分类

te=[];

for i=1:11

    tot=0;

    for j=1:209

        if y(j)==i

            tot=tot+1;

        end

    end

    te=[te tot];

end

c1=[];c2=[];c3=[];c4=[];c5=[];c6=[];c7=[];c8=[];c9=[];c10=[];c11=[];

for i=1:fl

    if y(i)==1 c1=[c1 i];

    elseif y(i)==2 c2=[c2 i];

    elseif y(i)==3 c3=[c3 i];

    elseif y(i)==4 c4=[c4 i];

    elseif y(i)==5 c5=[c5 i];

    elseif y(i)==6 c6=[c6 i];

    elseif y(i)==7 c7=[c7 i];

    elseif y(i)==8 c8=[c8 i];

    elseif y(i)==9 c9=[c9 i];

    elseif y(i)==10 c10=[c10 i];

    else c11=[c11 i];

    end

end

c1=sort(c1);c2=sort(c2);c3=sort(c3);c4=sort(c4);c5=sort(c5);c6=sort(c6);c7=sort(c7);c8=sort(c8);c9=sort(c9);c10=sort(c10);c11=sort(c11);

%人工干预

photo2_1_b:

% 求最左边的图片

leftque=[];

cntnum=zeros(fl,1);

maxnum=0;

tmax=0;tpos=0;

for i=1:fl

    imi=imread([path,files(i).name]);

    [ro,co]=size(imi);

    thresh=graythresh(imi);%图像二值化;

    imj=im2bw(imi,thresh);

    tot=0;

    for j=1:co

        if min(imj(:,j))==1

            tot=tot+1;

        else

            break;

        end

    end

    if i==8

        maxnum=tot;

    end

    if tmax<tot

        tmax=tot;tpos=i;

    end

    cntnum(i)=tot;

end

for i=1:fl

    if cntnum(i)>=maxnum

        leftque=[leftque i];

    end

end

leftque;

%构建209*209相似度矩阵

x2smat=ones(209,209);

for i=1:fl

    imi=imread([path,files(i).name]);

%     thresh=graythresh(imi);

%     tmi=im2bw(imi,thresh);

    for j=1:fl

        imj=imread([path,files(j).name]);

%         thresh=graythresh(imj);

%         tmj=im2bw(imj,thresh);

        timi=im2double(imi);

        timj=im2double(imj);

        x2ssum=sum(abs(timi(:,co)-timj(:,1)));

        x2smat(i,j)=x2ssum;

    end

end

photo2_1_c:

%%组合每行碎片

c=[ 8    1  33    46    54    57    69    71    94   127 138   139   154   159   167   175   176   197   209;

    169    2    19    24    27    31    42    51    63    77  87    88   101   121   143   148  180   192  196;

    50    3    12    23    29    55    58    66    92  96   119   130   142   144   179   187   189   191    193;

    15    4    13    32    40    52    74    83   108  116   129   135   136   160   161   170   177   200  204;

    90    5    41   102   103   109   114   115   118  120   124      141   147   152   155   156   186  195   208;

    126    14    17    22    67   107   110   111  140   146   151     158   174   182   183   185   188   198   205 ;

    30     6    11    38    45    49    56    60    65    76   93    99   105   112   172   173   181   202   207;

    62    7    20    21    37    53    64    68    70    73  79    80    97   100   117   132   163   164   178;

    39    9    10    25    26    36   47    75    82   89   104   106   123   131   149   162   168   190 194;

    72    16    18    28    34    61   81    84    86   133   134   153   157   166   171   199   201   203   206;

    95    35    43    44    48    59    78    85    91     98   113   122   125   128   137   145   150   165   184];

p1=c(1,:);p2=c(2,:);p3=c(3,:);p4=c(4,:);p5=c(5,:);p6=c(6,:);p7=c(7,:);p8=c(8,:);

p9=c(9,:);p10=c(10,:);p11=c(11,:);

q1=cocom(p1);%根据相似度求每列的顺序

q2=cocom(p2);q3=cocom(p3);q4=cocom(p4);q5=cocom(p5);q6=cocom(p6);q7=cocom(p7);

q8=cocom(p8);q9=cocom(p9);q10=cocom(p10);q11=cocom(p11);

f1=[8 209 139 159 127  69 176  46 175 1 138 54  57  94 154 71 167  33 197 ];

f2=[169  101  77  63 143  31  42  24 148 192 51   180   121    87 196  27 2    88   19];

f3=[50 55  66 144 187  3  58 193 179 119   191    96    12    23   130    29    92   189   142];

f4=[15   129     4   160    83   200   136    13    74 161   204   170   135    40    32    52   108   116 177];

f5=[90   147   103   155   115    41   152   208   156 141   186   109   118     5   102   114   195   120 124];

f6=[126    14   183   110   198    17   185   111   188  67   107   151    22   174   158   182   205   140  146];

f7=p7(q7);

f8=[62    20    79    68    70   100   163    97   132 80    64   117   164    73     7   178    21    53 37];

f9=[39   149    47   162    25    36    82   190   123 104   131   194    89   168    26     9    10   106  75];

f10=[72   157    84   133   201    18    81    34   203   199    16   134 171   206    86   153   166    28    61];

f11=[95    35    85   184    91    48   122    43   125   145    78   113   150    98   137   165   128    59    44];

tm1=[];tm2=[];tm3=[];tm4=[];tm5=[];tm6=[];tm7=[];tm8=[];tm9=[];tm10=[];tm11=[];

for i=1:19

    tpim=imread([path,files(f1(i)).name]);

    tm1=[tm1 tpim];

    tpim1=imread([path,files(f2(i)).name]);

    tm2=[tm2 tpim1];

    tpim2=imread([path,files(f3(i)).name]);

    tm3=[tm3 tpim2];

    tpim3=imread([path,files(f4(i)).name]);

    tm4=[tm4 tpim3];

    tpim4=imread([path,files(f5(i)).name]);

    tm5=[tm5 tpim4];

    tpim5=imread([path,files(f6(i)).name]);

    tm6=[tm6 tpim5];

    tpim6=imread([path,files(f7(i)).name]);

    tm7=[tm7 tpim6];

    tpim7=imread([path,files(f8(i)).name]);

    tm8=[tm8 tpim7];

    tpim8=imread([path,files(f9(i)).name]);

    tm9=[tm9 tpim8];

    tpim9=imread([path,files(f10(i)).name]);

    tm10=[tm10 tpim9];

    tpim10=imread([path,files(f11(i)).name]);

    tm11=[tm11 tpim10];

end

%%组合11张横向碎片

%first 50 tm3;

[a,b]=size(tm3);

eu=[];ed=[];%tm矩阵的第一行和最后一行;

eu=[eu tm1(1,:);tm2(1,:);tm3(1,:);tm4(1,:);tm5(1,:);tm6(1,:);tm7(1,:);tm8(1,:);tm9(1,:);tm10(1,:);tm11(1,:)];

ed=[ed tm1(a,:);tm2(a,:);tm3(a,:);tm4(a,:);tm5(a,:);tm6(a,:);tm7(a,:);tm8(a,:);tm9(a,:);tm10(a,:);tm11(a,:)];

eu=im2double(eu);

ed=im2double(ed);

ei=zeros(11,11);

for i=1:11

    for j=1:11

        eimat(i,j)=sum(abs(ed(i,:)-eu(j,:)));

    end

end

first=3;

q2ue=[];ff=ones(11,1);

q2ue=[q2ue first];ff(first)=-1;

flag=1;

pos=1;

tf=first;

while(flag==1)

    t=255*255*254;

    tp=0;

    for i=1:11

        if i~=tf && ff(i)~=-1;

            if t>eimat(tf,i)

                t=eimat(tf,i);

                tp=i;

            end

        end

    end

    tf=tp;ff(tf)=-1;

    pos=pos+1;

    q2ue=[q2ue tp];

    if pos==11

        flag=0;

    end

end

% imshow(tm4)

fim=[tm3;tm8;tm2;tm9;tm10;tm4;tm11;tm6;tm7;tm1;tm5;];

imshow(fim);

cocom函数：

function fc=cocom(p)%根据相似度求行的排列顺序

load x2smat.mat

p1=p;tp1=p1;

fc=[1];

lo=1;tp1(1)=-1;

flag=1;tot=0;

while(flag==1)

    tot=tot+1;

    tpo=0;

    min=255*255*255;

    for i=1:19

        if i~=lo && tp1(i)~=-1

            if min>x2smat(p1(lo),p1(i))

                min=x2smat(p1(lo),p1(i));

                tpo=i;

            end

        end

    end

    fc=[fc tpo];

    lo=tpo;

    tp1(lo)=-1;

    if tot==18

        flag=0;

    end

end

% la=0;

% for i=1:19

%     if tp1(i)~=-1

%         la=i;

%     end

% end

% fc=[fc la];

　　Photo2_2_a:

clear;clc;

path='E:\B\附件4\';

files=dir('E:\B\附件4\*.bmp');

fl=length(files);

eig=[];%特征向量矩阵

for i=1:fl

    imi=imread([path,files(i).name]);

    [ro,co]=size(imi);

    thresh=graythresh(imi);%图像二值化;

    imj=im2bw(imi,thresh);

    t=[];

    for j=1:ro

        if sum(imj(j,:))<56

            t=[t 1];

        else

            t=[t 0];

        end

    end

    eig=[eig t'];

end

P=eig;

net=newc(minmax(P),11,0.1);

net.trainParam.epochs=250;

net=init(net);

% net=newsom(P,[11]);

net=train(net,P);

Y=sim(net,P);

y=vec2ind(Y);

% y=kmeans(P',11);

% 测试分类

te=[];

for i=1:11

    tot=0;

    for j=1:209

        if y(j)==i

            tot=tot+1;

        end

    end

    te=[te tot];

end

c1=[];c2=[];c3=[];c4=[];c5=[];c6=[];c7=[];c8=[];c9=[];c10=[];c11=[];

for i=1:fl

    if y(i)==1 c1=[c1 i];

    elseif y(i)==2 c2=[c2 i];

    elseif y(i)==3 c3=[c3 i];

    elseif y(i)==4 c4=[c4 i];

    elseif y(i)==5 c5=[c5 i];

    elseif y(i)==6 c6=[c6 i];

    elseif y(i)==7 c7=[c7 i];

    elseif y(i)==8 c8=[c8 i];

    elseif y(i)==9 c9=[c9 i];

    elseif y(i)==10 c10=[c10 i];

    else c11=[c11 i];

    end

end

c1=sort(c1);c2=sort(c2);c3=sort(c3);c4=sort(c4);c5=sort(c5);c6=sort(c6);c7=sort(c7);c8=sort(c8);c9=sort(c9);c10=sort(c10);c11=sort(c11);

%人工干预

　　Photo2_2_b:

% 求最左边的图片

leftque=[];

cntnum=zeros(fl,1);

maxnum=0;

tmax=0;tpos=0;

for i=1:fl

    imi=imread([path,files(i).name]);

    [ro,co]=size(imi);

    thresh=graythresh(imi);%图像二值化;

    imj=im2bw(imi,thresh);

    tot=0;

    for j=1:co

        if min(imj(:,j))==1

            tot=tot+1;

        else

            break;

        end

    end

    if i==172

        maxnum=tot;

    end

    if tmax<tot

        tmax=tot;tpos=i;

    end

    cntnum(i)=tot;

end

for i=1:fl

    if cntnum(i)>0

        leftque=[leftque i];

    end

end

leftque=[19    20    70    81  86   132    159   171   191   201   208];leftque=leftque+1;

%构建209*209相似度矩阵

x2smat=ones(209,209);

for i=1:fl

    imi=imread([path,files(i).name]);

%     thresh=graythresh(imi);

%     tmi=im2bw(imi,thresh);

    for j=1:fl

        imj=imread([path,files(j).name]);

%         thresh=graythresh(imj);

%         tmj=im2bw(imj,thresh);

        timi=im2double(imi);

        timj=im2double(imj);

        x2ssum=sum(abs(timi(:,co)-timj(:,1)));

        x2smat(i,j)=x2ssum;

    end

end

　　Photo2_2_c:

%%组合每行碎片

c=[

10	11	17	19	36	43	45	56	57	67	75	84	135	146	153	158	172	184	206;

6	25	30	31	38	41	47	52	59	60	87	93	99	108	118	128	151	159	187;

16	21	37	42	44	46	74	77	80	109	117	136	137	144	162	174	180	200	208;

2	32	39	51	54	64	86	98	121	124	130	139	140	154	160	161	176	188	204;

8	22	34	50	55	62	63	113	119	120	134	143	163	169	170	190	193	198	209;

7	18	27	29	79	81	92	95	101	102	104	114	147	149	165	171	197	199	202;

20	23	58	72	83	89	94	106	115	122	127	142	152	156	166	177	183	195	203;

4	14	26	28	35	70	96	111	112	131	133	145	164	167	168	179	182	189	207;

9	15	24	48	61	69	71	85	91	97	100	110	123	138	157	173	175	186	196;

1	13	49	53	73	78	82	88	90	103	116	125	126	129	132	141	178	194	201;

3	5	12	33	40	65	66	68	76	105	107	148	150	155	181	185	191	192	205];

p1=c(1,:);p2=c(2,:);p3=c(3,:);p4=c(4,:);p5=c(5,:);p6=c(6,:);p7=c(7,:);p8=c(8,:);

p9=c(9,:);p10=c(10,:);p11=c(11,:);

q1=cocom(p1);%根据相似度求每列的顺序

q2=cocom(p2);q3=cocom(p3);q4=cocom(p4);q5=cocom(p5);q6=cocom(p6);q7=cocom(p7);

q8=cocom(p8);q9=cocom(p9);q10=cocom(p10);q11=cocom(p11);

f1=[172	43	67	206	11	158	75	146	84	135	56	19	57	36	17	10	184	153	45];

f2=[87	52	108	30	41	159	187	99	25	118	151	6	60	59	93	31	38	47	128];

f3=[21	42	109	117	137	74	37	208	136	16	77	44	200	46	174	80	162	180	144];

f4=[160	140	2	130	64	139	154	54	39	124	121	176	86	51	161	188	98	204	32];

f5=[209	22	8	50	62	120	34	143	169	63	170	55	193	134	119	190	163	198	113];

f6=[202	149	171	197	199	95	114	165	79	104	92	81	102	27	101	7	18	29	147];

f7=[20	195	94	142	89	122	127	106	156	115	177	183	152	23	58	203	72	166	83];

f8=[133	182	96	70	168	164	167	189	112	145	207	4	131	35	14	111	26	28	179];

f9=[71	85	61	15	69	175	138	196	9	48	173	157	97	24	100	123	91	186	110];

f10=[82	78	129	201	132	53	126	141	194	88	90	49	73	13	178	125	1	103	116];

f11=[192	76	12	155	191	185	3	105	181	65	107	5	150	33	205	66	40	68	148];

tm1=[];tm2=[];tm3=[];tm4=[];tm5=[];tm6=[];tm7=[];tm8=[];tm9=[];tm10=[];tm11=[];

for i=1:19

    tpim=imread([path,files(f1(i)).name]);

    tm1=[tm1 tpim];

    tpim1=imread([path,files(f2(i)).name]);

    tm2=[tm2 tpim1];

    tpim2=imread([path,files(f3(i)).name]);

    tm3=[tm3 tpim2];

    tpim3=imread([path,files(f4(i)).name]);

    tm4=[tm4 tpim3];

    tpim4=imread([path,files(f5(i)).name]);

    tm5=[tm5 tpim4];

    tpim5=imread([path,files(f6(i)).name]);

    tm6=[tm6 tpim5];

    tpim6=imread([path,files(f7(i)).name]);

    tm7=[tm7 tpim6];

    tpim7=imread([path,files(f8(i)).name]);

    tm8=[tm8 tpim7];

    tpim8=imread([path,files(f9(i)).name]);

    tm9=[tm9 tpim8];

    tpim9=imread([path,files(f10(i)).name]);

    tm10=[tm10 tpim9];

    tpim10=imread([path,files(f11(i)).name]);

    tm11=[tm11 tpim10];

end

%组合11张横向碎片

%first 50 tm3;

[a,b]=size(tm3);

eu=[];ed=[];%tm矩阵的第一行和最后一行;

eu=[eu tm1(1,:);tm2(1,:);tm3(1,:);tm4(1,:);tm5(1,:);tm6(1,:);tm7(1,:);tm8(1,:);tm9(1,:);tm10(1,:);tm11(1,:)];

ed=[ed tm1(a,:);tm2(a,:);tm3(a,:);tm4(a,:);tm5(a,:);tm6(a,:);tm7(a,:);tm8(a,:);tm9(a,:);tm10(a,:);tm11(a,:)];

eu=im2double(eu);

ed=im2double(ed);

ei=zeros(11,11);

for i=1:11

    for j=1:11

        eimat(i,j)=sum(abs(ed(i,:)-eu(j,:)));

    end

end

first=3;

q2ue=[];ff=ones(11,1);

q2ue=[q2ue first];ff(first)=-1;

flag=1;

pos=1;

tf=first;

while(flag==1)

    t=255*255*254;

    tp=0;

    for i=1:11

        if i~=tf && ff(i)~=-1;

            if t>eimat(tf,i)

                t=eimat(tf,i);

                tp=i;

            end

        end

    end

    tf=tp;ff(tf)=-1;

    pos=pos+1;

    q2ue=[q2ue tp];

    if pos==11

        flag=0;

    end

end

% imshow(tm7)

fim=[tm11;tm6;tm2;tm7;tm4;tm3;tm5;tm9;tm8;tm1;tm10];

imshow(fim);

巴特西

2013B题碎纸片拼接

最新文章

热门文章