第07章图像复原

第7章

图像复原内容提要:

7.1图像退化原因与复原技术分类化的数学模型

7.2逆滤波复原

7.3约束复原

7.4非线性复原方法

7.5盲图像复原

7.7几何失真校正

7.7图像修复技术简介

7.8实验：图像复原本章小结

知识要点重点了解图像处理的任务、基本的图像处理系统、微机图像处理系统、数字图像的表示、MATLAB图像处理工具箱的初步使用。图像退化的常见原因图像退化模型图像退化与图像增强的关系线性代数复原非线性复原几何失真校正盲图像复原MATLAB图像处理工具箱去模糊函数

教学建议重点了解数字图像复原的基本任务、图像退化的各种原因、图像复原的常用方法，能够用MATLAB图像处理工具箱解决简单的图像退化问题。先修知识包括：线性代数（循环矩阵的表示；矩阵的广义逆等）信号与线性系统数字信号处理（圆周卷积、离散卷积定理等）随机过程（平稳随机过程等）优化理论。注意本章与“图像增强”一章的联系与区别。

7.1图像退化原因与复原技术分类图像在形成、传输和记录过程中，由于受到多方面的影响，造成图像质量的退化（degradation）。（1）射线辐射、大气湍流等造成的照片畸变。（2）A/D过程会损失部分细节，造成图像质量下降。（3）镜头聚焦不准产生的散焦模糊。（4）成像系统中始终存在的噪声干扰。（5）相机与景物之间的相对运动产生的运动模糊。（7）底片感光、图像显示时会造成记录显示失真。（7）成像系统的像差、非线性畸变、有限带宽。（8）携带遥感仪器的飞机或卫星运动的不稳定，以及地球自转等因素引起的照片几何失真。目的和任务目的在研究图像退化原因的基础上，以退化图像为依据，根据一定的先验知识，建立一个退化模型，然后用相反的运算，以恢复原始景物图像。图像复原要明确规定质量准则衡量接近原始景物图像的程度。图像复原模型可以用连续数学或离散数学处理。图像复原根据退化的数学模型对退化图像进行处理，其实现可在空间域卷积或在频域相乘。在初级视觉处理中的地位在航空航天、国防公安、生物医学、文物修复等领域具有广泛的应用。传统的复原方法基于平稳图像、线性空间不变的退化系统、图像和噪声统计特性的先验知识已知等条件下讨论的现代的复原方法对非平稳图像处理（如卡尔曼滤波）、非线性方法（如神经网络）、信号与噪声的先验知识未知（如盲图像复原）等前提下开展工作。7.1.1连续图像退化的数学模型连续图像退化的一般模型如图7.1所示。输入图像f(x,y)经过一个退化系统或退化算子H(x,y)后产生的退化图像g(x,y)可以表示为：

g(x,y)=H[f(x,y)]

（7.1）

如果仅考虑加性噪声的影响，则退化图像可表示为：g(x,y)=H[f(x,y)]+n(x,y)

（7.2）f(x,y)的最佳估计退化的图像是由成像系统的退化加上额外的系统噪声而形成的。若已知H(x,y)和n(x,y)，图像复原是在退化图像的基础上，作逆运算，得到f(x,y)的一个最佳估计。“最佳估计”而非“真实估计”。由于存在可能导致图像复原的病态性。

导致图像复原的病态性的原因（1）最佳估计问题不一定有解。由于图像复原中可能遇到奇异问题；（2）逆问题可能存在多个解。点扩展函数PSF（Point-spreadFunction）在退化算子H表示线性和空间不变系统的情况下，输入图像f(x,y)经退化后的输出为g(x,y)：

h(x,y)称为退化系统的冲激响应函数。在图像形成的光学过程中，冲激为一光点。又被称为退化系统的点扩展函数PSF。

空间域分析与频率分析退化系统的输出就是输入图像f(x,y)与点扩展函数h(x,y)的卷积，考虑到噪声的影响，即

在频域上可以写成（7.6）G(u,v)、F(u,v)、N(u,v)分别是g(x,y)、f(x,y)、n(x,y)的傅立叶变换

H(u,v)是h(x,y)的傅立叶变换，为系统的传递函数。（7.7）7.1.2离散图像退化的数学模型设f(x,y)大小为A×B，h(x,y)被均匀采样为C×D大小。为避免交叠误差，采用添零延拓的方法，将它们扩展成M=A+C-1和N=B+D-1个元素的周期函数。（7.8a）（7.8b）则输出的降质数字图像为（7.9）二维离散退化模型可以用矩阵形式

