图像复原之逆滤波

人生第一个播客试试水~

实现一个简单的图像逆滤波,网上看到的和自己的想法有点出入,自己实现一下。

选用lena女神图像做测试

图像复原之逆滤波_第1张图片

1.添加运动模糊(不加噪声)

PSF = fspecial('motion',20,32);%len theta
g = imfilter(I, PSF, 'circular');%replicate circular

图像复原之逆滤波_第2张图片

直接逆滤波复原,采用频域相除

   cFT = psf2otf(PSF, size(g));

   gFilt = abs(ifft2(fft2(g)./cFT));  

figure,imshow(gFilt )

图像复原之逆滤波_第3张图片

结论:采用cicular形式的图像逆滤波复原效果较好。

2.添加运动模糊(加噪声)

g_noised=imnoise(g,'gaussian',0,0.0001);

图像复原之逆滤波_第4张图片

cFT = psf2otf(PSF, size(g));
g_noiedFilt = abs(ifft2(fft2(g_noised)./cFT));
g_noiedFilt=mat2gray(g_noiedFilt);%归一化
figure,imshow(g_noiedFilt)

图像复原之逆滤波_第5张图片

结论:可以看出逆滤波对噪声即为敏感,有噪声存在时效果很差。、

3.将带噪声的频谱选择最佳逆滤波半径滤波,优化逆滤波算法

to be continued

你可能感兴趣的:(图像复原,计算机视觉)