5.4 逆滤波
一、无约束滤波
二、逆滤波
频域上图像退化与恢复过程
H(u,v)=0或很小,N(u,v)不为0难以计算或者比F(u,v)大得多
问题:恢复出来的结果与预期结果相差很大,甚至面目全非。
改进1
令逆滤波器的转移函数为M(u,v)
其中k和d均为小于1的常数
改进2
令逆滤波器的转移函数为M(u,v)
式中D0是逆滤波器的空间截止频率一般选择D0位于H(u,v)通带内某一适当位置使复原图像的信噪比较大。
三、消除匀速运动模糊
求模糊图像的傅里叶变换;
观察图像中感兴趣的物体或目标,分别估计水平方向与垂直
方向的移动距离a和b,按上述公式确定退化转移函数;
计算复原图像的傅里叶变换;
执行傅里叶反变换,得复原图像;
例
[MF,map]=imread('image3-MF.jpg'); %装入运动模糊图像
figure(1);
imshow(MF); %显示模糊图像
LEN=30;
THETA=45;
INITPSF=fspecial('motion',LEN,THETA); %建立复原点扩散函数
[J P]= deconvblind(MF,INITPSF,30); %去卷积
figure(2); imshow(J); %显示结果图像如图
figure(3); imshow(P,[],'notruesize'); %显示复原点扩散函数