g=Hf

（7.10）（7.11）（7.12）给定了退化图像g(x,y)、退化系统的点扩展函数h(x,y)和噪声分布n(x,y)，就可以得到原始图像f的估计。实际计算的工作量十分庞大。（7.13）g=Hf+n

（7.14）通常有两种解决上述问题的途径：假设图像大小M=N，则H的大小为N4，要解出f(x,y)需要解N2个联立方程组。（1）通过对角化简化分块循环矩阵，再利用FFT快速算法可以大大地降低计算量且能极大地节省存储空间。（2）分析退化的具体原因，找出H的具体简化形式。匀速运动造成模糊的PSF就可以用简单的形式表示，这样使复原问题变得简单。各种代数复原方法可能是通过无约束条件而得到原始图像f的估计。也可能是约束复原f。7.2逆滤波复原非约束复原根据对退化系统H和噪声n的了解，已知退化图像g的情况下，在一定的最小误差准则下，得到原始图像f的估计。逆滤波是最早使用的一种无约束复原方法成功地应用于航天器传来的退化图像。n=g－Hf

（7.15）

当对n的统计特性不确定时，希望对原始图像f的估计满足这样的条件，使H在最小二乘意义上近似于g。使得噪声项的范数希望找到一个最小。即目标函数为最小。

在M=N的情况下，H为方阵且H有逆阵H-1，则（7.20）（7.21）当系统H逆作用于退化图像g时，可以得到最小平方意义上的非约束估计。对式（7.20）进行傅立叶变换，则逆滤波法的特点优点：形式简单适用于极高信噪比条件下的图像复原问题，且降质系统的传递函数H不存在病态性质。缺点：具体求解的计算量很大，需要根据循环分块矩阵条件进行简化。当H等于0或接近于0时，还原的图像将变得无意义。这时需要人为对传递函数进行修正，以降低由于传递函数病态而造成的恢复不稳定性。7.3约束复原约束复原除了对降质系统的PSF有所了解外，还需要对原始图像和外加噪声的特性有先验知识。根据不同领域的要求，有时需要对f和n作一些特殊的规定，使处理得到的图像满足某些条件。7.3.1约束复原的基本原理约束最小二乘法复原问题令Q为f的线性算子，要设法寻找一个最优估计，使形式为

的、服从约束条件的函数最小化。7.3.1约束复原的基本原理7.3.2维纳滤波方法7.3.3平滑度约束最小平方滤波7.3.1约束复原的基本原理约束最小二乘法复原问题令Q为f的线性算子，要设法寻找一个最优估计使下面的目标函数为最小

（7.23）(7.24)

式中，γ=α-1

式中，α为拉格朗日乘子,f的最佳估值7.3.2维纳滤波方法最小二乘法滤波复原的核心是如何选择一个合适的变换矩阵Q。Q的形式不同，可得到不同类型的复原方法。选用图像f和噪声n的自相关矩阵Rf和Rn表示Q就可得到维纳滤波复原方法。将f和n近似地看成是平稳随机过程。假设Rf和Rn为f和n的自相关矩阵。Rf

=E{ffT}（7.25a）

Rn

=E{nnT}（7.25b）定义QTQ=R-1fRn，代入式(7.23)，得

（7.26）

假设M=N，Sf和Sn分别为图像和噪声的功率谱，则

（7.27）

分3种情况对式(7.27)作分析：（1）如果γ=1，系统函数Hw(u,v)是维纳滤波器的传递函数。（2）如果γ=0，系统变成单纯的去卷积滤波器，系统的传递函数即为。尽管γ≠0但无噪声影响，Sn(u,v)=0，复原系统亦为理想的逆滤波器，可以看成是维纳滤波器的一种特殊情况。（3）若γ为可调其他参数，此时为参数化维纳滤波器。图7.2去卷积维纳滤波复原方法（a）加性噪声污染和运动模糊图像（b）复原后的图像【例7.1】采用维纳滤波复原的具体实现方法【例7.1】原始图像如图7.2（a），使用函数DECONVWNR对图7.2（b）所示的无噪声模糊图像进行复原重建。

（a）原始图像（b）无噪声模糊图像图7.3原始图像及无噪声模糊图像

