sparse

基于压缩感知的图像融合程序 matlab编写

%  本程序实现基于压缩传感的图像图像融合
%  程序作者：赵昊，南京航空航天大学信息工程系
%  输入变量：A是待融合的图像矩阵B是待融合的图像矩阵
%  输出变量：X3是融合后的图像矩阵，大小于AB相同
%  参考程序：沙威，WaveletOMP.m

function X1=CSfusion（AB）

%对原始图像进行预处理
[a0b0]=size（A）;
a01=nextpow2（a0）;
b01=nextpow2（b0）;
if a01>b01
    M1=2^a01;
else
    M1=2^b01;
end
Am=[Arandn（a0M1-b0）;randn（M1-a0M1）];
Bm=[Brandn（a0M1-b0）;randn（M1-a0M1）];


Am=double（Am）;
Bm=double（Bm）;
[ab]=size（Am）;

%  小波变换矩阵生成
ww=DWT（a）;
size（ww）
%  小波变换让图像稀疏化（注意该步骤会耗费时间，但是会增大稀疏度）
A1=ww*sparse（Am）*ww‘;
size（sparse（Am））%skjfskf;lsa;lf
A1=full（A1）;
B1=ww*sparse（Bm）*ww‘;
B1=full（B1）;

%  随机矩阵生成
M=190;
R=randn（Ma）;

%  测量并融合
Y1=R*A1;
Y2=R*B1;
Y=（Y1+Y2）/2;

%  OMP算法
X2=zeros（ab）;  %  恢复矩阵
for i=1:b  %  列循环
    rec=omp（Y（:i）Ra）;
    X2（:i）=rec;
end

%  压缩传感恢

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----

     文件       2376  2010-04-28 08:16  CSfusion.m

     文件       1065  2009-10-14 18:27  DWT.m

----------- ---------  ---------- -----  ----

                 3441                    2