510 likes | 604 Views
4.1 图像增强引言 4.2 空域图像增强 4.3 频域图像增强 4.3.1 频域滤波器 4.3.2 从频域规范产生空域模板. 为什么要用频率域进行图像增强呢?. 图像从空间域变换到频率域后, 1 、低频分量对应图像中灰度变化比较缓慢的区域; 2 、高频分量表征了图像中物体的边缘和随机噪声等信息。 3 、频率域上计算得到“简化”。. 频域增强原理. u. 边缘、噪音、 变化陡峭部分. 变化平缓部分. v. 4.3.1 频域滤波器. 频域滤波器 1 )低通滤波:保留低频分量,抑制高频分量 2 )高通滤波:保留高频分量,抑制低频分量.
E N D
4.1 图像增强引言 4.2 空域图像增强 4.3 频域图像增强 4.3.1 频域滤波器 4.3.2 从频域规范产生空域模板
为什么要用频率域进行图像增强呢? 图像从空间域变换到频率域后, 1、低频分量对应图像中灰度变化比较缓慢的区域; 2、高频分量表征了图像中物体的边缘和随机噪声等信息。 3、频率域上计算得到“简化”。
频域增强原理 u 边缘、噪音、 变化陡峭部分 变化平缓部分 v
4.3.1 频域滤波器 • 频域滤波器 1)低通滤波:保留低频分量,抑制高频分量 2)高通滤波:保留高频分量,抑制低频分量
频域滤波的理论基础 • 图像函数f(x,y)与线性位移不变算子h(x,y)的卷积结果是g(x,y),即g(x,y)= f(x,y)* h(x,y) • 由卷积定理可得:G(u,v)=F(u,v) H(u,v) • H(u,v)称为传递函数或滤波器函数。
频域滤波器的主要步骤: 1、对原始图像f(x,y)进行傅里叶变换得到F(u,v); 2、将F(u,v)与传递函数H(u,v)相乘得到G(u,v); 3、将G(u,v)进行傅里叶反变换得到增强图像g(x,y)。
4.3.1 频域滤波器:低通滤波 1)低通滤波 • 理想低通滤波器 • Butterworth低通滤波器 • 指数低通滤波器 • 梯形低通滤波器
4.3.1 频域滤波器:低通滤波 • 理想低通滤波器 • 理想低通滤波器的定义 • 理想低通滤波器截止频率的设计 • 理想低通滤波器的分析
4.3.1 频域滤波器:低通滤波 • 理想低通滤波器的定义 • 一个二维的理想低通滤波器(ILPF)的转换函数满足(是一个分段函数) 其中:D0 为截止频率 D(u,v)为距离函数 D(u,v)=(u2+v2)1/2
4.3.1 频域滤波器:低通滤波 • 理想低通滤波器的截面图 H(u,v) H(u,v)作为距离函数D(u,v)的函数的截面图 1 0 D0 D(u,v)
H(u,v) u 4.3.1 频域滤波器:低通滤波 • 理想低通滤波器的三维透视图 v H(u,v)作为u、v的函数的三维透视图
4.3.1 频域滤波器:低通滤波 • 理想低通滤波器的截止频率的设计 • 先求出总的信号能量PT: 其中: p(u,v) = |F(u,v)|2 = R2(u,v) + I2(u,v) 是能量模
4.3.1 频域滤波器:低通滤波 • 理想低通滤波器的截止频率的设计 • 如果将变换作中心平移,则一个以频域中心为原点,r为半径的圆就包含了百分之β的能量 其中: (u2 + v2) 1/2 <r
4.3.1 频域滤波器:低通滤波 • 理想低通滤波器的截止频率的设计 • 求出相应的D0 r = D0 =(u2 + v2)1/2
4.3.1 频域滤波器:低通滤波 • 理想低通滤波器的计算 • 对于给定的β • 用下面的公式计算出截止频率D0 r = D0 =(u2 + v2)1/2 3. 用频域理想低通滤波器H(u,v)与F(u,v)相乘
4.3.1 频域滤波器:低通滤波 • 理想低通滤波器的分析 • 例如:整个能量的90%被一个半径为8的小圆周包含,大部分尖锐的细节信息都存在于被去掉的10%的能量中 • 被钝化的图像被一种非常严重的振铃效果——理想低通滤波器的一种特性所影响
(b) (a) 理想低通滤波结果 半径分别为15,30,80,滤去的能量为5.4%、3.6%、2%。 (c)
4.3.1 频域滤波器:低通滤波 • 理想低通滤波器的分析 • 振铃效果——理想低通滤波器的一种特性
振铃效应 f(t)是一个输入信号,g(t)是f(t)通过一冲激响应为h(t)的线性系统后的输出。由于h(t)的两个负边带的存在,导致h(t)和f(t)卷积后,在输出图像上信号的两侧出现过冲现象,即称为振铃效应。
振铃效应 G(u,v)=H(u,v)F(u,v) g(x,y)=h(x,y)*f(x,y) (a)半径为5的脉冲图像 (b)相应的空间滤波器 (c)空域的5个脉冲 (d)滤波结果
4.3.1 频域滤波器:低通滤波 • 梯形低通滤波器 性能介于理想低通滤波器和Butterworth滤波器之间,有一定的模糊和振铃效应 调整D1,平滑图像,保持足够的清晰度。
= n H ( u , v ) exp{[ln( 1 2 )][ D ( u , v ) D ] } 0 = - n exp{ 0 . 347 [ D ( u , v ) D ] } 0 4.3.1 频域滤波器:低通滤波 • 指数低通滤波器 图像模糊程度比Butterworth严重,无振铃效应。
4.3.1 频域滤波器:低通滤波 • Butterworth低通滤波器 • Butterworth低通滤波器的定义 • Butterworth低通滤波器截止频率的设计 • Butterworth低通滤波器的分析
4.3.1 频域滤波器:低通滤波 • Butterworth低通滤波器的定义 • 一个截止频率在与原点距离为D0的n阶Butterworth低通滤波器(BLPF)的变换函数如下:
4.3.1 频域滤波器:低通滤波 • Butterworth低通滤波器的截面图 H(u,v) 由于H(u,v)在通过频率和滤去频率之间没有明显的不连续,更无阶跃或突变,而是存在一个平滑的过渡带,因而平滑图像不存在振铃现象,较理想低通滤波器好。 1 0.5 D(u,v)/D0 0 1 2 3 4
4.3.1 频域滤波器:低通滤波 • Butterworth滤波器截止频率的设计 • 变换函数中不存在一个不连续点作为一个通过的和被滤波掉的截止频率的明显划分 • 通常把H(u,v)开始小于其最大值的一定比例的点当作其截止频率点
4.3.1 频域滤波器:低通滤波 有两种选择: 1、选择H(u,v) = 0.5时的 D(u,v) 为截止频率D0,其传递函数为: 2、选择H(u,v) = 1/2 时的 D(u,v)为截止频率 D0 ,其传递函数为:
ELPF 特性曲线 BLPF 特性曲线
2阶BLPF滤波的结果 (a)原图像(b)半径15 (b)半径30 (d)半径80 2阶ELPF滤波的结果 (a)原图像(b)半径15 (b)半径30 (d)半径80
4.3.1 频域滤波器:低通滤波 • Butterworth低通滤波器的分析 • 在任何经BLPF处理过的图像中都没有明显的振铃效果,这是滤波器在低频和高频之间的平滑过渡的结果 • 低通滤波是一个以牺牲图像清晰度为代价来减少干扰效果的修饰过程
4.3.1 频域滤波器:低通滤波 • Butterworth低通滤波器的分析 • BLPF处理过的图像中都没有振铃效果
4.3.1 频域滤波器:低通滤波 • Butterworth低通滤波器的分析
4.3.1 频域滤波器:高通滤波 2)高通滤波 • 频域高通滤波的基本思想 • 理想高通滤波器 • Butterworth高通滤波器
4.3.1 频域滤波器:高通滤波 • 频域高通滤波的基本思想 • G(u,v)=F(u,v)H(u,v) • F(u,v)是需要锐化图像的傅立叶变换形式。 • 目标是选取一个滤波器变换函数H(u,v),通过它减少F(u,v)的低频部分来得到G(u,v)。 • 运用傅立叶逆变换得到锐化后的图像。
4.3.1 频域滤波器:高通滤波 • 理想高通滤波器 • 理想高通滤波器的定义 • 理想高通滤波器截止频率的设计 • 理想高通滤波器的分析
4.3.1 频域滤波器:高通滤波 • 理想高通滤波器的定义 • 一个二维的理想高通滤波器(ILPF)的转换函数满足(是一个分段函数) 其中:D0 为截止频率 D(u,v)为距离函数 D(u,v)=(u2+v2)1/2
4.3.1 频域滤波器:高通滤波 • 理想高通滤波器的截面图 H(u,v) H(u,v)作为距离函数D(u,v)的函数的截面图 1 0 D0 D(u,v)
H(u,v) u v 4.3.1 频域滤波器:高通滤波 • 理想高通滤波器的三维透视图 H(u,v)作为u、v的函数的三维透视图
4.3.1 频域滤波器:高通滤波 • 梯形高通滤波器
4.3.1 频域滤波器:高通滤波 • Butterworth高通滤波器 • Butterworth高通滤波器的定义 • Butterworth高通滤波器截止频率设计 • Butterworth高通滤波器的分析
4.3.1 频域滤波器:高通滤波 • Butterworth高通滤波器的定义 • 一个截止频率在与原点距离为D0的n阶Butterworth高通滤波器(BHPF)的变换函数如下:
4.3.1 频域滤波器:高通滤波 • Butterworth高通滤波器的截面图 H(u,v) H(u,v)作为D(u,v)/D0的函数的截面图 1 0.5 D(u,v)/D0 0 1 2 3
4.3.1 频域滤波器:高通滤波 • Butterworth高通滤波器截止频率设计 • 变换函数中不存在一个不连续点作为一个通过的和被滤波掉的截止频率的明显划分 • 通常把H(u,v)开始小于其最大值的一定比例的点当作其截止频率点
4.3.1 频域滤波器:高通滤波 有两种选择: 1、选择H(u,v) = 0.5时的 D(u,v) 为截止频率D0 ,其 传递函数为: 2、选择H(u,v) = 1/2 时的 D(u,v)为截止频率 D0 ,传递函数为:
4.3.1 频域滤波器:高通滤波 • Butterworth高通滤波器的分析 • 问题:低频成分被严重地消弱了,使图像失去层次 • 改进措施: • 加一个常数到变换函数 H(u,v) + A 这种方法被称为高频强调 • 为了解决变暗的趋势,在变换结果图像上再进行一次直方图均衡化。这种方法被称为后滤波处理
BHPF IHPF EHPF