运动模糊图像盲恢复程序matlab版

function [out_impatch_location]=automatic_patch_selector（impatch_sizeweightsat_mask）
  
% Author: Rob Fergus
% Version: 1.0 distribution code.
% Project: Removing Camera Shake from a Single Image SIGGRAPH 2006 paper
% Copyright 2006 Massachusetts Institute of Technology

  SMOOTH_SIGMA = 3;
  
  %% Assume input image is 0:255
  %% patch_size is integer （odd）
  %% sat_mask is binary & same size as im
    
  %% weight is tuning parameter btw. variance and non-saturated pixels
    
  %%% Get size of input im
  [IIJJ] = size（im）;
  
  %%% Compute centre weighting mask
  [xxyy] = meshgrid（[1:JJ]-round（JJ/2）[1:II]-round（II/2））;
  centre_weight_mask = exp（-weight/（JJ^2）*（xx.^2+yy.^2））;
  II = II*2; JJ = JJ*2;
  %% shift by patch_size
  centre_weight_mask = real（ ifft2（ fft2（centre_weight_maskIIJJ） .* fft2（delta_kernel（patch_size）IIJJ） ） ）;
  
  %%% Get patch mask
  pmask = ones（patch_size）/patch_size.^2;
  
  %%% Find patch with largest variance 
  ei2 = real（ ifft2（ fft2（im.^2IIJJ） .* fft2（pmaskIIJJ） ） ）;
  mu2 = real（ ifft2（ fft2（imIIJJ） .* fft2（pmaskIIJJ） ） ）.^2;
  w = ei2 - mu2;
  
  %%% Compute convolution with sautration mask
  q = real（ ifft2（ fft2（（sat_mask）IIJJ） .* fft2（pmaskIIJJ） ） ）;
  %% q is small if more pixels are available for use 


  combined = （centre_weight_mask）.*w./（q*mean（im（:））.^2+1）; %% more variance less saturation

  %% now find stable maximum （smooth resonse image）
  f = fspecial（‘gaussian‘[8 8]SMOOTH_SIGMA）;
  combined = real（ ifft2（ fft2（combinedIIJJ） .* fft2（fIIJJ） ） ）;

  %%% crop to avoid edge effects
  combined = combined（patch_size:II/2patch_size:JJ/2）;
  
  %%% find max
  [tmpmm] = max（combined（:））;
  [sysx] = ind2sub（size（combined）mm）;

  %%% get coords for axis
  patch_location = [sx sy] -1;
  
  %% chop out patch
  out_im = im（sy-1:sy-2+patch_sizesx-1:sx-2+patch_size）;
  
  
  

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

     文件        358  2007-02-22 16:41  blind_deblur（FergusR）
esultsdelta_kernel.m

     文件       6055  2007-02-22 16:42  blind_deblur（FergusR）
esultsian1.m

     文件   26656662  2010-05-15 13:00  blind_deblur（FergusR）
esultsian1.mat

     文件     254346  2010-05-15 13:00  blind_deblur（FergusR）
esultsian1_blurry.jpg

     文件     447159  2010-05-15 13:00  blind_deblur（FergusR）
esultsian1_final.jpg

     文件     236278  2010-05-15 14:34  blind_deblur（FergusR）
esultsian1模糊核.bmp

     文件       1624  2007-02-22 16:42  blind_deblur（FergusR）priorslinear_street_4.mat

     文件       1624  2007-02-22 16:42  blind_deblur（FergusR）priorslinear_whiteboard_4.mat

     文件     441727  2007-02-22 16:42  blind_deblur（FergusR）imagesian1.jpg

     文件     853977  2007-02-22 16:42  blind_deblur（FergusR）imageslyndsey2.jpg

     文件       1879  2007-02-22 16:41  blind_deblur（FergusR）codeautomatic_patch_selector.m

     文件        556  2007-02-22 16:41  blind_deblur（FergusR）codeclip_image.m

     文件      25600  2011-03-21 10:24  blind_deblur（FergusR）codecode笔记文档.doc

     文件       2527  2007-02-22 16:41  blind_deblur（FergusR）codecreate_greenspan_settings.m

     文件      31970  2011-02-20 18:34  blind_deblur（FergusR）codedeblur.m

     文件        320  2007-02-22 16:47  blind_deblur（FergusR）codedeconvlucy_intens.m

     文件        358  2007-02-22 16:41  blind_deblur（FergusR）codedelta_kernel.m

     文件       2870  2007-02-22 16:41  blind_deblur（FergusR）codeestimate_priors2.m

     文件       1472  2007-02-22 16:41  blind_deblur（FergusR）codeExportFig.m

     文件       5126  2011-02-19 21:44  blind_deblur（FergusR）codefiddle_lucy3.m

     文件       5276  2007-02-22 16:48  blind_deblur（FergusR）codefiddle_lucy4.m

     文件        869  2007-02-22 16:41  blind_deblur（FergusR）codefix_image.m

     文件       5685  2007-02-22 16:41  blind_deblur（FergusR）codeGaussianMixtures1D.m

     文件       1088  2007-02-22 16:41  blind_deblur（FergusR）codegreenspan.m

     文件       1559  2007-02-22 16:41  blind_deblur（FergusR）codehistmatch.m

     文件       5863  2010-05-21 09:25  blind_deblur（FergusR）codeinitialize_parameters2.m

     文件        545  2007-02-22 16:41  blind_deblur（FergusR）codeinvDel2.m

     文件        495  2007-02-22 16:41  blind_deblur（FergusR）codemix_exponentials.m

     文件       2641  2007-02-22 16:41  blind_deblur（FergusR）codemove_level.m

     文件        571  2007-02-22 16:41  blind_deblur（FergusR）code
ormMDpdf.m

............此处省略23个文件信息