基于压缩感知的图像重构

文件中包含多种压缩感知图像重构方法，CoSaMp,omp,sp等，能实现图像重构

function CoSa

clc;clear

%  读文件
X=imread（‘lena256.bmp‘）;
X=double（X）;
[ab]=size（X）;

%  小波变换矩阵生成
ww=DWT（a）;

%  小波变换让图像稀疏化（注意该步骤会耗费时间，但是会增大稀疏度）
X1=ww*sparse（X）*ww‘;
X1=full（X1）;

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

%  测量
Y=R*X1;

%  OMP算法
X2=zeros（ab）;  %  恢复矩阵

tic
for i=1:b  %  列循环       
    rec=CS_CoSaMP（Y（:i）Ra）;
    X2（:i）=rec;
end
toc

%  原始图像
figure（1）;
imshow（uint8（X））;
title（‘原始图像‘）;
%  变换图像
figure（2）;
imshow（uint8（X1））;
title（‘小波变换后的图像‘）;
%  压缩传感恢复的图像
figure（3）;
X3=ww‘*sparse（X2）*ww;  %  小波反变换
X3=full（X3）;
imshow（uint8（X3））;
title（‘恢复的图像‘）;
%  误差（PSNR）
errorx=sum（sum（abs（X3-X）.^2））;        %  MSE误差
%psnr=10*log10（255*255/（errorx/a/b））   %  PSNR

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

     文件        830  2018-03-01 17:03  cs\CoSa.m

     文件       1220  2018-03-12 15:11  cs\cs.m

     文件        917  2018-03-12 15:51  cs\cs64.m

     文件       2170  2018-03-01 15:34  cs\CS_CoSaMP.m

     文件       1115  2018-02-28 20:11  cs\DWT.m

     文件       3652  2018-03-08 09:40  cs\lena128.jpg

     文件      66614  2009-06-08 18:00  cs\lena256.bmp

     文件       1388  2018-03-08 14:37  cs\lena64.jpg

     文件       1346  2018-03-02 10:19  cs\omp.m

     文件        842  2018-03-05 15:55  cs\PartHadamardMtx.m

     文件        444  2018-03-07 16:26  cs\SparseRandom.m

     文件        444  2018-03-07 16:30  cs\SparseRandomMtx.m

     目录          0  2018-03-12 17:03  cs

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

                80982                    13