资源描述：
《数字图像处理第六次作业.doc》由会员上传分享，免费在线阅读，更多相关内容在行业资料-天天文库。
1、...数字图像处理第六次作业.......摘要本次报告主要记录第六次作业中的各项任务完成情况。本次作业以Matlab2013为平台，结合matlab函数编程实现对所给图像文件的相关处理：1.在测试图像上产生高斯噪声lena图-需能指定均值和方差；并用滤波器（自选）恢复图像；2.实现下边要求；(a)实现模糊滤波器如方程Eq.(5.6-11)；(b)模糊lena图像：45度方向，T=1；(c)在模糊的lena图像中增加高斯噪声，均值=0，方差=10pixels以产生退化图像；(d)分别利用方程Eq.(5.8-6)和(5.9-4)，恢复图像。以上任务完成后均得到
2、了预期的结果。1.在测试图像上产生高斯噪声lena图-需能指定均值和方差；并用滤波器（自选）恢复图像；（1）实验原理与方法图像复原处理是建立在图像退化的数学模型基础上的，这个退化数学模型能够反映图像退化的原因。图像的退化过程可以理解为施加于原图像上的运算和噪声两者联合作用的结果，图像退化模型如图1所示，可以表示为：图1图像退化模型高斯噪声是指它的概率密度函数服从高斯分布（即正态分布）的一类噪声。一个高斯随机变量z的PDF可表示为：其中z代表灰度，u是z的均值，是z的标准差。高斯噪声的灰度值多集中在均值附近。本文采用5×5模板的中值滤波器和高斯滤波器（σ=1
3、.5）作业四中已经介绍过，中值滤波器是使用一个像素邻域中灰度级的中值来替代该像素值，即，。高斯滤波是一种根据高斯函数的形状来选择模板权值的线性平滑滤波方法，具体操作是：用一个模板（或称卷积）扫描图像中的每一个像素，用模板确定的邻域像素的加权平均灰度值去替代模板中心像素点的值。利用matlab中imnoise函数加入高斯噪声：g=imnoise(f,type,parameters)调用格式：g=imnoise(I,type)g=imnoise(I,type,parameters)参数Type对应的噪声类型如下：'gaussian'高斯白噪声'localvar
4、'0均值白噪声'poisson'泊松噪声.......'salt&pepper'盐椒噪声'speckle'乘性噪声滤波程序同作业四（2）处理结果..............（3）结果分析通过imnoise函数产生了被均值和方差可选的高斯噪声污染的图像。当高斯噪声均值不变为0时，随着方差增加，图像噪声越严重；当高斯噪声方差不变时，均值会影响到整个图像的灰度值，使整个图像变亮。与理论上均值和方差对图像的影响一致。分别使用高斯滤波器和中值滤波器对加噪图像进行恢复。两种方法在一定程度上都可以降低噪声。高斯滤波器降低噪声的同时保存的图像细节更丰富，亮度比原噪声图像和
5、中值滤波后图像暗更接近原始图像，中值滤波后图像亮度基本与原噪声图像相同。2．实现下边要求：(a)实现模糊滤波器如方程Eq.(5.6-11).(b)模糊lena图像：45度方向，T=1；(c)在模糊的lena图像中增加高斯噪声，均值=0，方差=10pixels以产生退化图像；(d)分别利用方程Eq.(5.8-6)和(5.9-4)，恢复图像；（5.6-11）利用上式模糊lena图像：45度方向，T=1，即使a=b=0.1，T=1。对原始图像的图像矩阵做傅里叶变换并移至图像中心得到频域矩阵F，使H与F相乘后反傅里叶变换到空域得到变换后图像。维纳滤波综合了退化函数
6、和噪声统计特性两个方面进行复原处理，其目标是寻找一个滤波器，使得复原后图像与原始图像的均方误差最小：.......因此维纳滤波器又称为最小均方误差滤波器，在频率中用下式表达：其中G(u,v)是退化图像的傅里叶变换，H(u,v)是退化函数。为噪声功率谱，为未退化图像的功率谱。式5.8-6为属于维纳滤波式5.9-4为其中，γ是一个参数，必须对它进行调整以满足的条件，属于约束最小二层方滤波利用以上两式恢复图像的流程与对模糊原始图像的流程相似：对原始图像的图像矩阵做傅里叶变换并移至图像中心得到频域矩阵F，通过H得到，使与F相乘后反傅里叶变换到空域得到变换后图像。在
7、实现5.9-4时借助matlab工具包以得到更好的效果。用fspecial和imnoise函数得到45度方向，T=1的模糊lena图像，并在此图像上实现维纳滤波和约束最小二乘方滤波。1）imfilter功能：对任意类型数组或多维图像进行滤波。用法：B=imfilter(A,H)B=imfilter(A,H,option1,option2,...)或写做g=imfilter(f,w,filtering_mode,boundary_options,size_options)其中，f为输入图像，w为滤波掩模，g为滤波后图像。filtering_mode用于指定在
8、滤波过程中是使用“相关”还是“卷积”。boundary_optio
