输入图像经过经一层小波分解后，被分成4个子
图：
LL1—逼近子图，它代表输入图像水平和垂直
两个方向的低频成分； HL1—细节子图，它代表输入图像水平方向的 高频成分和垂直方向的低频成分；
LH1—细节子图，它代表输入图像水平方向的
低频成分和垂直方向的高频成分； HH1—细节子图，它代表输入图像水平和垂直 方向高频成分。
特征空间的分类方法可降低数据维数，降低计算复 杂性，但问题相关性较强，与特征提取的方法和效 果有很大关系。
3. 小波多尺度统计特征抽取及图像分类
图 像 特 征 提 取 及 分 类 方 法
图像的小波特征提取首先对输入图像做J层二维
小波分解；
因为小波变换具有很好的时频局部化特性，所
以可以将图像的不同底层特征变换为不同的小波系 数；
二维金字塔分解算法
令I(x,y)表示大小为MN的原始图像，l(i)表示相对于分析小 波的低通滤波器系数，i=0,1,2,…,Nl-1， Nl表示滤波器L的支 撑长度； h(i)表示相对于分析小波的高通滤波器系数， i=0,1,2,…,Nh-1， Nh表示滤波器H的支撑长度，则
I L x, y I H x, y x 0 ,1, , 1 Nl 1 Nh M 2
在逼近子图 LL1 上重复二维小波分解过程，进行
二层小波分解，如此继续分解，得到子图序列{LLJ， [HLk，LHk，HHk](k=1,2,…,J)}。
小波基与分解层次的选取是非常重要的，目前
还没有一个统一的标准。
小波基的选取一般考虑下列因素：
线性相位：如果小波具有线性相位或至少具有广义线性
N l 1
l i M , y M , y
l j I 2 x j mod
j0
1; y 0 ,1, , N 1
I LL x , y I LH x , y I HL x , y I HH x , y x 0 ,1, ,
i0 N h 1 H
N

N
h j I x , 2 x j mod
j0

1; y 0 ,1, ,
N 2
1
对逼近子图重复此过程，直到确定的分解水平，下 图是二层小波分解的示意图。
图6 图像多尺度分解，(a)一层分解，(b)二层分解
2. 图像分类问题现状
其他
分解层次
分解层次一般2-5层均可，要视具体应用而定，我们取为3， 即作3层小波分解，共得到10个子图，如图7。
特征抽取
每个子图抽取四个特征： 最大的小波系数； 最小的小波系数； 小波系数均值； 小波系数均方差。 这样对于一幅图像，可得到一个40(36)维的向量。
均值和方差的计算公式：
相位，则可以避免小波分解和重构时的图像失真，尤其是 图像在边缘处的失真；
紧支性和衰减性：紧支性和衰减性是小波的重要性质，
紧支宽度越窄或衰减越快，小波的局部化特性越好。计算 复杂度越低，便于快速实现；
正交性：用正交小波基对图像做多尺度分解，可得一正
交的镜像滤波器。低通子带数据和高通子带数据分别落在 相互正交的L2(R2)的子空间中，使个子带数据相关性减少；
LL x x y ;
x x y ; HL x x y ; HH x x y

LH
图像的二维小波变换包括沿行向(水平方向)和列向(垂直 方向)滤波和2-下采样，如图所示：
图1 小波例1
图2 小波例2
不是小波的例
图3
图4
x , x 构成Vj+1的正交基。
j k j k
x 和 x 满足下列关系式(二尺度方程)：
x x
2 l n 2 x n
n Z
2 h n 2 x n
n Z
其中 l n 称为低通滤波器， 且 h n ＝ 1 l 1 n
n
h n 称为高通滤波器。
信号的多尺度分解：
f x
j k

cn x n
0 j k
n Z

k
ck k x
J J

j 1 k
J
d k
j
J k
x
1 Nl 1 Nh 1 Nl 1 Nh M 2
N l 1
l i I x , 2 x i mod
L i0 N h 1 L
N

N
h j I x , 2 x j mod
j0 H

