第 36 卷 第 9 期 自 动 化 学 报 Vol. 36, No. 9
2010 年 9 月 ACTA AUTOMATICA SINICA September, 2010
基于 Mean Shift 算法和 NMI
特征的目标跟踪算法研究
甘明刚
1, 2
陈 杰
1, 2
王亚楠
1, 2
金代中
1, 2
摘 要 针对传统 Mean shift 跟踪算法对空中运动目标跟踪效果不
理想的问题, 提出了基于 Mean shift 算法和归一化转动惯量 (Normal-
ized moment of inertia, NMI) 特征的目标跟踪算法. 算法中引入了
目标 NMI 特征, 建立了基于虚警概率最小原则和相似度二级判决门限
的跟踪策略, 对目标模型进行更新. 同时利用卡尔曼滤波, 在目标被遮挡
后进行估计预测. 实验表明该算法在空中运动目标存在较大形变、被遮
挡等情况下, 能够进行实时、稳定跟踪.
关键词 目标跟踪, 均值漂移, NMI 特征, 卡尔曼滤波
DOI 10.3724/SP.J.1004.2010.01332
A Target Tracking Algorithm
Based on Mean Shift and Normalized
Moment of Inertia Feature
GAN Ming-Gang
1, 2
CHEN Jie
1, 2
WANG Ya-Nan
1, 2
JIN Dai-Zhong
1, 2
Abstract The paper presents a target tracking algorithm
which is based on the mean shift algorithm and the normal-
ized moment of inertia (NMI) feature, because the result of the
moving target tracking in the air is not satisfactory by the tra-
ditional mean shift tracking algorithm. The NMI feature of the
target is introduced, and the studied tracking strategy based on
the minimum principle of the false alarm probability and the
two-stage decision threshold of the similarity is constructed in
the algorithm. The Kalman filter is used for the estimation and
prediction when the target is occluded. The experimental results
show that with the algorithm even the moving target in the air
is in large deformation and occlusion, the system can effectively
guarantee that the tracking is real-time and stable.
Key words Target tracking, mean shift, normalized moment
of inertia (NMI) feature, Kalman filter
基于图像信息的运动目标跟踪, 是以图像处理技术为核
心, 有机融合了计算机、传感器、模式识别、人工智能等多
种理论和方法的新型目标识别跟踪技术. Mean shift 算法作
为一种基于核函数密度梯度的无参估计方法, 由 Comaniciu
等
[1−2]
首先将该算法应用到目标跟踪领域中来. 相对于基于
边缘特征、变形模板、几何模型等原理的跟踪方法, 其计算量
小、实时性好、跟踪精度较高、适用性较强以及跟踪更加稳
定
[3]
, 近年来已被广泛应用在运动目标跟踪中.
常用的 Mean shift 目标跟踪算法多是以颜色直方图作
为特征模板的, 具有一定的鲁棒性. 然而, 在基本的 Mean
收稿日期 2009-11-11 录用日期 2010-03-30
Manuscript received November 11, 2009; accepted March 30, 2010
国家杰出青年科学基金 (60925011) 资助
Supported by National Science Fund for Distinguished Young
Scholars (60925011)
1. 北京理工大学自动化学院 北京 100081 2. 北京理工大学复杂系统智能控
制与决策教育部重点实验室 北京 100081
1. School of Automation, Beijing Institute of Technology, Beijing
100081 2. Key Laboratory of Complex Systems Intelligent Control
and Decision, Ministry of Education, Beijing Institute of Technology,
Beijing 100081
shift 算法中, 假设目标模型在跟踪过程中是基本上不变的,
在很多情况下跟踪效果不佳. 对此, 许多研究者提出一些解
决方法. 文献 [4] 研究了 Mean shift 算法中核窗宽自动选取
问题, 提供了一种目标模型更新方法; 文献 [5− 6] 选取了区
分度明显的不同目标特征来建立和更新模型; 文献 [7−8] 将
卡尔曼滤波器应用在跟踪中的模型更新; 文献 [9] 提出了基
于选择性子模型的目标更新. 但是, 当目标的运动状态较为
复杂时, 目前解决目标更新问题的方法不够完备, Mean shift
跟踪算法仍然会出现跟踪偏差或者完全跟不上目标的情况.
同时, 如果目标特征模型的建立和更新的计算复杂度很高,
将会很难满足实时跟踪的要求.
因此, 当前对目标模型更新的研究仍有许多不足, 特别
是当目标外形发生较大改变或被遮挡时, 将会很容易出现因
目标模型的更新问题而导致的跟踪失败
[10]
. 比如, 在对空中
运动目标的跟踪过程中, 由于目标自身的运动、非目标物干
扰、图像获取手段或其他外界因素的影响, 图像中目标的形
状、大小、明暗会不断变化, 目标模型发生较大改变, 基本的
Mean shift 跟踪算法不能满足跟踪要求. 因此, 在实际应用
中, 必须研究不同状态下目标模型的更新方法和跟踪策略,
以保证跟踪算法的稳健性. 针对以上问题, 本文提出了一种
基于 Mean shift 和归一化转动惯量 (Normalized moment of
inertia, NMI) 特征判断相结合的目标跟踪算法, 来实现对灰
度序列图像中的空中运动目标进行实时、准确、稳定的跟踪.
1 Mean shift 跟踪算法的基本原理
Mean shift 跟踪算法通常将颜色直方图作为目标的特征
描述, 建立概率密度函数, 用相似度函数度量目标模型和候
选目标之间的相似性, 通过求相似度函数的最大值得到关于
目标的均值漂移向量, 从而将跟踪问题转化为基于均值漂移
的模式匹配问题, 最终迭代搜索到目标的位置
[1−2]
.
假设 {x
∗
i
}
i=1,2,··· ,n
为目标模型的归一化像素位置, 其中
心位置为 0, 则目标模型特征可以表示为
b
q
q
q = {bq
u
}
u=1,2,··· ,m
,
且
P
m
u=1
bq
u
= 1, m 为直方图函数的特征空间大小, 目标模
型中特征值 u 出现的概率为
bq
u
= C
n
X
i=1
k(kx
∗
i
k
2
)δ[b(x
∗
i
) − u] (1)
其中, k(x) 为核函数的轮廓函数, δ 为 Kronecker delta 函数;
映射 b(x
∗
i
) 为 x
∗
i
处像素对应直方图中的颜色索引; C 为归
一化系数. 相应地, 设 {x
i
}
i=1,2,··· ,n
h
为候选目标模型的归
一化像素位置, 其中心位置为 y, 则候选目标模型特征表示为
b
p
p
p = {bp
u
(y)}
u=1,2,··· ,m
, 且
P
m
u=1
bp
u
(y) = 1, 候选目标模型中
特征值 u 出现的概率为
bp
u
(y) = C
h
n
h
X
i=1
k(
°
°
°
y − x
i
h
°
°
°
2
)δ[b(x
i
) − u] (2)
其中, C
h
为归一化系数,h 为核函数的窗宽. 由跟踪目标模型
与候选目标模型, 定义两者之间的距离为
d(y) =
p
1 − ρ[
b
p
p
p(y),
b
q
q
q] (3)
采用 Bhattacharyya 系数来度量两者的相似性, 表示为
bρ(y ) ≡ ρ[
b
p
p
p(y),
b
q
q
q] =
m
X
u=1
p
bp
u
(y)bq
u
(4)