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

压缩感知-单像素相机-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 = Phi\y;
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  cscameradata\Ball_128.mat
     文件        8429  2006-11-15 04:17  cscameradata\Ball_32.mat
     文件       30027  2006-11-15 04:17  cscameradata\Ball_64.mat
     文件       29716  2006-11-15 04:07  cscameradata\Dice_128.mat
     文件        8457  2006-11-15 04:17  cscameradata\Dice_32.mat
     文件       30658  2006-11-15 04:17  cscameradata\Dice_64.mat
     文件        1346  2006-11-15 04:07  cscameradata\Logo_128.mat
     文件        8241  2006-11-15 04:17  cscameradata\Logo_32.mat
     文件       29783  2006-11-15 04:17  cscameradata\Logo_64.mat
     文件       33542  2006-11-15 04:08  cscameradata\Mug_128.mat
     文件        8003  2006-11-15 04:17  cscameradata\Mug_32.mat
     文件       30163  2006-11-15 04:17  cscameradata\Mug_64.mat
     文件    13061523  2006-11-15 04:08  cscameradata\Phi_128.mat
     文件      237246  2006-11-15 04:08  cscameradata\Phi_32.mat
     文件     2936553  2006-11-15 04:08  cscameradata\Phi_64.mat
     文件       29775  2006-11-15 04:08  cscameradata\R_128.mat
     文件        8221  2006-11-15 04:17  cscameradata\R_32.mat
     文件       30768  2006-11-15 04:17  cscameradata\R_64.mat
     文件       32502  2014-06-22 18:00  cscameradata\ball64.fig
     文件        1606  2014-06-22 17:33  cscameradata\cscamera_test.m
     文件        1316  2014-06-22 17:13  cscameradata\cscamera_test_hang.asv
     文件        1316  2014-06-22 17:14  cscameradata\cscamera_test_hang.m
     文件        1219  2014-06-22 17:07  cscameradata\hang.jpg
     文件        2488  2014-06-22 15:52  cscameradata\hang64.jpg
     文件        9673  2006-11-15 04:51  cscameradata\pathdef.m
     文件        2880  2006-04-21 05:33  cscameradata\tvqc_logbarrier.m
     文件        5563  2006-04-21 05:33  cscameradata\tvqc_newton.m