dwt_ihs_fusion 基于离散小波变换和IHS变换的遥感图像融合代码

基于离散小波变换和IHS变换的遥感图像融合代码，主要用于多光谱和全色图像的融合

function dwt = fuse_dwt_his（EZ PAN wname level）
[M N n] = size（EZ）;
K = M * N;

X = reshape（EZ [K n]）;

RGB = X‘;
trans = [       1/3          1/3           1/3
         -sqrt（2）/6   -sqrt（2）/6   2*sqrt（2）/6
          1/sqrt（2）   -1/sqrt（2）             0];
Iv1v2 = trans * RGB;

I = Iv1v2（1 :）;  % H = tan-1（v2/v1） S = sqrt（v1^2+v2^2）. 
I = reshape（I [M N]）;

% 全色图像和I分量直方图匹配.
h = imhist（I 256）;
im2 = match（PAN h）;

[CI S] = wavedec2（I level wname）;
[CPAN S] = wavedec2（im2 level wname）;
 
r = prod（S（1 :））;
Cf = [CI（1:r） CPAN（（r+1）:end）];

temp = waverec2（Cf S wname）;

Iv1v2（1 :） = temp（:）‘;

itrans = [1  -1/sqrt（2）   1/sqrt（2）
          1  -1/sqrt（2）  -1/sqrt（2）
          1     sqrt（2）           0];
Y = itrans * Iv1v2;

dwt = reshape（Y

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----
     文件        1065  2011-06-27 19:49  match.m
     文件         468  2011-06-27 19:49  paramshow.m
     文件        1400  2011-06-27 19:49  psnr.m
     文件         304  2011-06-27 19:49  sdd.m
     文件         899  2011-06-27 19:49  fuse_dwt_his.m