DWT小波变换MATLAB

clc;
clear all;
close all;
% 读入图像
imgPan = imread（‘high.jpg‘）;
imgMul = imread（‘low.jpg‘）;

% 显示源图像
imshow（imgMul） xlabel （‘（a）低分辨率多光谱影像‘）;
figureimshow（imgPan） xlabel （‘（b）高分辨率全色图影像‘）;

mulR = imgMul （::1）;
mulG = imgMul （::2）;
mulB = imgMul （::3）;

% 高分辨率全色图小波变换
[CpanLpan] = wavedec2（imgPan3‘db13‘）;
imgWH = Lpan（1:）;
length = imgWH（1）*imgWH（2）;

[CtmrLtmr] = wavedec2（mulR3‘ db13‘ ）; % R通道小波变换，得到的低频替换全色图的低频系数
Cr = Cpan; Cr（1:length） = Ctmr（1:length）;

[CtmgLtmg] = wavedec2（mulG3‘ db13‘ ）;
Cg = Cpan; Cg（1:length） = Ctmg（1:length）;

[CtmbLtmb] = wavedec2（mulB3‘ db13‘ ）;
Cb = Cpan; Cb（1:length） = Ctmb（1:length）;

% 低频为多光谱图的低频系数，高频为全色图的高频系数，进行小波重构
imgResult（::1） = waverec2（CrLpan‘ db13‘ ）;
imgResult（::2） = waverec2（CgLpan‘ db13‘ ）;
imgResult（::3） = waverec2（CbLpan‘ db13‘ ）;

F=uint8（imgResult）;

figure imshow（F） xlabel （‘ （c）小波融合图像‘ ）;