添加椒盐噪声并用低通巴特沃斯滤波器设计步骤去噪matlab

数字图像处理实验1-9点击下列链接囿源码和链接:

      (1)噪声: 图像噪声是指存在于图像数据中的不必要的或多余的干扰信息

1.高斯噪声:高斯噪声是指它的概率密度函数服从高斯分布(即正态分布)的一类噪声。常见的高斯噪声包括起伏噪声、宇宙噪声、热噪声和散粒噪声等等除常用抑制噪声的方法外,对高斯噪声的抑制方法常常采用数理统计方法高斯噪声最常见最普通的噪声。正态分布(高斯分布)的公式(高中知识不详解了):

2.椒盐噪聲:椒盐噪声(salt-and-pepper noise)又称脉冲噪声它随机改变一些像素值,在二值图像上表现为使一些像素点变白一些像素点变黑(可以理解为随机概率嘚像素点变色)。

3.泊松噪声:泊松噪声就是符合泊松分布的噪声模型,泊松分布适合于描述单位时间内随机事件发生的次数的概率分布

泊松分布公式(高中知识不解释了):

(2)滤波:图像滤波,即在尽量保留图像细节特征的条件下对目标图像的噪声进行抑制是图像预处悝中不可缺少的操作,其处理效果的好坏将直接影响到后续图像处理和分析的有效性和可靠性

针对不同噪声,也有不同的滤波方式:

1.高斯滤波-------顾名思义解决高斯噪声问题

高斯滤波是一种线性平滑滤波,适用于消除高斯噪声广泛应用于图像处理的减噪过程。通俗的讲高斯滤波就是对整幅图像进行加权平均的过程,每一个像素点的值都由其本身和邻域内的其他像素值经过加权平均后得到。方法主要有兩种:

一是用离散化窗口滑窗卷积:就是选定一定大小的滑窗遍历整个图像。比如我取像素9x9的矩阵滑窗遍历整个图像进行卷积,一般峩们都这样处理高斯噪声

只有当离散化的窗口非常大,用滑窗计算量非常大(即使用可分离滤波器的实现)的情况下可能会考虑基于傅里叶变化的实现方法。

2.中值滤波-------解决椒盐噪声效果比较好

中值滤波法是一种非线性平滑技术,它将每一像素点的灰度值设置为该点某鄰域窗口内的所有像素点灰度值的中值中值滤波的基本原理是把数字图像或数字序列中一点的值用该点的一个邻域中各点值的中值代替,让周围的像素值接近的真实值从而消除孤立的噪声点。一般用于处理椒盐噪声

3.均值滤波-------解决泊松噪声效果比较好。

均值滤波是典型嘚线性滤波算法它是指在图像上对目标像素给一个模板,该模板包括了其周围的临近像素(以目标像素为中心的周围8个像素构成一个濾波模板,即去掉目标像素本身)再用模板中的全体像素的平均值来代替原来像素值。

回答题目的问题:平滑处理就是用滤波器去除噪聲!



3、平滑处理(图像滤波)

(1)对添加了高斯噪声的图像进行高斯滤波:

(2)使用中值滤波处理椒盐噪声

3)使用均值滤波处理泊松噪聲

}

    图像去噪是数字图像处理中的重偠环节和步骤去噪效果的好坏直接影响到后续的图像处理工作如图像分割、边缘检测等。图像信号在产生、传输过程中都可能会受到噪聲的污染一般数字图像系统中的常见噪声主要有:高斯噪声(主要由阻性元器件内部产生)、椒盐噪声(主要是图像切割引起的黑图像仩的白点噪声或光电转换过程中产生的泊松噪声)等;

目前比较经典的图像去噪算法主要有以下三种:

均值滤波算法:也称线性滤波,主偠思想为邻域平均法即用几个像素灰度的平均值来代替每个像素的灰度。有效抑制加性噪声但容易引起图像模糊,可以对其进行改进主要避开对景物边缘的平滑处理。

中值滤波:基于排序统计理论的一种能有效抑制噪声的非线性平滑滤波信号处理技术中值滤波的特點即是首先确定一个以某个像素为中心点的邻域,一般为方形邻域也可以为圆形、十字形等等,然后将邻域中各像素的灰度值排序取其中间值作为中心像素灰度的新值,这里领域被称为窗口当窗口移动时,利用中值滤波可以对图像进行平滑处理其算法简单,时间复雜度低但其对点、线和尖顶多的图像不宜采用中值滤波。很容易自适应化

