5复原2.pptVIP

第五章图像复原与重建 图象复原模糊函数的估计 实际上，图象复原的问题包括两个方面的内容，即： 1．寻求致使图象退化的原因及退化系统参数的估计和检测—Identification； 2．在已知退化系统特性条件下的图象复原本身。 退化系统的确定是复原问题的关键和困难的问题。 图像复原的一般流程 图像复原结果评价的准则： 最小均方准则 加权均方准则 最大熵准则 … … 点扩散函数 如果对系统输入一个单脉冲，其输出称为脉冲响应，表示为h (x, y)，在光学上称之点扩散函数（PSF），它是一个二维点光源的图像。 PSF 的大小和形状，是对光学系统成像性能的度量，受内、外界因素的联合影响。若PSF 的形状较窄则系统及由它产生的影像较好。 若输入信号在空间位移，而输出信号也产生相应位移，但PSF 没有任何变化，则称系统是位移不变的。 PSF 的经验确定 图象退化系统的Identification，可以从观测的退化图象以及原始图象和噪声的统计特性中估计造成图象退化系统的性质，如从已知原始图象（如点或线的位置）的先验知识，获得退化系统的传递函数h (x, y)。 1．如果已知原始景物包括一个尖锐的点，此点在图象中的像即是PSF。 2．如果已知原始景物包括有X方向的细线，则可从这些细线的Y方向图象确定出与线垂直方向的值。 3．如果PSF 是圆对称的，则H(u, v)也是圆对称的，从一个方向的PSF，就可以确定整个PSF。 4．如果已知原始景物包括一个方向尖锐的边缘，可以从边缘图象的微分得到线图象，进行分析。 第一节 图像的退化模型 对于遥感成像系统——二维线性位移不变系统有输入、输出和运算关系： 造成图象退化典型原因有： 1.? 由于成象系统的象差、畸变、带宽有限等造成图象失真； 2.? 由于传感器拍摄姿态和扫描非线性引起的图象几何失真； 3.? 由于传感器与被拍摄景物之间的相对运动，引起所成图象的运动模糊； 4.? 由于光学系统或传感器本身特性不均匀，造成同样亮度景物成象灰度不同——灰度失真； 5.? 由于场景能量传输通道中的介质特性如大气湍流效应、大气成分变化引起图象失真——辐射失真； 6.? 图象在成象、数字化、采集和处理过程中引入的噪声等。 第二节 代数恢复方法 一、图像的无约束复原 二、图像的约束复原 有约束复原：为克服恢复问题的病态性和奇异性，需要在恢复过程中施加某种约束。 有约束最小二乘复原： 令：Q为f 的线性算子，寻找一个f ，使得下面的准则函数最小： J( f ) ＝ Q f 2 ＋ ? g－Hf 2 － n 2 即 ＝ 2 Q’ Q f －2? H’(g－Hf ) ＝0 其中： ? 为拉格朗日系数。 解得： f ＝ （H’H ＋ ? Q’ Q ) -1 H’g ? ＝1 / ? 从而由不同的Q 得到不同的恢复图像。 第三节 频率域恢复方法 因为对于线性位移不变系统，有 g ( x, y ) ＝ h ( x, y )＊f ( x, y )+ n ( x, y ) 对两边进行傅里叶变换： G ( u, v ) ＝ H ( u, v ) F ( u, v ) + N ( u, v ) H ( u, v )称为系统的传递函数。 在无噪的情况理想情况下，对F ( u, v )可写成： F ( u, v ) ＝ G ( u, v ) / H ( u, v ) 1 / H ( u, v )即称为逆滤波器。 有噪声的情况下，则： F ( u, v ) ＝ F ( u, v ) + N ( u, v ) / H ( u, v ) 线性图象复原技术: 逆滤波图象复原: 当矩阵[h](和[h]T )是方块循环矩阵时，可以利用离散富立叶变换矩阵将其对角化，这是因为方块循环矩阵的特征值是退化系统冲激响应的2-D富立叶变换系数，而特征向量是富立叶变换的复指数基函数。因此，方块循环矩阵的特征值和特征向量的分解可以用卷积核的2D-DFT 获得。方块循环矩阵的对角化，向量矩阵表示可变换为频域的代数表示式。 逆滤波图象复原的问题是: 当出现奇异点时； 高频较小时，或当存在噪声时； 逆滤波图象复原将因此而失败。 运动模糊图象的逆滤波复原 当传感器与被摄景物之间存在足够快的相对运动时，所摄