blinddeconv盲卷积去模糊

在matlabR2017a上运行ok。 盲卷积去模糊，提供了几张图片以供测试，有比较好的效果。 直接运行test_blind_deconv.m文件即可。 在test_blind_deconv.m中修改输入图片路径，在ms_blind_deconv.m文件中输出图片。

function [x y k] = center_kernel_separate（x y k）
  
%
% Center the kernel by translation so that boundary issues are mitigated. Additionally
% if one shifts the kernel the the image must also be shifted in the
% opposite direction.
  
% get centre of mass
mu_y = sum（[1:size（k 1）] .* sum（k 2）‘）;
mu_x = sum（[1:size（k 2）] .* sum（k 1））;    
  
% get mean offset
offset_x = round（ floor（size（k 2） / 2） + 1 - mu_x ）;
offset_y = round（ floor（size（k 1） / 2） + 1 - mu_y ）;

fprintf（‘CenterKernel: weightedMean[%f %f] offset[%d %d]\n‘ mu_x-1 mu_y-1 offset_x offset_y）;

% make kernel to do translation
shift_kernel = zeros（abs（offset_y * 2） + 1 abs（offset_x * 2） + 1）;
shift_kernel（abs（offset_y） + 1 + offset_y abs（offset_x） + 1 + offset_x） = 1;
    
% shift both image and blur kernel
kshift = conv2（k shift_kernel ‘same‘）;
k = kshift;

xshift = conv2（x flipud（fliplr（shift_kernel）） ‘same‘）;
x = xshift;

yshift = conv2（y flipud（fliplr（shift_kernel）） ‘same‘）;
y = yshift;

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----
     目录           0  2018-11-27 20:21  blinddeconv_ok\
     文件      178213  2013-05-24 15:07  blinddeconv_ok\1.jpg
     文件       47109  2013-05-24 15:15  blinddeconv_ok\2.jpg
     文件         988  2011-06-28 00:41  blinddeconv_ok\center_kernel_separate.m
     文件        3169  2011-05-26 09:49  blinddeconv_ok\fast_deconv_bregman.m
     文件       59141  2011-05-26 03:09  blinddeconv_ok\fishes.jpg
     文件     3957468  2011-05-26 03:09  blinddeconv_ok\lyndsey.tif
     文件        5939  2019-01-04 10:25  blinddeconv_ok\ms_blind_deconv.m
     文件      142355  2011-05-26 03:09  blinddeconv_ok\mukta.jpg
     文件       50268  2018-11-22 10:14  blinddeconv_ok\original.PNG
     文件         448  2011-02-03 21:30  blinddeconv_ok\pcg_kernel_core_irls_conv.m
     文件        2194  2011-05-26 09:54  blinddeconv_ok\pcg_kernel_irls_conv.m
     文件     1489716  2011-05-26 03:09  blinddeconv_ok\pietro.tif
     文件         725  2011-06-28 01:41  blinddeconv_ok\README.txt
     文件        6074  2011-01-30 03:04  blinddeconv_ok\solve_image_bregman.m
     文件        3385  2011-06-28 00:49  blinddeconv_ok\ss_blind_deconv.m
     文件        1920  2018-11-27 20:21  blinddeconv_ok\test_blind_deconv.m
     文件        3700  2018-11-26 20:26  blinddeconv_ok\X.jpg
     文件        3687  2013-05-24 15:04  blinddeconv_ok\zhang.jpg