基于正则化的图像超分辨重建matlab代码

matlab代码，基于正则化的图像超分辨重建与处理，用PSNR值确定重建效果

function out = deconvtv（g H mu opts）
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
% out = deconvtvl1（g H mu opts）
% deconvolves image g by solving the following TV minimization problem
%
% min  mu  || Hf - g ||_1 + ||f||_TV
% min mu/2 || Hf - g ||^2 + ||f||_TV
%
% where ||f||_TV = sum_{xyt} sqrt（ a||Dxf||^2 + b||Dyf||^2 + c||Dtf||^2）
% Dxf = f（x+1y t） - f（xyt）
% Dyf = f（xy+1 t） - f（xyt）
% Dtf = f（xy t+1） - f（xyt）
%
% Input:      g      - the observed image can be gray scale color or images
%             H      - point spread function
%            mu      - regularization parameter
%     opts.method    - either ‘l1‘ or {‘l2‘}
%     opts.rho_r     - initial penalty parameter for ||u-Df||   {2}
%     opts.rho_o     - initial penalty parameter for ||Hf-g-r|| {50}
%     opts.beta      - regularization parameter [a b c] for weighted TV norm {[1 1 0]}
%     opts.gamma     - update constant for rho_r {2}
%     opts.max_itr   - maximum iteration {20}
%     opts.alpha     - constant that determines constraint violation {0.7}
%     opts.tol       - tolerance level on relative change {1e-3}
%     opts.print     - print screen option {false}
%     opts.f         - initial  f {g}
%     opts.y1        - initial y1 {0}
%     opts.y2        - initial y2 {0}
%     opts.y3        - initial y3 {0}
%     opts.z         - initial  z {0}
%     ** default values of opts are given in { }.
%
% Output: out.f      - output images
%         out.itr    - total number of iterations elapsed
%         out.relchg - final relative change
%         out.Df1    - Dxf f is the output images
%         out.Df2    - Dyf f is the output images
%         out.Df3    - Dtf f is the output images
%         out.y1     - Lagrange multiplier for Df1
%         out.y2     - Lagrange multiplier for Df2
%         out.y3     - Lagrange multiplier for Df3
%         out.rho_r  - final penalty parameter
% %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
path（pathgenpath（pwd））;

if nargin<3
    error（‘not enough inputs try again \n‘）;
elseif nargin==3
    opts = [];
end

if ~isnumeric（mu）
    error（‘mu must be a numeric value! \n‘）;
end

[rowscolsframes] = size（g）;
memory_condition = memory;
max_array_memory = memory_condition.MaxPossibleArrayBytes/16;
if rows*cols*frames>0.1*max_array_memory
    fprintf（‘Warning: possible memory issue \n‘）;
    reply = input（‘Do you want to continue? [y/n]: ‘ ‘s‘）;
    if isequal（reply ‘n‘）
        out.f = 0;
        return
    end
end

if ~isfield（opts‘method‘）
    method = ‘l2‘;
else
    method = opts.method;
end

switch method
    case ‘l2‘
        out = deconvtvl2（gHmuopts）;
    case ‘l1‘
        out = deconvtvl1（gHmuopts）;
    otherwise
        error（‘unknown method \n‘）;
end

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

     文件       6489  2018-04-26 16:26  deconvtvl2.m

     文件        917  2018-04-25 20:59  Example_image_deblur.m

     文件        962  2018-04-25 20:55  Example_image_denoise.m

     文件        412  2018-04-25 22:02  psnr.m

     文件      76309  2011-01-21 06:52  data\building.jpg

     文件      12708  2011-01-21 07:02  data\data0001.jpg

     文件      12855  2011-01-21 07:02  data\data0002.jpg

     文件      12599  2011-01-21 07:02  data\data0003.jpg

     文件      12720  2011-01-21 07:02  data\data0004.jpg

     文件      12508  2011-01-21 07:02  data\data0005.jpg

     文件      13025  2011-01-21 07:02  data\data0006.jpg

     文件      12927  2011-01-21 07:02  data\data0007.jpg

     文件      12748  2011-01-21 07:02  data\data0008.jpg

     文件      12901  2011-01-21 07:02  data\data0009.jpg

     文件      12601  2011-01-21 07:02  data\data0010.jpg

     文件     387443  2018-04-18 09:52  data\lena_color_512.png

     文件     153785  2016-11-29 00:34  data\lena_gray_512.png

     文件       8806  2011-01-21 06:53  data\wind.jpg

     文件     105816  2018-04-25 20:20  output\guassdenoise.jpg

     文件      45866  2018-04-25 21:00  output\output-001.jpg

     文件       2814  2018-04-25 21:42  deconvtv.m

     文件       6707  2018-04-26 16:26  deconvtvl1.m

     目录          0  2018-05-20 22:29  data

     目录          0  2018-05-20 22:29  output

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

               923918                    24