qiulp.m

拉平拉斯图像融合

figure(9)
m0=0;
n0=0;
E0l(1,1)=2*Lp0l(1,1)^2+Lp0l(1,2)^2+Lp0l(2,1)^2;
E0l(1,258)=2*Lp0l(1,258)^2+Lp0l(1,257)^2+Lp0l(2,258)^2;
E0l(140,1)=2*Lp0l(140,1)^2+Lp0l(139,1)^2+Lp0l(140,2)^2;
if m0==1
    for n0=2:257
        E0l(1,n0)=2*Lp0l(1,n0)^2+Lp0l(1,n0-1)^2+Lp0l(1,n0+1)^2+Lp0l(2,n0)^2;
    end
end
if n0==1
    for m0=2:139
        E0l(m0,1)=2*Lp0l(m0,1)^2+Lp0l(m0-1,1)^2+lp0l(m0+1,1)^2+Lp0l(m0,2)^2;
    end
end
for m0=2:140
    for n0=2:258
        if m0+1==141
            Lp0l(m0+1,n0)=0;
            Lp0h(m0+1,n0)=0;
        elseif n0+1==259
            Lp0l(m0,n0+1)=0;
            Lp0h(m0,n0+1)=0;
        end
        E0l(m0,n0)=2*Lp0l(m0,n0)^2+Lp0l(m0-1,n0)^2+Lp0l(m0+1,n0)^2+Lp0l(m0,n0-1)^2+Lp0l(m0,n0+1)^2;
        E0h(m0,n0)=2*Lp0h(m0,n0)^2+Lp0h(m0-1,n0)^2+Lp0h(m0+1,n0)^2+Lp0h(m0,n0-1)^2+Lp0h(m0,n0+1)^2;
        E0lh(m0,n0)=2*Lp0l(m0,n0)*Lp0h(m0,n0)+Lp0l(m0-1,n0)*Lp0h(m0-1,n0)*+Lp0l(m0+1,n0)*Lp0h(m0+1,n0)+Lp0l(m0,n0-1)*Lp0h(m0,n0-1)+Lp0l(m0,n0+1)*Lp0h(m0,n0+1);
        M0(m0,n0)=2*E0lh(m0,n0)/(E0l(m0,n0)+E0h(m0,n0));
        if M0(m0,n0)<0.75
            if E0l(m0,n0)>E0h(m0,n0)
                Lp0(m0,n0)=Lp0l(m0,n0);
            else
                Lp0(m0,n0)=Lp0h(m0,n0);
            end
        else
            Wmin(m0,n0)=1/2-1/2*(1-M0(m0,n0))/0.25;
            Wmax(m0,n0)=1-Wmin(m0,n0);
            if E0l(m0,n0)>E0h(m0,n0)
                Lp0(m0,n0)=Wmax(m0,n0)*Lp0l(m0,n0)+Wmin(m0,n0)*Lp0h(m0,n0);
            else
                Lp0(m0,n0)=Wmin(m0,n0)*Lp0l(m0,n0)+Wmax(m0,n0)*Lp0h(m0,n0);
            end
        end
        if isnan(Lp0(m0,n0));
            Lp0(m0,n0)=0;
        end
    end
end    
subplot(2,2,1)
image(Lp0);
colormap(map)
axis([0 258 0 140])
axis off
title('配准后第0层图像')
%求LP0140×258
m1=0;
n1=0;
E1l(1,1)=2*Lp1l(1,1)^2+Lp1l(1,2)^2+Lp1l(2,1)^2;
E1l(1,129)=2*Lp1l(1,129)^2+Lp1l(1,128)^2+Lp1l(2,129)^2;
E1l(70,1)=2*Lp1l(70,1)^2+Lp1l(69,1)^2+Lp1l(70,2)^2;
if m1==1
    for n1=2:128
        E1l(1,n1)=2*Lp1l(1,n1)^2+Lp1l(1,n1-1)^2+Lp1l(1,n1+1)^2+Lp1l(2,n1)^2;
    end
end
if n1==1
    for m1=2:69
        E0l(m1,1)=2*Lp1l(m1,1)^2+Lp1l(m1-1,1)^2+lp1l(m1+1,1)^2+Lp1l(m1,2)^2;
    end