图7.4不同PSF产生的复原效果比较（a）使用真实的PSF复原（b）使用较“长”的PSF复原（c）使用较“陡峭”的PSF复原7.3.3平滑度约束最小平方滤波使用逆滤波器一类的方法进行图像复原时，由于退化算子H的病态性质，导致在零点附近数值起伏过大，使复原后的图像产生了人为的噪声和边缘（振铃）。通过选择合理的Q（高通滤波器），并对‖Qf‖2进行优化，可将这种不平滑性降低至最小。使某个函数的二阶导数最小（如Q选用拉普拉斯算子形式表示），可以推导出以平滑度为基础的约束最小二乘方复原方法。7.4非线性复原方法经典复原滤波器的显著特点约束方程和准则函数中的表达式都可以改为矩阵乘法。非线性复原方法所采用的准则函数都不能进行对角化，因而线性代数的方法在这里是不适用的。设S是非线性函数，当考虑图像的非线性退化时，图像的退化模型可以表示成（7.34a）（7.34b）

方法简介7.4.1最大后验复原7.4.2最大熵复原1．正性约束条件2.最大熵复原原理3．Friend和Burg复原方法（1）Friend最大熵复原（2）Burg最大熵复原7.4.3投影复原7.4.4同态滤波复原（重点介绍）7.4.4同态滤波复原自然景物的图像是由照明函数和反射函数两个分量的乘积所组成。同态滤波法复原方法是基于图像的乘性结构理论而提出来的。当降质图像是由两个分量相乘得到时，可先对降质图像取对数，得到两个相加的分量，再进行滤波处理，最后通过指数变换得到复原图像。退化图像g(x,y)可以分为两部分乘积，即

g(x,y)=i(x,y)r(x,y)（7.44）取对数得

logg(x,y)=logi(x,y)+logr(x,y)（7.45）设同态滤波器冲激响应为l(x,y)，其复原结果。同态滤波复原过程可用图7.9所示。图7.9同态滤波器复原同态复原结果=10l(x,y)*[logi(x,y)+logr(x,y)](7.46)同态滤波器的传递函数（7.47）同态滤波技术也可以用于图像增强。采用同态滤波可以实现同态增晰能使图像的灰度动态范围压缩又能使感兴趣的物体图像灰度级扩展。7.5盲图像复原很多情况下难以确定退化的点扩展函数和噪声的统计特性。盲图像复原法是在没有图像退化必要的先验知识的情况下，对观察的图像以某种方式提取出退化信息，采用盲去卷积算法对图像进行复原。对具有加性噪声的模糊图像作盲图像复原的方法有两种：7.5.1直接测量法7.5.2间接估计法【例7.3】盲图像进行重建。（a）原始图像（b）有噪声模糊图像图7.10原始图像及其有噪声模糊化图像（a）较小的PSF（b）较大的PSF（c）真实的PSF图7.11采用不同大小的PSF进行图像复原的效果比较（a）较大的PSF（b）较小的PSF（c）真实的PSF

图7.12不同大小PSF的图像比较7.6几何失真校正图像在获取过程中，由于成像系统的非线性、飞行器的姿态变化等原因，成像后的图像与原景物图像相比，会产生比例失调，甚至扭曲。以上图像退化现象称之为几何失真。有几何畸变的图像不但视觉效果不好而且在对图像进行定量分析时提取的形状、距离、面积等数据也不准确。7.6.1典型的几何失真1．系统失真光学系统、电子扫描系统失真而引起的斜视畸变、枕形、桶形畸变等，都可能使图像产生几何特性失真。原图像；(b)梯形失真；(c)枕形失真；(d)桶形失真图7.13典型的系统几何失真2．非系统失真从飞行器上所获得的地面图像，由于飞行器的姿态、高度和速度变化引起的不稳定与不可预测的几何失真这类畸变一般要根据航天器的跟踪资料和地面设置控制点办法来进行校正。几何畸变校正一般分两步来做几何畸变校正要对失真图像进行精确的几何校正通常是先确定一幅图像为基准，然后去校正另一幅图像的几何形状。第一步：图像空间坐标的变换；第二步：重新确定在校正空间各像素点的取值。7.6.2空间几何坐标变换按照一幅标准图像f(x,y)或一组基准点去校正另一幅几何失真图像g(x‘,y’)。根据两幅图像的一些已知对应点对建立起函数关系式，将失真图像的x‘-y’坐标系变换到标准图像x-y坐标系，从而实现失真图像按标准图像的几何位置校正，使f(x,y)中的每一像点都可在g(x',y')中找到对应像点。（7.57）（a）x'−y'坐标系中的失真图像（b）x−y坐标系中的标准图像