压缩感知-单像素相机-RICE大学的源代码

% Compressive Imaging Reconstruction Tests
%
% Written by: Marco F. Duarte
%             Digital Signal Processing Group
%             Rice University
%             November 14 2006
%
% Requires l1-Magic available at http://www.l1-magic.org
% NOTE: The ‘Optimization‘ folder in the l1-Magic toolbox must be in the
% path

% 
% CONFIGURABLE PARAMETERS
%
clear all;
% Test picture name. Options include ‘Ball‘‘Mug‘‘Dice‘‘Logo‘ ‘R‘
picname = ‘Ball‘;
% Picture size （sidelength in pixels）. Options include 32 64 128
sidelength = 64;
% Subsampling rate: ratio between M （number of samples） and image size in
% pixels （picsize^2）
subrate = 0.5;
% Optimization parameters
epsilon = 0.01;

%
% script
%

% Load measurements - Column Vector
eval（[‘load ‘ picname ‘_‘ num2str（sidelength）]）;

% Image size
N = sidelength^2;

% Get desired number of samples
y = y（1:（round（N*subrate）））;

% Number of measurements
M = length（y）;

% Load measurement matrix and truncate
% Each row of Phi contains a measurement vector of length N
% and is a raster scan of the 0/1 pattern fed to the DMD
eval（[‘load Phi_‘ num2str（sidelength）]）
Phi = Phi（1:M:）;

% Set initial solution and run TV minimizer
% Output is in column format raster scan of reconstructed image
% Initial point: minimum energy
x0 = Phiy;
estIm = tvqc_logbarrier（x0 Phi [] y epsilon1e-42）;

% Reshape output; flip upside down
xtvqc = reshape（estImsqrt（N）sqrt（N））;
xtvqc = flipud（xtvqc）;

% Normalize output
xtvqc = （xtvqc-min（min（xtvqc）））/（max（max（xtvqc））-min（min（xtvqc）））;

% Show reconstruction
figureimagesc（xtvqc）colormap（gray）axis image

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----
     目录           0  2014-06-23 20:36  cscameradata
     文件       29656  2006-11-15 04:04  cscameradataBall_128.mat
     文件        8429  2006-11-15 04:17  cscameradataBall_32.mat
     文件       30027  2006-11-15 04:17  cscameradataBall_64.mat
     文件       29716  2006-11-15 04:07  cscameradataDice_128.mat
     文件        8457  2006-11-15 04:17  cscameradataDice_32.mat
     文件       30658  2006-11-15 04:17  cscameradataDice_64.mat
     文件        1346  2006-11-15 04:07  cscameradataLogo_128.mat
     文件        8241  2006-11-15 04:17  cscameradataLogo_32.mat
     文件       29783  2006-11-15 04:17  cscameradataLogo_64.mat
     文件       33542  2006-11-15 04:08  cscameradataMug_128.mat
     文件        8003  2006-11-15 04:17  cscameradataMug_32.mat
     文件       30163  2006-11-15 04:17  cscameradataMug_64.mat
     文件    13061523  2006-11-15 04:08  cscameradataPhi_128.mat
     文件      237246  2006-11-15 04:08  cscameradataPhi_32.mat
     文件     2936553  2006-11-15 04:08  cscameradataPhi_64.mat
     文件       29775  2006-11-15 04:08  cscameradataR_128.mat
     文件        8221  2006-11-15 04:17  cscameradataR_32.mat
     文件       30768  2006-11-15 04:17  cscameradataR_64.mat
     文件       32502  2014-06-22 18:00  cscameradataall64.fig
     文件        1606  2014-06-22 17:33  cscameradatacscamera_test.m
     文件        1316  2014-06-22 17:13  cscameradatacscamera_test_hang.asv
     文件        1316  2014-06-22 17:14  cscameradatacscamera_test_hang.m
     文件        1219  2014-06-22 17:07  cscameradatahang.jpg
     文件        2488  2014-06-22 15:52  cscameradatahang64.jpg
     文件        9673  2006-11-15 04:51  cscameradatapathdef.m
     文件        2880  2006-04-21 05:33  cscameradata vqc_logbarrier.m
     文件        5563  2006-04-21 05:33  cscameradata vqc_newton.m