%?2D最大似然法（空域），卷積函數(shù)：g=imfilter?（f?h?‘conv‘?‘replicate‘?）。

clear
%?f=imread（‘cell.tif‘）;?????%?讀入圖像f
f=imread（‘Lena256.tif‘）;?????%?讀入圖像f
f=im2double（f）;???????????%?從8?位型轉(zhuǎn)為雙精度以便后面的運(yùn)算
figure（1）subplot（131）imshow（f）;??
R=3;
h=fspecial（‘disk‘R）;?????%?創(chuàng)建離散點(diǎn)擴(kuò)散函數(shù)h
%?h=fspecial（‘gaussian‘51）;????%?創(chuàng)建離散點(diǎn)擴(kuò)散函數(shù)h
g=imfilter?（f?h?‘conv‘?‘replicate‘?）;?%?原圖像f與點(diǎn)擴(kuò)散函數(shù)h卷積，得散焦模糊圖像g。加邊界處理，所以無(wú)黑邊
figure（1）subplot（132）imshow（g[]）;???