absda.m

拉平拉斯图像融合

figure(6)
for m0=1:140
    for n0=1:258
        if abs(Lp0l(m0,n0))>abs(Lp0h(m0,n0))
            Lp0(m0,n0)=Lp0l(m0,n0);
        else
            Lp0(m0,n0)=Lp0h(m0,n0);
        end
    end
end
subplot(2,2,1)
image(Lp0)
colormap(map)
axis([0 258 0 140])
axis off
title('配准后第0层图像');
for m1=1:70
    for n1=1:129
        if abs(Lp1l(m1,n1))>abs(Lp1h(m1,n1))
            Lp1(m1,n1)=Lp1l(m1,n1);
        else
            Lp1(m1,n1)=Lp1h(m1,n1);
        end
    end
end
subplot(2,2,2)
image(Lp1)
colormap(map)
axis([0 258 0 140])
axis off
title('配准后第1层图像');
for m2=1:35
    for n2=1:65
        if abs(Lp2l(m2,n2))>abs(Lp2h(m2,n2))
            Lp2(m2,n2)=Lp2l(m2,n2);
        else
            Lp2(m2,n2)=Lp2h(m2,n2);
        end
    end
end
subplot(2,2,3)
image(Lp2)
colormap(map)
axis([0 258 0 140])
axis off
title('配准后第2层图像');
for m3=1:18
    for n3=1:33
        if abs(Lp3l(m3,n3))>abs(Lp3h(m3,n3))
            Lp3(m3,n3)=Lp3l(m3,n3);
        else
            Lp3(m3,n3)=Lp3h(m3,n3);
        end
    end
end
subplot(2,2,4)
image(Lp3)
colormap(map)
axis([0 258 0 140])
axis off
title('配准后第3层图像');