da.m

来自「拉平拉斯图像融合」· M 代码 · 共 60 行

for m0=1:140
    for n0=1:258
        if Lp0l(m0,n0)>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 Lp1l(m1,n1)>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 Lp3l(m3,n3)>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层图像');