模糊图像复原代码

I = imread（‘peppers.png‘）;%读取图像
I = I（10+[1:256]222+[1:256]:）;%选定图像的范围
figure（1）;  imshow（I）;%显示原图像
title（‘Original Image‘）;
LEN = 31; THETA = 11;
PSF = fspecial（‘motion‘LENTHETA）;%产生运动模糊函数
blurred = imfilter（IPSF‘circular‘‘conv‘）;%产生运动模糊图像
figure（2）; subplot（221） 
imshow（blurred）;%显示运动模糊图像
title（‘Blurred‘）;
wnr1 = deconvwnr（blurredPSF）;%使用维纳滤波器进行滤波
subplot（222）; imshow（wnr1）;%显示维纳滤波器后滤波的图像
title（‘Restored True PSF‘）;
wnr2 = deconvwnr（blurred...
   fspecial（‘motion‘2*LENTHETA））;%使用两倍长度的点扩散函数进行维纳滤波
subplot（223）; imshow（wnr2）; %显示恢复的图像
title（‘Restored “Long“ PSF‘）;
wnr3 = deconvwnr（blurred...
   fspecial（‘motion‘LEN2*THETA））;%使用两倍角度的点扩散函数进行维纳滤波
subplot（224）; imshow（wnr3）; %显示恢复的图像
title（‘Restored Steep‘）;
noise = 0.1*randn（size（I））;%生成随机噪声
blurredNoisy