第五章数字图像恢复要点分析.pptVIP

再把y带入，就可得到最终结果： 替换 由运动模糊引起的降质图像的恢复。 a为原图，b为降质图，c为恢复图 沿水平方向匀速运动造成的模糊图像的恢复处理例子。 （a）是模糊图像，（b）是恢复后的图像。 去除由匀速运动引起的模糊 这种由于在图像成像积分时间内摄象机与被摄景物之间产生相对运动，导致图像模糊的现象是我们在实际应用中常常遇到的问题。 对于这种模糊图像的恢复是十分重要的。 * 为 的范数的平方。 * A可用g图像估计 实际降质过程中，降质的另一个复杂因素是随机噪声，考虑有噪声的图像恢复，必需知道噪声统计特性以及噪声和图像信号的相关情况，这是非常复杂的。 ? 假设1：噪声与图像内容不相关 ? 假设2：退化图像可以表示为噪声与图像的叠加 ? 噪声频谱的特点 –随机噪声的频谱通常较宽 –图像内容的带宽总是有限的 由于离散卷积和离散傅里叶变换均是针对周期函数定义的，为了避免离散卷积的周期性序列之间发生相互重叠现象，必须对函数f(x)和h(x )进行周期性延拓，并取M=A+B-1，则 已知退化图像的傅里叶变换和系统的传递函数，就可以求原始图像的傅里叶变换，经傅里叶反变换就可以求得原始图像。 二维离散退化模型 考虑到噪声，离散图像退化模型可以写成如下形式： 写成矩阵形式为： 3.图像的离散退化模型 其中，每个 都是一个N×N 矩阵，由延拓函数 第j 行而来： 可见， 是一个循环矩阵，H是分块循环矩阵。 上述离散退化模型都是在线性空间不变的前提下得出的，这种退化模型已为许多恢复方法所采用，并有良好的恢复效果。 恢复任务：去除图像畸变、补偿图像模糊和减弱噪声效应，以使空间图像恢复。 基本思路：由给定的降质图像g和对退化模型H及噪声n的先验了解，寻找一个对原图像f的最优估计 ，而使事先确定的最优准则为最小。 5.3代数恢复方法 若图像的降质模型为： 对n没有先验知识， 对f的估计 使得： 取值最小，J称为准则或目标函数。 这里： 以上形成了一个对f无约束的复原问题。 5.3.1无约束最小二乘方恢复 常用的最小二乘方准则函数。 求极小值： 令M＝N，且H可逆： 上式表明在最小二乘方准则下寻找出的最优估计图像 可由降质图像 和降质系统冲激响应的逆矩阵 求得。 约束最小二乘法恢复：除了要求知道系统的冲激响应外，还需要知道噪声的统计特性和噪声、图像的相关情况。 5.3.2有约束最小二乘方恢复 1、有约束的最小二乘恢复的方法 若图像的降质模型为： 令Q为 f 的线性算子，在约束 下, 使 为最小的原图f的最佳估计 。 求极小值： 可利用拉格朗日乘数法进行处理。 构造辅助函数： 由此可得恢复步骤为： ① 选取一个 代入 ，把求得的 代入式 ② 当结果大于 ，减小 ，返回①； ③ 当结果小于 ，增加 ，返回①； ④ 重复上述迭代过程，直到式 满足为止。此时 为恢复图像。 之差小于某一阈值 问题的核心转化为求合适的矩阵Q，不同的选择导致不同的恢复方法。 如：Q＝I，能量约束恢复， 维纳滤波图像恢复。 2、维纳滤波图像恢复 设 分别为原始图像f和噪声n的相关矩阵，定义如下： 维纳滤波选择： 获得对原图像的均方误差最小的估计。这时： ＝1为标准维纳滤波， 不为 1时对应含参数的维纳滤波。 可以通过迭代的方法得到。 5.4频率域恢复方法 5.4.1逆滤波恢复法 对于线性移不变系统 对上式两边进行傅立叶变换得 H(u,v)称为系统的传递函数。从频率域角度看，它使图像退化，因而反映了成像系统的性能。 1、 通常在无噪声的理想情况下，上式可简化为 则 进行傅立叶反变换可得到 。 以上就是逆滤波恢复的基本原理。 1/H(u,v)称为逆滤波器。 2、 但实际获取的影像都有噪声，因而只能求F(u,v)的估计值 。 再作傅立叶逆变换得 结论： 1.若噪声为