图像的反卷积恢复

2007年ICIP国际会议论文的源码，用于实现对图像的反卷积恢复

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

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

     文件       1666  2007-07-31 14:37  TwIST_v1\conv2c.m

     文件       2593  2007-11-23 18:40  TwIST_v1\demo_l2_l0.m

     文件       4509  2007-11-23 18:40  TwIST_v1\demo_l2_l1.m

     文件       3434  2007-11-23 18:41  TwIST_v1\demo_l2_l1_debias.m

     文件       3560  2007-11-23 18:41  TwIST_v1\demo_l2_TV.m

     文件       3644  2007-11-23 18:42  TwIST_v1\demo_MRI.m

     文件       3644  2007-11-21 00:19  TwIST_v1\demo_Piecewise_cubic_polynomial.m

     文件       4700  2007-11-21 18:59  TwIST_v1\demo_sinusoids_l2_l1.m

     文件       3408  2007-11-23 18:40  TwIST_v1\demo_wave_DWT_deconv.m

     文件         52  2007-11-09 00:31  TwIST_v1\diffh.m

     文件         53  2007-11-09 00:31  TwIST_v1\diffv.m

     文件        226  2007-10-17 14:26  TwIST_v1\hard.m

     文件         96  2007-10-19 17:35  TwIST_v1\l0norm.m

     文件        693  2007-11-09 00:51  TwIST_v1\masked_FFT.m

     文件        435  2007-11-09 00:51  TwIST_v1\masked_FFT_t.m

     文件      98304  2007-11-09 00:31  TwIST_v1\mex_vartotale.dll

     文件       1106  2007-07-31 14:28  TwIST_v1\mex_vartotale.m

     文件        477  2007-11-09 00:52  TwIST_v1\MRImask.m

     文件       2355  2007-11-23 18:47  TwIST_v1\README.txt

     文件         73  2007-10-18 17:08  TwIST_v1\soft.m

     文件         70  2007-11-09 00:31  TwIST_v1\TVnorm.m

     文件      22642  2007-11-20 23:49  TwIST_v1\TwIST.m

     目录          0  2007-12-23 23:56  TwIST_v1

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

               157958                    24