《数字图像处理及应用 含 1CD 教学课件 ppt 作者 李朝晖 主编图像恢复》由会员分享，可在线阅读，更多相关《数字图像处理及应用 含 1CD 教学课件 ppt 作者 李朝晖 主编图像恢复（52页珍藏版）》请在金锄头文库上搜索。
1、第五讲 图像复原,图像退化及复原,什么是图像退化?,图像的质量变坏叫做退化。退化的形式有图像模糊、图像有干扰等,图像退化的原因,无论是由光学、光电或电子方法获得的图像都会有不同程度的退化；,退化的形式多种多样。如传感器噪声、摄像机未聚焦、物体与摄像设备之间的相对移动、随机大气湍流、光学系统的相差、成像光源或射线的散射等；,如果我们对退化的类型、机制和过程都十分清楚，那么就可以利用其反过程来复原图像。,用巴特沃思带阻滤波器复原受正弦噪声干扰的图像 a) 被正弦噪声干扰的图像 b) 滤波效果图,维纳滤波器应用 a) 受大气湍流的严重影响的图像 b) 用维纳滤波器恢复出来的图像,图像恢复：将降质了的
2、图像恢复成原来的图像，针对引起图像退化原因，以及降质过程某先验知识，建立退化模型，再针对降质过程采取相反的方法，恢复图像。,一般地讲，复原的好坏应有一个规定的客观标准，以能对复原的 结果作出某种最佳的估计。,图像还原与增强的区别： 1、图像退化原因决定还原方法,2、评价标准不同：a）突出感兴趣的那部分主观评估 b）利用退化的逆过程恢复原始图像， 客观评估： 接近原图像,无约束恢复 技术 有约束恢复 自动方法 图像恢复： 策略 交互方法 根据是否需要外来干预 空域 处理域 频域 图像一般模型：是线性移不变系统 标准：非线性恢复、线性恢复,5.1 图像的退化模型 5.2 离散图像退化模型 5.3
3、恢复图像的滤波法 5.3.1 非约束还原 5.3.2 有约束还原 5.3.3 非线性约束还原 5.4 图像恢复的相关指标 5.5 点扩散函数的估计 5.6 几何失真矫正,5.1 图像的退化模型,以后讨论中对降质模型H作以下假设： H是线性的（对两个输入函数之和的响应等于它对2个输入图像响应的和）, H是空间（或移位）不变的 对任一个f(x，y)和任一个常数 和都有： H f(x-,y-) = g(x-,y-) 就是说图像上任一点的运算结果只取决于该点的输入值，而与坐标位置无关。,f(i, j)：原始图像 g(i, j)：降质图像 H()： 成像系统的作用，则：,由于 函数的筛选性质（一幅图像可
4、以看作是由一系列点源函数组成的）,系统H的冲激响应，在光学中冲激为一个光点， 退化可 以理解为系统冲激响应造成图像的降质。,其中*表示卷积运算。如果H()是一个h可分离系统，即,则二维运算可以分解为列和行两次一维运算来代替。,在加性噪声情况下，图像退化模型可以表示为,其中n(x, y)为噪声图像,线性位移不变的图像退化模型则表示为：,重要结论：一个线性系统完全可以由它的点扩散函数 来表征。若系统的PSF已知，则系统在（x，y）点的输出响应可看成是不同坐标 处输入函数 所产生的脉冲响应在（x，y）处的叠加。,而在实际降质过程中，降质的另一个复杂因素是随机噪声，考虑有噪声的图像恢复，必经知道噪声统
5、计特性以及噪声和图像信号的相关情况，这是非常复杂的。,实际中假设是白噪声频谱密度为常数，且与图像不相关， （一般只要噪声带宽比图像带宽大得多时，此假设成立的）， 由此得出图像退化模型。,讨论的前提是假设H线性，下面一些恢复方法都是对上述模型 的近似估计。,两边进行付氏变换：,讨论恢复问题： 若略去噪音N，得：,反变换,可求 F f,若H有零点，G也有零点出现，0/0的不定值，这样模型不保证所有逆过程都有解？,由于引起退化的因素众多，而且性质不同，而目前又没有统一的恢复方法，许多人根据不同的物理模型，采用不同的退化模型、处理技巧和估计准则，从而导出了多种恢复方法。 有效方法：针对特定条件，用特定
6、模型处理。,5.2 离散图像退化模型,对于图像降质过程进行数学建模,f(i, j)：原始图像 y(i, j)：降质图像 h(i, j; k, l)：点扩散函数 图像为MN维,假设为空间移不变h(i, j; k, l)，则：,5.3 恢复图像的滤波法,寻找滤波传递函数，通过频域图像滤波得到复原图像的傅立叶变换， 再求反变换，得到复原图像,5.3.1 非约束还原 5.3.2 有约束还原 5.3.3 非线性约束还原,5.3.1 非约束还原 一、逆滤波复原 退化模型： 逆过程：复原图像：,当H(u,v)为0或很小时， ，原点附近： 图像完全被噪声淹没，造成噪声放大,解决方法：1、去除原点 2、设置原点
7、值 3、原点、邻域 均不计算,5.3.2 有约束还原,一、维纳滤波：,Wiener滤波恢复正是在假定图像信号可近似看作平稳随机过程的前 提下，按照使原图像f(x,y)与恢复后的图像 之间的均方误差e2达到 最小的准则，来实现图像恢复的。即：,满足这一要求的转移函数为：,发现： 1）H(u,v)=0，无病态现象，分母不为0 2）SNR高时，同逆滤波器 3）SNR低时，效果不满意。,原因：维纳滤波是基于平稳随机过程模型，且假设退化模型为线性空 间不变系统的原因，这与实际情况存在一定差距。另外，最小 均方误差准则与人的视觉准则不一定匹配,二、约束最小平方滤波图像恢复 准则：以函数平滑为基础 1）使函
8、数的二阶导数为最小。二阶导数是突出图像边缘、轮廓,约束最小二乘需反复迭代才能完成,约束条件：,2）用内积来考察函数f的平滑性,5.3.3 非线性约束还原,最大后验恢复 最大熵恢复,即确定：h(x,y)与n(x,y) g(x,y)=H f(x,y)+n(x,y),5.4 点扩散函数的估计,一、退化的原因为已知 对退化过程有先验知识，如希望能确定PSF和噪声特性,1根据导致模糊的物理过程（先验知识） 大气喘流造成的传递函数 PSF,光学系统散焦传递函数,光学系统散焦时，点光源的像将成圆盘。 从公式可看出，散焦系统的传递函数在以原点为中心，d的倍数为半径处存在零点，形成一些同心的暗环，由散焦图像的频
9、谱上估计出这些同心圆的半径，可得到H（u,v）,均匀聚焦不准-模糊,相机聚焦不准确引起，（不聚焦由许多参数决定：如相机的焦距、相机孔的大小、形状、物体和相机之间的距离等）,在研究中为了简单起见，用下列函数表示聚焦不准引起的模糊：,通常模糊算子相当于一个低通滤波器，因此当模糊算子作用于原始图像时，会引起图像中边缘和轮廓的模糊。77均匀二维模糊算子作用于图像Camera的结果如下图所示,原始图像,77均匀二维模糊算子作用后的结果,模糊后图像任意点的值 ：,特点：图像的频谱在垂直于该方向上存在暗直线，可估出 的大小，运动方向 也可由图像的频谱估计出来,已知：设相机不动，对象运动，运动分量x，y分别为
10、x0(t)，y0(t）相机快 门速度是理想的，快门开启时间（暴光）T。,匀速直线运动模糊下的PSF 相机与景物之间相对运动造成图像降质， H(u,v)运动模糊：,运动模糊图像的恢复处理 a) 原始图像 b) 模糊图像 c) 复原图像,总结：后两种模糊中，可从频谱的零点来估计其有关参数,二、由图像中的点或线估计（后验知识） 1）原始景物中有一清晰的点或点光源。由所成的像得到退化系统的PSF,2）原始景物中确定一条线，成像，由直线产生模糊，根据模糊可以测 定在于边缘垂直方向上的PSF断面曲线，得出一维PSF，如果PSF对称， 旋转一维PSF得到二维PSF,三、由功率谱估计PSF,四、噪声的确定：n
11、(x,y) 分为与图像相关、与图像不相关两类。 还需了解噪声统计性质及其与图像的相关性质。,一般假设： 白色噪声，与图像无关。指图像上各点噪声不相关，其频谱密度为常数。只要噪声带宽远大于图像带宽即可作白噪声处理。,从退化图像大块平坦区中估计，一般不具备噪声先验知识。不同方法要不同特征参数方差，频谱。,5.5 图像恢复的相关指标,模糊信噪比(BSNR, the Blurred Signal-to-Noise Ratio ),表示由模糊和叠加噪声引起的降质程度,ISNR（the Improvement in SNR）,ISNR只是评价图像恢复算法好坏的一个客观指标，ISNR高并不一定主观视觉效果好
12、。,5.6 几何畸变的复原,例如： 从太空中宇航器拍摄的地球上的等距平行线，图像会变为歪斜或不等距； 用光学和电子扫描仪摄取的图像常会有桶形畸变和枕形畸变； 用普通的光学摄影与测试雷达拍摄的同一地区的景物二者在几何形状上有较大的差异。,以一副图像为基准，去校正另一种方式摄入的图像，以校正其几何 畸变，就叫做图像的几何畸变复原或者几何畸变校正。,2 steps: (1)空间变换：对图像平面上的像素进行重新排列以恢复原空间关系 (2)灰度插值：对空间变换后的像素赋予相应的灰度值以恢复原位置的灰度值,事实上，几何校正就是一种几何变换，是图像的几何畸变的反运算 与几何变换类似，几何校正是由输出图像像素
13、坐标反算输入图像坐标，然后通过灰度再采样求出输出像素灰度值。,5.6.1 几何矫正 5.6.2 灰度再采样 5.6.3 几何变形,5.6.1 几何校正,几何基准图像的坐标系统用(x, y)来表示 需要校正的图像的坐标系统用(x, y)表示,设两个图像坐标系统之间的关系用解析式表示,通常h1(x,y)和h2(x,y)用多项式来表示：,通常用线性畸变来近似较小的几何畸变,更精确一些可以用二次型来近似,若基准图像为f(x,y)，畸变图像为g(x,y)，对于景物上的同一个点，假定其灰度不变，则,h1(x,y)和h2(x,y)未知,通常用已知的多对对应点来确定系数a, b,线性畸变,可由基准图找出三个点
14、(x1, y1), (x2, y2), (x3, y3)与畸变图像上三个点(x1, y1), (x2, y2), (x3, y3)一一对应。,将对应点代入，有：,解联立方程组，得出6个系数。,这时,二次型畸变,有12个未知量，需要6对已知对应点,代入上式,记做向量矩阵形式,同样有,解方程组，得到ai, bi. 12个系数。,5.6.2 灰度再采样,Resampling,f(x,y),g(x, y),几何变换是由输出图像像素坐标反算出输入图像坐标，但该坐标并 非整数，需要进行灰度再采样。,再采样是通过灰度插值来完成的,最近邻插值 双线性插值 Nearest Neighbor Bilinear,双
15、线性插值,5.6.3 几何变形,几何变形也是一种几何变换，是在正常图像上人为地生成各种各样的几何畸变，以产生不同的视觉效果。,常见几何变形 简单类 歪斜skew 透视perspective 扭曲distort 复杂类 球面 极坐标 旋涡 波纹,作业：,问题： 1、图像退化的原因是什么？什么是图像恢复？图像恢复与 图像增强的区别？ 2、a）图像的退化模型是什么？对降质模型有哪两个假设？ b) 点扩散函数是如何由冲击函数推导出来的？ 3、离散图像退化模型是什么？ 4、a）逆滤波复原、维纳滤波、约束最小平方滤波图像恢复的 约束条件及表示式是怎样的？ b) 非线性约束还原的两种方法？ 5、大气喘流造成的传递函数 PSF、光学系统散焦传递函数 匀速直线运动模糊下的PSF 6、图像恢复的相关指标有哪些？ 7、几何矫正的两个步骤是什么？为何要进行灰度重采样？ 有哪几种重采样方法？,
