第五讲图像复原

第五讲图像复原第1页，课件共72页，创作于2023年2月图像复原的目标：改善给定的图像。图像复原方法是利用退化的某种先验知识来重建或复原被退化的图像，因而图像复原技术是把退化模型化，并采用相反的过程进行处理，以便复原出原图像。本讲主要介绍退化噪声模型，空间域和频率域复原，退化函数的估计以及几种不同的滤波（逆滤波、维纳滤波（最小均放误差滤波）和几何均值滤波等）。第2页，课件共72页，创作于2023年2月图像复原技术是通过去模糊函数去除图像模糊，获得原始图像的最佳近似估计。图像复原与增强的区别和联系：二者都是在某一个最佳准则下，通过特定的处理而产生期望的最佳结果。复原主要是一种客观的过程，而增强主要是主观的过程。增强技术基本是一个探索的过程，是为了人类视觉系统的生理接受特点而设计一种改善图像的方法。如对比度拉伸提供一种可接受的视觉图像。第3页，课件共72页，创作于2023年2月

噪声退化函数H复原滤波+退化复原图像退化/复原过程模型一、图像退化/复原过程的模型退化过程可以模型化为一个退化函数和一个加性噪声项。第4页，课件共72页，创作于2023年2月退化图像的空间域描述：频率域描述：本讲中图像均是基于此退化模型而言的！第5页，课件共72页，创作于2023年2月二、噪声模型噪声来源：图像的获取和传输过程中。如传感器的工作环境和自身的质量：如CCD摄像机的光照和传感器的温度；传输过程中受到的干扰等。几种典型噪声：高斯噪声：瑞利噪声：第6页，课件共72页，创作于2023年2月伽马（爱尔兰）噪声：指数分布噪声：脉冲噪声（椒盐噪声）：第7页，课件共72页，创作于2023年2月几种典型概率密度函数示意图高斯瑞利伽马均匀指数脉冲第8页，课件共72页，创作于2023年2月直方图直方图直方图样本噪声图像和它们的直方图高斯瑞利伽马第9页，课件共72页，创作于2023年2月直方图直方图直方图唯一视觉可见的噪声类型指数均匀椒盐第10页，课件共72页，创作于2023年2月周期噪声：图形获取中从电力或机电干扰中产生的正弦噪声污染的图像正弦噪声污染的图像傅里叶谱（与一个正弦波对应的每一对共轭脉冲）第11页，课件共72页，创作于2023年2月噪声参数的估计：周期噪声参数是通过检测傅里叶谱来估计的。噪声PDF参数的估计：第12页，课件共72页，创作于2023年2月二、噪声存在下的唯一空间滤波复原

1、均值滤波算术均值滤波几何均值滤波谐波均值滤波逆谐波均值滤波第13页，课件共72页，创作于2023年2月例一、均值滤波对高斯噪声的作用电路板的X射线图高斯(均值=0、方差=400)噪声污染后的电路板X射线图第14页，课件共72页，创作于2023年2月3x3算术均值滤波器滤波结果3x3几何均值滤波器滤波结果清晰模糊第15页，课件共72页，创作于2023年2月例二、逆谐波滤波器对“椒盐”噪声的作用以0.1的概率被“胡椒”污染的图像以0.1的概率被“盐”污染点的图像第16页，课件共72页，创作于2023年2月用3x3，阶数为1.5的逆谐波滤波器滤波图用3x3，阶数为-1.5的逆谐波滤波器滤波图暗区模糊背景清晰背景模糊暗区清晰去噪效果很好第17页，课件共72页，创作于2023年2月例三、在逆谐波均值滤波器中错误选择符号的结果用3x3Q=-1.5逆滤波器滤波的结果图（对p=0.1胡椒噪声污染的图像的处理）用3x3Q=1.5逆滤波器滤波的结果图（对P=0.1盐噪声污染的图像的处理）第18页，课件共72页，创作于2023年2月2、顺序统计滤波器中值滤波最大值和最小值滤波器中点滤波修正后的阿尔法均值滤波第19页，课件共72页，创作于2023年2月例四、中值滤波器对“椒盐”噪声的作用用概率Pa=Pb=0.1椒盐噪声污染的图像用3x3中值滤波器滤波的图像效果好一些噪声第20页，课件共72页，创作于2023年2月用3x3进行二次中值滤波器滤波的图像用3x3进行三次中值滤波器滤波的图像微少噪声完全滤除多次使用中值滤波器会带来图像的模糊化第21页，课件共72页，创作于2023年2月例五、最大值和最小值滤波器对“椒盐”噪声的作用用3x3最大值滤波器滤波的图像用3x3最小值滤波器滤波的图像亮色变小暗色变大亮色变大暗色变小第22页，课件共72页，创作于2023年2月例六、几种滤波器的滤波性能比较加性均匀噪声(均值为0方差为800)污染的图像在加性均匀噪声基础上叠加Pa=Pb=0.1的椒盐噪声污染的图像第23页，课件共72页，创作于2023年2月5x5算术均值滤波器5x5几何均值滤波器5x5中值滤波器5x5，d=5的修正后的阿尔法均值滤波器作用不明显更平滑第24页，课件共72页，创作于2023年2月3、自适应滤波器：噪声图像在（x，y）点的值：干扰f（x，y）以形成g（x，y）的噪声方差：在Sxy上象素点的局部方差：在Sxy上象素点的局部均值自适应、局部噪声消除滤波器第25页，课件共72页，创作于2023年2月例七、自适应局部噪声消除滤波零均值方差1000的加性高斯噪声污染的图像7x7算术均值滤波器噪声平滑图像模糊第26页，课件共72页，创作于2023年2月7x7几何均值滤波器7x7自适应噪声削减滤波（噪声方差为1000）噪声平滑模糊减小噪声平滑图像锐化第27页，课件共72页，创作于2023年2月自适应中值滤波器：处理更大概率的冲激噪声，平滑非冲激噪声时可保留细节。B层：B1=zxy–zminB2=zxy–zmax

