目录 l 傅立叶变换 l 低通高通滤波器 l 带通和带阻滤波器 l 同态滤波器

Size: px

Start display at page:

Download "目录 l 傅立叶变换 l 低通高通滤波器 l 带通和带阻滤波器 l 同态滤波器"

宽贺
5 years ago
Views:

1 图像处理与分析频域内滤波授课教师 : 孙剑 jiansun@mail.xjtu.edu.cn 西安交通大学数学与统计学院

2 目录 l 傅立叶变换 l 低通高通滤波器 l 带通和带阻滤波器 l 同态滤波器

3 傅里叶变换 2-D 傅里叶变换傅里叶变换定理快速傅里叶变换

4 2-D 傅里叶变换 1-D 正变换 N 1 F {f ( x )}= F ( u)= 1 N f ( x )exp[ j2 πux /N ] x= 0 对 1 个连续函数 f (x) 等间隔采样 f ( x) x

5 2-D 傅里叶变换 1-D 反变换 N 1 F 1 {F ( u) }= f ( x)= F ( u) exp[j2πux/n ] u= 0 变换表达 F ( u)= R( u)+ ji ( u)= F ( u) exp[j φ( u)] 频谱 ( 幅度 ) F (u) = [R 2 (u)+ I 2 ( u)] 1/2 相位角 φ(u)= arctan[ I (u) /R (u)]

6 2-D 傅里叶变换 2D 傅里叶变换公式 F ( u, v )= 1 N x= 0 f ( x, y )= 1 N u= 0 N 1 N 1 y= 0 N 1 N 1 v= 0 f ( x, y) exp[ j2 π( ux + vy )/N ] F ( u, v) exp[j2π( ux +vy )/N ] 频谱 ( 幅度 ) 相位角功率谱 F (u, v ) = [R 2 ( u, v)+ I 2 (u, v )] 1/2 φ(u, v )= arctan[ I ( u, v )/R(u, v)] P(u, v)= F ( u, v ) 2 = R 2 ( u, v)+ I 2 ( u, v)

7 傅里叶变换定理分离性质 F ( u, v ) = N 1 F ( x, v)= N[ 1 N y= 0 f ( x, y) exp[ j2πvy/n ]] N 1 N 1 F ( x, v ) exp[ j2πux /N ] x= 0 1 次 2-D 通过 2 次 1-D 计算 O(N 4 ) 减为 O(N 2 ) (N-1) (0,0) Y (0,0) V 列变换行变换 f (x, y ) F(x, v ) F(u, v) 乘以 N X X - (N 1) (N-1) (N-1) (N-1) (0,0) V (N-1) U

8 傅里叶变换定理 1 平移定理 f ( x, y ) F ( u, v) f ( x, y )exp[ j2π( cx+ dy)/n ] F (u c, v d ) F{ f(x,y)exp[ j2π(cx +dy) / N] } = 1 N = 1 N N 1 N 1 x=0 y=0 N 1 N 1 x=0 y=0 = F[(u c),(v d)] f(x,y)exp[ j2π(cx +dy) / N]exp[ j2π(ux +vy) / N] f(x,y)exp{ j2π[(u c)x +(v d)y] / N}

9 傅里叶变换定理 4 剪切定理 ( 水平方向 ) 纯剪切 x'= x+ by y '= y ( 垂直方向 ) 纯剪切 f ( x+by, y) F (u, v bu) x'= x y '= dx+ y f ( x,dx+ y) F (u du, v )

10 傅里叶变换定理 5 组合剪切定理平移 + 旋转 + 尺度 f ( x+by, dx+ y ) 1 1 bd F ( u dv 1 bd, bu+ v ) 1 bd x ' [ = 1 d b 1] x 水平剪切 [ 1 b 1] 垂直剪切 [ 1 0 1] 0 d

11 傅里叶变换定理 6 卷积定理 f ( x) g( x) F (u)g (u) 2-D f ( x)g( x) F (u) G (u) f ( x, y ) g( x, y )= f ( p,q) g( x p, y q)dpdq f ( x, y ) g( x, y ) F (u, v)g (u, v) f ( x, y )g( x, y ) F (u, v) G (u, v)

