第3章图象变换讲义.doc

第章图象变换 3.1 傅里叶变换 1. 1-D 傅里叶变换 对一个连续函数f (x)等间隔采样可得到一个离散序列。设共采了N个样，则这个离散序列可表示为{f (0)，f (1)，f (2)，…，f (N – 1)}。借助这种表达，并令n为离散实变量，u为离散频率变量，可将离散傅里叶变换对定义为： 每个u值都确定所对应的正弦和余弦对的频率，所以称为频率变量。 2. 2-D傅里叶变换 设图象的大小为M(N，2-D傅里叶变换为： 其反变换为： 例 实际图象的傅里叶频谱 下图给出两幅实际图象和他们的傅里叶频谱图。图(a)的图象反差比较柔和，反映在傅里叶频谱上低频分量较多，频谱图中心值较大（中心为频域原点）。图(b)的图象中有较规则的线状物，反映在傅里叶频谱上也有比较明显的射线状条带。 ? (a) (b) 3. 卷积定理 两个1-D连续函数的卷积定义为： 如果f (x)的傅里叶变换是F(u)，g(x)的傅里叶变换是G(u)，那么： 推广到2-D，两个2-D函数的卷积定义为： 2-D时的卷积定理为： 设f(x,y)与g(x,y)分别是大小为A(B、C(D 的两个图象，将其周期性扩展为M(N，M=A+C-1，N=B+D-1。则f(x,y)与g(x,y)的卷积的2-D傅里叶变换为f(x,y)与g(x,y)的2-D傅里叶变换的乘积，即： 这样，图象的线性滤波结果可以方便地用图象与滤波器单位响应的2-D傅里叶变换的乘积快速计算。设函数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(x, y)，h(x, y)，f (x, y)的傅里叶变换。 然后，通过2-D傅里叶反变换求得滤波器的输出结果： 一维Gabor函数的富里叶变换为： 当时，，此时带宽，等效带宽。 二维Gabor函数二进制尺度γ伸缩均匀旋转角度θ富里叶变换为： 可见，经过伸缩和旋转的Gabor函数相当于带宽和中心频率都随尺度因子变化的带通滤波器，但带宽和中心频率比不变ω0=π，δ=1。当γ=0，1，2，3；θ=0时Gabor滤波器的纵切面如图γ=0，1，2，3；θ=00，，，，，如图。 带宽和中心频率比：采用γθ=00，，，γ=0，1，2，3的Gabor滤波器对蝴蝶斑纹图象的亮度滤波效果如图（左上角为花斑图像原图）。 =00 =450 =900 =1350 =0 =1 =2 =3 =4 3-1 不同尺度和角度的Gabor滤波器滤波效果3.2 离散余弦变换 2-D离散余弦变换（DCT）和其反变换由以下两式定义： 其中a(u)由下式定义： 2-D DCT的正变换核表达式为： DCT的变换核具有可分离性和对称性，即 例 离散余弦变换示例 下两图给出离散余弦变换的一个示例，其中左图是一幅原始图象，右图是对左图的离散余弦变换结果（变换幅值）。右图左上角对应低频分量，由图可见，左图中的大部分能量在低频部分。 ? 3.3 小波变换 1．连续小波变换 小波变换是利用基本小波的尺度伸缩和位移对信号进行变换。基本小波ψ(x)是具有震荡衰减特性的实值函数，且满足下列条件： （3-9） 其频谱满足： （3-10） 小波基函数为： （3-11） 其中，尺度参数a为实数，指示小波基函数的伸缩尺度。位移参数b为实数，指示沿x轴的平移位置。 函数f(x)连续小波变换定义为： （3-12） 连续小波逆变换为： （3-13） ?? 小波变换用卷积表示： ? 重构f(x)同样可以由滤波器族的输出经再次滤波后组合而成，即： 因而，不同尺度的小波变换相当于用不同带宽和中心频率的滤波器对信号进行滤波，小波变换系数相当于滤波器的输出。二维函数f(x,y)的连续小波变换定义为： （3-14） （3-15） 二维连续小波逆变换为： （3-16） ，其中为整数。离散化的方式： 特别地，取，称二进小波变化。 二进小波族： ? （3-17） 二进小波变换系数： （3-18） 任何都可以表示为： （3-19）小波函数可以由尺度函数的伸缩和平移的线性组合生成，而尺度函数也满足双尺度方程，即某一