图像复原之逆滤波

人生第一个播客试试水~

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

选用lena女神图像做测试

1.添加运动模糊（不加噪声）

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

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

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

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

figure,imshow(gFilt )

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

2.添加运动模糊（加噪声）

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

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

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

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

to be continued