N l 1
l i I x , 2 x i mod
小波变换的一个重要性质是它在时域和频域均 具有很好的局部化特征，它能够提供目标信号各 个频率子段的频率信息。这种信息对于信号分类 是非常有用的。 小波变换一个信号为一个小波级数，这样一个 信号可由小波系数来刻画。
1.1 一维小波变换（一维多尺度分析） 设有L2(R )空间的子空间序列：
V 0 V1 V 2
每次实验从中选取视觉相似度较高的两类图像，在每一
类中随机选取40幅作为训练集，另32幅作为测试集，所以训 练集包含80幅图像，测试集包含64幅图像。
采用Db4小波对实验图像做三层小波分解，共有10个不同
频率字段的子图，每个子图抽取出4个特征，这样共有40个 特征.SVM分类器，高斯核函数。实验结果列于表1。
图5 图像滤波采样
说明：如图所示，首先对原图像I(x,y)沿行向(水平 方向)进行滤波和2-下采样，得到系数矩阵IL(x,y)和 IH(x,y)，然后再对IL(x,y)和IH(x,y)分别沿列向(垂直方 向)滤波和2-下采样，最后得到一层小波分解的4个 子图： ILL (x,y)—I(x,y)的（粗）逼近子图 IHL(x,y) — I(x,y)的水平方向细节子图 ILH (x,y) — I(x,y)的垂直方向细节子图 IHH (x,y) — I(x,y)的对角线方向细节子图
计算：
c 称为尺度系数， ck
j
d 称为小波系数，它们的
d
j k
j 1 c k l n 2 k n Z 一维 MALLAT 算法 j 1 d k h n 2 k n Z
1.2 二维小波变换（二维多尺度分析）
二维小波变换是由一维小波变换扩展而来的，二维尺度 函数和二维小波函数可由一维尺度函数和小波函数张量 积得到，即：
NEXT
RETURN
图8 COIL-20图像
RETURN
图9 7个位置的图像
RETURN 表1 两种方法的实验结果比较
4. 下一步的工作
从图像小波系数中抽取其它特征，如多尺度熵
特征； 小波函数逼近与径向基函数逼近的联系； 完善实验设计。
6. 参考文献
[1] Olivier Chapelle, Patrick Haffner, and Vladimir N. Vapnik. Support Vector Machines for Histogram-Based Image Classification. IEEE TRANSACTIONS ON NEURAL NETWORKS, VOL. 10, NO. 5, SEPTEMBER 1999. [2] Kwang In Kim, Keechul Jung, Se Hyun, and Hang Joon Kim, Support Vector Machine for Texture classification, IEEE TRANSACTIONS ON PATTERN ANALYSIS AND MACHINE INTELLIGENCE, VOL. 24, NO. 11 NOVEMBER 2002, pp. 1542-1550. [3] 万华林，M. U. Chowdhury. 基于支持向量机的图像语义 分类. 软机学报，2003，VOL.14 NO.11,PP. 1892-1899.
i i
1 M N 1 M N

M N
M
N
I i ( x, y ) Ii ( x, y ) i
2
x 1 y 1

x 1 y 1
NEXT
I(x,y) [128128]
RETURN
I1(x,y) [6464]
I1H(x,y) [6464]
I1V(x,y) [6464]
基于小波多尺度统计特征的图像分类
报告内容



1. 小波变换 2. 图像分类问题现状 3. 小波多尺度统计特征抽取及图像分类 4. 实验比较 5. 下一步工作 6. 参考文献
1. 小波变换
小波变换是强有力的时频分析(处理)工具，是在 克服傅立叶变换缺点的基础上发展而来的。已成 功应用于很多领域，如信号处理、图像处理、模 式识别等。
Vj 的正交基函数是由一个称为尺度函数的函数(x)经伸缩
平移得到的
k x 2 x k
j j
设Wj 是Vj 相对于Vj+1的正交补空间， Wj 的正交基函数是 由一个称为小波函数的函数(x)经伸缩平移得到的

j k
x 2
j
xk

小波函数必须满足以下两个条件的函数： (1) 小波必须是振荡的； (2) 小波的振幅只能在一个很短的一段区间上非零，即是局 部化的。如：
图像分类就是利用计算机对图像进行定量分析，
把图像或图像中的每个像元或区域划归为若干个类 别中的某一种，以代替人的视觉判读。 图像分类方法可分为： 图像空间的分类方法—利用图像的灰度、颜 色、纹理、形状、位置等底层特征对图像进行 分类；例如： 文献[1]利用灰度直方图特征对图像进行 分类； 文献[2]利用纹理特征对图像进行分类； 文献[3]采用纹理、边缘和颜色直方图混 合特征对图像进行分类 ；