《图像的恢复与重构》由会员分享，可在线阅读，更多相关《图像的恢复与重构（51页珍藏版）》请在金锄头文库上搜索。
1、第五章 图像的恢复与重构,什么是图像退化：图像的质量变坏叫做退化。退化的形式有图像模糊、图像有干扰等。 图像退化的处理方法：无论是由光学、光电或电子方法获得的图像都会有不同程度的退化；退化的形式多种多样。如传感器噪声、摄像机未聚焦、物体与摄像设备之间的相对移动、随机大气湍流、光学系统的相差、成像光源或射线的散射等；如果我们对退化的类型、机制和过程都十分清楚，那么就可以利用其反过程来复原图像。典型的图像复原方法是根据图像退化的先验知识建立一个退化模型，以此模型为基础，采用滤波等手段进行处理，使得复原后的图像符合一定的准则，达到改善图像质量的目的。,一、图像退化模型,f(i, j)：原始图像 g(
2、i, j)：降质图像 T()：成像系统的作用，则： g(x,y)=Tf(x,y) 设T是线性移不变的。 一幅连续的图像f(x,y)可以用抽样函数的二维卷积表示：,因此， 令h(x,;y,) =T (x-,y-) ，则有：,定义于不在原点的二维函数,由于f(, )与x，y没有关系,称 h(x,;y,) 为点扩散函数(PSF)或系统冲击响应。多数情况下它表现为时不变的，反映在图像中为位移不变的，则 h(x,;y,) 可以表示为h(x-,y-),其中*表示卷积运算。如果T()是一个h可分离系统，即： h(x,;y,)=h1 (x,) h2 (y, ) 则二维运算可以分解为列和行两次一维运算来代替。在
3、加性噪声情况下，图像退化模型可以表示为： g(x,y)=f(x,y)*h(x,y)+n(n,y) 其中n(x, y)为噪声图像。,二、离散图像退化模型,对于图像降质过程进行数学建模，设： f(i, j)为原始图像；y(i, j)为降质图像；h(i, j; k, l)为点扩散函数；图像为MN维。有,假设为空间不移变h(i, j; k, l)，则：,线性位移不变的图像退化模型则表示为：g(x,y)=f(x,y)*h(x,y)+n(x,y),结论：如果已知g(x,y)、 n(x,y)、 h(x,y)，则 f(x,y)可以计算出来。对等式两端取傅立叶变换有：G(u,v)=F(u,v)H(u,v)+N(
4、u,v)F(u,v)= (G(u,v)- N(u,v)/ H(u,v)f(x,y)= F-1F(u,v),三、循环矩阵及傅立叶化,一个一维离散序列通过一个系统发生失真的过程可用下图表示,f= H-1g- n,用矩阵表示，可以写成,如果考虑加性噪声，根据离散序列的卷积定理，有,扩展为周期为M的序列,由于离散卷积的周期性，有he(x)=he(x+M)，H可以写成,H是一个循环阵。 结论：离散卷积都可以写成：输入矩阵循环矩阵！对数字图象的二维离散函数也是如此。 对图像退化模型而言，有,A=5,B=5,M=9,M=9,用矩阵形式表示上式：g=Hf+n g、f和n分别表示MN的函数矩阵ge(i, j)、
5、fe(i, j)和ne(i, j)的各行前后相连而成的列矢量(堆叠矢量)。如果假设原始图像是MN维矩阵，则H是MNMN循环矩阵，且H是一个分块(MM个)循环矩阵：,每一个子矩阵Hi自身也是循环矩阵NN：,1、一维信号序列循环矩阵的对角化和傅立叶化,解矩阵方程： f= H-1g- n 最简单的计算方法就是对角化，H H-1也是对角阵。 对角化H的方法求取其特征值和特征矢量。 对循环矩阵而言，设：其有M个特征值和特征矢量。,由于w(k)是由傅立叶系数构成的，因此w(k)彼此是正交的。所以，由w(k)构成的变换矩阵是可逆的。,用特征矢量组成的矩阵：W=w(0),w(1)w(M-1) 生成对角矩阵D：
6、D=W -1HW；且D(k,k)=(k)。而,根据周期性， M i = - i,凑一个常数,所以 H=WDW -1,对g=Hf+n而言，可以写成g=WDW-1f+n ，有W1g=DW-1f+ W1n；其中对 W-1g列矢量的每一行G(k)而言，有,对W-1f列矢量的每一行F(k)而言,对W-1n而言有同样的结果。所以对W1g=DW-1f+ W1n而言，G(u)=MH(u)F(u)+N(u)。上面的过程称之为循环矩阵的傅立叶化。,如果图像的g、f、n采用堆叠矢量的方法构成， g=Hf+n。同一维的情况类似，不同的地方是H为块循环矩阵，以及其中的傅立叶变换是二维的，但最后结论是一样的。 G(u)=
7、MH(u)F(u)+N(u) 一维情况下的结论 G(u,v)=MNH(u,v)F(u,v)+N(u,v)二维情况下的结论 在实际应用中，认为f(x,y)、h(x,y)、g(x,y)、n(x,y)的维数是相等的。 F(u,v) =G(u,v)-N(u,v)/ MNH(u,v) f(x,y) =F-1F(u,v),2、二维信号序列块循环矩阵的对角化和傅立叶化,3、H(u,v)的获取,要知道一个图象降质系统的H(u,v)是一件非常困难的事情。但因为f(x,y)*h(x,y)=g(x,y)，有F(u,v)H(u,v)=G(u,v)如果f(x,y)=(x,y)，F(x,y)=1则H(u,v)=G(u,v
8、) 所以，可以用实验的方法得到h(x,y)和H(u,v)；H(u,v)可用点源的输出图像的傅立叶变换来近似。 另外，有一些图象降质系统的H(u,v)有固定的或近似的数学模型。,四、常见的线性移不变降质算子,运动模糊：通常在拍摄过程中，相机或物体移动造成的运动模糊可以用一维均匀邻域像素灰度的平均值来表示：大气扰动模糊：这种模糊经常出现在遥感和航空摄影中，由于曝光时间过长引起的模糊可用高斯点扩散函数来表示： 式中K是一个归一化常数，保证模糊的大小为单位值，2可以决定模糊的程度。,均匀不聚焦模糊 这是由于相机聚焦不准确引起的，虽然不聚焦由许多参数决定，如相机的焦距、相机光圈的大小、形状、物体和相机之
9、间的距离等，但在研究中为了简单起见，我们用下列函数表示聚焦不准引起的模糊：均匀二维模糊 这是最常见的一种模糊，可以用来近似聚焦不准引起的模糊： 其中L是奇数。,五、无约束恢复,逆滤波 对于图像退化模型： g(x,y)=f(x,y)*h(x,y)+n(x,y) 两边取傅立叶变换：G(u,v)= F(u,v) H(u,v) +N(u,v)H(u,v)又称为系统的转移函数(或滤波函数)，它使图像退化。 在无噪声的情况下，上式可以简化为：G(u,v)= F(u,v) H(u,v)F(u,v) =G(u,v)/H(u,v) 这种1/H(u,v)的形式称为逆滤波。再进行傅立叶逆变换就可以得到f(x,y)。
10、,什么是无约束恢复,当对噪声一无所知时，使n无约束的小。 由于g=Hf+n，假设通过恢复可以得到一个不错的f的估计f。显然，f应满足关系gHfn，我们希望n尽可能的小。于是问题转化为f在什么情况下n最小对矩阵而言就是它的迹(对角线之和)的平方最小。,实际情况中，噪声是不可避免的，因而只能求F(u,v)的估计值：,如果H(u,v)有许多零点，必然使得复原的结果受到极大影响。 或者如果H(u,v)不为零但是有非常小的值，也即病态条件，也会使复原效果受到影响。 解决这个问题的方法是避开H(u,v)的零点。 幸好一般的H(u,v)在低频附近的有限区域内不为零。 因此逆滤波可以在原点附近进行，相当于在频
11、域乘上一低通窗口函数W(u,v)。,为了防止随着u、v的增大H(u,v)的迅速减小而增设一些条件由于截断地原因，被恢复的图象振铃较大。一种改进的方法是取(?),其中，d、k均为小于1的常数。逆滤波的应用条件：退化图像g(x,y)是信噪比较高的图像。,设：,对降质图像g=Hf+n而言，所谓恢复可以看成是对g通过某种线性变换L而得到原图像的估计值f。其中，L是H的某种逆过程。如果f与f之间的差满足某些条件约束条件，例如：满足给定的均方误差，则称该图像被恢复了。,六、有约束恢复(某种准则)的一般公式,噪声是独立的，与图像不相关。求其期望应为0。,期望值越小，其迹也越小,因此,对求迹而言，式中第2项和
12、第3项的相等，有,问题转化为：L等于什么时，均方误差最小。,由此可得最佳L：,0,0,将Rf，Rn看成是一个循环块矩阵，对其进行对角化，有：,图像的相关度约为2030个象素,合理的假设,由最佳L：,可得,两边乘以W和W-1：,又因为当D为对角阵时有,所以,当L为分块循环矩阵时，WL W-1也是一对角阵。其对角线上元素的值为循环序列(“主值”序列)的傅立叶变换值。同理，D，D*，A，B都为对角阵，其对角线上的元素的值为对应循环序列(“主值”序列)的傅立叶变换值。于是上式 可以写成标量形式,合理的假设f(x,y)=g(x,y)*l(x,y),七、维纳滤波恢复,基本原理研究发现，逆滤波复原方法对噪声
13、极为敏感，要求信噪比较高(100以上)，通常不满足该条件。因此希望找到一种方法，在有噪声条件下，从退化图像g(x,y)复原出f(x,y)的估计值f(x,y) ，该估计值应符合一定的准则。用向量f, g, n来表示f(x,y), g(x, y), n(x,y)，Q为对f的线性算子，在约束条件|g-Hf|2=|n|2 下求Qf的最小化|Qf|2而得到f的最佳估计f。,用拉格朗日法求解,设Rf和Rn为f和n的相关矩阵，根据前述分析结果：,若QTQ用Rf1Rn来代替，有：,拉格朗日系数,当D为对角阵，分块循环矩阵为,因此： 写成标量形式为：,统计量的得到是非常困难的一件事。,其中Sf(u,v), Sn
14、(u,v)分别是f(x,y)和n(x,y)的功率谱。,在实际应用中，当噪声情况未知时,八、最小平方(最大平滑)恢复 与维纳滤波恢复的不同在于QTQ的选择。维纳滤波中QTQ用Rf1Rn ，根据前面的推导，其依据的准则为均方误差最小。而最小平方恢复依据的准则为,换言之：在约束条件|g-Hf|2=|n|2条件下，恢复出来的图像的梯度最小(或称最大平滑)。根据此条件选择Q。,由图像增强一章可知，对图像求梯度可用如下操作,如果f(x,y)尺寸是AB，则图像与p(x,y)相卷的结果后的图像为MN，MA+3-1， NB+3-1。构造p(x,y)、f(x,y)的扩展周期图像pe(x,y)与fe(x,y)做卷积
15、。,* * * * * * * * * * * * * * * * * * * * * * * * *,0 -1 0 -1 4 -1 0 -1 0,如果图像的g、f采用堆叠矢量的方法构成， 有g=Cf,对C进行对角化，有 E=W-1CW 在约束条件|g-Hf|2=|n|2条件下，恢复出来的图像的梯度最小(或称最大平滑)。根据此条件选择Q。,令QC根据 有：,左乘W-1，有,根据对角化的讨论，有,与维纳滤波不同此值为已知，不是统计量。,剩下的问题是对s的估测。构造残差矢量r=g-Hf，调节s使其满足|r|2=|n|2a(a是一个准确度系数)。 赋给s某个初始值；计算f和|r|2。 如果满足|r|2=|n|2a，停止计算；如果|r|2|n|2+a，减少s，继续计算。,点对点计算,九、几何失真校正,在不同条件下拍摄的图像，一个物体的图像常会发生几何失真，出现歪斜变形的现象。 从太空中宇航器拍摄的地球上的等距平行线，其图像会变为歪斜或不等距； 用光学和电子扫描仪摄取的图像常会有桶形失真和枕形失真； 用普通的光学摄影与测试雷达拍摄的同一地区的景物二者在几何形状上有较大的差异。 以一幅图像为基准，去校正另一种方式摄入的图像，以期校正其几何失真，就叫做图像的几何失真复原或者几何失真校正。,1、空间变换,
