2. 大连东软信息学院 计算机科学与技术系, 辽宁 大连 116023
2. Department of Computer Science and Technology, Dalian Neusoft University of Information, Dalian 116023, China
视觉跟踪是计算机视觉领域的热点研究问题,这主要是由于它在智能交通、人机交互、视频监控和军事等领域有着广泛的应用[1].传统的解决视觉跟踪问题的方法可以分为两类,即:基于随机采样的方法和确定性方法.在基于随机采样的跟踪算法中,粒子滤波算法是典型的代表[2, 3],它在处理非线性、非高斯多模态的跟踪问题中显示出较强的适应性.粒子滤波算法通过一个加权的粒子集合近似地表示后验滤波分布,其中,粒子是从后验滤波分布的状态空间中抽取的,后验的更新则依赖于粒子的序贯传播;粒子的权值则通过观测似然进行更新.在多目标跟踪问题中,一般采用多个随机变量描述目标状态,导致状态空间维数过高,而马氏链蒙特卡罗方法(Markov chain Monte Carlo,简称MCMC)[4]可以针对这类高维状态空间降低计算耗费.在确定性跟踪方法中,均值漂移(mean shift,简称MS)算法[5, 6]最具代表性,该算法一般通过计算目标区域与模板区域的距离,如巴特查里亚距离,并使其最大化来计算目标的位置.以上两种方法都受到研究人员的重视并取得了成功的应用,但它们都是基于平滑运动假设,即被跟踪目标在运动中不发生状态突变.在现实应用中,由于低帧率视频、目标动力学特征突变以及摄像机频繁的镜头切换所引起的突变运动大量存在,使得在传统的平滑运动限制背景下的跟踪算法很难有效地实现这些场景下的目标跟踪.
处理突变运动最直接的方法就是在跟踪过程中搜索整个状态空间,但这在实际应用中难以实现,因为搜索整个状态空间会使得计算时间急剧增长.文献[7]基于轮廓表示法,提出一种基于粒子群优化和水平集评估的突变运动跟踪框架,但目标的轮廓表示法过度依赖于图像边缘的准确定位,在低分辨率图像中很难提取.目前,针对突变运动跟踪,较多采用的方法是基于模板的矩形框目标表示.文献[8]针对低帧率视频引起的突变运动跟踪问题,提出了一种级联粒子滤波算法.在该算法中,以粒子滤波算法框架为基础,引入识别方法来检测被跟踪目标,在人脸跟踪中取得了较好的效果.但是这种跟踪方法需要较为可靠的先验目标观测模型以及离线学习过程.Kwon等人[9]提出一种基于Wang-Landau蒙特卡罗采样(WLMC)的突变运动跟踪算法,该算法以贝叶斯滤波框架为基础,使用基于态密度(density of states,简称DOS)的先验分布,将每一帧图像划分为多个相同大小的子空间,然后使用WLMC算法更新每一个子区域的态密度,并运用态密度引导跟踪器进行目标状态的更新.WLMC跟踪算法可以有效地避免传统的MCMC方法易陷入局部最优的问题,这在统计物理领域具有非常重要的意义.但是研究表明:WLMC算法尚未有严格的收敛理论作为支持,在某些应用中该算法也只能达到有限的统计精度[10].Kwon等人[11]将WLMC跟踪算法进行了一定的改进,提出了一种N-fold WLMC(NWLMC)突变运动跟踪算法,该算法能够处理目标位置突变和尺度突变的问题.NWLMC跟踪算法可以分为两个阶段,即:开发(exploitation)和探索(exploration).在开发阶段,采用边缘似然项的估计进行平滑运动的跟踪.在探索阶段,采用DOS项的估计跟踪突变运动.NWLMC算法能够以较少的样本数获得相对较好的跟踪结果.以Kwon等人提出的DOS思想为基础,周修庄等人[12, 13]针对突变运动跟踪难题提出了一种自适应随机逼近蒙特卡罗跟踪算法(adaptive stochastic approximation Monte Carlo,简称ASAMC).ASAMC算法是一种基于粒子滤波算法框架和随机逼近蒙特卡罗采样(SAMC)的方法,在粒子滤波算法框架内运用SAMC采样和DOS动态更新粒子权值,在突变运动跟踪中取得了优于WLMC算法的跟踪结果.但WLMC算法和ASAMC算法都采用了DOS的思想,需要将包含目标的图像切分成大小相等的子空间,在将每一帧图像进行切分时,如果图像尺寸过大,则会降低跟踪算法的效率,而且基于子空间DOS的目标状态转移不够鲁棒.此外,由于MCMC算法大都基于随机游动方法,使得从复杂、高维目标分布进行采样时非常缓慢,易陷入局部最优,为取得一个相对理想的目标状态,通常需要较长的迭代过程.
本文针对以上问题,提出一种基于有序超松弛Hamiltonian MCMC方法的突变运动跟踪算法.Hamiltonian马氏链蒙特卡罗方法,又称为混合蒙特卡罗方法(hybrid Monte Carlo,简称HMC),由Duane等人[14]提出,并被应用于诸如神经网络、机器学习等统计问题[15-17].Alfaki等人[18]将HMC算法应用于参数估计与不确定性量化问题. Choo等人[19]首次将该算法应用于高维空间的三维人体运动跟踪问题,实验结果表明,在高达28维的状态空间中,该算法比传统的粒子滤波算法快数千倍.HMC算法基于Hamiltonian动力学(Hamiltonian dynamics,简称HD),由于HD具有可逆性、保护Hamiltonian特性、容量保持以及对偶性等特征,使得HD可以与MCMC算法集成[20].HMC方法集成了动态状态序列的梯度信息,可以在一定程度上抑制传统MCMC方法的随机游动,确保了算法的快速混频以及更快速的收敛,并改进了马氏链的有效性[21].目前,多数文献中HD的实现采用蛙跳技术,HMC算法分为两步:在提议阶段,首先采用Gibbs采样方法抽取动量项,然后采用蛙跳方法迭代L次,每次步长为e,迭代完成后得到当前时刻的提议联合状态.在接受阶段,判断是否接受提议联合状态.文献[22]采用基于蛙跳方法实现的HMC算法跟踪突变运动,取得了优于WLMC跟踪算法的结果.然而,基于HMC的突变运动跟踪算法仍然不可避免地受到Gibbs采样阶段的随机游动行为的影响.同时,步长参数e的值难以确定:若e的值较小,则能保证更好的搜索状态空间,但算法时间耗费过高;如果值太大,则会使得算法提议状态的接受率偏低.因此,本文将通过两种策略对HMC跟踪算法进行改进:首先,采用有序超松弛(ordered over-relaxation,简称OR)方法抑制提议阶段Gibbs采样中的随机游动行为,OR方法适用于任何存在Gibbs采样的系统[23];其次,采用自适应步长的HD实现方法,在跟踪过程中动态调整步长.实验结果表明:基于以上两种策略的跟踪算法,能够在不增加时间耗费的前提下提高算法跟踪不同类型的突变运动的准确性.
1 基本HMC算法与有序超松弛方法 1.1 HMC算法在HMC算法中,将原始的状态变量XÎRk扩张为一个新的联合状态变量(X,V),VÎRk表示动量,k表示维数.然后,构造一个Hamiltonian函数如下:
H(X,V)=r(X)+K(V) (1)
r(X)=-log(q(X)) (2)
(3)
其中,q(X)表示目标概率分布,r(X)和K(V)分别表示动能与势能.基于该函数,可以从新的与exp(-H)成比例的概率密度函数中抽取随机样本.根据统计力学中正则分布的概念,Hamiltonian函数是关于联合状态(X,V)的能量函数,因此定义如下形式的联合概率分布:
(4)
其中,T表示系统的温度,Z是归一化常数.在Gibbs采样阶段,从高斯分布中抽取动量项V,然后经过L次蛙跳迭代,提议新的状态(X¢,V¢).在接受阶段计算接受概率a,然后判断该状态是否可以作为马氏链的下一个状态.
a=min(1,exp(-H(X¢,V¢)+H(X,V))) (5)
如果提议状态被拒绝,则当前时刻的状态(X,V)将会作为马氏链的下一个状态.
HMC算法是一种典型的具有各态历经性的算法[20],因此在搜索过程中不会陷入状态空间的某个子空间中,从而避免了传统MCMC算法易陷入局部最优的缺陷.
1.2 有序超松弛方法有序超松弛方法是一种用于抑制Gibbs采样中的随机游动行为的方法[23],在实际应用中可用于使用Gibbs采样方法的所有系统中.在Gibbs采样方法的一次迭代中,从条件分布中抽取样本值,有序超松弛方法需要从条件分布中抽取k个样本值.
假设函数CDF(V)为条件分布的累积分布函数,CDF-1(V)为其逆函数,对Vj进行有序超松弛的操作过程如下所示:
步骤1. 计算uj=CDF(Vj).
步骤2. 对uj进行有序超松弛变换,得到
步骤3. 计算即为Vj超松弛之后的结果.
2 基于有序超松弛Hamiltonian马氏链蒙特卡罗方法的跟踪算法 2.1 贝叶斯目标跟踪框架定义在t时刻的目标状态Xt包括目标的位置信息和尺度信息,即这里,分别表示目标的x,y坐标,表示目标的尺度.目标跟踪问题可以描述为贝叶斯滤波问题:在t时刻,根据观测值序列Y1:t= {Y1,Y2,…,Yt}递归地估计目标的状态Xt.贝叶斯滤波通过式(6)递归地更新目标状态的后验概率密度:
(6)
其中,p(Xt|Xt-1)表示系统状态转移模型,p(Yt|Xt)表示观测模型,c是一个归一化常数.计算出后验概率p(Xt|Y1:t)后,在数目为N的样本集上,通过执行最大后验估计获得当前时刻的目标状态估计值:(7)
即为当前时刻的目标状态估计值,表示在当前时刻的观测值序列下,最优的目标状态估计值.
然而,在实际运用贝叶斯滤波算法进行状态估计时,计算公式(6)所示的复杂积分是很难实现的,尤其是对于高维状态空间.为解决该问题,有研究者提出了运用Metropolis Hastings(MH)算法从提议分布密度Q中采样的思想[24],运用采样得到的样本集合近似目标的后验概率分布密度p(Xt|Y1:t).这样,公式(6)所示的复杂积分运算就被近似为求和运算.
2.2 基于有序超松弛Hamiltonian马氏链蒙特卡罗方法的跟踪算法如前所述,传统的MCMC方法通过模拟收敛于稳态分布p(Xt|Y1:t)的马氏链来逼近后验概率分布,但该类算法受随机游动的影响极易陷入局部最优.如果状态变量之间的相关性很高,采样分布就会变得非常狭窄而导致需要较长的迭代过程访问分布空间.图 1所示为二维状态空间中的随机游动行为示例:在相邻的两帧图像中,目标位置发生了突变,为了捕获被跟踪目标,基于随机游动的跟踪算法需要多次迭代来搜索状态空间,以找到关于目标的最有希望的区域.
Neal等人[23]提出采用有序超松弛方法可以抑制Gibbs采样的随机游动行为,本文首次将该方法引入到突变运动跟踪问题中,在Hamiltonian马氏链蒙特卡罗方法的跟踪框架下,使用该方法抑制提议阶段的Gibbs采样方法中的随机游动行为,改善算法的跟踪性能.在该跟踪框架中,首先将目标状态Xt扩张为一个新的联合状态(Xt,V),其中,V表示Xt的动量.在提议阶段,首先使用有序超松弛方法产生动量的值V¢,其步骤如下:
步骤1. 假设累积分布函数CDF(×)为高斯分布,计算V在该函数上的值u=CDF(V).
步骤2. 定义松弛度k,以k和u为参数构建二项分布,并从中随机抽取一个整数r~binomial(k,u).
步骤3. 如果r>k-r,构建贝塔分布,并生成随机数w~Beta(k-r+1,2r-k),令u¢=uxw; 否则,如果r<k-r,生成随机数w~Beta(r+1,k-2r),令u¢=1-(1-u)xw; 否则,如果r=k-r,令u¢=u.
步骤4. 根据累积分布函数的逆函数计算动量的值V¢=CDF-1(u¢).
在得到动量值之后,执行蛙跳迭代,并最终得到提议联合状态(X(*),V(*)).
在接受阶段,计算接受概率a*,并判断是否接受该提议联合状态,计算公式如下:
(8)
其中,(X(0),V(0))为当前时刻的初始化状态.算法执行完成后,得到当前时刻的样本集合基于有序超松弛Hamiltonian马氏链蒙特卡罗方法的跟踪算法(ordered over-relaxation Hamiltonian Markov chain Monte Carlo method,简称ORHMC)如算法1所示.
算法1. 基于ORHMC的跟踪算法.
初始化:
1. 初始化参数L=30,e=0.02.
2. 初始化松弛度k=10.
3. 初始化
FOR i=1,…,N,
1. 计算
2. 生成服从二项分布的随机数:r=Binomial(k,u).
3. IF r>k-r,
(1) 生成随机数w~Beta(k-r+1,2r-k).
(2) u¢=uxw.
ELSE IF r<k-r,
(1) 生成w~Beta(r+1,k-2r)
(2) u¢=1-(1-u)xw.
ELSE IF r=k-r
u¢=u
END IF
4. 计算V¢=CDF-1(u¢).
5. 令:
6. FOR j=1,…,L,
(1)
(2)
(3)
END FOR
7. 提议新的联合状态(X(*),V(*))=(X(L),V(L)).
8. 从均匀分布中抽取aa~U(0,1).
9. 根据式(8)计算接受概率a*.
10. IF aa<a*,接受提议状态,
ELSE 拒绝提议状态,
END FOR
得到样本集合,根据公式(7)计算目标状态估计值.
END
算法中的ÑV表示V的梯度,e为迭代步长.
ORHMC跟踪算法的关键是选择合适的参数e:如果其值过大,则会导致提议状态的拒绝率过高;如果其值偏小,则会导致计算量的增加.但通常很难确定一个最优的值,在跟踪过程中自适应地调整步长值,是解决该问题的一种有效方法.
2.3 自适应步长方法采用蛙跳方法实现Hamiltonian动力学比较普遍.如前所述,为了获得较好的跟踪结果,需要选择一个最优的步长值.Huang等人[25]提出了一种修正的Stömer-Verlet方法来实现Hamiltonian动力学.在该方法中,引入一个虚构变量g和尺度函数S来调整步长的大小.调整的过程如下所示:
(9)
(10)
(11)
(12)
(13)
其中,n表示第n次迭代.尺度函数S定义如下:
(14)
其中,V¢(X)=[ÑV(X)]T,V²(X)=[ÑV¢(X)]T.
如果尺度函数S满足S(X,V)=S(X,-V),那么该方法是完全明确、对称和时间可逆的[26].虚构参数g可以初始化为g0=S(X0,V0).与算法1中基于蛙跳方法的实现相比,该方法通过公式(10)中的项自适应地调整步长值,将该过程集成到算法1中的蛙跳步骤中,就得到自适应ORHMC算法(adaptive ORHMC,简称A-ORHMC),算法步骤如算法2所示.
算法2. 自适应ORHMC跟踪算法.
初始化:如算法1.
FOR i=1,…,N,
1. 执行有序超松弛变换得到V¢,如算法1.
2. 令
3. g0=S(X(0),V(0)).
4. FOR j=1,…,L,
(1)
(2)
(3)
(4)
(5)
END FOR
5. 提议步骤,如算法1.
6. 接受步骤,如算法1.
ENDFOR
得到样本集合,根据公式(7)计算目标状态估计值.
END
3 算法实现与实验结果分析本文采用11个视频序列来测试A-ORHMC跟踪算法的跟踪性能,见表 1.为对比分析实验结果,将A- ORHMC算法与其他几种算法进行比较,包括WLMC算法[9]、ASAMC算法[12]、HMC算法[22]、MCMC-PF算法[27]、自适应MCMC(AMCMC)算法[28].
在算法实现时,目标状态转移模型采用标准的二阶恒加速度动态模型,跟踪模型采用文献[29]中的基于颜色的跟踪模型,基于HSV颜色直方图的相似度定义滤波分布的似然函数如下:
(15)
其中,HSr是参考目标模型;HS(Xt)为t时刻的候选目标模型,是预先定义的调节参数;D(×)用于计算HSV颜色直方图的巴氏距离.本文采用的视频序列均为国外研究机构的开放视频序列.
为了评估本文提出算法的性能,我们采用信息检索领域中较为常用的f-measure作为度量指标,该指标也常用于视觉跟踪领域[9].召回率f(recall)也称查全率,准确率z(precision)又称为正确率,用来度量不同跟踪算法的跟踪性能,定义如下:
(16)
(17)
其中,为t时刻估计的目标区域位置,为标定的t时刻的目标区域位置.f-measure定义为(18)
如果估计目标位置与真实目标位置重合,则f-measure为1.0;其值越大,则跟踪结果越好.
3.1 定性比较为了定性地比较A-ORHMC跟踪算法的性能,实验中针对不同情形的突变运动进行了测试,分别是:摄像机镜头切换、低帧率视频、快速移动目标、运动学突变.同时,针对遮挡问题进行测试.若不作特殊说明,各算法使用的样本(粒子)数默认为300个.
· 目标快速运动
目标的快速移动和低帧率视频能够引起运动突变,本实验采用Face序列(数据来自于Clemson大学工程科学学院[30])和Tennis序列测试A-ORHMC算法处理快速运动目标的能力.实验中采用的样本数目为:A-ORHMC算法50个,WLMC,ASAMC,HMC算法100个,MCMC-PF和AMCMC算法各150个.跟踪结果如图 2所示:A-ORHMC算法能够以较少的样本数目准确跟踪快速移动的人脸和网球运动员,而其他5种算法不能连续地捕获快速运动带来的不确定性,整体跟踪结果比A-ORHMC算法要差.该实验结果表明:A-ORHMC算法在处理目标快速移动带来的突变时,性能优于其他算法.
· 镜头切换
在本实验中,实验数据采用文献[9]中的ChoiHongMan和YoungKi序列,跟踪场景中的行人和拳击运动员.实验结果如图 3所示.
从结果中可以看出:目标不发生镜头切换时,A-ORHMC算法可以有效地跟踪目标.
· 在ChoiHongMan序列中,在第236帧~第248帧发生镜头切换时,本算法仍然可以有效地捕获拳击运动员;
· YoungKi序列中的突变在第546帧和第547帧,A-ORHMC算法能够准确地跟踪行人.WLMC算法、ASAMC算法、HMC算法在部分发生突变的帧可以捕获目标,但无法实现连续的有效跟踪;MCMC和AMCMC算法则频繁地失去目标.在处理镜头切换引起的运动突变时,A-ORHMC算法的跟踪结果优于其他算法.
· 目标动力学突变
该实验采用Pingpong序列,跟踪一个在球拍上来回运动的白色乒乓球.当乒乓球敲击球拍时,球的运动发生突然的变化,在视频序列的最后几帧发生了摄像机镜头的变化,使得跟踪小球的难度加大.结果表明:A-ORHMC算法能够在发生运动突变以及摄像机镜头发生变化时准确地跟踪小球,如图 4所示;而其他5种算法在特定帧上可以取得可接受的跟踪结果,但整体性能比A-ORHMC算法要差.
· 目标频繁遮挡
采用SeqMS视频序列测试不同算法在处理频繁遮挡时的能力.场景中,人脸被双手频繁地全部或部分遮挡.实验中,A-ORHMC算法使用100个样本,WLMC,ASAMC,HMC算法各300个,其他算法600个,跟踪结果如图 5所示.A-ORHMC算法能够在遮挡前后准确地捕获人脸,在人脸被完全遮挡时仍然能够准确定位其位置;而其他5种算法无法持续地捕获人脸,跟踪的准确性比A-ORHMC算法要差.
3.2 定量比较
在定量比较实验中,主要考察各算法的跟踪成功率、降采样区间对成功率的影响以及算法的计算时间等指标.
· 算法跟踪成功率
对视频序列的某一帧来说,如果f-measure的值大于0.5,则认为在当前帧目标被准确跟踪.成功率定义为被准确跟踪到目标的帧数和视频序列总帧数的比值.
图 6所示为各算法在不同视频序列上的成功率对比.A-ORHMC算法与ASAMC,WLMC,HMC算法使用300个样本,MCMC和AMCMC算法各600个样本.
从结果中可以看出:除在ChoiHongMan序列上WLMC算法比A-ORHMC算法的成功率稍高之外,在其他序列上,A-ORHMC算法的跟踪成功率明显高于其他算法.ASAMC算法、WLMC算法及HMC算法在部分序列上的成功率相差不多,总体优于MCMC和AMCMC算法,A-ORHMC算法的成功率则高于其他算法.
· 降采样区间的影响
在处理低帧率视频引起的突变运动时,降采样区间的大小对算法的跟踪结果具有非常重要的影响.实验中,采用Tennis10,Tennis15,Tennis20,Tennis30,Tennis35序列,针对6种算法在不同降采样区间的低帧率视频上进行跟踪测试,样本数为:A-ORHMC算法、ASAMC算法、WLMC算法、HMC算法各300个,其他算法各600个.图 7所示为各算法在不同降采样区间上的成功率曲线,本文算法的成功率明显高于其他几种算法,A-ORHMC算法、WLMC算法和ASAMC算法不易受降采样区间的影响,而MCMC算法和AMCMC算法则受降采样区间的影响比较明显.本实验结果表明,A-ORHMC算法对不同类型的突变运动跟踪具有鲁棒性.
· 算法计算时间
表 2所示为不同跟踪算法在Tennis序列上的运行时间(秒/帧).实验中,各算法均采用300个样本.结果显示: A-ORHMC,HMC,WLMC,ASAMC这4种算法的运行时间几乎一样,而AMCMC与MCMC算法的运行时间与以上4种算法也相差不多.这表明,A-ORHMC算法与ASAMC,WLMC,HMC算法相比没有引入额外的计算时间.
4 结 论
本文提出一种基于有序超松弛Hamiltonian马氏链蒙特卡罗方法的跟踪算法,解决突变运动跟踪中的难点问题.为抑制Gibbs采样过程中的随机游动行为,采用有序超松弛方法抽取动量值,该方法可以有效地消除样本之间的依赖性,提高跟踪算法的跟踪性能.同时,采用自适应步长方法实现Hamiltonian动力学,在跟踪过程中自适应地调整步长.在贝叶斯滤波框架下集成以上两种策略,得到A-ORHMC跟踪算法.实验结果表明:该算法在处理不同类型的突变运动时,与传统的基于随机游动的MCMC跟踪算法相比,表现出较好的性能,同时不需要额外的计算时间.本文提出的跟踪算法可以与复杂的目标外观模型结合,来处理目标外观剧变和尺度剧变的跟踪问题,这也是我们下一步的研究目标.
致谢 感谢韩国首尔大学计算机视觉实验室Junseok Kwon博士和首都师范大学信息工程学院周修庄博士在此研究中提供的帮助.
[1] | Yilmaz A, Javed O, Shah M. Object tracking: A survey. ACM Computing Surveys, 2006,38(4):1-45 . |
[2] | Li PH. A novel color based particle filter algorithm for object tracking. Chinese Journal of Computers, 2009,32(12): 2454-2463 (in Chinese with English abstract) . |
[3] | Martinez-del-Rincon J, Orrite C, Medrano C. Rao-Blackwellised particle filter for color-based tracking. Pattern Recognition Letters, 2011,32(2):210-220 . |
[4] | Smith K, Perez, DG, Odobez J. Using particles to track varying numbers of interacting people. In: Proc. of the Int’l Conf. on Computer Vision and Pattern Recognition. San Diego: IEEE Press, 2005. 962-969 . |
[5] | Comaniciu D, Ramesh V, Meer P. Real-Time tracking of non-rigid objects using mean shift. In: Proc. of the Int’l Conf. on Computer Vision and Pattern Recognition. Hilton Head: IEEE Computer Society, 2000. 142-149 . |
[6] | LI PH. An improved Mean Shift algorithm for object tracking. Acta Automatica Sinica, 2007,33(4):347-354 (in Chinese with English abstract). |
[7] | Li W, Zhang XQ, Hu WM. Contour tracking with abrupt motion. In: Proc. of the Int’l Conf. on Image Processing. Egypt: IEEE Press, 2009. 3593-3596 . |
[8] | Li Y, Ai HZ, Yamashita T, Lao S, Kawade M. Tracking in low frame rate video: A cascade particle filter with discriminative observers of different life spans. In: Proc. of the Int’l Conf. on Computer Vision and Pattern Recognition. Minneapolis: IEEE Press, 2007. 1-8 . |
[9] | Kwon JS, Lee KM. Tracking of abrupt motion using Wang-Landau Monte Carlo estimation. In: Proc. of the 10th European Conf. on Computer Vision, Vol.5302. Berlin, Heidelberg: Springer-Verlag, 2008. 387-400 . |
[10] | Liang FM, Liu CH, Carroll RJ. Stochastic approximation in Monte Carlo computation. Journal of the American Statistical Association, 2007,102(477):305-320 . |
[11] | Kwon JS, Lee KM. Wang-Landau Monte Carlo-based tracking methods for abrupt motions. IEEE Trans. on Pattern Analysis and Machine Intelligence, 2013,35(4):1011-1024 . |
[12] | Zhou XZ, Lu Y. Abrupt motion tracking via adaptive stochastic approximation Monte Carlo sampling. In: Proc. of the Int’l Conf. on Computer Vision and Pattern Recognition. San Francisco: IEEE Press, 2010. 1847-1854 . |
[13] | Zhou XZ, Lu Y, Lu JW, Zhou J. Abrupt motion tracking via intensively adaptive Markov-chain Monte Carlo sampling. IEEE Trans. on Image Processing, 2012,21(2):789-801 . |
[14] | Duane S, Kennedy AD, Pendleton BJ, Roweth D. Hybrid Monte Carlo. Physics Letters B, 1987,195:216-222 . |
[15] | Ishwaran H. Applications of hybrid Monte Carlo to Bayesian generalized linear models: Quasicomplete separation and neural networks. Journal of Computational and Graphical Statistics, 1999,8(4):779-799 . |
[16] | Akhmatskaya E, Reich S. New hybrid Monte Carlo methods for efficient sampling: From physics to biology and statistics. Progress in Nuclear Science And Technology, 2011,2:447-462. http://www.aesj.or.jp/publication/pnst002/index.html |
[17] | Schimidt MN. Function factorization using warped gaussian processes. In: Proc. of the Int’l Conf. on Machine Learning. Montreal: ACM Press, 2009. 921-928 . |
[18] | Mohammed A. Improving efficiency in parameter estimation using the Hamiltonian Monte Carlo algorithm [MS. Thesis]. Bergen: University of Bergen, 2008. |
[19] | Choo K, Fleet DJ. People tracking using hybrid Monte Carlo filtering. In: Proc. of the Int’l Conf. on Computer Vision. Vancouver: IEEE Press, 2001. 321-328 . |
[20] | Brooks S, Gelman A, Jones G, Meng X. Handbook of Markov Chain Monte Carlo. Chapman and Hall/CRC Press, 2011. 113-162. http://mcmchandbook.net/index.html |
[21] | Hanson KM. Markov chain Monte Carlo posterior sampling with the Hamiltonian method. Proc. of the SPIE, 2001,4322:456-567 . |
[22] | Wang F, Lu M. Hamiltonian Monte Carlo estimator for abrupt motion tracking. In: Proc. of the Int’l Conf. on Pattern Recognition. Tsukuba: IEEE Press, 2012. 3066-3069. |
[23] | Neal RM. Suppressing random walks in Markov chain Monte Carlo using ordered overrexlation. Technical Report, No.9508, Department of Statistics, University of Toronto, 1995. [http://arxiv.org/abs/bayes-an/9506004] |
[24] | Brooks S, Gelman A, Jones G, Meng X. Handbook of Markov Chain Monte Carlo. Chapman and Hall/CRC Press, 2011. 3-48. |
[25] | Huang WZ, Leimkuhler B. The adaptive verlet method. SIAM Journal on Scientific Computing, 1997,18(1):239-256 . |
[26] | Holder T, Leimkuhler B, Reish S. Explicit variable step-size and time reversible integration. Applied Numerical Mathematics, 2001, 39(3-4):367-277 . |
[27] | Khan Z, Balch T, Dellaert F. MCMC-Based particle filter for tracking a variable number of interacting targets. IEEE Trans. on Pattern Analysis and Machine Intelligence, 2005,27(11):1805-1819 . |
[28] | Roberts GO, Rosenthal JS. Examples of adaptive MCMC. Journal of Computational and Graphical Statistics, 2009,18(2):349-367 . |
[29] | Perez P, Hue C, Vermaak J, Gangnet M. Color-Based probabilistic tracking. In: Proc. of the 7th European Conf. on Computer Vision, Vol.2350. Berlin, Heidelberg: Springer-Verlag, 2002. 661-675 . |
[30] | Birchfield S. BMP image sequences for elliptical head tracking. 1998. http://www.ces.clemson.edu/~stb/research/headtracker/seq/ headtracker_sequences.zip |
[2] | 李培华.一种新颖的基于颜色信息的粒子滤波器跟踪算法.计算机学报,2009,32(12):2454-2463 . |
[6] | 李培华.一种改进的Mean Shift跟踪算法.自动化学报,2007,33(4):347-354. |