end
for m1=2:70
    for n1=2:129
        if m1+1==71
            Lp1l(m1+1,n1)=0;
            Lp1h(m1+1,n1)=0;
        elseif n1+1==130
            Lp1l(m1,n1+1)=0;
            Lp1h(m1,n1+1)=0;
        end
        E1l(m1,n1)=2*Lp1l(m1,n1)^2+Lp1l(m1-1,n1)^2+Lp1l(m1+1,n1)^2+Lp1l(m1,n1-1)^2+Lp1l(m1,n1+1)^2;
        E1h(m1,n1)=2*Lp1h(m1,n1)^2+Lp1h(m1-1,n1)^2+Lp1h(m1+1,n1)^2+Lp1h(m1,n1-1)^2+Lp1h(m1,n1+1)^2;
        E1lh(m1,n1)=2*Lp1l(m1,n1)*Lp1h(m1,n1)+Lp1l(m1-1,n1)*Lp1h(m1-1,n1)*+Lp1l(m1+1,n1)*Lp1h(m1+1,n1)+Lp1l(m1,n1-1)*Lp1h(m1,n1-1)+Lp1l(m1,n1+1)*Lp1h(m1,n1+1);
        M1(m1,n1)=2*E1lh(m1,n1)/(E1l(m1,n1)+E1h(m1,n1));
        if M1(m1,n1)<0.75
            if E1l(m1,n1)>E1h(m1,n1)
                Lp1(m1,n1)=Lp1l(m1,n1);
            else
                Lp1(m1,n1)=Lp1h(m1,n1);
            end
        else
            Wmin(m1,n1)=1/2-1/2*(1-M1(m1,n1))/0.25;
            Wmax(m1,n1)=1-Wmin(m1,n1);
            if E1l(m1,n1)>E1h(m1,n1)
                Lp1(m1,n1)=Wmax(m1,n1)*Lp1l(m1,n1)+Wmin(m1,n1)*Lp1h(m1,n1);
            else
                Lp1(m1,n1)=Wmin(m1,n1)*Lp1l(m1,n1)+Wmax(m1,n1)*Lp1h(m1,n1);
            end
        end
         if isnan(Lp1(m1,n1));
            Lp1(m1,n1)=0;
        end
    end
end
subplot(2,2,2)
image(Lp1);
colormap(map)
axis([0 258 0 140])
axis off
title('配准后第1层图像')
%求LP170×129
m2=0;
n2=0;
E2l(1,1)=2*Lp2l(1,1)^2+Lp2l(1,2)^2+Lp2l(2,1)^2;
E2l(1,65)=2*Lp2l(1,65)^2+Lp2l(1,64)^2+Lp2l(2,65)^2;
E2l(35,1)=2*Lp2l(35,1)^2+Lp2l(34,1)^2+Lp2l(35,2)^2;
if m2==1
    for n2=2:64
        E2l(1,n2)=2*Lp2l(1,n2)^2+Lp2l(1,n2-1)^2+Lp2l(1,n2+1)^2+Lp2l(2,n2)^2;
    end
end
if n2==1
    for m2=2:34
        E2l(m2,1)=2*Lp2l(m2,1)^2+Lp2l(m2-1,1)^2+lp2l(m2+1,1)^2+Lp2l(m2,2)^2;
    end
end
for m2=2:35
    for n2=2:65
        if m2+1==36
            Lp2l(m2+1,n2)=0;
            Lp2h(m2+1,n2)=0;
        elseif n2+1==66
            Lp2l(m2,n2+1)=0;
            Lp2h(m2,n2+1)=0;
        end
        E2l(m2,n2)=2*Lp2l(m2,n2)^2+Lp2l(m2-1,n2)^2+Lp2l(m2+1,n2)^2+Lp2l(m2,n2-1)^2+Lp2l(m2,n2+1)^2;
        E2h(m2,n2)=2*Lp2h(m2,n2)^2+Lp2h(m2-1,n2)^2+Lp2h(m2+1,n2)^2+Lp2h(m2,n2-1)^2+Lp2h(m2,n2+1)^2;
        E2lh(m2,n2)=2*Lp2l(m2,n2)*Lp2h(m2,n2)+Lp2l(m2-1,n2)*Lp2h(m2-1,n2)*+Lp2l(m2+1,n2)*Lp2h(m2+1,n2)+Lp2l(m2,n2-1)*Lp2h(m2,n2-1)+Lp2l(m2,n2+1)*Lp2h(m2,n2+1);
        M2(m2,n2)=2*E2lh(m2,n2)/(E2l(m2,n2)+E2h(m2,n2));
        if M2(m2,n2)<0.75
            if E2l(m2,n2)>E2h(m2,n2)
                Lp2(m2,n2)=Lp2l(m2,n2);
            else
                Lp2(m2,n2)=Lp2h(m2,n2);
            end
        else
            Wmin(m2,n2)=1/2-1/2*(1-M2(m2,n2))/0.25;
            Wmax(m2,n2)=1-Wmin(m2,n2);
            if E2l(m2,n2)>E2h(m2,n2)
                Lp2(m2,n2)=Wmax(m2,n2)*Lp2l(m2,n2)+Wmin(m2,n2)*Lp2h(m2,n2);
            else
                Lp2(m2,n2)=Wmin(m2,n2)*Lp2l(m2,n2)+Wmax(m2,n2)*Lp2h(m2,n2);
            end
        end
        if isnan(Lp2(m2,n2));
            Lp2(m2,n2)=0;
        end 
    end
