1、内容简介
略
522-可以交流、咨询、答疑
2、内容说明
略
3、仿真分析
Image = im2double(rgb2gray(imread( 'flower.jpg')));
[n,m] = size(Image);
LEN = 21; THETA = 11;
PSF= fspecial( 'motion', LEN, THETA);
BlurredI = conv2(PSF, Image);
figure,imshow(BlurredI) ,title('运动模糊图像') ;
[nh, mh] = size(PSF) ;
n=n+ nh-1; m=m + mh- 1;
noise = imnoise( zeros(n, m),'salt & pepper', 0.001);
BlurandnoiseI =BlurredI + noise;
figure,imshow(BlurandnoiseI), title('运动模糊加噪声图像');
h1 = zeros(n,m);
h1 (1:nh, 1: mh) = PSF;
H= fftshift(fft2(h1));
K= sum( noise(:).^2)/sum( Image(:).^2) ;
M= (1./(abs(H).^2 +K) ).^0.5;
G=fftshift(fft2(BlurandnoiseI));
f =ifft2(ifftshift(G.*M));
result = f(1:n- nh+ 1, 1:m- mh+1);
figure, imshow(abs(result)),title( '等功率谱滤波复原图像') ;
4、参考论文
略