您好,欢迎来到纷纭教育。
搜索
您的当前位置:首页基于压缩感知的粒子滤波跟踪算法

基于压缩感知的粒子滤波跟踪算法

来源:纷纭教育
第37卷第11期 系统工程与电子技术 Systems Engineering and Electronics Vo1.37 No.11 November 2O15 网址:www.sys—ele.corn 2015年11月 文章编号:1001—506X(2015)11—2617—06 基于压缩感知的粒子滤波跟踪算法 吴晓雨,吴凌琳,杨 磊 (中国传媒大学信息工程学院,北京100024) 摘 要:针对运动目标跟踪存在的目标遮挡和光照变化问题,提出一种基于压缩感知的粒子滤渡跟踪算法。 将改进的压缩感知跟踪算法提取的特征融合到粒子滤波跟踪框架中,并对压缩感知提取的特征和原始粒子滤波 中的颜色特征进行可信度判定,能够较好地处理图像序列中由于目标遮挡和光照变化所带来的影响。此算法在 公开数据库中进行测试,实验结果表明,提出的算法与已有改进压缩感知跟踪算法和粒子滤波跟踪算法相比,鲁 棒性更好,能准确实时地对目标进行跟踪。 关键词:运动目标跟踪;粒子滤波;压缩感知;可信度判定 中图分类号:TN 911.73 文献标志码:A DOI:10.3969/j.issn.1001~506X.2015.11.30 Particle filtering tracking based on compressive sensing WU Xiao—yu,WU Ling~lin。YANG Lei (School of Information Engineering,Communication University of China,Beijing 100024,China) Abstract:To deal with the target occlusion problem and illumination changes in moving target tracking,a particle filtering algorithm based on compressive sensing is proposed.The extracted features are added by eom— pressive sense of the improved compressive tracking(CT)algorithm into the framework of particle filtering tracking.The credibility of extracted features including the color features of original particle filtering and com— pressive sensing features is judged,which deals with the target occlusion effects and illumination changes.The algorithm is tested in public database and experimental results show that the proposed algorithm brings about better robustness and tracks targets accurately in real time in comparison with the improved CT algorithm and particle filtering algorithm. Keywords:moving target tracking;particle filter;compressive sensing;credibility j udge 0 引 言 运动目标跟踪作为机器视觉研究领域中的一个研究热 的跟踪,粒子重采样策略较好应对了目标被遮挡问题,但在 目标环境的光照发生变化时,粒子滤波的鲁棒性很差。文 献[7]基于文献E6]的压缩感知跟踪(compressive tracking, CT)算法提出了改进的算法,在原始的特征提取矩阵的基 础上,生成一个提取灰度的特征提取矩阵,该算法能在目标 形态和光照发生剧变等情况下准确实时地跟踪目标。但当 目标被遮挡时,改进的CT算法则由于正负样本不断更新, 点,已经被广泛应用到智能人机交互、交通安全监控、医疗 成像、军事指导、天文探测、电视制作等领域。作为一个研 究热点,运动目标跟踪针对不同的环境场景存在不同的问 题,如目标旋转、大小改变、光照变化、目标遮挡等,这些问 题都给目标跟踪带来不同方面的难题。 近年来,针对这些问题,学者们提出了不同的跟踪算 正样本被遮挡物取代,导致遮挡后无法跟上。 视频序列中,经常会同时出现目标遮挡和光照变化 的情况,而上述算法都不能对其进行有效的跟踪。针对 上述存在的跟踪问题,本文首先剖析了粒子滤波跟踪和 压缩感知跟踪算法的优缺点,而后提出了一种基于压缩 感知的粒子滤波跟踪算法,即在粒子滤波框架下结合改 进的CT算法中的压缩感知提取特征。在公开测试库上 的实验结果表明,提出的跟踪算法同时解决光照变化和 法。文献[1]提取了目标区域的SIFT特征并结合mean— shift的跟踪算法,较好地解决了跟踪目标的尺度及旋转变 化的问题。文献[2]提出的局部无序跟踪(1ocally orderless tracking,LOT)算法,对于运动目标大小变化,都能有效地 跟踪。但文献[1—2]对于光照变化和有遮挡的运动跟踪鲁 棒性有待提高。文献[3—5]基于粒子滤波框架实现对目标 收稿日期:2014~12—29;修回日期:2015—04—24;网络优先出版日期:2015—07—06。 网络优先出版地址:http://www.cnki.net/kcms/detail/11.2422.TN.20150706.1705.015.html 基金项目:国家科技支撑计划(2014BAHlOF00,2012BAH01F01—01)资助课题 ・2618 ・ 系统工程与电子技术 第37卷 目标遮挡时跟踪失败的问题,在复杂场景下具有较好的 跟踪效果。 l相关算法介绍 1.1粒子滤波跟踪算法 粒子滤波器是一种通过蒙特卡罗仿真实现递推贝叶斯 滤波的方法,实质是随机采样运算将积分转化为有限样本 点的求和运算,用一个带权重值的随机样本集(即粒子集 (z 一 ))来表示后验密度函数,以此来计算系统的状态估 计。样本的数量足够多时可较好地描述目标状态的后验概 率分布,得到目标状态的最优贝叶斯估计。k时刻的后验 密度可以近似表示为 N p( l 21 :)≈2 议J (z 一z ) (1) 一1 式中,占(・)为狄克拉函数;N为粒子总数;k为序列帧数;i 为粒子标志。 总体来说,粒子滤波跟踪[8 主要步骤总结如下:首先 随机采样一个粒子集并提取粒子特征;其次通过特征计算 粒子权重来表示后验密度函数,找到最相似粒子作为跟踪 目标;最后对粒子进行重采样。 1.2压缩感知跟踪算法 文献E63提出了基于压缩感知的跟踪算法,利用图像的 稀疏性和压缩感知 算法生成了一个稀疏的随机测量矩 阵,然后图像原始特征通过向此随机测量矩阵投影,去除了 大量的冗余信息,并得到压缩的低维特征 对这些低维特 征采用朴素贝叶斯分类学习的跟踪算法框架进行目标分类 跟踪。上述压缩感知提取特征的公式为 V—RX (2) 式中,Xq 为目标候选区域转换成1维后的信号;R∈ R (n《m)为特征提取矩阵;V∈R 为特征。而式(2)中 不同的R会提取出不同的特征。文献[6]对R定义如下:  I1 概率为 S 一 ×{。 概率为1 (3) S 1 l一 t,概率为 2S 式中,s通过平均概率在2~ 中随机选 取,文献[7]指出实 际特征提取为 roRectso一 (4) √NR 式中,Rects为在候选区域随机选取的图像块(即R中为√5 或 s的区块);NR为图像块数目(NR—s);j表示第J个 图像块IrI,在1与一1中等概率随机选取,R的值在同一图 像块中相同。 2 基于压缩感知的粒子滤波算法 由于粒子滤波跟踪算法中粒子随机采样和重采样的核 心思想,使得该算法对于目标有遮挡的运动有很好的跟踪 效果;但粒子滤波跟踪算法中仅用到颜色特征,对于光照的 变化鲁棒性很差。而CT算法可以快速提取到低维且有效 的灰度和纹理特征,有效解决粒子滤波存在的问题,然而由 于对正样本的不断更新,导致在目标遮挡后正样本偏离原 始目标而使跟踪丢失。因此,本文提出将CT算法中的压 缩感知提取特征融合到粒子滤波框架中,能够同时解决视 频跟踪过程中遮挡和光照变化的问题。本文算法流程图如 图1所示,首先在目标区域周围采样两组粒子集,利用直方 图和压缩感知分别提取颜色特征和灰度纹理特征;其次通 过直方图距离和贝叶斯分类分别计算两组粒子的权重;再 以一定的决策方法对两组粒子的跟踪结果进行可信度判定 来确定当前的跟踪结果;最后进行重采样来防止粒子退化。 I初始化跟踪目标I l 采样两组粒子 ,P2 I采样 -。—。_ … 一十…………… 提gZP,粒子的l I压缩感知提取 I特征 颜色直方图 f l粒子的灰度和纹理I提取 + +  II巴特沃斯系数Il确定粒子权熏I I斯分类一l朴素贝叶l、1分类器I更新r权l计算 重  ll。 ● I粒子跟踪结果决策I确定日标 一 Y 。 I 粒子重采样 I重采样 图1本文算法流程图 2.1粒子滤波框架下的特征提取 第一组粒子P 利用直方图来提取粒子特征,颜色直方 图在目标遮挡前后变化不大,并且在目标区域周围选用了 高斯分布随机采样,让采样粒子尽可能覆盖目标的各种状 态。另外利用重采样降低粒子退化现象,保持粒子的有效 性和多样性。因此,粒子滤波能够有效地解决目标遮挡问 题。但对于跟踪过程中的光照变化,颜色直方图鲁棒性 很差。 为了减少光照的影响,需要加入图像灰度和纹理等特 征,而这些特征的提取比较又会大大增加算法的计算量。 而压缩感知中心思想是采样和压缩共同进行,得到的压缩 样本数据量远远小于传统采样所得的数据量。因此,第二 组粒子P 采用压缩感知来提取图像灰度和纹理特征。 文献[7]指出特征提取矩阵(见式3)每行中1和一1同 时存在的概率为0.71,提取的特征更多地表现为图像的纹 理特征,纹理特征在目标形态和光照发生剧变等情况下并 不稳定。为了增强跟踪的稳定性,文献[7]中利用式(3)生 成R ,而后根据R 生成新的特征提取矩阵Rz,从而利用 R 和R 提取具有互补的纹理特征和灰度两种特征。新的 特征提取矩阵生成公式如下: I I rh I, 中同时存在1和一1 f I(-1) f rh l,r¨中存在1或一1 第11期 吴晓雨等:基于压缩感知的粒子滤波跟踪算法 2.3粒子跟踪结果决策 以两组粒子中权重最大粒子位置作为最相似粒子位 置,P 和P 的最相似粒子位置分别为P 和C ,可以确定 式中,i表示特征提取矩阵R 和R。的第i行;k表示R 中 第i行中第k个非零值。因此本文参考文献E7],将两个特 征矩阵融入到粒子滤波跟踪框架中,对采样的粒子集P。提 取纹理特征和灰度均值特征。 2.2粒子权重 2.2.1提取颜色特征的粒子权重计算 目标的最优估计位置: 一ap +(1一a)C 。这里 是融合 的权重。由于两组粒子的特征都不能同时解决目标遮挡和 光照变化的情况,因此需要有个判断条件来调整粒子融合 权值a。 本文通过计算两组粒子的权重来表示两个后验密度函 数。其中针对粒子集P ,本文采用经典的巴特沃斯距离来表 示粒子与目标模板之间颜色直方图的相似度(存放的是粒子的 对于粒子集P 采用颜色直方图特征,粒子集越聚合, 估计的目标位置方差就会越小,故基于颜色直方图特征估 颜色直方图),希望找到与目标颜色分布最接近的粒子。因此, 巴特沃斯距离越大相似度越高,而粒子的权重也就越大。 2.2.2提取灰度和纹理特征的粒子权重计算 针对粒子集尸2,本文通过朴素贝叶斯分类器对压缩感知 提取的特征进行相似性判别,寻找与标准目标最相似的粒子。 对每个样本x∈R- ,它的低维表示是y∈ ( 《m),假定y 中的各元素是分布的,构建朴素贝叶斯分类器为 H(V)一ln Ⅱp(v I Y一1)p(y一1)  !.:::. ................................... ........................... .....................一 Ⅱp(v。I Y—O)p(y—o) i=1 In(麓— ) 其中,假设两个类的先验概率相等 ( 一1)一 (y—O),Y∈ {0,1}代表二值样本标签, 一0表示负样本, 一1表示正样 本。文献[143指出高维向量随机投影向量几乎满足高斯分 布。所以,分类器H(V)中的条件概率p(v lj,一1)和P(73 l —O)被假定属于高斯分布,那么 fp(v l Y一1)一户( l Y一1)夕( l Y一1) l ( l Y一1)~N( , 1 ) j 夕(口 l Y一1)~N( 。, 1。 (7) l p(v I y—o)一p(口 j Y—o) ( l Y—o) I p(v l —o)~N(/2 , 0 ) 【p(v {y一0)~N( , ) 式中, 和 分别表示两种特征矩阵提取的特征; , 和 , 分别表示R 和R2提取特征的正样本的均值和标 准差; 0 , 和 ? , 分别表示R 和R。提取特征的负样 本的均值和标准差。 H(y)为后验概率,其值越大则与标准目标特征区域越相 似,因此第二组粒子以H(y)值作为粒子权重。选取H(y) 值最大的粒子作为新一帧的目标区域,并对相关系数 和 进行更新 t— +(1一 ) 一 A ( ) +(1一 )( ) + (1一 )( 一 ) (8) 式中,z一1代表目标样本,z一0对应背景样本;r一1时利用 R 求取特征,r一2时利用R 计算特征; 对应着更新率,其 值越大更新越慢,保留之前特征越多。 计的目标位置被赋予高的可信度。当在画面出现光照变化 时,粒子空间位置会从原来较为集中变得相对分散,其位置 方差随之变大,因此当方差大于某一阈值T 时,可判定为 粒子滤波不可信(本文取T 一45),则a一0。 对于粒子集P。采用压缩感知提取特征,当目标被遮挡 时,采样的图像片与负样本更为相似,因此H(y)将突然变 小,gosh于某一阈值Tz(本文取Tz一2),以此来判断压缩感 知提取的特征是否可信,则a一1。 当第一组粒子空间位置方差D和第二组压缩感知概 率值H( )均在设定的阈值内,则代表两组粒子都可信。 对两组粒子的结果进行归一化来确定a取值。对于P 粒 子,归一化值为a一(T 一D)/T ;对于P 粒子,归一化值 为6一(H(V)一T2)/ , 为设定的一个压缩感知概率中 可能存在的最大值,本文取 一100,当6>1时,将b置为 1。那么两组粒子的融合权重为:a—n/(n+6)。 2.4粒子重采样 粒子滤波器在历经几轮迭代后,大多数粒子权重都变 得很小,即出现粒子退化现象,此时的粒子集将无法准确描 述实际的目标后验概率分布。为解决粒子退化问题,重采 样方法思想在保证粒子采样总数不变的基础下,通过复制 权重较大的粒子,增多了权值大的粒子数量,消除了权值小 的粒子,提高了粒子的多样性。本文对两组粒子都采用此 重采样方法来防止粒子的退化。 2.5本文算法流程 概括起来本文算法主要步骤如下: 步骤1初始化跟踪目标。手动框选初始化跟踪区 域,确定跟踪目标。 步骤2采样粒子。在目标区域周围按照高斯分布采 样两组粒子集P 和P 。 步骤3粒子特征提取。粒子集P 以颜色直方图提 取颜色特征;粒子集P 通过压缩感知提取灰度和纹理 特征。 步骤4粒子权重计算。采用巴特沃斯距离计算P 粒子与目标模板之间颜色直方图的相似度,以确定第是帧 粒子集P 各个粒子的权重 ;以贝叶斯分类器计算P 粒 子的权重,在目标位置周围进行采样正负样本,用第2.1节 中提到的R 和R ,对这些样本进行压缩感知特征提取,得 到特征向量y 和y。,使用式(6)对P 粒子集的特征向量 y】和V。进行分类,以H(V)值大小作为 粒子的权重。 ・2622 ・ 系统工程与电子技术 第37卷 仅比粒子滤波算法慢2~3 ms,这是由于加入的压缩感知 提取的特征数据量极少,整个算法运算速度几乎不受影响。 本文还将CT算法的每帧采样图像片用粒子滤波的每帧重 采样替换,运算量比CT算法降低很多,因此每帧图像的运 算速度比CT算法提高了一倍。 表1算法运行速度比较 ms/fp 4 结 论 针对目标跟踪存在目标遮挡和环境光照变化的情况, 本文提出将改进的CT算法中的压缩感知特征融合入粒子 滤波框架中。由于利用压缩感知提取灰度和纹理特征,使 得在图像信号压缩后数据量大大减少的情况下,可以有效 地缓解光照变化对目标跟踪的影响。而粒子滤波中采用 了颜色直方图特征,以及随机采样和重采样的核心思想, 从而能在目标有遮挡的情况下,对其进行有效的跟踪。因 此,融入压缩感知特征的粒子滤波算法在运算速度并没有 降低的情况下,能同时解决目标遮挡和光照变化两个问 题,从而提高了算法的鲁棒性。未来,如何白适应地调节 判定算法可信度的阈值来使跟踪算法更为稳定是本文的 研究重点。 参考文献: Eli Zhao Q J,Zhang Y,Zhao D B.A scale and rotation adaptive algo— rithm for object tracking[J].Journal of Central South University- Science and Technology,2013,44(6):2354—2360.(赵钦君,张勇, 赵东标.一种尺度和旋转自适应的目标跟踪算法EJ].中南大学学 报(自然科学版),2013,44(6):2354—2360.) E2]Oron S,Bar-Hillel A,Levi D,et a1.Locally orderless tracking[C]// Proc.ofthe n EE Confefence on Computer Vision and Pattern Recog~ nition,2012:1940—1947. [33 Nummiaro K,Koller-Meier E,Van Gool L.An adaptive color based particle filter[J].Image and Vision Computing,2003,21 (1):99—110. [4]Mei X,Ling H.Robust visual tracking and vehicle classification via sparse representation[J]. IEEE Trans.on Pattern Analysis andMachine Intelligence,20U,33(11):2259—22 72. [53 Hess R,Fern A.Discriminatively trained particle filters for complex muhi-obieet tracking[el}{Proc.of the IEEE Con{erence Com purer Vision and Pattern Recognition,2009:240—247. [6]Zhang K,Zhang L,Yang M H.Real-time compressive tracking[c]∥ Proc.ofthe 12th European Conference on Computer Vision,2012:864 877. [7]Zhu Q P,Yan J,Zhang H,et a1.Real—time tracking using mul— tiple features based on compressive sensing[J].Optics and Pre cision Engineering,2013,21(2):437—443.(朱秋平,颜佳,张虎, 等.基于压缩感知的多特征实时跟踪[J].光学精密工程,2013, 21(2):437—443.) [8]Hu S Q,Jing Z L.Overview of particle filter algorithm[J]. Control and Decision,2005,20(4):361—365.(胡士强,敬忠良. 粒子滤波算法综述[J].控制与决策,2005,2O(4):361—365.) [9]Sanjeev A M,Maskell S,Gordon N,et a1.A tutorial on particle filters for online nonlinear/non—Gaussian Bayesian tracking[J]. IEEE Trans.on Signal Processing,2002,50(2):174—188. ElO]Djuric P M,Kotecha J H,Zhang J,et a1.Particle filtering[J]. IEEETrans.onSignal ProcessingMagazine,2003,20(5):19—38. [11]Donoho D.Compressed sensing[J].IEEE Trans.on Informa tion Theory,2006,52(4):1289—1306. [1z]Candes E,Tao T.Near optimal signal recovery from random projections and universal encoding strategies[J].IEEE Trans. on Information Theory,2006,52(12):5406—5425. [13]Candes E,Wakin M.An introduction to compressive sampling[J]. IEEETrans.onSignal ProcessingMagazine,2008,25(2):21—30. [14]Diaconis P,Freedman D.Asymptotics of graphical projection pursuit[J].T^PAnnals of Statistics,1984,12(3):793—81 5. [15]wuY,um J,YangMH.Online objecttracking:a benchmark[(;]∥ Proc.oftheIEEEConfer ̄q ̄e onComputerVision armPattererRec ognition,2013:2411—2418. 作者简介: 吴晓雨(1979一),女,副教授,博士,主要研究方向为图像处理、视频 智能分析。 E—mail:wuxiaoyu@cue.edu.cn 吴凌琳(1990一),通信作者,女,硕士研究生,主要研究方向为图像处 理、视频智能分析。 E—mail:carinel@163.corn 杨磊(1960一),男,教授,主要研究方向为安防领域、视频分析。 E—mail:yanglei@CUC.edu.cn 

因篇幅问题不能全部显示,请点此查看更多更全内容

Copyright © 2019- fenyunshixun.cn 版权所有 湘ICP备2023022495号-9

违法及侵权请联系:TEL:199 18 7713 E-MAIL:2724546146@qq.com

本站由北京市万商天勤律师事务所王兴未律师提供法律服务