mlmlm.m

融合遥感图像

load bust
X1=X;
map1=map;
subplot(221);image(X1);
colormap(map1);title('原始图像1');
axis square
load mask
X2=X;
map2=map;
for i=1:256
for j=1:256
if(X2(i,j)>100)
X2(i,j)=1.2*X2(i,j);
else
X2(i,j)=0.5*X2(i,j);
end
end
end
subplot(222);
image(X2);colormap(map);title('原图像2');
axis square
[c1,s1]=wavedec2(X1,2,'sym4');
sizec1=size(c1);
for i=1:sizec1
c1(i)=1.2*c1(i);
end
[c2,s2]=wavedec2(X2,2,'sym4');
c=c1+c2;
c=0.5*c;
s=s1+s2;
s=0.5*s;
XX=waverec2(c,s,'sym4');
subplot(223);image(XX);title('融合图像');
axis square
%%%%%或者命令 XFUS = wfusimg(X1,X2,'db2',5,'max','max','plot');