資源簡介
文件中包含多種壓縮感知圖像重構(gòu)方法,CoSaMp,omp,sp等,能實現(xiàn)圖像重構(gòu)
代碼片段和文件信息
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);??%??恢復(fù)矩陣
tic
for?i=1:b??%??列循環(huán)???????
????rec=CS_CoSaMP(Y(:i)Ra);
????X2(:i)=rec;
end
toc
%??原始圖像
figure(1);
imshow(uint8(X));
title(‘原始圖像‘);
%??變換圖像
figure(2);
imshow(uint8(X1));
title(‘小波變換后的圖像‘);
%??壓縮傳感恢復(fù)的圖像
figure(3);
X3=ww‘*sparse(X2)*ww;??%??小波反變換
X3=full(X3);
imshow(uint8(X3));
title(‘恢復(fù)的圖像‘);
%??誤差(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
評論
共有 條評論