[c0,s0] = wavedec2(M1, N, wtype)；

采用小波分解行

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);