资源描述：
《图像恢复(数字图像处理实验报告)》由会员上传分享，免费在线阅读，更多相关内容在教育资源-天天文库。
1、数字图像处理作业——图像恢复摘要数字图像恢复是数字图像处理的一个基本的和重要的课题，它是后期图像处理（分析和理解）的前提。图像在摄取、传输、储存的过程中不可避免地引起图像质量的下降（图像退化），图像恢复就是试图利用退化过程的先验知识使已退化的图像恢复本来面貌，即根据退化的原因，分析引起退化的环境因素，建立相应的数学模型，并沿着使图像降质的逆过程恢复图像。本文首先对测试图像进行模糊及加噪处理，然后用不同的图像恢复方法，如维纳滤波恢复、约束最小二乘滤波进行图像恢复，并比较它们的处理效果。发现维纳滤波较约束最小二乘法滤波效果要好
2、，这是因为前者利用了原图像的统计信息，采用了真实的PSF函数来恢复。无论何种算法，它们都要依据获取的相关信息才能有效地实施，算法利用的信息越多，信息的准确性越高，复原图像的质量也就越高。实验原理：图像复原处理是建立在图像退化的数学模型基础上的，这个退化数学模型能够反映图像退化的原因。图像的退化过程可以理解为施加于原图像上的运算和噪声两者联合作用的结果，图像退化模型如图1所示，可以表示为：g(x,y)=H[f(x,y)]+n(x,y)=f(x,y)*h(x,y)+n(x,y)（1）图1图像退化模型(1)在测试图像上产生高斯噪
3、声lena图-需能指定均值和方差；并用滤波器（自选）恢复图像；实验原理：噪声是最常见的退化因素之一，也是图像恢复中重点研究的内容，图像中的噪声可定义为图像中不希望有的部分。噪声是一种随机过程，它的波形和瞬时振幅以及相位都随时间无规则变化，因此无法精确测量，所以不能当做具体的处理对象，而只能用概率统计的理论和方法进行分析和处理。本文中研究高斯噪声对图像的影响及其去噪过程。①高斯噪声的产生：所谓高斯噪声是指它的概率密度函数服从高斯分布（即正态分布）的一类噪声。一个高斯随机变量z的PDF可表示为：P（z）=（2）其中z代表灰度，
4、u是z的均值，是z的标准差。高斯噪声的灰度值多集中在均值附近。图2高斯函数可以通过不同的算法用matlab来产生高斯噪声。②高斯噪声对信号的影响噪声影响图像处理的输入、采集、处理的各个环节以及输出结果的全过程，在图像中加高斯噪声通常会使图像变得模糊并且会出现细小的斑点，使图像变得不清晰。11③去除高斯噪声的一些方法去除高斯噪声的方法有直方图变换，低通滤波，高通滤波，逆滤波，维纳滤波，中值滤波等。本文应用高斯平滑滤波进行去噪处理。处理结果如下图：维纳滤波对高斯白噪声的图像滤波效果较好，具有比较好的选择性，可以更好地保存图像的
5、边缘和高频细节信息。所以，维纳滤波在大多数情况下都可以获得满意的结果，尤其对含有高斯噪声的图像。(1)推导维纳滤波器并实现下边要求；实验原理：维纳滤波综合了退化函数和噪声统计特性两个方面进行复原处理，其目标是寻找一个滤波器，使得复原后图像与原始图像的均方误差最小：11因此维纳滤波器又称为最小均方误差滤波器。在频率中用下式表达：其中，是退化图像的傅立叶变换，是退化函数。，是的复共轭。，为噪声的功率谱。，为未退化图像的功率谱。维纳滤波器的推导：合理假设要估计的信号f(x,y)为0均值平稳随机过程，噪声为0均值平稳随机过程。g(
6、x,y)=h(x,y)*f(x,y)+(x,y)根据已退化图像g(x,y),利用线性估计器来估计原始图像。先证：g(x,y)=f(x,y)+(x,y)，最小化均方误差：=,e(x,y)=f(x,y)−fˆ(x,y),fˆ(x,y)=w(x,y)*g(x,y)最小化均方误差可以由正交原则求解，即最优求解的误差与观测的信号值不相关。得：E[e(m,n)g*(x,y)]=E[f(m,n)−fˆ(m,n)]g*(x,y)]=0进一步分解得：E[f(m,n)g*(x,y)]=E[fˆ(m,n)g*(x,y)]==E{w(k1,k2)
7、[g(m−k1,n−k2)g*(x,y)]}=E{w(k1,k2)(Rg(m−k1-x,n−k2-y)}=w(m-x,n-y)*Rg(m-x,n-y)11记：Rfg(m-x,n-y)=E[f(m,n)g*(x,y)]因此，Rfg(m,n)=w(m,n)*Rg(m,n)上式作Flouier变换得：Pfg(w1,w2)=W(w1,w2)*Pg(w1,w2)即：W(w1,w2)=Pfg(w1,w2)/Pg(w1,w2).（3）又有：Rfg(m-x,n-y)=E[f(m,n)g*(x,y)]=E[f(m,n).[h(k1,k2)f
8、(x-k1,y-k2)+(x,y)]*=h*(k1,k2)Rf(m-x+k1,n-y+k2)=h*(-(m-x),-(n-y))*Rf(m-x,n-y)因此，Rfg(m,n)=h*(-m,-n)*Rf(m,n)对上式作Flouier变换得：Pfg(w1,w2)=H*(w1,w2)Pf(w1,w2)（4）
