采用小波分解行
clear ;
X1 = imread('20170901204704.jpg'); imshow(X1); figure;
M1 = double(X1) / 256; imshow(M1);
N = 4;
wtype = 'sym4';
[c0,s0] = wavedec2(M1, N, wtype);
M1导入后为809*803*3，numel(M1) = 1948881;
而length(c0) = 2011320;
两者并非是按照小波变换的定义数获得的。
下例中，作者认为原来程序的低频部分仅为 c0(1:s1(1,1))和c1(1:s1(1,1)是理解错误，这部分是缺少全部的低频信息。他认为低频应该包括c1(1:s1(1,1)*s1(1,2))和c2(1:s1(1,1)*s1(1,2)。
按照小波变换的定义
c0(1:s1(1,1)) 包含的是cAn的部分信息。
c1(1:s1(1,1)*s1(1,2))部分则包含的是真正的信息。这个与
说说wavedec2函数_hugebawu的博客-CSDN博客_wavedec2函数
说明的一致。
使用MATLAB工具wfusimg函数进行图像的融合 - it610.com
x1=imread('cs1.jpg');
x2=imread('cs2.jpg');
x1=double(x1)/256;
x2=double(x2)/256;
xfus=wfusimg(x1,x2,'sym4',5,'max','max'); %使用wfusimg函数进行图像的融合，其实就是把对低频和高频的处理封装了。
figure;
imshow(xfus);
和上一个方法相比效果要好得多，这是因为分别对低频个高频进行了处理
wfusimg(x1,x2,'sym4',5,'max','max')说明
x1，x2为需要融合的图像，采用小波：sym4，分解为5层，近似信号取两幅图中绝对值最大的值，细节信号取两幅图中绝对值最大值，
可选：max,min,mean,img1,img2,rand.
以下是转载的具体方式：
当然，我们也可以不使用小波工具箱中给出的现成函数，而选择自己动手实现基于小波的图像融合算法。这对我们实际了解该算法是如何进行的大有裨益。下面这段示例程序实现了与前面代码段相同的作用，请读者留意注释说明的部分。另需说明的是，下面这段程序仅仅是为了演示算法实现而编写的，因此并没做异常处理的考虑，我们默认待处理的两幅原始图像的尺寸是一样的。
X1 = imread(‘cathe1.bmp’);
X2 = imread(‘cathe2.bmp’);
M1 = double(X1) / 256;
M2 = double(X2) / 256;
N = 4;
wtype = ‘sym4’;
[c0,s0] = wavedec2(M1, N, wtype);
[c1,s1] = wavedec2(M2, N, wtype);
length = size(c1);
Coef_Fusion = zeros(1,length(2));
%低频系数的处理，取平均值
Coef_Fusion(1:s1(1,1)) = (c0(1:s1(1,1))+c1(1:s1(1,1)))/2;
%处理高频系数，取绝对值大者，这里用到了矩阵乘法
MM1 = c0(s1(1,1)+1:length(2));
MM2 = c1(s1(1,1)+1:length(2));
mm = (abs(MM1)) > (abs(MM2));
Y = (mm.*MM1) + ((~mm).*MM2);
Coef_Fusion(s1(1,1)+1:length(2)) = Y;
%重构
Y = waverec2(Coef_Fusion,s0,wtype);
imshow(Y,[]);
上面这段代码中用到了小波重构函数waverec2()，函数waverec2()其实就是wavedec2()的相反过程，该函数的常用语法形式如下：
X = waverec2(C,S,’wname’)
X = waverec2(C,S,Lo_R,Hi_R)
其中，参数C、S、N和’wname’的意义都与二维小波分解函数wavedec2()中定义的一样，这里不再赘言。Lo_R是重构低通滤波器，Hi_R是重构高通滤波器。
自己修改后的代码：（主要是感觉上面的代码没有把低频部分全部计算进去）
clc;
clear all;
p=imread('12.jpg');
q=imread('21.jpg');
p=double(p)/256;
q=double(q)/256;
imshow(p);
figure;
imshow(q);
figure;
[c1,s1]=wavedec2(p,4,'sym4');
[c2,s2]=wavedec2(q,4,'sym4');
length=length(c1);
hecheng=zeros(1,length);
hecheng(1:s1(1,1)*s1(1,2))=(c1(1:s1(1,1)*s1(1,2))+c2(1:s1(1,1)*s1(1,2))/2;
MM1=c1(s1(1,1)*s1(1,2)+1:length);
MM2=c2(s1(1,1)*s1(1,2)+1:length);
mm=(abs(MM1)>abs(MM2));
Y=(mm.*MM1)+(~mm.*MM2);
hecheng(s1(1,1)*s1(1,2)+1:length)=Y;
Y=waverec2(hecheng,s2,'sym4');
imshow(Y);
