一种基于旋翼无人机平台的目标跟踪方法
【专利摘要】本发明公开了一种基于旋翼无人机平台的目标跟踪方法。本发明目标跟踪方法通过多尺度的样本采集并结合实时更新的分类器，有效实现了旋翼无人机平台对移动目标的快速准确地在线实时跟踪；本发明在跟踪过程中，根据当前帧的最大分类器响应值以及其相比前一帧的最大分类器响应值的变化情况，并结合前一帧的跟踪性能判定结果，来判定当前帧的跟踪性能是否稳定，当跟踪性能不稳定时，及时对分类器的分类结果进行修正，可有效防止目标跟踪过程中由于遮挡而导致的跟踪目标丢失。相比现有技术，本发明有效实现了旋翼无人机平台对移动目标的快速准确地在线实时跟踪，且具备多尺度跟踪特性，跟踪性能更好，算法复杂度更低，实时性更好，对硬件资源的要求较低。
【专利说明】
-种基于旋翼无人机平台的目标跟踪方法
技术领域
[0001] 本发明设及图像处理及模式识别技术领域，尤其设及一种基于旋翼无人机平台的 目柄跟踪方法。
【背景技术】
[0002] 近年来，随着科学技术水平的日益提高，无人机事业的发展蒸蒸日上。由于其本身 具有结构简单、成本较低、灵活性高、环境适应性强、便于装载各类传感器及装备等特点，加 之数字通信、传感器、微电子、虚拟现实等技术的快速进步，无人机在空中作战、城市武装反 恐、海上捜救等军用领域W及民用领域中都具有广阔的应用前景和发展空间，而利用无人 机所拍摄的视频对目标实现跟踪则是其中一项重要应用。
[0003] 根据背景的不同，运动目标的识别跟踪可W分为静态背景和动态背景下的目标识 别跟踪。静态背景下运动目标识别是指摄像机在整个过程中静止;而动态背景下运动目标 识别是指摄像机在识别过程中发生了诸如平动、旋转或多自由度运动等移动的情况。显然， 基于无人机运样旋翼无人机平台的目标跟踪属于上述的第二类范畴。
[0004] 基于计算机视觉的目标跟踪算法主要分为四类，分别是基于目标区域的跟踪、基 于目标特征的跟踪、基于目标变形模板的跟踪和基于目标模型的跟踪。基于区域的跟踪算 法首先通过图像分割或先验知识进行运动区域的预先提取，并得到相应的模板，然后运用 相关算法对序列图像中的目标进行跟踪。相对其它跟踪算法来说，由基于区域的跟踪算法 提取出的目标模板比较完整，因而得到的图像信息更为丰富，可稳定跟踪目标，具有较高的 跟踪精度。但它也存在一些缺点，如计算量大、对图像变形问题模板匹配困难、目标遮挡或 形变时效果不佳、复杂环境下模板漂移等问题。基于目标特征的跟踪算法利用匹配算法在 图像序列中寻找与目标特征相符的运动物体。其主要步骤为特征提取、特征匹配W及运动 信息计算。该方法对遮挡、照明、视角等问题造成的图像变化表现出很好的鲁棒性。在W往 的文献中，学者将其与Adaboost、卡尔曼滤波等预测算法进行结合，取得了令人满意的跟踪 效果。如何选择合适的特征是该算法应用过程中的关键点和难点。变形模板是纹理或边缘 可W按照一定的条件变形的面板或曲线。基于目标变形模板的跟踪算法对可变性目标十分 有效。基于目标模型的跟踪算法中，运动物体的表达方式分为=种：线图模型、2D模型和3D 模型。在实际情况中，采用3D模型的跟踪具有更广的应用范围。一般来说，算法首先根据先 验知识获得目标的立体结构模型和运动模型，然后结合实际场景对目标物体进行跟踪。基 于3D模型的目标跟踪算法性能可靠，及时目标运动状态发生变化也能取得很好的效果。但 在实际应用中很难获得运动目标的精确几何模型，且运算量大，算法实时性差。
[0005] 无人机搭载摄像头后，运动目标检测是基于动态背景的，虽然摄像头的视野得到 扩大，但是目标和摄像机之间存在复杂的相对运动，因而对目标跟踪算法提出了更高的要 求。
【发明内容】
[0006] 本发明所要解决的技术问题在于克服现有技术不足，提供一种基于旋翼无人机平 台的目标跟踪方法，可快速准确地实现对移动目标的在线跟踪，且可有效防止目标被遮挡 导致跟踪目标丢失的情况。
[0007] 本发明具体采用W下技术方案解决上述技术问题：
[0008] -种基于旋翼无人机平台的目标跟踪方法，包括W下步骤：
[0009] Sl、对旋翼无人机平台所拍摄的视频的每一帖进行自动目标识别，W初次识别出 目标的视频帖作为目标跟踪的初始帖，并在初始帖中目标周围划定矩形的初始目标跟踪 框；
[0010] S2、W初始目标跟踪框为基准，在初始帖中采集一组与初始目标跟踪框同尺度的 正、负样本，提取所述样本的尺度不变特征对分类器进行训练，并W初始目标跟踪框中的图 像作为正样本模板；
[0011] S3、在当前帖中，W前一帖的目标跟踪框位置为中屯、的一定区域内采集一组不同 尺度的矩形区域作为检测样本，并在每个检测样本上叠加均值为0的高斯白噪声;提取叠加 噪声后的各检测样本的尺度不变特征并分别作为所述分类器的输入，得到各检测样本的分 类器响应值；
[0012] S4、判断当前帖的跟踪性能是否稳定，如性能不稳定，则转至S5,否则，转至S8;如 W下条件之一得到满足，则当前帖的跟踪性能不稳定：
[0013] (1)当前帖检测样本的最大分类器响应值小于0;
[0014] (2)当前帖检测样本的最大分类器响应值大于等于0,同时当前帖检测样本的最大 分类器响应值小于前一帖检测样本的最大分类器响应值，且前一帖的跟踪性能不稳定；
[0015] (3)当前帖检测样本的最大分类器响应值在0和一个预设的大于0的阔值之间，同 时当前帖检测样本的最大分类器响应值小于前一帖检测样本的最大分类器响应值，且前一 帖的跟踪性能稳定；
[0016] S5、计算当前帖各检测样本与正样本模板之间的归一化相似度，并W归一化相似 度作为修正系数对当前帖相应检测样本的分类器响应值进行修正；W修正后的分类器响应 值最大的检测样本所对应的矩形区域作为当前帖的目标跟踪框；
[0017] S6、判断连续跟踪不稳定的帖数是否达至顺设上限值，如是，则将所述高斯白噪声 的标准差增大后转SlO;否则，转至S7;
[0018] S7、W当前帖的目标跟踪框为基准，在当前帖中采集一组与当前帖的目标跟踪框 同尺度的正、负样本，从所有正样本中选出与正样本模板之间相似度最大的一部分;然后提 取所选出的正样本W及负样本的尺度不变特征并对分类器进行重新训练后转至S10;
[0019] S8、W分类器响应值最大的检测样本所对应的矩形区域作为当前帖的目标跟踪 框，并用前一帖正样本模板与初始帖正样本模板的加权和来更新正样本模板；
[0020] S9、W当前帖的目标跟踪框为基准，在当前帖中采集一组与当前帖的目标跟踪框 同尺度的正、负样本，提取所述正、负样本的尺度不变特征并对分类器进行重新训练；
[0021] S10、如视频已结束，则退出目标跟踪;否则，W下一帖作为当前帖，转至S3。
[0022] 优选地，所述分类器为朴素贝叶斯分类器。
[0023] 优选地，所述尺度不变特征为尺度不变压缩特征。
[0024] 进一步地，所述分类器的学习速率A在目标跟踪过程中按照下式进行自适应更新：
[0025]
[0026] 式中，X为当前帖的目标跟踪框中图像高斯分布均值与前一帖的目标跟踪框中图 像高斯分布均值之间差值的绝对值；
[0027] 优选地，检测样本与正样本模板之间的相似度使用W下方式度量:检测样本与正 样本模板的各维特征值的差值绝对值之和。
[0028] 优选地，利用W下方法对旋翼无人机平台所拍摄的视频的每一帖进行自动目标识 另IJ:首先提取当前帖中的感兴趣区域;然后将感兴趣区域与目标模板先后进行主颜色匹配 和轮廓匹配，如存在颜色匹配和轮廓匹配均匹配成功的感兴趣区域，则该感兴趣区域即为 在当前帖中识别出的目标，否则，继续对下一帖进行目标识别。
[0029] 相比现有技术，本发明技术方案具有W下有益效果：
[0030] 本发明有效实现了旋翼无人机平台对移动目标的快速准确地在线实时跟踪，且具 备多尺度跟踪特性；
[0031 ]本发明可有效防止目标跟踪过程中由于遮挡而导致的跟踪目标丢失，跟踪效果更 好；
[0032] 本发明算法复杂度低，实时性好，对硬件资源的要求较低。
【附图说明】
[0033] 图1是本发明目标跟踪方法流程图；
[0034] 图2是【具体实施方式】中所使用目标识别方法的流程图；
[0035] 图3是选取多尺度待检测样本方案的示意图；
[0036] 图4是尺度不变压缩特征向量提取示意图；
[0037] 图5是分类器更新阶段正负样本采集区域示意图；
[0038] 图6是学习速率自适应更新函数示意图。
【具体实施方式】
[0039] 下面结合附图对本发明的技术方案进行详细说明：
[0040] 本发明基于旋翼无人机平台的目标跟踪方法的基本流程如图1所示，具体包括W 下步骤：
[0041 ] Sl、对旋翼无人机平台所拍摄的视频的每一帖进行自动目标识别，W初次识别出 目标的视频帖作为目标跟踪的初始帖，并在初始帖中目标周围划定矩形的初始目标跟踪 框；
[0042]本发明可采用现有的各种自动目标识别方法，例如基于轮廓的目标识别方法、基 于主颜色的目标识别方法等，为了提高目标识别的准确性，本发明将基于轮廓的目标识别 方法与基于主颜色的目标识别方法相结合，即：首先提取当前帖中的感兴趣区域;然后将感 兴趣区域与目标模板先后进行主颜色匹配和轮廓匹配，如存在颜色匹配和轮廓匹配均匹配 成功的感兴趣区域，则该感兴趣区域即为在当前帖中识别出的目标，否则，继续对下一帖进 行目标识别。如图2所示，目标识别的具体步骤如下：
[0043] (1)图像灰度化
[0044] 本文中，使用加权平均法将彩色图像转换成灰度图像，具体公式为：
[0045] 1 = 0.3XR+0.59XG+0.IlXB
[0046] 其中R，G，B为彩色图像中对应的红、绿、蓝=个通道的灰度分量，I为转换后的单通 道灰度图像。
[0047] (2)图像二值化，用最大亮度值1来表示前景，用最小亮度值0来表示背景 [004引（3)提取轮廓
[0049] 物体形状的边缘轮廓借助边缘检测进行提取，然后对物体形状的边缘轮廓进行跟 踪，并得到边界点的相应坐标，按照顺时针或者逆时针的顺序，用n个像素点坐标CO，Cl... Cn 来描述运些边界点的位置，W此得到物体的形状。本文默认C日=Cn,即目标物体的轮廓是封 闭的。
[0050] (4)使当前帖的原彩色图像仅保留彩色前景部分；
[0051] (5)从前景中提取主颜色：
[0052] 使用K均值聚类算法提取主颜色：
[0053] (a)将图像的颜色空间由RGB转化到服V;
[0054] (b)确定主颜色数目k和初始聚类中屯、Cl;
[0055] (C)循环步骤(cO到(e)，直到每个聚类中屯、都不再发生变化，转至(f);
[0056] (d)计算每个像素点与聚类中屯、。的距离，然后将其归进相似度最高的聚类中。若 用Ci=化i, Si ,Vi)，Cj =化j ,Sj ,Vj)分别表示HSV空间中的两种颜色，则二者的相似性计算公 式如下：
[0化7]
[0化引其中34在0至1之间取值，服V空间是单位高度的圆柱体。
[0059] (e)重新计算每个类的聚类中屯、，其数值取该聚类中所有像素的平均值；
[0060] (f)得到图像的主颜色直方图Ci W及对应的主颜色频率Pi;
[0061] (6)对当前帖和模板图像进行主颜色匹配；
[0062] 假设待检测图像对应的前景图像主色集合为主色频率为Pq,主色共M个;模板的 前景图像主色集合为Cl，主色频率为Pl，主色共N个。构造主色集，使两个直方图的颜色数量 和颜色集合都相同，C = CQUCL，主色共有M+N个，颜色直方图为P'Q=(Pql，Pq2，...，Pqm，0， 0, . . .0)，P'L=(0,0, . . .0，pii，pi2, . . .pin)。将主色集合的主色由RGB空间转至HSV空间，并 计算C中任意两种颜色Ci，Cj的相似性Sij(i，j = 1，2，. . .M+N)，并W
[0063] 此构造颜色集合的相似矩阵S，两个任意主色集之间的距离为：
[0064] d2=(PVP'L)S(PVP'L)T
[0065] (7)对当前帖和模板图像进行轮廓匹配。采用化矩匹配度量两个轮廓的相似
[0066] 度，化矩计算公式如下：
[0067] hi=ri2〇 巧 02
[006引 h2= (ri20-n02)2+化 A
[0069] h3= (ri3〇-3rU2)2+(3ri2i-n〇3)2
[0070] h4=(ri30 巧 12)2+(1121 巧 03)2
[007。 hs= (ri30-3rU2) (1130巧12) [ (ri30巧i2)2-3(ri2i-n03)2]
[007^ +(31121-11。3)(1121巧。3)[3(化21巧。3)2-(1121巧03)2]
[0073] h6=(n2〇-n〇2)[(ri30 巧 12)2-(1121 巧 03)2]+化 11(1130 巧 12)(1121 巧 03)
[0074] h?=(化21-11日3 )(物巧日3 ) [ 3 (ri3日巧22 ) 2- (1121 巧03 ) 2 ]
[0075] -(ri3〇-ni2) (n2i+n〇3) [ 3 (n30+ni2) (n2i+n03) ^ ]
[0076] 其中rip,q为图像的归一化矩。相似度指标采用如下公式：
[0077]
[007引 "和皆分别是A和B的Hu矩。
[0079] 当主颜色和轮廓均匹配时，该感兴趣区域即为当前帖中识别出的目标，否则，继续 进行下一帖的自动目标识别。
[0080] S2、W初始目标跟踪框为基准，在初始帖中采集一组与初始目标跟踪框同尺度的 正、负样本，提取所述样本的尺度不变特征对分类器进行训练，并W初始目标跟踪框中的图 像作为正样本模板；
[0081 ]所述尺度不变特征可W是SIFT特征、SURF特征、BRISK特征等，为了降低算法复杂 度，提高目标跟踪实时性，本发明优选采用尺度不变压缩特征，并优选朴素贝叶斯分类器。 具体地，对于每一个样本图像ZGlfT，它在压缩域中的压缩向量表达为1' = 0，,，'''，vf 6?"，其 中n?m。压缩跟踪算法中假设低维特征向量V中的所有元素都是独立分布的，算法采用一个 朴素贝叶斯分类器对它们进行建模：
[0082]
[0083] 其中y是表达样本标签的二值变量，只有0和1两个取值。我们假设正样本和负样本 的先验概率是相等的，即P(y = 1) =P(y=0)。
[0084] W往的理论中指出，由随机稀疏投影矩阵得到的稀疏投影分布几乎总是可W看作 是高斯的。所W，朴素贝叶斯分类器中的两个样本分布P(ViIy = I)和P(ViIy = O)可W看作高 斯分布，并且它们的分布特征可由四个参数如，武片,。，《)表征：
[0085] Q
[0086] S3、在当前帖中，W前一帖的目标跟踪框位置为中屯、的一定区域内采集一组不同 尺度的矩形区域作为检测样本，并在每个检测样本上叠加均值为0的高斯白噪声；
[0087] 参照图3,在上一帖最佳位置附近选择待检测样本时，将上一帖跟踪框的位置在原 有基础上分别向上下左右4个方向增大或缩小1/10的长度，并在每个样本上叠加一定的高 斯白噪声noise,则可W得到81个待检测样本，定义S为当前帖跟踪框长度与上一帖最佳跟 踪框宽之比，S=Wp/w。
[008引 wp=w±w/10+noisew
[0089] hp = h±h/10+noiseh
[0090] Wp，hp分别为当前帖待检测矩形样本的宽和高，w，h是上一帖跟踪框位置的宽和高， noisew,noiseh均为人工叠加的均值为O的局斯白噪声，初始状态下令二者标准差为Stdw = 5, Stdh= 2.5。
[0091] 提取叠加噪声后的各检测样本的尺度不变压缩特征向量并分别作为所述分类器 的输入，得到各检测样本的分类器响应值；
[0092] 参照图4,尺度不变压缩特征向量提取的具体步骤如下：
[0093] (1)保持初始随机测量矩阵Ro中所有非零元数值不变，非零元对应的矩形参数fx (i，t)，fy(i，t)，fw(i，t)，扣(i，t)分别与S相乘，并四舍五入取整：
[0094]
[0095]
[0096]
[0097]
[0098] 其中1 = 1，2-，1，*=1，2-，9，1为压缩特征向量维数，9为初始随机测量矩阵第1行 非零元素的个数；
[0099] (2)计算尺度不变压缩特征向量第i维的特征值：
[0100]
[0101]其中，P(i，t)是初始随机测量矩阵Ro中非零元素的值，PsumQ，t)是修正后的随机 测量矩阵化中第i行第t列个非零元素所对应的矩形内像素值的总和，其计算公式如下：
[01 02] Psum( i , t) = I (Xmin , Ymin) +1 (Xmax , Ymax) ~1 (Xmax , Ymin) ~1 (Xmin , Ymax)
[0103] Xmin = y+fys(i ,t)-2
[0104] Xmax = y+fys(i ,t)+fhs(i ,t)-2
[0105] Ymin = X+fXs(i ,t)-2
[0106] Ymax = X+fXs(i ,t)+fWs(i ,t)-2
[0107] (3)最终可W得到M维尺度不变压缩特征向量V。
[0108] S4、判断当前帖的跟踪性能是否稳定，如性能不稳定，则转至S5,否则，转至S8;
[0109] 本发明根据当前帖的最大分类器响应值W及其相比前一帖的最大分类器响应值 的变化情况，并结合前一帖的跟踪性能判定结果，来判定当前帖的跟踪性能是否稳定;具体 为出现W下任一种情况时，判定当前帖的跟踪性能不稳定：
[0110] (1)当前帖检测样本的最大分类器响应值小于0;
[0111] (2)当前帖检测样本的最大分类器响应值大于等于0，同时当前帖检测样本的最大 分类器响应值小于前一帖检测样本的最大分类器响应值，且前一帖的跟踪性能不稳定；
[0112] (3)当前帖检测样本的最大分类器响应值在0和一个预设的大于0的阔值之间，同 时当前帖检测样本的最大分类器响应值小于前一帖检测样本的最大分类器响应值，且前一 帖的跟踪性能稳定。
[0113] 用Pt表示当前帖的跟踪性能，Pt = O表示跟踪性能稳定，Pt=I表示跟踪性能不稳 定。Pt值的确定方法如下:假设Ct表示当前的第t帖中最优检测样本的分类器响应值(即最大 分类器响应值），且A Ct = Ct-Ct-ID为Ct分别设定两个阔值th和ti，且A Ct的值被限定在一个 常数丫内。则我们判定如下几种情况为跟踪性能不稳定(Pt=I):
[0114] (l)Ct<0:当前帖所有样本的分类器响应值都小于0,即分类器得到的目标样本其 实是负样本，运种情况下目标可能遭遇遮挡或者跟踪框已经漂移甚至丢失目标；
[011引（2)Ct〉0，A Ct<0 ,Pt-I = 1:在上一帖跟踪性能不稳定的情况下，若当前帖最佳样本 的分类器响应值大于0,但正样本特征的质量比起上一帖有所下降，此时我们仍然认为跟踪 性能不稳定；
[0116] (3)0<Ct<ti，ACt<0，Pt-i = 0:在上一帖跟踪性能稳定的情况下，若当前帖最佳样本 响应值小于阔值，且比上一帖的跟踪性能差，则认为跟踪性能不稳定。S5、计算当前帖各检 测样本与正样本模板之间的归一化相似度，并W归一化相似度作为修正系数对当前帖相应 检测样本的分类器响应值进行修正；W修正后的分类器响应值最大的检测样本所对应的矩 形区域作为当前帖的目标跟踪框；
[0117] 检测样本与正样本模板之间的相似性度量可采用欧式距离、汉明距离、相关系数 等度量方式，为了简化计算复杂度，本【具体实施方式】中采用绝对差(AD)来衡量两个向量之 间的相似度：
[011 引
[0119] 其中，Vci为检测样本特征向量第i维特征值，Vp;为正样本模板特征向量第i维特征 值。对于第j个样本，得到一个样本差值特征向量A V后对其各维特征值进行求和，结果为一 个样本与正样本模板之间的特征差值绝对值之和A Vsumj。对所有待检测样本完成上述计算 后，将所有结果归一化，并与朴素贝叶斯分类器结果相乘，即可得到所有样本的修正后的分 类器响应值向量。
[0120]
[0121]
[0122] C(V)=H(V)L(V)
[0123] 其中，y是表达样本标签的二值变量，只有0和1两个取值。我们假设正样本和负样 本具有相同的先验概率，即P(y = l)=P(y = 〇) emax( A Vsimj)为所有待检测样本对应的A Vsumj中的最大特征值项。当Pt = O，即跟踪性能稳定时，算法中令正样本相似度函数为全1项， 从而跳过上述相似度衡量步骤，减小计算量。
[0124] S6、判断连续跟踪不稳定的帖数是否达到预设上限值，如是，则将所述高斯白噪声 的标准差增大后转SlO;否则，转至S7;
[0125] 本【具体实施方式】中设定:如果跟踪不稳定帖数大于等于5，则将高斯白噪声的标准 差扩大至初始值的2倍W扩大样本框捜索范围。
[0126] S7、W当前帖的目标跟踪框为基准，在当前帖中采集一组与当前帖的目标跟踪框 同尺度的正、负样本，从所有正样本中选出与正样本模板之间相似度最大的一部分;然后提 取所选出的正样本W及负样本的尺度不变特征并对分类器进行重新训练后转至S10;
[0127] 参照图5,得到当前帖最佳样本后，W目标框中屯、为基准，在附近的一定半径内利 用滑窗法采集与之大小一致的正样本和负样本
[012引 D?s= 口 I I I I(Z)-Io I I <p0srad}
[0129] 〇NEG= {z I p〇srad+4< I I I(Z)-Io I I <negrad}
[0130] W区间（0，posrad)为半径采集正样本，W区间（posrad+4 ,negrad)为半径采集负 样本，本【具体实施方式】中posrad的值设为4，negrad的值设为30,并在跟踪过程中保持不变。 二者的值与原始的压缩跟踪算法中相同，在实际应用中，可W根据目标尺寸、运动速度、图 像分辨率等因素对posrad和negrad的数值进行调整，W取得不同情况下较为理想的跟踪效 果。
[0131] 用滑窗法采集正负样本更新分类器相关步骤如下：
[0132] rowmin=max(0,y-posrad)；
[0133] roWmax=min(row-1 ,y+posrad)；
[0134] colmin=max(0,x-posrad)；
[0135] colmax=min(col_l ,x+posrad);
[0136] 其中，；r〇W = ；r〇Wimage-h-l ; col = Colimage-W-I ;r〇Wimage 是当前帖图像一列的像素数， CO 1 image是当前帖图像一行的像素数，X，y，W，h分别是当前帖得分最高的跟踪框所对应的左 上顶点横纵坐标、宽和高。posrad为正样本采集半径。
[0137] W样本的左上顶点为中屯、，行捜索Wrowmin为起点，rowmax为结束点，列捜索W COlmin为起点，COlmax为结束点，将最高得分跟踪框窗日依次滑动过上述像素所在位置，选取 正负样本。
[0138] 在此过程中4
NUMmax为算法所需要的最大样本 数量，本算法中采集正样本时将其设为1000,采集负样本时将其设为100。若当前帖图像所 能提供的样本数量能够满足算法需要的样本数量，则此时b<l;若当前帖图像所能提供的样 本数量小于算法需要的样本数量，则此时b> 1。设置一个随机数发生器，可W随机产生一个0 至1之间的数C。本【具体实施方式】中，当b〉c且dist<posrad2时采集正样本，当b〉c且dist〉 negrad2时采集负样本，运里posrad为正样本采集半径，negrad为负样本采集半径，本具体 实施方式中分别定为4和30。
[0139] 对于采集到的正样本，将其依次与当前正样本模板进行比较，计算低维特征向量 相似度，选择绝对差之和最小的前1/2更新朴素贝叶斯分类器：
O[0142] 如图6所示，本【具体实施方式】中还采用自适应的方式更新学习速率入：
[0140] ，
[0141]
[0143]
[0144] 其中，^=|4-化|，4为当前帖目标跟踪框中图像高斯分布均值，化为前帖目标跟踪 框中图像高斯分布均值；当x<t时，r = Ti，x〉t时，r = K JiJU、t均为预设参数，本具体实施 方式中 t = 0.4, ri = 0.4, K = 0.5。
[0145] 结合图6可知，若当前帖样本与已有目标高斯分布样本均值相差较大，则学习速率 上升较快，从而分类器更新变慢;若二者高斯分布均值相差不大，说明跟踪性能良好，则学 习速率较低，从而分类器更新速度维持在较快水平。
[0146] S8、W分类器响应值最大的检测样本所对应的矩形区域作为当前帖的目标跟踪 框，并用前一帖正样本模板与初始帖正样本模板的加权和来更新正样本模板；
[0147] 正样本模板建立在样本对应的低维特征向量基础上，本【具体实施方式】中中正样本 模板的罕就公古责.
[014 引
[0149] 其中，9表示当前帖正样本模板，它由上一时刻正样本模板巧和初始时刻目标框低 维特征向量口 /按照一定比重共同构成，本算法中y取0.9。
[0150] S9、W当前帖的目标跟踪框为基准，在当前帖中采集一组与当前帖的目标跟踪框 同尺度的正、负样本，提取所述正、负样本的尺度不变特征并对分类器进行重新训练；
[0151] S10、如视频已结束，则退出目标跟踪;否则，W下一帖作为当前帖，转至S3。
【主权项】
1. 一种基于旋翼无人机平台的目标跟踪方法，其特征在于，包括以下步骤： 51、 对旋翼无人机平台所拍摄的视频的每一帧进行自动目标识别，以初次识别出目标 的视频帧作为目标跟踪的初始帧，并在初始帧中目标周围划定矩形的初始目标跟踪框； 52、 以初始目标跟踪框为基准，在初始帧中采集一组与初始目标跟踪框同尺度的正、负 样本，提取所述样本的尺度不变特征对分类器进行训练，并以初始目标跟踪框中的图像作 为正样本模板； 53、 在当前帧中，以前一帧的目标跟踪框位置为中心的一定区域内采集一组不同尺度 的矩形区域作为检测样本，并在每个检测样本上叠加均值为O的高斯白噪声;提取叠加噪声 后的各检测样本的尺度不变特征并分别作为所述分类器的输入，得到各检测样本的分类器 响应值； 54、 判断当前帧的跟踪性能是否稳定，如性能不稳定，则转至S5,否则，转至S8;如以下 条件之一得到满足，则当前帧的跟踪性能不稳定： (1) 当前帧检测样本的最大分类器响应值小于〇; (2) 当前帧检测样本的最大分类器响应值大于等于0,同时当前帧检测样本的最大分类 器响应值小于前一帧检测样本的最大分类器响应值，且前一帧的跟踪性能不稳定； (3) 当前帧检测样本的最大分类器响应值在0和一个预设的大于0的阈值之间，同时当 前帧检测样本的最大分类器响应值小于前一帧检测样本的最大分类器响应值，且前一帧的 跟踪性能稳定； 55、 计算当前帧各检测样本与正样本模板之间的归一化相似度，并以归一化相似度作 为修正系数对当前帧相应检测样本的分类器响应值进行修正；以修正后的分类器响应值最 大的检测样本所对应的矩形区域作为当前帧的目标跟踪框； 56、 判断连续跟踪不稳定的帧数是否达到预设上限值，如是，则将所述高斯白噪声的标 准差增大后转SlO;否则，转至S7; 57、 以当前帧的目标跟踪框为基准，在当前帧中采集一组与当前帧的目标跟踪框同尺 度的正、负样本，从所有正样本中选出与正样本模板之间相似度最大的一部分;然后提取所 选出的正样本以及负样本的尺度不变特征并对分类器进行重新训练后转至S10; 58、 以分类器响应值最大的检测样本所对应的矩形区域作为当前帧的目标跟踪框，并 用前一帧正样本模板与初始帧正样本模板的加权和来更新正样本模板； 59、 以当前帧的目标跟踪框为基准，在当前帧中采集一组与当前帧的目标跟踪框同尺 度的正、负样本，提取所述正、负样本的尺度不变特征并对分类器进行重新训练； S10、如视频已结束，则退出目标跟踪;否则，以下一帧作为当前帧，转至S3。2. 如权利要求1所述目标跟踪方法，其特征在于，所述分类器为朴素贝叶斯分类器。3. 如权利要求1所述目标跟踪方法，其特征在于，所述尺度不变特征为尺度不变压缩特 征。4. 如权利要求1所述目标跟踪方法，其特征在于，所述分类器的学习速率λ在目标跟踪 过程中按照下式进行自适应更新·式中，X为当前帧的目标跟踪框中图像高斯分布均值与前一帧的目标跟踪框中图像高 斯分布均值之间差值的绝对I均为预设参数。5. 如权利要求1所述目标跟踪方法，其特征在于，检测样本与正样本模板之间的相似度 使用以下方式度量:检测样本与正样本模板的各维特征值的差值绝对值之和。6. 如权利要求1所述目标跟踪方法，其特征在于，利用以下方法对旋翼无人机平台所拍 摄的视频的每一帧进行自动目标识别：首先提取当前帧中的感兴趣区域;然后将感兴趣区 域与目标模板先后进行主颜色匹配和轮廓匹配，如存在颜色匹配和轮廓匹配均匹配成功的 感兴趣区域，则该感兴趣区域即为在当前帧中识别出的目标，否则，继续对下一帧进行目标 识别。7. 如权利要求6所述目标跟踪方法，其特征在于，所述轮廓匹配中所使用的轮廓特征为 Hu不变矩特征。
【文档编号】G06K9/62GK106023257SQ201610355179
【公开日】2016年10月12日
【申请日】2016年5月26日
【发明人】贾春叶, 甄子洋, 马坤
【申请人】南京航空航天大学