基于区域能量的小波图像融合

%图象融合（基于小波变化与局部能量）
clc;
clear all;
Y=imread（‘clockA-t.jpg‘）; %装入原图像
figure（1）
subplot（221）
imshow（Y）;
title（‘clockA-t.jpg‘）
X1=double（Y）/255;
subplot（222）
imhist（X1）;
Y1=imread（‘clockB-t.jpg‘）; %装入原图像
subplot（223）;
imshow（Y1）;
title（‘clockB-t.jpg‘）
X2=double（Y1）/255;
subplot（224）;
imhist（X2）;
%对上述二图像进行三层小波分解
%X1=double（X1）/255;
%X2=double（X2）/255;
[c1s1]=wavedec2（X13‘sym4‘）;
[c2s2]=wavedec2（X23‘sym4‘）;

cA31=appcoef2（c1s1‘sym4‘3）;
 cH31=detcoef2（‘h‘c1s13）;
 cV31=detcoef2（‘v‘c1s13）;
 cD31=detcoef2（‘d‘c1s13）;
 cH21=detcoef2（‘h‘c1s12）;    
 cV21=detcoef2（‘v‘c1s12）;
 cD21=detcoef2（‘d‘c1s12）;
cH11=detcoef2（‘h‘c1s11）;
cV11=detcoef2（‘v‘c1s11）;
cD11=detcoef2（‘d‘c1s11）;


 cA32=appcoef2（c2s2‘sym4‘3）;
 cH32=detcoef2（‘h‘c2s23）;
 cV32=detcoef2（‘v‘c2s23）;
 cD32=detcoef2（‘d‘c2s23）;
 cH22=detcoef2（‘h‘c2s22）;
 cV22=detcoef2（‘v‘c2s22）;
 cD22=detcoef2（‘d‘c2s22）;
cH12=detcoef2（‘h‘c2s21）;
cV12=detcoef2（‘v‘c2s21）;
cD12=detcoef2（‘d‘c2s21）;


 cA33=cA31.*cA32;
 
 Ea=nlfilter（cA31[3 3]@myfunction1）;
 Eb=nlfilter（cA32[3 3]@myfunction1）;
 E=nlfilter（cA33[3 3]@myfunction1）;
 T=1.5;
 for i=1:numel（Ea）
 M（i）=2*E（i）/（Ea（i）+Eb（i））;
 Wmin（i）=0.5-0.5*（（1-M（i））/（1-T））;
 end
Wmax=1-Wmin;
for k=1:n