12 傅里叶变换定理 7 相关定理 1D 相关运算 f ( x) g( x)= f ( z)g( x+ z) dz f ( x, y ) g( x, y )= f ( p, q)g( x+ p, y+q) dpdq f ( x, y ) g( x, y ) F (u, v)g (u, v) f ( x, y )g( x, y ) F (u, v) G (u, v)

13 快速傅里叶变换 (fft) 直接进行一个 N N 的 2-D 傅里叶变换需要 N 4 次复数乘法运算和 N 2 (N 2 1) 次复数加法运算 N 1 F {f ( x )}= F ( u)= 1 N x= 0 f ( x )exp[ j2πux /N ] u= 0, 1,, N 1 1-D: 复数乘法和加法的次数都正比于 N 2 快速傅里叶变换 (FFT): 将复数乘法和加法的次数减少为正比于 N log 2 N

14 傅里叶变换例子 im = double(imread('lena.jpg')); imfft = fft2(im(:,:,1)); imfftshift = fftshift(imfft); R = real(imfftshift); I = imag(imfftshift); A1 = sqrt(r.^2 + I ^ 2); A1=(A1-min(min(A1)))/(max(max(A1))- min(min(a1)))*225; figure,imshow(a1)

15 离散余弦变换一种可分离正交对称的变换 1-D 离散余弦变换 (DCT) C( u) N = ( ) a u x= 1 0 f (2x + 1) uπ ( x)cos 2N u = 0,1,, N 1 N 1 ( 2x+ 1) uπ f ( x)= a( u)c ( u) cos[ ] x= 0, 1,, N 1 u= 0 2N a( u)= { 1/N 当 u= 0 2/N 当 u= 1, 2,, N 1

16 离散余弦变换 2-D 离散余弦变换 (DCT) N 1 N 1 ( 2x+ 1) uπ C ( u, v )= a( u) a( v ) f ( x, y) cos[ ] x= 0 y= 0 2N cos [ ( 2y +1) vπ ] 2N N 1 N 1 ( 2x+ 1) uπ f ( x, y )= a( u) a( v )C ( u, v) cos[ ] u= 0 v= 0 2N cos [ ( 2y +1) vπ ] 2N + 讨论可分离性和对称性 h( x, y, u, v )= h 1 ( x, u) h 2 ( y, v) h( x, y, u, v )= h 1 ( x, u) h 1 ( y, v )

17 离散余弦变换例子 im = imread('lena.jpg'); gray = rgb2gray(im); D = dct2(gray); figure(2);imshow(log(abs(d)),[]);

18 离散余弦变换例子变换域内信息压缩 D(abs(D) < 10) = 0 逆变换 D(abs(D)<10)=0; I=idct2(D)/255; figure(3);imshow(i)

19 频域图像增强频域增强原理低通滤波高通滤波带通和带阻滤波同态滤波频域技术与空域技术

20 频域增强原理卷积理论是频域技术的基础设函数 f (x, y) 与线性位不变算子 h(x, y) 的卷积结果是 g(x, y), 即 g(x, y) = h(x, y) * f (x, y), 那么根据卷积定理在频域有 : G (u, v)= H (u, v) F ( u, v) 其中 G(u, v),h(u, v),f(u, v) 分别是 g(x, y),h(x, y),f (x, y) 的傅里叶变换用线性系统理论的话来说,H(u, v) 是转移函数

21 频域增强原理在具体增强应用中,f (x, y) 是给定的 ( 所以 F(u, v) 可利用变换得到 ), 需要确定的是 H(u, v), 这样具有所需特性的 g(x, y) 就可由算出 G(u, v) 而得到 : g( x, y)= T 1 {E H [T [f ( x, y)]]} 步骤 : (1) 转换到频域 (2) 在频域增强 (3) 转换回空域

22 频域增强原理卷积定理增强图 G (u, v)= H (u, v) F ( u, v) g( x, y)= T 1 [H (u, v) F (u, v )] 步骤 (1) 计算图象的变换 (2) 在频域滤波 (3) 反变换回图象空间频域滤波低通, 高通, 带通 / 带阻, 同态

23 低通滤波低通滤波器图象中的边缘和噪声都对应图象傅里叶变换中的高频部分, 所以如要在频域中消弱其影响就要设法减弱这部分频率的分量 ; 根据频域增强技术的原理, 需要选择一个合适的 H(u, v) 以得到消弱 F(u, v) 高频分量的 G(u, v); 以下讨论对 F(u, v) 的实部和虚部影响完全相同的滤波转移函数具有这种特性的滤波器称为零相移滤波器

24 低通滤波 1 理想低通滤波器理想是指小于 D 0 的频率可以完全不受影响地通过滤波器, 而大于 D 0 的频率则完全通不过

25 低通滤波 1 理想低通滤波器 H(u, v): 转移 / 滤波函数 H(u,v) = 1 if D(u,v) D 0 0 if D(u,v) > D 0 D 0 : 截断频率 ( 非负整数 ) D(u, v) 是从点 (u, v) 到频率平面原点的距离 D(u, v) = (u 2 +v 2 ) 1/2

26 低通滤波 2 理想低通滤波器的模糊理想低通滤波产生振铃现象

27 低通滤波 2 理想低通滤波器的模糊理想低通滤波所产生的振铃现象在 2-D 图象上表现为一系列同心圆环圆环半径反比于截断频率理想低通滤波产生模糊效应 B: 能量百分比,R: 圆周半径,P(u, v): 功率谱 B= 100 [ u R v R N 1 N 1 P(u, v)/ u= 0 v= 0 P(u, v)]

28 低通滤波 3 巴特沃斯低通滤波器物理上可实现 ( 理想低通滤波器在数学上定义得很清楚, 在计算机模拟中也可实现, 但在截断频率处直上直下的理想低通滤波器是不能用实际的电子器件实现的 ) 减少振铃效应, 高低频率间的过渡比较光滑阶为 n H ( u, v)= 1 1+ [D ( u, v)/d 0 ] 2n

29 低通滤波 3 巴特沃斯低通滤波器截断频率使 H 最大值降到某个百分比的频率在 D(u, v) = D 0 时 H(u, v) = 1/2 H ( u, v)= 1 1+ [D ( u, v)/d 0 ] 2n

30 低通滤波 3 巴特沃斯低通滤波器图象由于量化不足产生虚假轮廓时常可用低通滤波进行平滑以改进图象质量效果比较 ( 相同截断频率 ): 教材图理想低通滤波器阶数为 1 的巴特沃斯低通滤波器

31 低通滤波 4 其他低通滤波器梯形指数

32 高通滤波 1 理想高通滤波器形状与低通滤波器的形状正好相反 H(u,v) = 0 if D(u,v) D 0 1 if D(u,v) > D 0

33 高通滤波 2 巴特沃斯高通滤波器形状与巴特沃斯低通滤波器的形状正好相反截断频率使 H 值上升到最大值某个百分比的频率 H(u, v) = 1/2 H(u, v) = 1/2 1/2 H ( u, v)= 1 1+ [D 0 /D ( u, v)] 2n

34 高通滤波 3 高频增强滤波器傅里叶变换 :G(u, v) = H(u, v)f(u,v) 高频增强转移函数 :H e (u, v) = k H(u, v) + c 高频增强输出图的傅里叶变换 : G e (u, v) = k G(u, v) + c F(u, v) 反变换回去 : g e (x, y) = k g(x, y) + c f (x, y)

35 高通滤波 4 高频提升滤波器用原始图减去低通图得到高通滤波器的效果把原始图乘以一个放大系数 A 再减去低通图就可构成高频提升 (high-boost) 滤波器高通滤波器 :A = 1 高频增强滤波器 :?

36 带通和带阻滤波带阻滤波器阻止一定频率范围 ( 允许其它频率范围 ) 1 H ( u,v) W D ( u,v) 0 D 0 章毓晋 (TH-EE-IE)

37 带通和带阻滤波带阻滤波器傅里叶变换的对称性 > 两两工作 H(u,v) = 0 if D 1 (u,v) D 0 or D 2 (u,v) D 0 1 else D 1 ( u, v )= [( u u 0 ) 2 +( v v 0 ) 2 ] 1/2 D 2 ( u, v )= [( u+u 0 ) 2 +(v+ v 0 ) 2 ] 1/2

38 带通和带阻滤波放射对称的带阻滤波器 H(u,v) = 1 D(u,v) < D 0 W / 2 0 D 0 W / 2 D(u,v) D 0 +W / 2 1 D(u,v) > D 0 +W / 2 H ( u, v)= 1 D ( u, v )W 2]2n 1+[ D 2 ( u, v) D 0

39 带通和带阻滤波带通滤波器与带阻滤波器互补允许一定频率范围 ( 阻止其它频率范围 ) H P (u, v )= [H R (u, v) 1]= 1 H R (u, v) 1 H ( u,v) W D ( u,v) 0 D 0

40 带通和带阻滤波放射对称的带通滤波器 H(u,v) = 0 if D(u,v) < D 0 W / 2 1 if D 0 W / 2 D(u,v) D 0 +W / 2 0 if D(u,v) > D 0 +W / 2 H ( u, v)= [ D ( u, v) D 0 D ( u, v) W ]2n 1

41 频域技术与空域技术空间滤波器的工作原理可借助频域进行分析空间平滑滤波器消除或减弱图象中灰度值具有较大较快变化部分的影响, 这些部分对应频域中的高频分量, 所以可用频域低通滤波来实现空间锐化滤波器消除或减弱图象中灰度值缓慢变化的部分, 这些部分对应频域中的低频分量, 所以可用频域高通滤波来实现

42 频域技术与空域技术空域中的平滑滤波器在频域里对应低通滤波器频域越宽, 空域越窄, 平滑作用越弱频域越窄, 空域越宽, 模糊作用越强 F ( u ) f ( x ) 0 0 u x

43 频域技术与空域技术空域中的锐化滤波器在频域里对应高通滤波器空域有正负值, 一旦变为负数不再变为正数频域 / 空域的宽窄有什么关系和含义?

44 作业和练习 1. 基于 Matlab, 尝试将图像进行 fft 变换 ; 在变换域内进行低通高通滤波等操作, 并将滤波后信号用逆 fft 变换获得图像 2. 课后题 4-1

数字带通带阻高通滤波器的设计把一个归一化原型模拟低通滤波器变换成另一个所需类型的模拟滤波器, 再将其数字化直接从模拟滤波器通过一定的频率变换关系完成所需类型数字滤波器的设计先设计低通型的数字滤波器, 再用数字频率变化方法将其转换成所需类型数字滤波器

数字带通带阻高通滤波器的设计把一个归一化原型模拟低通滤波器变换成另一个所需类型的模拟滤波器, 再将其数字化直接从模拟滤波器通过一定的频率变换关系完成所需类型数字滤波器的设计先设计低通型的数字滤波器, 再用数字频率变化方法将其转换成所需类型数字滤波器数字带通带阻高通滤波器的设计把一个归一化原型模拟低通滤波器变换成另一个所需类型的模拟滤波器, 再将其数字化直接从模拟滤波器通过一定的频率变换关系完成所需类型数字滤波器的设计先设计低通型的数字滤波器, 再用数字频率变化方法将其转换成所需类型数字滤波器模拟原型方法 : 模拟低通 - 模拟带通 H ( j) H ( j) 3 3 3 模拟原型方法 : 模拟低通 - 模拟带通 H ( j) 模拟低通

目录 l 傅立叶变换 l 低通 高通滤波器 l 带通和带阻滤波器 l 同态滤波器

目录 l 傅立叶变换 l 低通高通滤波器 l 带通和带阻滤波器 l 同态滤波器