%??本程序實現基于壓縮傳感的圖像圖像融合
%??程序作者：趙昊，南京航空航天大學信息工程系
%??輸入變量：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