gailapulasitaxingfenjie.m

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

figure(4)
Lp3l=D2l;
subplot(4,1,1)
image(Lp3l)
colormap(map)
axis([0 258 0 140])
axis off
title('3层：18×33')
F1l=interp2(D2l,'spline');
Lp2l=C2l-F1l;
subplot(4,1,2)
image(Lp2l)
colormap(map)
axis([0 258 0 140])
axis off
title('2层：35×65')
F2l=interp2(C2l,'spline');
F2l(70,129)=0;
Lp1l=B2l-F2l;
subplot(4,1,3)
image(Lp1l)
colormap(map)
axis([0 258 0 140])
axis off
title('1层：70×129')
F3l=interp2(B2l,'spline');
F3l(140,258)=0;
Lp0l=Bl-F3l;
subplot(4,1,4)
image(Lp0l)
colormap(map)
axis([0 258 0 140])
axis off
title('0层：140×258')
figure(5)
Lp3h=D2h;
subplot(4,1,1)
image(Lp3h)
colormap(map)
axis([0 258 0 140])
axis off
title('3层：18×33')
F1h=interp2(D2h,'spline');
Lp2h=C2h-F1h;
subplot(4,1,2)
image(Lp2h)
colormap(map)
axis([0 258 0 140])
axis off
title('2层：35×65')
F2h=interp2(C2h,'spline');
F2h(70,129)=0;
Lp1h=B2h-F2h;
subplot(4,1,3)
image(Lp1h)
colormap(map)
axis([0 258 0 140])
axis off
title('1层：70×129')
F3h=interp2(B2h,'spline');
F3h(140,258)=0;
Lp0h=Bh-F3h;
subplot(4,1,4)
image(Lp0h)
colormap(map)
axis([0 258 0 140])
axis off
title('0层：140×258')