matlab置乱图像复原

f=checkerboard(8)

%读图像

psf=fspecial('gaussian',7,10)

%高斯核

sd=0.01%噪声标准差

g=imnoise(imfilter(f,psf),'gaussian',0,sd^2)%图像中加入模糊核噪声

subplot(3,2,1)%三子图第一张

imshow(f),title('(a)')%画原图

subplot(3,2,2)%三子图第二张

imshow(g),title('(b)')%画降质图像

dampar=10*sd%去模糊参数

lim=ceil(size(psf,1)/2)%坐标缩影

weight=zeros(size(g))%去模糊权重初始化

weight(lim+1:end-lim,lim+1:end-lim)=1%去和族厅唤隐模糊权重赋值穗信

numit=5%迭代次数

f5=deconvlucy(g,psf,numit,dampar,weight)%去模糊

subplot(3,2,3)%三子图第三张

imshow(f5),title('(c)')%画复原图像

PSF

=

fspecial('motion',len,ang)

%建立扩散子颤宴，兄洞粗其中len是模糊长度，ang是模糊角度

img2=deconvlucy(img,PSF,n)

%用lucy-richardson方法复原图像，其中img是运动模羡镇糊图像，PSF是扩散子，n是迭代次数，img2是复原图像

A = imread(FILENAME,FMT)

读取图像到A，如果宏亏文件是包含一灰度图像，A是一二维矩阵，

如果文件是包含一真彩色图像（RGB），A是一三维矩阵（轿世M-by-N-by-3）。

FILENAME ：图像文件名；闭绝肢

FMT：图像文件格式；

---