均值漂移算法在无人机侦察目标图像跟踪中的研究

Ｅｑｕｉｐｍｅｎｔ　Ｍａｎｕｆａｃｔｕｒｉｎｇ　Ｔｅｃｈｎｏｌｏｇｙ　Ｎｏ．６，２０１　３　均值漂移算法在无人机侦察目标图像跟踪中的研究　杨勇，吴洋　（中国航天空气动力技术研究院，北京１０００７４）　摘要：通过分析均值漂移算法及改进算法在元人机侦察目标图像跟踪中的应用，采用公式法及文献法对其进行研究。　研究结果表明基本的均值漂移算法对比较简单的侦察目标图像跟踪有较多的优点，但对于比较复杂的侦察目标图像跟　踪需要使用改进的均值漂移算法。该研究对无人机侦察目标图像跟踪有参考依据。　关键字：均值漂移算法；核函数；目标图像跟踪；改进算法　中图分类号：ＴＭ９１１．７　文献标识码：Ａ　文章编号：１６７２—５４５Ｘ（２０１３）０６－００５２—０３　均值漂移算法（ｍｅａｎ　ｓｈｉｆｔ）在无人机侦察目标图　中我们规定，如果ｌ　Ｉｅ（ｒｘ）一　ＩＩ＜占，则循环结束；否　像跟踪中应用广泛，可以将要跟踪目标的特征信息　则，把ｍ（ｘ）的值赋予　，继续循环。在上述迭代过程　与其环境信息有效的结合起来，并根据计算、分析得　中，样本戈经过的一系列位置称为样本　的轨迹，可　　，ｍ（ｘ），ｍ（ｍ（ｘ）），…｝。　‘　到的数据，做到及时、准确的对无人机侦察目标图像　表示为｛跟踪进行定位。本文针对比较复杂的侦察目标图像　一通过公式（２）可以看出：均值漂移算法总是沿着　跟踪需要使用改进的均值漂移算法进行研究，以期　概率密度的梯度方向移动，这就是说均值漂移算法　定会收敛到样本　处附近的峰值。　对无人机侦察目标图像跟踪提供参考依据。　１　均值漂移算法　２　基于均值漂移算法的无人机目标图像跟踪　１无人机目标图像跟踪的相关介绍　均值漂移（ｍｅａｎ　ｓｈｉｆｔ）算法最早来源于一篇估计　２．所谓的无人机目标图像跟踪是指通过分析目标　概率密度梯度函数的重要文献，后来Ｙｉｚｏｎｇ　Ｃｈｅｎｇ对　其做了两个方面的推广：第一，定义了一族核函数；第　图像序列中各帧图像的特征，锁定并找出被跟踪的　目标图像。图像的坐标空间（位形空间）我们用　来表　围进一步扩大，被广泛应用到到图像分割、特征空间分　示；在图像中出现的可能是跟踪目标的区域我们用　跟踪窗口来表示；被跟踪的目标用目标模型表示；与　析、模式分类、目标跟踪以及图像平滑等相关领域。　均值漂移算法的计算步骤如下：概率密度函数　之相关的变量用上带短波浪线的字符表示。跟踪窗　二，设定了一个权重系数，使得均值漂移算法的适用范　即满足一定统计约束条件的核函数Ｈ（ｘ），若由密度　Ｉ＝１的信息与目标模型的信息可分别表示为：　，　一－函数ｇ（　）经过　次抽样得到的样本集　）：　，那么得　和　，ｕ　ＪＮ　１，其中　，　∈　；　，ｕ　∈Ｔ；Ｎ表示跟踪窗口　　．　～　到的密度估计值为：　ｇ（　）＝∑　Ｉ－Ｉ（ｘ一　）　（１）　的个数，＾嚷示目标模型中像素点的个数。　由于被跟踪目标的外部颜色及形状往往会比较　其中，　为权重系数，且满足约束条件∑　＝１。　复杂，所以对跟踪目标进行有效的建模变得十分困　如果核函数日（　）是某核函数　（　）的影子集，那么其　难。为有效地对目标进行建模我们引入了基于均值　均值漂移向量为：　ｍ（　）一　：乓　厶ｆ　１－１　一Ｘｉ／　一　（２）　漂移算法的目标跟踪。　２．２基于均值漂移算法的无人机目标图像跟踪的基　本算法　基本的均值漂移算法对所关心的图像区域通常　其中，ｍ（　）是样本　处的样本均值。在上述过程　收稿日期：２０１３—０３—１１　作者简介：杨勇（１９８４一），男，山东青州人，工程师，硕士学历，研究方向为无人机载荷应用与情报处理；吴洋（１９８７一），男，新疆　乌鲁木齐人，工科硕士，助理工程师，研究方向为无人机机载设备与测控系统。　５２　《装备制造技术）２０１３年第６期　使用统计直方图进行建模，并且对应图像区域内的　概率密度函数也可以近似的用直方图表示。用；表示　目标模型的直方图；用Ｊ；（　）表示另外一个变量的直方　图。如果用　（　）＝ｐ（　（　），　）作为描述以上两个变量　的相似性函数，那么使　（　）达到最大值（局部最大　Ⅳ０。＝∑∑ｘｌ（ｘ，　），ＮＩ。＝∑∑ｚ，（　）（７）　＾，　＾，　计算出质心：　＝　。上，），ｃ＝　』Ｖ００　』Ｖ００　（３）对跟搜索窗口的大小重新设置，使之满足条　ｒ　值）的观测窗口的坐标就是我们所要求的目标位置。　为方便计算，通常对搜集到的图像特征进行ｍ级量　件：ｓ：ｓ＝２、／Ｖ　３０　　。　（４）循环过程（２）和（３），直到质心的变化小于已　化，其中　和　的量化结果分别用ｂ　和ｂ　表示。　由于初始的直方图没能抓住自身在位形空间中　的分布情况，因此，需要用统计直方图的特征点与目　标中心之间的距离来对初始特征进行修正：　ｇ　＝　∑　（　）６（ｂ—　）　（３）　Ｐ　＝Ｃ。∑　Ｈｏ＇（ｘ　一　）６（６一　）　（４）　其中，　和　为归一化常数；Ｏｒ代表窗宽，具体　表示跟踪窗口的大小，经过修正的直方图信息则比　较全面，即包含了其特征在位形空间中的分布情况，　又包含了自身及目标的特征信息。在相似性方面使　用巴氏系数，即使用公式　（ｚ）＝∑　、／　（　）　进行　计算，　（　）进行最大值求解，我们就可以得出被跟　踪目标位置的计算公式（迭代算式）：　＝　器　（５）　其中，系数　＝∑　Ｖ　（　）　（６　一ｕ）。　２．３均值漂移算法的改进　２．３．１自动调整跟踪窗口尺度的改进算法　由于在跟踪过程中，目标的尺寸会发生变化，所　以，自动调整跟踪窗口的尺度显得尤为重要，跟踪窗　口过大或者过小都有可能导致跟踪失败。基于以上　原因，我们可以根据目标区域内颜色直方图或者灰　度质心的计算结果，并依据窗口尺寸与质心已有的　数量关系，自动地调整跟踪窗口的尺度。具体的改进　过程如下：　（１）在给定的、初始的目标图像中，对目标的大　小和位置进行人为的初始化，以固定搜索窗口　的　尺寸大小（其中用ｓ来表示搜索窗口的大小）。　（２）令搜索窗口　中任一坐标（　，　）的像素颜　色概率分布函数为Ｉ（ｘ，　），并计算零阶矩阵：　＝∑∑如，　）　（６）　继续计算　的一阶矩阵和　的一阶矩阵：　经设置的迭代次数或阈值，停止计算。　２．３．２有效解决目标遮挡问题的改进算法　当被跟踪目标出现绝大部分或者全部图像被遮　挡的情况时，当前帧中的目标图像已经找不到被跟　踪的目标区域，但是按照要求，基本的均值漂移算法　必须找到目标区域，由于自身的跟踪能力有限，当前　找到的目标区域的特征已经发生了较大的变化，在　遮挡结束以后，基本的均值漂移算法仍然把已经发　生较大变化的目标特征作为跟踪的目标，从而造成　跟踪失败。跟踪失败根本的原因是基本的均值漂移　算法只有对确定的数据进行估计的能力，无法解决　概率密度分布函数出现多个峰值的情况，可通过加　入粒子滤波算法来解决上述问题。　粒子滤波跟踪算法的主要过程为：　（１）随机选取一组初始粒子的权值和位置，当被　跟踪目标的位置确定后，便可通过粒子集得到被跟　踪目标的初始化颜色直方图　（２）在目标图像的变化过程中，粒子按照预先设　定的经验传播，借此估计被跟踪目标的下一个状态。　尽可能多的使用之前所有的观测值是一个好的传播　模型所必须具有特点，二阶自回归模型就具有这一　特点：　；　一　＝　一　ｉ＋ｌ一　＝　ｆ一　ｌ＋　＋’’１￣　ｉ（８）　其中，　表示粒子在ｉ时刻所具有的状态；Ⅳｉ代　表随机产生的高斯噪音。从（８）中可以明显的看出：　假设目标始终保持相同的运动方向和运动速度。　（３）对比每个粒子此时和初始值的颜色直方图，　把得到的二者的相似度赋予每个粒子，并以此作为　每个粒子的权值。把传播后粒子和参考粒子的颜色　直方图（分别用ｊ；　（　）和　来表示）进行对比，并利用　Ｂｈａｔｔａｅｈａｒｙｙａ系数来测量二者的相似度：　！ｒ—————一　（　）＝ｐ（　（　），；）＝∑、／＝ｌ　　（ｚ）；　（９）　假设Ｂｈａｔｔａｃｈａｒｙｙａ系数与粒子的相似度成正比　５３　Ｅｑｕｉｐｍｅｎｔ　Ｍａｎｕｆａｃｔｕｒｉｎｇ　Ｔｅｃｈｎｏｌｏｇｙ　Ｎｏ．６，２０１　３　例关系，对于粒子赋予的权重，相似度越大赋予的权　当重要，若阈值选取的较小，则会导致产生较大的误　重就越大；相似度越小赋予的权重就越小。　差；如若阈值选取的较小，则很难收敛，在实际的应　（４）筛选经过上述传播过程之后得到的新的粒　用中，阈值一般预设在０．８—０．９之间。在第（６）步中，　子，为避免把计算时间浪费在无用的粒子上，最大可　通常只用简化的粒子滤波算法。　能的保留权重大的粒子，丢掉权重较小的粒子。　通过多次的实验证明：此改进算法（有效解决目　能够有效地解决最复杂的　为保证能在更大的区域内搜索被跟踪目标，粒　标遮挡问题的改进算法），子滤波算法通过增加预测目标的数目来实现，与此　多目标互相遮挡的情况，实现对目标的实时、准确地　　同时，根据相似度更新粒子，使得跟踪算法不断收　跟踪。敛，最终收敛到与被跟踪目标最相似的粒子上。因　次，粒子滤波算法可以在发生遮挡时和发生遮挡后，　３结束语　重新在目标图像中搜索被跟踪目标，大大提高了遮　挡结束后成功跟踪目标的可能性。所以，在均值漂移　运用均值漂移算法对无人机目标图像跟踪，其　算法中加入粒子滤波算法，可解决目标被遮挡后出　算法简单、快捷、易操作、实用功能强大，能够对简单　现的目标丢失问题。　环境中的视频目标进行实时跟踪，而对于较复杂的　２．３．３改进算法的计算步骤　情况下（如目标的尺寸不断发生变化或者目标的绝　（１）选择被跟踪的目标区域，创建被跟踪的目标　大部分或全部被遮挡）的无人机目标图像跟踪则需　要利用改进的算法。目前，很多的专家、学者正进一　模型；　；　（２）在图像不断更新的过程中，计算预测目标区　域的直方图　，　（ｚ）；　步研究该算法在无人机目标图像发生较大形变或者　旋转的情况下如何进行有效地跟踪等课题，期待他　们早Ｅｔ取得丰厚的科研成果。　（３）通过使用均值漂移算法中的迭代计算方法　计算（预测）当前帧中的跟踪目标位置；　参考文献：　（４）利用公式来计算当前帧中跟踪目标模型与　［１】李乡儒，吴福朝，胡占义．均值漂移算法的收敛性叨．软件学　跟踪目标区域的相似度系数；　报，２００９，ｌ６（３）：３６５．　２】卞丽．基于改进的ＩＭＭ目标跟踪算法研究卞丽叽．科技通　（５）若（４）中的相似度系数小于预设的阈值，则　【报，２０１２（０４）：４９—５１．　转到第（６）步；否则，转到第（７）步；　（６）启用粒子滤波目标跟踪算法；　［３］Ｃｏｍａｎｉｃｉｕ　Ｄ，Ｒａｍｅｓｈ　Ｖ，Ｍｅｅｒ　Ｐ．Ｋｅｒｎｅｌ－ｂａｓｅｄ　Ｏｂｊｅｃｔ　Ｔｒａｃｋ－　ｉｎｇ［￣．ＩＥＥＥ　Ｔｒａｎｓ．ｏｎ　ＰａＲｅｒｎ　Ａｎａ１．Ｍａｃｈｉｎｅ　ＩｎｔｅＵ．，２００９，２５　（７）把计算出的结果用作下一帧的目标模型，然　（５）：５６４．　后返回到上述过程中的第（２）步。　［４］彭宁嵩，杨杰，刘志，等．Ｍｅａｎ—ｓｈｉｆｔ跟踪算法中核函数　在上述的过程中，第（５）步中预设阈值的选取相　窗宽的自动选取叨．软件学报，２０１０，１６（９）：１５４２．　Ａｖｅｒａｇｅ　Ｖａｌｕｅ　Ｄｒｉｆｔｉｎｇ　Ａｌｇｏｒｉｔｈｍ　ｉｎ　Ｕｎｍａｎｎｅｄ　Ａｅｒｉａｌ　Ｖｅｈｉｃｌｅ　Ｒｅｃｏｎｎａｉｓｓａｎｃｅ　Ｏｂｊｅｃｔｉｖｅ　Ｉｍａｇｅ　Ｔｒａｃｋ　Ｒｅｓｅａｒｃｈ　ＹＡＮＧ　Ｙｏｎｇ．ＷＵ　Ｙａｎｇ　（Ａｅｒｏｓｐａｃｅ　Ｃｈｉｎａ　Ａｅｒｏｄｙｎａｍｉｃ　Ｆｏｒｃｅ　Ｔｅｃｈｎｉｃａｌ　Ｉｎｓｔｉｔｕｔｅ，Ｂｅｉｊｉｎｇ　１０００７４，Ｃｈｉｎａ）　Ａｂｓｔｒａｃｔ：Ｔｈｒｏｕｇｈ　ｔｈｅ　ａｎａｌｙｓｉｓ　ｏｆ　ｍｅａｎ　ｓｈｉｆｔ　ａｌｇｏｒｉｔｈｍ　ａｎｄ　ａｐｐｌｉｃａｔｉｏｎ　ｏｆ　ｉｍｐｒｏｖｅｄ　ａｌｇｏｒｉｔｈｍ　ｆｏｒ　ＵＡＶ　ｒｅｃｏｎｎａｉｓｓａｎｃｅ　ｉｍａｇｅ　ｔｒａｃｋｉｎｇ，ｓｔｕｄｉｅｄ　ｕｓｉｎｇ　ｆｏｒｍｕｌａｓ　ａｎｄ　ｌｉｔｅｒａｔｕｒｅ　ｌａｗ；ｓｔｕｄｉｅｓ　ｈａｖｅ　ｓｈｏｗｎ　ｂａｓｉｃ　ｃｏｍｐａｒｉｓｏｎ　ｏｆ　ｔｈｅ　ｍｅａｎ　ｓｈｉｆｔ　ａｌｇｏｒｉｔｈｍ　ｆｏｒ　ｓｉｍｐｌｅ　ｒｅｃｏｎｎａｉｓｓａｎｃｅ　ｔｒａｃｋｅｄ　ｔｈｅ　ｔａｒｇｅｔ　ｉｍａｇｅ　ｈａｓ　ｍｏｒｅ　ａｄｖａｎｔｇｅｓ，ｂｕｔａ　ｏｒ　ｆｍｏｒｅ　ｃｏｍｐｌｅｘ　ｒｅｃｏｎｎａｉｓｓａｎｃｅ　ａｎｄ　ｔｒｇｅｔａ　ｔｒａｃｋｉｎｇ　ｕｓｉｎｇ　ｉｍｐｒｏｖｅｄ　ｍｅａｌｌ　ｓｈｉｈ　ａｌｇｏｒｉｔｈｍ　ｆｏｒ　ｉｍａｇｅ．Ｔｈｅ　ｓｔｕｄｙ　ｏｎ　ＵＡ￣ｒｅｃｏｎｎａｉｓｓａｎｃｅ　ｉｍａｇｅ　ｔｒａｃｋｉｎｇ　ｒｅｆｅｒｅｎｃｅ．　Ｋｅｙ　ｗｏｒｄｓ．＂ｍｅａｎ　ｓｈｉｔｆ　ｌｇｏｒａｉｔｈｍ；ｆｏｒ　ｋｅｒｎｅｌ　ｆｕｎｃｔｉｏｎ；ｉｍａｇｅ　ｒｔａｃｋｉｎｇ；ｉｍｐｒｏｖｅｄ　ｌｇｏｒａｉｔｈｍ