wavelet.m

clc;
clear all;
X1=imread（‘1.tif‘）;
X1=double（X1）/256;  %这里转化成double类型，否则使用小波变换输出的会有大量大于1的存在，会导致图像显示有问题
figure（1）;
imshow（X1）title（‘左焦距‘）;
axis square;
X2=imread（‘2.tif‘）;
X2=double（X2）/256;
imshow（X2）title（‘右焦距‘）;
axis square;
[c1s1]=wavedec2（X13‘db3‘）; %将x1进行2维，使用‘sym4’进行变换
sizec1=size（c1）;
for I=1:sizec1（2）
    c1（I）=1.2*c1（I）; %将分解后的值都扩大1.2倍
end
[c2s2]=wavedec2（X23‘db3‘）;
c=c1+c2; %计算平均值
c=0.5*c;
s=s1+s2;
s=0.5*s;
xx=waverec2（cs‘db3‘）; %进行重构
figure;
imshow（xx）title（‘融合后的‘）; 
axis square;
% [cH1cV1cD1]=detcoef2（‘all‘cs1）;%尺度1的所有方向的高频系数
% [cH2cV2cD2]=detcoef2（‘