图像处理(七).docVIP

上次: 一、引言 二、图像运算 三、图像变换 四、图像增强 五、二值形态学及其应用 六、图像复原 图像复原又称图像恢复，是图像处理中的重要技术。所谓的图像复原就是在研究图像退化原因的基础上，以退化图像为依据，根据一定的先验知识设计一种算子，从而估计出理想场景的操作。为了刻画成像系统的特征，通常将成像系统看成是一个线性平移不变系统（因为线性平移不变系统已经有了完整的理论体系，能使用F-变换，便于谱域处理，便于计算机处理，而且能够反映出系统的主要特性），从中推导出物体输入与图像的输出关系的通用数学表达式，从而建立成像系统的退化模型，在此基础上研究图像的复原技术。图像复原在电子监测、医疗摄像等领域具有重要用途。 景物在成像过程中可能出现畸变、模糊、失真或混入噪声，使所得成像图像降质，这种现象称为图像退化。图像复原的关键是在于建立图像退化模型。 1．图像退化模型 图像复原一定是建立在图像退化的数学模型基础上的，这个退化数学模型应该能够反映退化的原因。 f(x,y) f*(x,y) g(x,y) n(x,y) 输入图像图经过一个退化系统，输出图像为，即 叠加系统噪声，就构成了退化图像，即 而 （1.0） 表示的是一副退化图像，式中一般假设系统噪声为加性噪声，且为白噪声。上图表示的就是图像的退化模型。 从上图中可以看出，图像的退化就是成像系统的退化加上额外的系统噪声而形成的，根据这个模型可知，图像复原就是在退化图像的基础上已知成像系统和系统噪声，然后进行反演运算，得到一个的最佳估计。或图像复原就是在退化图像经过复原滤波器，得到恢复的图像，即 由于数字图像都是离散的，所以在实际应用中都是采用的离散形式进行计算，退化图像其表达形式如下： （1.1） 式中，，函数，分别是周期为M和N的函数。是与，具有相同周期的函数，将函数矩阵，和各行堆叠成维列向量分别记为： ，， 则（1.1）式写成矩阵形式为: （1.2） 其中H为维矩阵。H包括个子矩阵，每个子矩阵的大小为，排列顺序如下： 式中每一个子矩阵都是由的第j行元素构成的 上述离散图像退化模型是在线性平移不变性的条件下得到的。图像复原就是根据给定的g，利用H和n来估计f。但是对于实际应用，直接估计出是一项十分繁杂的计算工作，必须寻找简化实现的方法。 （讲） 上面为时域模型，而频域模型为： ，而 （1.3） （讲） 因为点光源可近似为一个单位脉冲，而单位脉冲响应即点冲激响应的F-变换为，可表示为点冲激响应函数的退化模型。点冲激响应的退化函数也称为点扩展函数简称PSF，PSF决定了系统的成像质量。 （讲） 图像复原过程中的主要困难有两点： 其一，要确定引起退化的点扩散函数h(x,y)。分两种情况： 1）如果对于图像缺乏先验知识，则退化过程（模糊和加噪）建立一个模型，然后寻找产生退化的原因或削弱其影响的逆过程。在这种情况下建立的模型只可能是估计性的，所以其复原工作很困难； 2）如果已经具有足够的先验知识，则可以对原始图像建立一个数学模型，并将它与退化图像进行比较，从而找出退化的因素来确定点扩散函数。 其二，如果退化模型是一个病态积分方程。这就是说，在频域中，当H(u,v)（h(x,y)的频域）很小或等于零时，噪声将被放大。换句话说，退化模型的这种病态性质意味着退化图像中小的干扰在H(u,v)取值小的那些频谱上会对图像的恢复产生很大的影响，这给图像复原带来不少麻烦。因此，任何一种图像复原的方法都要考虑当存在病态性时如何控制噪声对复原结果的干扰。 评价估计效果的标准有很多，不同的复原方法实际上就是针对不同的估计评价准则而言的。 例如最小二乘滤波方法原理： (1) 无约束复原方程 ，，求，使 称得到的复原图像为原图像的无约束最小二乘滤波方法复原。 解： 若H可逆， （讲） (2) 有约束复原方程 令Q为的线性操作数向量，求，使 且满足约束条件。称得到的复原图像为原图像的有约束最小二乘滤波方法复原。 解：用Lagrange乘数法解条件极值 取准则函数： 其中为Lagrange乘数。 实际上来调节参数和Q，设法得到最佳的。 （讲） 在matlab中：有四种算法对图像复原处理： 带约束的最小二乘滤波算法 维纳滤波算法 盲去卷积算法 Lucy-Richardson迭代算法 可以通过调用deconvreg函