CASSI基于压缩感知的光谱成像源代码

本代码是基于压缩感知的光谱成像代码，与杜克大学文章想匹配

function y = conv2c（xh）
% Circular 2D convolution
x=wraparound（xh）;
y=conv2（xh‘valid‘）;


function y = wraparound（x m）
% Extend x so as to wrap around on both axes sufficient to allow a
% “valid“ convolution with m to return the cyclical convolution.
% We assume mask origin near centre of mask for compatibility with
% “same“ option.
[mx nx] = size（x）;
[mm nm] = size（m）;
if mm > mx | nm > nx
    error（‘Mask does not fit inside array‘）
end

mo = floor（（1+mm）/2）; no = floor（（1+nm）/2）;  % reflected mask origin
ml = mo-1;            nl = no-1;             % mask left/above origin
mr = mm-mo;           nr = nm-no;            % mask right/below origin
me = mx-ml+1;         ne = nx-nl+1;          % reflected margin in input
mt = mx+ml;           nt = nx+nl;            % top of image in output
my = mx+mm-1;         ny = nx+nm-1;          % output size

y = zeros（my ny）;
y（mo:mt no:nt） = x;      % central region
if ml > 0
    y（1:ml no:nt） = x（me:mx :）;                   % top side
    if nl > 0
        y（1:ml 1:nl） = x（me:mx ne:nx）;            % top left corner
    end
    if nr > 0
        y（1:ml nt+1:ny） = x（me:mx 1:nr）;          % top right corner
    end
end
if mr > 0
    y（mt+1:my no:nt） = x（1:mr :）;                 % bottom side
    if nl > 0
        y（mt+1:my 1:nl） = x（1:mr ne:nx）;          % bottom left corner
    end
    if nr > 0
        y（mt+1:my nt+1:ny） = x（1:mr 1:nr）;        % bottom right corner
    end
end
if nl > 0
    y（mo:mt 1:nl） = x（: ne:nx）;                   % left side
end
if nr > 0
    y（mo:mt nt+1:ny） = x（: 1:nr）;                 % right side
end

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----
     文件          52  2007-11-09 00:31  diffh.m
     文件          53  2007-11-09 00:31  diffv.m
     文件      199356  2008-10-22 16:59  digitalcamerapicture2.jpg
     文件         474  2008-06-04 10:19  mycalltoTVnew.m
     文件         246  2008-10-22 16:13  mywl.mat
     文件        1276  2008-06-04 08:46  projk.m
     文件        1294  2008-10-29 16:12  README.txt
     文件         355  2009-06-02 13:14  Rfuntwist.m
     文件         165  2009-06-02 13:13  RTfuntwist.m
     文件        2276  2009-06-02 13:20  RUNME.m
     文件         200  2008-06-19 16:08  shiftCube.m
     文件      308278  2008-06-09 14:26  staticscene1x150.BMP
     文件      308278  2008-06-09 14:36  staticscenex150xdark.BMP
     文件         143  2008-06-04 10:17  TVnormspectralimaging.m
     文件       22767  2009-06-02 13:01  TwISTmod.m
     文件       98708  2008-01-15 21:41  datacubeplotter\colorMatchFcn.m
     文件         271  2008-06-27 13:59  datacubeplotter\datacubeplottingTest.m
     文件         312  2009-06-02 13:20  datacubeplotter\dispCube.m
     文件        1090  2008-06-27 13:57  datacubeplotter\imagedatacube2.m
     文件        1095  2008-08-06 09:45  datacubeplotter\imagedatacube2gray.m
     文件         586  2008-01-15 21:40  datacubeplotter\ind2rgb8.m
     文件        9832  2008-01-15 21:40  datacubeplotter\ind2rgb8c.c
     文件       11430  2008-02-24 14:19  datacubeplotter\ind2rgb8c.mexa64
     文件        1530  2008-01-15 21:41  datacubeplotter\spectrumRGB.m
     文件        7649  2008-01-15 21:40  datacubeplotter\subimage2.m
     文件        6137  2008-06-13 11:50  datacubeplotter\subplot1.m
     目录           0  2009-06-02 13:14  datacubeplotter\
     文件        1666  2007-07-31 14:37  conv2c.m
     文件    20402793  2008-07-01 16:50  Cu10x1x1x1.mat