end
subplot(2,2,3)
image(Lp2);
colormap(map)
axis([0 258 0 140])
axis off
title('配准后第2层图像')
%求LP235×65
m3=0;
n3=0;
E3l(1,1)=2*Lp3l(1,1)^2+Lp3l(1,2)^2+Lp3l(2,1)^2;
E3l(1,33)=2*Lp3l(1,33)^2+Lp3l(1,32)^2+Lp3l(2,33)^2;
E3l(18,1)=2*Lp3l(18,1)^2+Lp3l(17,1)^2+Lp3l(18,2)^2;
if m3==1
    for n3=2:32
        E3l(1,n3)=2*Lp3l(1,n3)^2+Lp3l(1,n3-1)^2+Lp3l(1,n3+1)^2+Lp3l(2,n3)^2;
    end
end
if n3==1
    for m3=2:17
        E3l(m3,1)=2*Lp3l(m3,1)^2+Lp3l(m3-1,1)^2+lp3l(m3+1,1)^2+Lp3l(m3,2)^2;
    end
end
for m3=2:18
    for n3=2:33
        if m3+1==19
            Lp3l(m3+1,n3)=0;
            Lp3h(m3+1,n3)=0;
        elseif n3+1==34
            Lp3l(m3,n3+1)=0;
            Lp3h(m3,n3+1)=0;
        end
        E3l(m3,n3)=2*Lp3l(m3,n3)^2+Lp3l(m3-1,n3)^2+Lp3l(m3+1,n3)^2+Lp3l(m3,n3-1)^2+Lp3l(m3,n3+1)^2;
        E3h(m3,n3)=2*Lp3h(m3,n3)^2+Lp3h(m3-1,n3)^2+Lp3h(m3+1,n3)^2+Lp3h(m3,n3-1)^2+Lp3h(m3,n3+1)^2;
        E3lh(m3,n3)=2*Lp3l(m3,n3)*Lp3h(m3,n3)+Lp3l(m3-1,n3)*Lp3h(m3-1,n3)*+Lp3l(m3+1,n3)*Lp3h(m3+1,n3)+Lp3l(m3,n3-1)*Lp3h(m3,n3-1)+Lp3l(m3,n3+1)*Lp3h(m3,n3+1);
        M3(m3,n3)=2*E3lh(m3,n3)/(E3l(m3,n3)+E3h(m3,n3));
        if M3(m3,n3)<0.75
            if E3l(m3,n3)>E3h(m3,n3)
                Lp3(m3,n3)=Lp3l(m3,n3);
            else
                Lp3(m3,n3)=Lp3h(m3,n3);
            end
        else
            Wmin(m3,n3)=1/2-1/2*(1-M3(m3,n3))/0.25;
            Wmax(m3,n3)=1-Wmin(m3,n3);
            if E3l(m3,n3)>E3h(m3,n3)
                Lp3(m3,n3)=Wmax(m3,n3)*Lp3l(m3,n3)+Wmin(m3,n3)*Lp3h(m3,n3);
            else
                Lp3(m3,n3)=Wmin(m3,n3)*Lp3l(m3,n3)+Wmax(m3,n3)*Lp3h(m3,n3);
            end
        end
        if isnan(Lp3(m3,n3));
            Lp3(m3,n3)=0;
        end 
    end
end
subplot(2,2,4)
image(Lp3);
colormap(map)
axis([0 258 0 140])
axis off
title('配准后第3层图像')
%求LP318×33