Wiener维纳滤波:使原始图像和其恢复图像之间的均方误差最小的複原方法,是一种自适应滤波器根据局部方差来调整滤波器效果。对于去除高斯噪声效果明显

实验一:均值滤波对高斯噪声的效果

 

% 加叺均值为0,方差为0.005的高斯噪声

' 加入高斯噪声之后的图像 ' ); % 采用MATLAB中的函数filter2对受噪声干扰的图像进行均值滤波

fspecial函数用于创建预定义的滤波算子其语法格式为:
type='average' ,为均值滤波参数为n,代表模版尺寸用向量表示,默认值为[3,3]
type= 'gaussian' ,为高斯低通滤波器参数有两个,n表示模版尺寸默認值为[3,3]sigma表示滤波器的标准差单位为像素,默认值为

type= 'laplacian' 为拉普拉斯算子,参数为alpha用于控制拉普拉斯算子的形状,取值范围为[0,1]默认值為0.2
type= 'log' 为拉普拉斯高斯算子,参数有两个n表示模版尺寸,默认值为[3,3]sigma为滤波器的标准差,单位为像素默认值为 0.5
type= 'unsharp' ,为对比度增强滤波器参数alpha用于控制滤波器的形状,范围为[0,1]默认值为0.2

据我目测使用均值滤波去噪(高斯噪声)效果选用的邻域半径越大效果越好,当然其代价也会更大另外确切的去噪效果的好坏还需要用SNR等数据来度量。

实验二:二维自适应维纳滤波对高斯噪声的滤除效果
 

% 加入均值为0方差为0.005的高斯噪声

PS 维纳滤波的两个函数 wiener2deconvwnr都能够完成维纳滤波的功能,   deconvwnr强调图象复原方面 wiener2强调图象空间域锐化的作用,其中

实验三:對加入椒盐噪声的图像分别作均值、中值和维纳滤波
 

; % 均值滤波模板尺寸为3

' 均值滤波后的图像 ' ); ' 中值滤波后的图像 ' ); ' 维纳滤波后的图像

通过图3我們也可得出结论即中值滤波对于去除椒盐噪声效果最好,而维纳滤波去除效果则较差中值滤波对于去除椒盐噪声效果明显,是因为椒鹽噪声只在画面上的部分点随机出现而中值滤波根据数据排序,将未被污染的点代替噪声点的值的概率较大所以抑制效果好。对点、線和尖顶较多的图像不宜采用中值滤波因为一些细节点可能被当成噪声点。

实验四: 分别使用二维统计滤波对椒盐噪声和高斯噪声进行滤波
 
' 加椒盐噪声后的图像 ' ); ' 椒盐噪声滤波后的图像 ' ); ' 加高斯噪声后的图像
PSMATLAB小波分析工具箱提供的用于图像去噪的函数有wrcoef2wpdencmp其语法格式分别为:

其中,X=wrcoef2(‘type’,C,S,’wname’)返回基于小波分解结构[CS]的小波重构图像X。参数“type”等于a表示重构近似系数;等于h表示重构水平细节系数;等于v表示重构垂直细节系数等于d表示重构对角细节系数。

}

一)空间域滤波与频率域滤波

空間域滤波是指在图像空间中借助模板对图像领域进行操作处理图像每一个像素值。主要分为线性滤波和非线性滤波两类根据功能可分為平滑滤波器和锐化滤波器。平滑可通过低通来实现平滑的目的有两类,一是模糊目的是在提取较大的目标前去除太小的细节或将目標内的小尖端连接起来;二是去噪。锐化则可用高通滤波来实现锐化的目的是为了增强被模糊的细节。

      在matlab中实现空间域滤波有很多类型,如均值、中值、索贝尔、高斯、拉普拉斯、高斯-拉普拉斯等但各有差异。下面是用matlab实现的代码:

}

我要回帖

更多关于 低通巴特沃斯滤波器设计步骤 的文章

更多推荐

版权声明:文章内容来源于网络,版权归原作者所有,如有侵权请点击这里与我们联系,我们将及时删除。

点击添加站长微信