若B1>0且B2<0，输出zxy否则输出zmedA层：A1=zmed–zminA2=zmed–zmax

若A1>0且A2<0，转到B层否则增大窗口尺寸，若窗口尺寸≤Smax重复A层，否则输出zxy分为A层和B层：第28页，课件共72页，创作于2023年2月例八、自适应中值滤波器的作用被概率Pa=Pb=0.25的椒盐噪声污染的图像7x7中值滤波器的滤波效果Smax=7的自适应中值滤波器滤波效果噪声消除细节损失噪声消除细节保留第29页，课件共72页，创作于2023年2月四、频域滤波削减周期噪声削减或消除周期噪声的专用滤波器：带阻、带通和陷波滤波器。1、带阻滤波器。n阶巴特沃思带阻滤波器理想带阻滤波器第30页，课件共72页，创作于2023年2月几种带阻滤波器的透视图理想带阻滤波器一阶巴特沃思带阻滤波器高斯带阻滤波器高斯带阻滤波器第31页，课件共72页，创作于2023年2月例九、用带阻滤波器消除周期噪声被正弦噪声污染的图像被正弦噪声污染图像的傅里叶频谱巴特沃思带阻滤波器（白色为1）巴特沃思带阻滤波器滤波效果加性周期噪声近似二维正弦函数噪声成分对称亮点4阶第32页，课件共72页，创作于2023年2月2、带通滤波器。带通滤波器的传递函数是根据相应的带阻滤波器的传递函数按照下式：带通滤波器的操作与带阻滤波器相反。第33页，课件共72页，创作于2023年2月例十、带通滤波器提取噪声模型被正弦噪声污染的图像带通滤波器获得的噪声模式第34页，课件共72页，创作于2023年2月3、陷波滤波器。作用：阻止或通过预先定义的中心频率邻域内的频率。陷波滤波器必须以原点对称的形式出现，原因是傅里叶变换是对称的。特例陷波滤波器位于原点处时是以自身形式出现的。半径D0，中心在（u0，v0）且在（-u0，-v0）对称得理想陷波带阻滤波器函数为：第35页，课件共72页，创作于2023年2月n阶巴特沃思陷波带阻滤波器传递函数：高斯陷波带阻滤波器传递函数：陷波带通滤波器传递函数：第36页，课件共72页，创作于2023年2月几种陷波带阻滤波器的透视图理想陷波滤波器2阶巴特沃思陷波滤波器高斯陷波滤波器第37页，课件共72页，创作于2023年2月例十一、用陷波滤波器消除周期性噪声佛罗里达和墨西哥湾的卫星图(美国国家海洋和大气局提供)对应的傅里叶频谱叠加在频谱上的陷波带通滤波器滤波后图像的反傅里叶变换陷波带阻滤波器的效果水平扫描噪声扫描线消除图像未模糊第38页，课件共72页，创作于2023年2月4、最佳陷波滤波器。在多种噪声模式下，可在傅里叶频谱的每一尖峰下设置一个陷波带通滤波器H(u,v)，则干扰噪声为：则加性噪声下，未污染的图像为：噪声是真实噪声模式的近似，因此引入加权w(x,y)，得到图像的估计为：第39页，课件共72页，创作于2023年2月选取权函数，使得在某一邻域内，图像的方差最小。邻域内固定权函数第40页，课件共72页，创作于2023年2月滤波图像第41页，课件共72页，创作于2023年2月例十二、最佳陷波滤波器滤除噪声“水手6号”飞船拍摄的火星地形图（512x512美航空和宇宙航行局提供）周期干扰的傅里叶频谱（星形表明干扰模式包含不止一个正弦分量）干扰引起星形分量难以检测干扰精细第42页，课件共72页，创作于2023年2月水手6号火星地形图傅里叶频谱（美航空和宇宙航行局提供）第43页，课件共72页，创作于2023年2月噪声频谱和模式N（u，v）的傅里叶谱相应的噪声干扰模型η（u，v）（美航空和宇宙航行局提供）第44页，课件共72页，创作于2023年2月处理后的图像原始图像（美航空和宇宙航行局提供）第45页，课件共72页，创作于2023年2月五、线性、位置不变的退化由退化/复原模型有：若：则H具有加性和均匀性。若：则称H是线性、位置不变的系统函数。相应系统的冲激相应为：第46页，课件共72页，创作于2023年2月在有加性噪声条件下，线性退化模型表示为：在噪声随机，并与位置无关时，有：结论：具有加性噪声的线性空间不变退化系统，可表示为退化函数（点扩散函数）与图像的卷积。第47页，课件共72页，创作于2023年2月六、估计退化函数三种方法：试验法、观察法、数学建模法。1、观察法（仅有图像，对H一无所知）核心：收集图像信息：观察的子图像：构建的子图像2、试验估计法核心：用获得退化图像的装置，理论上可以得到准确的退化估计第48页，课件共72页，创作于2023年2月例十三、试验估计法的应用一个亮脉冲图像化的（退化的）冲激试验估计法的应用第49页，课件共72页，创作于2023年2月3、模型估计法模型化：大气湍流模型（Hufnagel和Stanley[1964]提出）数学建模：运动引起模糊图像的模型：水平运动的方程为：第50页，课件共72页，创作于2023年2月例十四、大气湍流模型可忽略的湍流剧烈湍流，k=0.0025第51页，课件共72页，创作于2023年2月中等湍流，k=0.001轻微湍流，k=0.00025第52页，课件共72页，创作于2023年2月例十五、运动引起的图像模糊原始图像a=b=0.1，T=1的模糊结果第53页，课件共72页，创作于2023年2月七、逆滤波逆滤波：用退化函数除退化图像的傅里叶频谱G(u，v)来计算原始图像的傅里叶估计，即：1、退化函数已知时，F(u，v)也是未知的，原因是N(u，v)是未知的。该式表明：2、若H(u，v)为零或者很小，N(u，v)/H(u，v)将决定原始图像的傅里叶估计。第54页，课件共72页，创作于2023年2月剧烈湍流，k=0.0025（480X480）（美航空航天局提供）例十六、逆滤波对图像进行复原退化函数为：k=0.0025第55页，课件共72页，创作于2023年2月全滤波复原的结果半径为40时截止H的结果退化值很小全逆滤波失效全逆滤波器失效，采用阶数为10的巴特沃思低通滤波器对进行截止。图像模糊第56页，课件共72页，创作于2023年2月半径为85时截止H的结果半径为70时截止H的结果图像清晰图像模糊第57页，课件共72页，创作于2023年2月八、最小均方误差滤波（维纳滤波）综合退化函数和噪声统计特征进行图像复原。目标是寻找一个未污染图像的估计，使得：前提条件：图像和噪声不相关；其中一个有零均值；估计的灰度级是退化图像灰度级的线性函数。第58页，课件共72页，创作于2023年2月此即维纳滤波。对白噪声，有：第59页，课件共72页，创作于2023年2月例十七、逆滤波和维纳滤波的比较剧烈湍流k=0.0025480X480半径受限逆滤波图全逆滤波的结果维纳滤波第60页，课件共72页，创作于2023年2月噪声幅度方差小1个量级噪声幅度方差小4个量级逆滤波的图像维纳滤波的图像例十八、维纳滤波和逆滤波的比较运动模糊及加性噪声污染的图像去模糊第61页，课件共72页，创作于2023年2月九、约束最小二乘方滤波器约束最小二乘方滤波器只需要噪声方差和均值。退化/复原模型有：设g(x，y)的尺寸为MXN，则g是由g(x，y)第一行元素构成向量的第一组N个元素，下一组N个元素用第二行构成，以此类推。结果向量为MNx1。同样方法构成f和η。H

矩阵为MNxMN，其元素由f(x，y)*h(x，y)的元素给出。第62页，课件共72页，创作于2023年2月上述方法的核心是H对噪声的敏感性问题。减小噪声敏感性的一种方法是采用平滑的最佳复原，为此引入最小准则函数C：上述方法的核心是H对噪声的敏感性问题。约束为：第63页，课件共72页，创作于2023年2月问题解为：P(u，v)是函数