森之张卫东 发表于 2015-8-6 09:50

注意初始化存储矩阵,否则,Matlab默认赋0,找最小值会出错!!!

clc
clear
%load 'C:\Users\Administrator\Desktop\碎纸拼接\程序\汉语\matlab_A.mat'
load 'H:\我们团队做过的题目\训练2\程序\二小问---程序\汉语\新发法2\bi_jiao_ju_zhen_A'
global A
x1=;
x2=;
x3=;
x4=;
x5=;
x6=;
x7=;
x8=;
x9=;
x10=;
x11=;
x12=;
x13=;
x14=;
x1=x1+1;x2=x2+1;x3=x3+1;x4=x4+1;x5=x5+1;x6=x6+1;x7=x7+1;x8=x8+1;
x9=x9+1;x10=x10+1;x11=x11+1;x12=x12+1;x13=x13+1;x14=x14+1;
%为了放到矩阵x中
x6(19)=0;x7(19)=0;x8(19)=0;x9(19)=0;x10(19)=0;x11(19)=0;x12(19)=0;
x13(19)=0;x14(19)=0;
global x
x=;

%构造比较矩阵
line=ones(180,1)*255;

%筛选左边的全白
n_zuo(14,2)=0;
for i=1:14
    t=0;
    for j=1:19
        if x(i,j)~=0
            c=isequal(A(:,1,x(i,j)),line);
            if  c==1
                t=t+1;
                n_zuo(i,t)=x(i,j);
            end
        end
    end
end
n_j_zuo=n_zuo;

%筛选右边的全白
n_you(14,2)=0;
for i=1:14
    t=0;
    for j=1:19
        if x(i,j)~=0
            c=isequal(A(:,72,x(i,j)),line);
            if  c==1
                t=t+1;
                n_you(i,t)=x(i,j);
            end
        end
    end
end
n_j_you=n_you;

%左边分类编程
zuo_liang_ge(1,2)=0;%左边有两个全白。
k_0=1;
zuo_yi_ge(1,1)=0;%左边只有一个是全白的。
k_1=1;
zuo_ling_ge(1,2)=0;%左边无全白的!
k_2=1;
for i=1:14
    k=0;
    for j=1:2
        if n_j_zuo(i,j)==0
            k=k+1; %k=0表示:左边有两个全白
        end        %k=1表示:左边有1个全白
    end
    if k==0%左边有两个全白。
        zuo_liang_ge(k_0,1:2)=n_j_zuo(i,:);
        zuo_liang_ge(k_0,4)=i;
        k_0=k_0+1;
    end
    if k==1%左边只有一个是全白的。
        zuo_yi_ge(k_1,1)=n_j_zuo(i,1);
        zuo_yi_ge(k_1,4)=i;
        k_1=k_1+1;
    end
    if k==2%左边无全白的!
        zuo_ling_ge(k_2,1:2)=n_j_zuo(i,:);
        zuo_ling_ge(k_2,4)=i;
        k_2=k_2+1;
    end
end

%右边分类编程
you_2_ge(1,2)=0;%右边有两个全白。
k_2=1;
you_1_ge(1,1)=0;%右边只有一个是全白的。
k_1=1;
you_0_ge(1,2)=0;%右边无全白的!
k_0=1;
for i=1:14
    k=0;
    for j=1:2
        if n_j_you(i,j)==0
            k=k+1; %k=0表示:右边有两个全白
        end         %k=1表示:右边有1个全白
    end
    if k==0%右边有两个全白。
        you_2_ge(k_2,1:2)=n_j_you(i,:);
        you_2_ge(k_2,4)=i;
        k_2=k_2+1;
    end
    if k==1%右边只有一个是全白的。
        you_1_ge(k_1,1)=n_j_you(i,1);
        you_1_ge(k_1,4)=i;
        k_1=k_1+1;
    end
    if k==2%右边无全白的!
        you_0_ge(k_0,1:2)=n_j_you(i,:);
        you_0_ge(k_0,4)=i;
        k_0=k_0+1;
    end
end



%左边全白为一个的排序
suo_ying=zuo_yi_ge;
bai=suo_ying(:,1);          %左边全白所对应的页数。
bai_row=suo_ying(:,4);  %左边全白所对应的 x 所在的行。
n_pai_xu=length(bai);  %要排几行。
pai_xu(1:n_pai_xu,1)=bai;  %排好之后的矩阵。
%----------------进行排序-----------------------------
flag=72+1;  %做标志位。
A(1,flag,:)=0;
for hang=1:n_pai_xu
    ye=bai(hang,1);          %左边全白所对应的页数。
    x_row=bai_row(hang,1);   %左边全白所对应的 x 所在的行。
    A(1,flag,ye)=1;    %首先,将最左边为全白的做标记。
    B=A(:,72,ye);
    %比较与最左边的最小距离。
   
    for i=1:18  %该循环安排新的最左边。
        d(1,1:19)=inf;  %注意初始化存储矩阵,否则,Matlab默认赋0,找最小值会出错!!!
        for j=1:19
            if x(x_row,j) ~=0
                if A( 1,flag, x(x_row,j) )==1
                    d(1,j)=inf;
                elseif A( 1,flag, x(x_row,j) )==0
                    R(j)= pdist2( double( B )',double( A( :,1, x(x_row,j) ) )' );
                    d(1,j)=sum( R(j) );
                end
            end
        end
        =min( d );
        A(1,flag, x(x_row,wei_zi) )=1;  %找到相邻图片后,做标志位。
        if  x(x_row,i+1)~=0      %判断 x 矩阵的元素是否为零。
            B=A(:,72,x(x_row,wei_zi) );   %并且更新最左边图片位置。
            pai_xu(hang,i+1)=x(x_row,wei_zi);
        elseif  x(x_row,i+1)==0    %如果为零,则将其位置赋inf
            pai_xu(hang,i+1)=inf;
        end
    end
end

%显示图片
for i=1:n_pai_xu
    for j=1:19
        if pai_xu(i,j)~=inf
            B(:, 72*(j-1)+1:72*j ,i)=A(:,1:72,pai_xu(i,j));
        end
    end
    %显示图片
   % figure
   % imshow(B(:,:,i));
end
pai_xu=pai_xu-1;

pai_xu_2_quan_bai=quan_bai_2(zuo_liang_ge);  %全部用第一个
%用 zuo_ling_ge,you_1_ge 先找出最右边为全白的图片
pai_xu_2_wu_bai=quan_bai_0(zuo_ling_ge,you_1_ge);  %全部用第一个




森之张卫东 发表于 2015-8-6 09:51

我的宝贵经验呀!!!
希望大家借鉴。

数学中国YY主管 发表于 2015-8-6 10:39

加油不错。。我到时候可以帮你收录总结一下。宝贵的经验啊

数学中国YY主管 发表于 2015-8-6 10:53

就是文中的符号被写成QQ表情了。。。

页: [1]
查看完整版本: 注意初始化存储矩阵,否则,Matlab默认赋0,找最小值会出错!!!