最大后验实现图像盲复原

% 2D最大后验概率法（空域），卷积函数：g=imfilter （f h ‘conv‘ ‘replicate‘ ）。

clear;close all
f=imread（‘cell.tif‘）;     % 读入图像f
% f=imread（‘boy.bmp‘）;     % 读入图像f
f=im2double（f）;           % 从8 位型转为双精度以便后面的运算
figure（1）imshow（f）;  
R=2;
h=fspecial（‘disk‘R）;     % 创建离散点扩散函数h
g=imfilter（fh‘conv‘‘replicate‘）;  % 原图像f与点扩散函数h卷积，得散焦模糊图像g。加边界处理，所以无黑边
figure（2）imshow（g[]）;   

% *******************************************************
% 按2D最大后验