DWT小波变换MATLAB

clc; clear all; close all; % 读入图像 imgPan = imread('high.jpg'); imgMul = imread('low.jpg'); % 显示源图像 imshow(imgMul), xlabel ('(a)低分辨率多光谱影像'); figure,imshow(imgPan), xlabel ('(b)高分辨率全色图影像'); mulR = imgMul (:,:,1); mulG = imgMul (:,:,2); mulB = imgMul (:,:,3); % 高分辨率全色图小波变换 [Cpan,Lpan] = wavedec2(imgPan,3,'db13'); imgWH = Lpan(1,:); length = imgWH(1)*imgWH(2); [Ctmr,Ltmr] = wavedec2(mulR,3,' db13' ); % R通道小波变换，得到的低频替换全色图的低频系数 Cr = Cpan; Cr(1:length) = Ctmr(1:length); [Ctmg,Ltmg] = wavedec2(mulG,3,' db13' ); Cg = Cpan; Cg(1:length) = Ctmg(1:length); [Ctmb,Ltmb] = wavedec2(mulB,3,' db13' ); Cb = Cpan; Cb(1:length) = Ctmb(1:length); % 低频为多光谱图的低频系数，高频为全色图的高频系数，进行小波重构 imgResult(:,:,1) = waverec2(Cr,Lpan,' db13' ); imgResult(:,:,2) = waverec2(Cg,Lpan,' db13' ); imgResult(:,:,3) = waverec2(Cb,Lpan,' db13' ); F=uint8(imgResult); figure, imshow(F), xlabel (' (c)小波融合图像' );