asd.txt

基于matlab的小波变换图像融合的程序代码源

imgPan = imread('D:\Program Files\MATLAB.V7.1\matlab7.0\work\小波1.BMP');
imgMul = imread('D:\Program Files\MATLAB.V7.1\matlab7.0\work\小波3.BMP');
subplot(1,3,1);imshow(imgMul); xlabel ('(a)低分辨率多光谱图像');
subplot(1,3,2); imshow(imgPan); xlabel ('(b)高分辨率全色图像');
%分离TM图像R、G、B分量
mulR = imgMul (:,:,1);
mulG = imgMul (:,:,2);
mulB = imgMul (:,:,3);
% 对SPOT全色图像进行db13小波三层分解
[Cpan,Lpan]=wavedec2(imgPan,3,'db13')
imgWH = Lpan(1,:);
length = imgWH(1)*imgWH(2);
% 对TM图像三分量分别进行db13小波三层分解
[Ctmr,Ltmr] = wavedec2(imgR,3,'db13');
Cr = Cpan; Cr(1:length) = Ctmr(1:length);
[Ctmg,Ltmg] = wavedec2(imgG,3,'db13');
Cg = Cpan; Cg(1:length) = Ctmg(1:length);
[Ctmb,Ltmb] = wavedec2(imgB,3,'db13');
Cb = Cpan; Cb(1:length) = Ctmb(1:length);
% 对TM图像各分量进行db13小波逆变换重构
imgResult(:,:,1) = waverec2(Cr,Lpan,'db13');
imgResult(:,:,2) = waverec2(Cg,Lpan,'db13');
imgResult(:,:,3) = waverec2(Cb,Lpan,'db13');
%写入并显示融合后图像文件
imwrite(uint8(imgResult), ' Merge.BMP');
subplot(1,3,3);imshow(uint8(imgResult), xlabel ('(c)小波融合图像'));