• 四川郎酒股份有限公司获第十二届人民企业社会责任奖年度环保奖 2019-05-13
  • 银保监会新规剑指大企业多头融资和过度融资 2019-05-12
  • 韩国再提4国联合申办世界杯 中国网友无视:我们自己来 2019-05-11
  • 中国人为什么一定要买房? 2019-05-11
  • 十九大精神进校园:风正扬帆当有为 勇做时代弄潮儿 2019-05-10
  • 粽叶飘香幸福邻里——廊坊市举办“我们的节日·端午”主题活动 2019-05-09
  • 太原设禁鸣路段 设备在测试中 2019-05-09
  • 拜耳医药保健有限公司获第十二届人民企业社会责任奖年度企业奖 2019-05-08
  • “港独”没出路!“梁天琦们”该醒醒了 2019-05-07
  • 陈卫平:中国文化内涵包含三方面 文化复兴表现在其中 2019-05-06
  • 人民日报客户端辟谣:“合成军装照”产品请放心使用 2019-05-05
  • 【十九大·理论新视野】为什么要“建设现代化经济体系”?   2019-05-04
  • 聚焦2017年乌鲁木齐市老城区改造提升工程 2019-05-04
  • 【专家谈】上合组织——构建区域命运共同体的有力实践者 2019-05-03
  • 【华商侃车NO.192】 亲!楼市火爆,别忘了买车位啊! 2019-05-03
    • / 17
    • 下载费用:30 金币  

    重庆时时彩怎样看冷热球: 一种基于分层模型的无人机目标跟踪方法.pdf

    关 键 词:
    一种 基于 分层 模型 无人机 目标 跟踪 方法
      专利查询网所有资源均是用户自行上传分享,仅供网友学习交流,未经上传用户书面授权,请勿作他用。
    摘要
    申请专利号:

    CN201611262156.2

    申请日:

    2016.12.30

    公开号:

    CN106683118A

    公开日:

    2017.05.17

    当前法律状态:

    实审

    有效性:

    审中

    法律详情: 实质审查的生效IPC(主分类):G06T 7/246申请日:20161230|||公开
    IPC分类号: G06T7/246(2017.01)I; G06T7/277(2017.01)I; G06N3/00 主分类号: G06T7/246
    申请人: 北京航空航天大学
    发明人: 丁文锐; 刘春蕾; 李红光
    地址: 100191 北京市海淀区学院路37号
    优先权:
    专利代理机构: 北京永创新实专利事务所 11121 代理人: 赵文颖
    PDF完整版下载: PDF下载
    法律状态
    申请(专利)号:

    CN201611262156.2

    授权公告号:

    |||

    法律状态公告日:

    2017.06.09|||2017.05.17

    法律状态类型:

    实质审查的生效|||公开

    摘要

    本发明公开了一种基于分层模型的无人机目标跟踪方法,属于无人机遥感图像处理领域,包括以下步骤:第一步,建立初始帧目标的特征表达;第二步,建立以均值漂移为基础的粗跟踪模型;第三步,建立以改进的粒子群优化为基础的精跟踪模型。本算法通过引入均值漂移与改进的粒子群优化结合的分层模型,有效克服了因无人机成像特殊性造成视频跟踪难度大,不准确且易丢失的问题,有效提高了目标跟踪的精度。

    权利要求书

    1.一种基于分层模型的无人机目标跟踪方法,包括以下步骤:
    第一步,建立初始帧目标的特征表达
    读入视频序列,选取要跟踪的目标区域,对目标区域建立基于核函数加权的归一化空
    间颜色直方图模型,并将该模型的特征向量作为目标区域的特征表达;
    第二步,建立以均值漂移为基础的粗跟踪模型
    1)建立均值漂移在跟踪中的数学表达
    设均值漂移偏移量为:
    <mrow> <mi>m</mi> <mrow> <mo>(</mo> <mi>y</mi> <mo>)</mo> </mrow> <mo>=</mo> <mi>y</mi> <mo>-</mo> <msub> <mi>y</mi> <mn>0</mn> </msub> <mo>=</mo> <mfrac> <mrow> <msubsup> <mi>&Sigma;</mi> <mrow> <mi>i</mi> <mo>=</mo> <mn>1</mn> </mrow> <mi>n</mi> </msubsup> <msub> <mi>x</mi> <mi>i</mi> </msub> <msub> <mi>w</mi> <mi>i</mi> </msub> <mi>g</mi> <mrow> <mo>(</mo> <mo>|</mo> <mo>|</mo> <mfrac> <mrow> <mi>y</mi> <mo>-</mo> <msub> <mi>x</mi> <mi>i</mi> </msub> </mrow> <mi>h</mi> </mfrac> <mo>|</mo> <msup> <mo>|</mo> <mn>2</mn> </msup> <mo>)</mo> </mrow> </mrow> <mrow> <msubsup> <mi>&Sigma;</mi> <mrow> <mi>i</mi> <mo>=</mo> <mn>1</mn> </mrow> <mi>n</mi> </msubsup> <msub> <mi>w</mi> <mi>i</mi> </msub> <mi>g</mi> <mrow> <mo>(</mo> <mo>|</mo> <mo>|</mo> <mfrac> <mrow> <mi>y</mi> <mo>-</mo> <msub> <mi>x</mi> <mi>i</mi> </msub> </mrow> <mi>h</mi> </mfrac> <mo>|</mo> <msup> <mo>|</mo> <mn>2</mn> </msup> <mo>)</mo> </mrow> </mrow> </mfrac> <mo>-</mo> <msub> <mi>y</mi> <mn>0</mn> </msub> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>15</mn> <mo>)</mo> </mrow> </mrow>
    其中:
    y为当前帧候选区域中心位置,y0为上一帧目标的中心位置,n为样本点的个数,即上一
    帧目标位置处的像素个数,xi为当前帧候选区域的第i个像素点的位置,h为带宽,k(·)为
    Epannechnikov核函数,k'(·)表示核函数k(·)的导数,fu表示目标模型中第u个特征值估
    计的概率密度,u为直方图中的颜色索引值,δ[·]为单位冲激响应,b(xi)为位置xi处的颜色
    量化等级,pu(y0)表示当前帧中位于y0处的候选目标的颜色概率密度;
    2)以均值漂移为基础的粗跟踪
    在无人机视频跟踪过程中,将步骤一中所建立的基于核函数加权的归一化空间颜色直
    方图模型作为均值漂移过程中目标的特征值,进行迭代,得到跟踪匹配结果,设置迭代次数
    N,当候选区域和目标区域的相似度大于阈值T时,跳出迭代;
    第三步,建立以改进的粒子群优化为基础的精跟踪模型
    判断目标是否丢失,如果目标未丢失,利用精跟踪模型对位置进行进一步精确定位;若
    目标丢失,利用精跟踪模型对目标重新进行搜索定位;若精跟踪后仍未搜索到目标,则再次
    对目标进行选取。
    2.根据权利要求1所述的一种基于分层模型的无人机目标跟踪方法,所述的第一步中:
    核函数采用Epanechnikov核函数,如式(1):
    <mrow> <mi>K</mi> <mrow> <mo>(</mo> <mo>|</mo> <mo>|</mo> <mi>r</mi> <mo>|</mo> <mo>|</mo> <mo>)</mo> </mrow> <mo>=</mo> <mfenced open = "{" close = ""> <mtable> <mtr> <mtd> <mrow> <mn>1</mn> <mo>-</mo> <mo>|</mo> <mo>|</mo> <mi>r</mi> <mo>|</mo> <msup> <mo>|</mo> <mn>2</mn> </msup> </mrow> </mtd> <mtd> <mrow> <mo>|</mo> <mo>|</mo> <mi>r</mi> <mo>|</mo> <mo>|</mo> <mo>&le;</mo> <mn>1</mn> </mrow> </mtd> </mtr> <mtr> <mtd> <mn>0</mn> </mtd> <mtd> <mrow> <mi>e</mi> <mi>l</mi> <mi>s</mi> <mi>e</mi> </mrow> </mtd> </mtr> </mtable> </mfenced> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>1</mn> <mo>)</mo> </mrow> </mrow>
    其中,K(·)为Epanechnikov核函数,||r||为像素距离目标中心的距离。
    3.根据权利要求1所述的一种基于分层模型的无人机目标跟踪方法,所述的第一步中:
    空间颜色直方图模型为RGB颜色空间,将R、G、B三通道中[0,255]的取值范围转换到[0,
    16],将三维向量映射为直方图上的一维向量,(R,G,B)=256×R+16×G+B,并将其做归一化
    处理,该归一化后的直方图即为初始化阶段中视频初始帧目标的特征表达。
    4.根据权利要求1所述的一种基于分层模型的无人机目标跟踪方法,所述的第三步中:
    改进的粒子群优化精跟踪模型包括:第一,对粒子群优化中的参数进行初始化;第二,
    建立在线更新的运动模型;第三,建立基于粒子群优化算法的自适应性参数调整跟踪模型;
    具体的:
    1)参数初始化
    在粒子群优化跟踪模型中:
    粒子群种群数目:N;
    迭代次数:Tpso;
    粒子0时刻的初始飞行速度:
    粒子0时刻的初始位置:
    惯性参数:w1;
    学习因子:c1、c2;
    随机数:r1、r2;
    2)建立运动模型
    假设目标在(x,y)处转弯,转弯时的速率为v,速度与x轴方向所成夹角为θ,将拐弯处的
    状态描述为(x,x',y,y'),其中x',y'表示x,y方向上的速度大??;假设转弯角速率为ω,得
    到其与时间t的关系;
    <mrow> <mi>&omega;</mi> <mrow> <mo>(</mo> <mi>t</mi> <mo>)</mo> </mrow> <mo>=</mo> <mfrac> <mrow> <mi>d</mi> <mi>&theta;</mi> </mrow> <mrow> <mi>d</mi> <mi>t</mi> </mrow> </mfrac> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>16</mn> <mo>)</mo> </mrow> </mrow>
    x方向上速度的大?。?br />x'=vx=vcosθ (17)
    y方向上速度的大?。?br />y'=vy=vsinθ (18)
    x方向上加速度的大?。?br />x”=v'x=-ω(t)vsinθ (19)
    y方向上加速度的大?。?br />y”=v'y=ω(t)vcosθ (20)
    建立转弯模型的状态描述方程为:
    <mrow> <mfenced open = "[" close = "]"> <mtable> <mtr> <mtd> <msup> <mi>x</mi> <mo>&prime;</mo> </msup> </mtd> </mtr> <mtr> <mtd> <msup> <mi>x</mi> <mrow> <mo>&prime;</mo> <mo>&prime;</mo> </mrow> </msup> </mtd> </mtr> <mtr> <mtd> <msup> <mi>y</mi> <mo>&prime;</mo> </msup> </mtd> </mtr> <mtr> <mtd> <msup> <mi>y</mi> <mrow> <mo>&prime;</mo> <mo>&prime;</mo> </mrow> </msup> </mtd> </mtr> </mtable> </mfenced> <mo>=</mo> <mfenced open = "[" close = "]"> <mtable> <mtr> <mtd> <mn>0</mn> </mtd> <mtd> <mn>1</mn> </mtd> <mtd> <mn>0</mn> </mtd> <mtd> <mn>0</mn> </mtd> </mtr> <mtr> <mtd> <mn>0</mn> </mtd> <mtd> <mn>0</mn> </mtd> <mtd> <mn>0</mn> </mtd> <mtd> <mrow> <mo>-</mo> <mi>&omega;</mi> <mrow> <mo>(</mo> <mi>t</mi> <mo>)</mo> </mrow> </mrow> </mtd> </mtr> <mtr> <mtd> <mn>0</mn> </mtd> <mtd> <mn>0</mn> </mtd> <mtd> <mn>0</mn> </mtd> <mtd> <mn>1</mn> </mtd> </mtr> <mtr> <mtd> <mn>0</mn> </mtd> <mtd> <mrow> <mi>&omega;</mi> <mrow> <mo>(</mo> <mi>t</mi> <mo>)</mo> </mrow> </mrow> </mtd> <mtd> <mn>0</mn> </mtd> <mtd> <mn>0</mn> </mtd> </mtr> </mtable> </mfenced> <mfenced open = "[" close = "]"> <mtable> <mtr> <mtd> <mi>x</mi> </mtd> </mtr> <mtr> <mtd> <msup> <mi>x</mi> <mo>&prime;</mo> </msup> </mtd> </mtr> <mtr> <mtd> <mi>y</mi> </mtd> </mtr> <mtr> <mtd> <msup> <mi>y</mi> <mo>&prime;</mo> </msup> </mtd> </mtr> </mtable> </mfenced> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>21</mn> <mo>)</mo> </mrow> </mrow>
    状态转移方程的形式为:
    Xk+1=FkXk+Wk (22)
    其中,Fk为转弯模型的系统状态转移矩阵,Xk+1,Xk表示k+1时刻和k时刻目标的状态,即
    [x,x',y,y']T,Wk为转弯时的系统状态噪声;
    将式(21)代入式(22),得到:
    <mrow> <mfrac> <mrow> <msub> <mi>dF</mi> <mi>k</mi> </msub> </mrow> <mrow> <mi>d</mi> <mi>t</mi> </mrow> </mfrac> <mo>=</mo> <mi>A</mi> <mrow> <mo>(</mo> <mi>t</mi> <mo>)</mo> </mrow> <msub> <mi>F</mi> <mi>k</mi> </msub> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>23</mn> <mo>)</mo> </mrow> </mrow>
    其中A(t)为
    <mrow> <mi>A</mi> <mrow> <mo>(</mo> <mi>t</mi> <mo>)</mo> </mrow> <mo>=</mo> <mfenced open = "[" close = "]"> <mtable> <mtr> <mtd> <mn>0</mn> </mtd> <mtd> <mn>1</mn> </mtd> <mtd> <mn>0</mn> </mtd> <mtd> <mn>0</mn> </mtd> </mtr> <mtr> <mtd> <mn>0</mn> </mtd> <mtd> <mn>0</mn> </mtd> <mtd> <mn>0</mn> </mtd> <mtd> <mrow> <mo>-</mo> <mi>&omega;</mi> <mrow> <mo>(</mo> <mi>t</mi> <mo>)</mo> </mrow> </mrow> </mtd> </mtr> <mtr> <mtd> <mn>0</mn> </mtd> <mtd> <mn>0</mn> </mtd> <mtd> <mn>0</mn> </mtd> <mtd> <mn>1</mn> </mtd> </mtr> <mtr> <mtd> <mn>0</mn> </mtd> <mtd> <mrow> <mi>&omega;</mi> <mrow> <mo>(</mo> <mi>t</mi> <mo>)</mo> </mrow> </mrow> </mtd> <mtd> <mn>0</mn> </mtd> <mtd> <mn>0</mn> </mtd> </mtr> </mtable> </mfenced> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>24</mn> <mo>)</mo> </mrow> </mrow>
    首先将A(t)离散化,得到
    <mrow> <mi>A</mi> <mrow> <mo>(</mo> <mi>T</mi> <mo>)</mo> </mrow> <mo>=</mo> <mfenced open = "[" close = "]"> <mtable> <mtr> <mtd> <mn>0</mn> </mtd> <mtd> <mn>1</mn> </mtd> <mtd> <mn>0</mn> </mtd> <mtd> <mn>0</mn> </mtd> </mtr> <mtr> <mtd> <mn>0</mn> </mtd> <mtd> <mn>0</mn> </mtd> <mtd> <mn>0</mn> </mtd> <mtd> <mrow> <mo>-</mo> <mi>&omega;</mi> <mrow> <mo>(</mo> <mi>T</mi> <mo>)</mo> </mrow> </mrow> </mtd> </mtr> <mtr> <mtd> <mn>0</mn> </mtd> <mtd> <mn>0</mn> </mtd> <mtd> <mn>0</mn> </mtd> <mtd> <mn>1</mn> </mtd> </mtr> <mtr> <mtd> <mn>0</mn> </mtd> <mtd> <mrow> <mi>&omega;</mi> <mrow> <mo>(</mo> <mi>T</mi> <mo>)</mo> </mrow> </mrow> </mtd> <mtd> <mn>0</mn> </mtd> <mtd> <mn>0</mn> </mtd> </mtr> </mtable> </mfenced> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>25</mn> <mo>)</mo> </mrow> </mrow>
    对式(23)进行拉普拉斯变换
    sL[Fk(s)]-Fk(0)=AFk(s) (26)
    合并同类项为
    (sI-A)Fk(s)=Fk(0) (27)
    其中I表示单位矩阵;
    已知
    Fk(0)=(f1(0),f2(0),f3(0),...,fn(0))T=I (28)

    <mrow> <msub> <mi>F</mi> <mi>k</mi> </msub> <mrow> <mo>(</mo> <mi>s</mi> <mo>)</mo> </mrow> <mo>=</mo> <msup> <mrow> <mo>(</mo> <mi>s</mi> <mi>I</mi> <mo>-</mo> <mi>A</mi> <mo>)</mo> </mrow> <mrow> <mo>-</mo> <mn>1</mn> </mrow> </msup> <msub> <mi>F</mi> <mi>k</mi> </msub> <mrow> <mo>(</mo> <mn>0</mn> <mo>)</mo> </mrow> <mo>=</mo> <mfrac> <msup> <mrow> <mo>(</mo> <mi>s</mi> <mi>I</mi> <mo>-</mo> <mi>A</mi> <mo>)</mo> </mrow> <mo>*</mo> </msup> <mrow> <mo>|</mo> <mi>s</mi> <mi>I</mi> <mo>-</mo> <mi>A</mi> <mo>|</mo> </mrow> </mfrac> <mi>I</mi> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>29</mn> <mo>)</mo> </mrow> </mrow>
    其中
    <mrow> <msup> <mrow> <mo>(</mo> <mi>s</mi> <mi>I</mi> <mo>-</mo> <mi>A</mi> <mo>)</mo> </mrow> <mrow> <mo>-</mo> <mn>1</mn> </mrow> </msup> <mo>=</mo> <mfenced open = "[" close = "]"> <mtable> <mtr> <mtd> <mfrac> <mn>1</mn> <mi>s</mi> </mfrac> </mtd> <mtd> <mn>0</mn> </mtd> <mtd> <mn>0</mn> </mtd> <mtd> <mn>0</mn> </mtd> </mtr> <mtr> <mtd> <mfrac> <mn>1</mn> <mrow> <msup> <mi>s</mi> <mn>2</mn> </msup> <mo>+</mo> <msup> <mi>&omega;</mi> <mn>2</mn> </msup> </mrow> </mfrac> </mtd> <mtd> <mfrac> <mi>s</mi> <mrow> <msup> <mi>s</mi> <mn>2</mn> </msup> <mo>+</mo> <msup> <mi>&omega;</mi> <mn>2</mn> </msup> </mrow> </mfrac> </mtd> <mtd> <mfrac> <mi>&omega;</mi> <mrow> <mi>s</mi> <mrow> <mo>(</mo> <msup> <mi>s</mi> <mn>2</mn> </msup> <mo>+</mo> <msup> <mi>&omega;</mi> <mn>2</mn> </msup> <mo>)</mo> </mrow> </mrow> </mfrac> </mtd> <mtd> <mfrac> <mi>&omega;</mi> <mrow> <msup> <mi>s</mi> <mn>2</mn> </msup> <mo>+</mo> <msup> <mi>&omega;</mi> <mn>2</mn> </msup> </mrow> </mfrac> </mtd> </mtr> <mtr> <mtd> <mn>0</mn> </mtd> <mtd> <mn>0</mn> </mtd> <mtd> <mfrac> <mn>1</mn> <mi>s</mi> </mfrac> </mtd> <mtd> <mn>0</mn> </mtd> </mtr> <mtr> <mtd> <mfrac> <mrow> <mo>-</mo> <mi>&omega;</mi> </mrow> <mrow> <mi>s</mi> <mrow> <mo>(</mo> <msup> <mi>s</mi> <mn>2</mn> </msup> <mo>+</mo> <msup> <mi>&omega;</mi> <mn>2</mn> </msup> <mo>)</mo> </mrow> </mrow> </mfrac> </mtd> <mtd> <mfrac> <mrow> <mo>-</mo> <mi>&omega;</mi> </mrow> <mrow> <msup> <mi>s</mi> <mn>2</mn> </msup> <mo>+</mo> <msup> <mi>&omega;</mi> <mn>2</mn> </msup> </mrow> </mfrac> </mtd> <mtd> <mfrac> <mn>1</mn> <mrow> <msup> <mi>s</mi> <mn>2</mn> </msup> <mo>+</mo> <msup> <mi>&omega;</mi> <mn>2</mn> </msup> </mrow> </mfrac> </mtd> <mtd> <mfrac> <mi>s</mi> <mrow> <msup> <mi>s</mi> <mn>2</mn> </msup> <mo>+</mo> <msup> <mi>&omega;</mi> <mn>2</mn> </msup> </mrow> </mfrac> </mtd> </mtr> </mtable> </mfenced> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>30</mn> <mo>)</mo> </mrow> </mrow>
    系统状态矩阵Fk Fk=L-1((sI-A)-1) (31)
    通过拉普拉斯逆变换,得:
    <mrow> <msub> <mi>F</mi> <mi>k</mi> </msub> <mo>=</mo> <mfenced open = "[" close = "]"> <mtable> <mtr> <mtd> <mn>1</mn> </mtd> <mtd> <mfrac> <mrow> <mi>sin</mi> <mrow> <mo>(</mo> <mrow> <mi>&omega;</mi> <mi>T</mi> </mrow> <mo>)</mo> </mrow> </mrow> <mi>&omega;</mi> </mfrac> </mtd> <mtd> <mn>0</mn> </mtd> <mtd> <mrow> <mo>-</mo> <mfrac> <mrow> <mn>1</mn> <mo>-</mo> <mi>cos</mi> <mrow> <mo>(</mo> <mrow> <mi>&omega;</mi> <mi>T</mi> </mrow> <mo>)</mo> </mrow> </mrow> <mi>&omega;</mi> </mfrac> </mrow> </mtd> </mtr> <mtr> <mtd> <mn>0</mn> </mtd> <mtd> <mrow> <mi>cos</mi> <mrow> <mo>(</mo> <mrow> <mi>&omega;</mi> <mi>T</mi> </mrow> <mo>)</mo> </mrow> </mrow> </mtd> <mtd> <mn>0</mn> </mtd> <mtd> <mrow> <mo>-</mo> <mi>sin</mi> <mrow> <mo>(</mo> <mrow> <mi>&omega;</mi> <mi>T</mi> </mrow> <mo>)</mo> </mrow> </mrow> </mtd> </mtr> <mtr> <mtd> <mn>0</mn> </mtd> <mtd> <mfrac> <mrow> <mn>1</mn> <mo>-</mo> <mi>cos</mi> <mrow> <mo>(</mo> <mrow> <mi>&omega;</mi> <mi>T</mi> </mrow> <mo>)</mo> </mrow> </mrow> <mi>&omega;</mi> </mfrac> </mtd> <mtd> <mn>1</mn> </mtd> <mtd> <mfrac> <mrow> <mi>sin</mi> <mrow> <mo>(</mo> <mrow> <mi>&omega;</mi> <mi>T</mi> <mo>)</mo> </mrow> <mo>)</mo> </mrow> </mrow> <mi>&omega;</mi> </mfrac> </mtd> </mtr> <mtr> <mtd> <mn>0</mn> </mtd> <mtd> <mrow> <mi>sin</mi> <mrow> <mo>(</mo> <mrow> <mi>&omega;</mi> <mi>T</mi> </mrow> <mo>)</mo> </mrow> </mrow> </mtd> <mtd> <mn>0</mn> </mtd> <mtd> <mrow> <mi>cos</mi> <mrow> <mo>(</mo> <mrow> <mi>&omega;</mi> <mi>T</mi> </mrow> <mo>)</mo> </mrow> </mrow> </mtd> </mtr> </mtable> </mfenced> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>32</mn> <mo>)</mo> </mrow> </mrow>
    得到转弯下状态矩阵的表达式;当ω=0时,状态矩阵Fk的表达式化简为:
    <mrow> <msub> <mi>F</mi> <mi>k</mi> </msub> <mo>=</mo> <mfenced open = "[" close = "]"> <mtable> <mtr> <mtd> <mn>1</mn> </mtd> <mtd> <mn>1</mn> </mtd> <mtd> <mn>0</mn> </mtd> <mtd> <mn>0</mn> </mtd> </mtr> <mtr> <mtd> <mn>0</mn> </mtd> <mtd> <mn>1</mn> </mtd> <mtd> <mn>0</mn> </mtd> <mtd> <mn>0</mn> </mtd> </mtr> <mtr> <mtd> <mn>0</mn> </mtd> <mtd> <mn>0</mn> </mtd> <mtd> <mn>1</mn> </mtd> <mtd> <mn>1</mn> </mtd> </mtr> <mtr> <mtd> <mn>0</mn> </mtd> <mtd> <mn>0</mn> </mtd> <mtd> <mn>0</mn> </mtd> <mtd> <mn>1</mn> </mtd> </mtr> </mtable> </mfenced> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>33</mn> <mo>)</mo> </mrow> </mrow>
    通过上述步骤,建立了直线运动和转弯运动情况下的状态矩阵,通过第二步均值漂移
    粗跟踪后,再利用粒子群优化的搜索方法进行精跟踪;
    无人机实际飞行转弯运动模型中,需要实时在线调整角速度,具体为:
    首先将路径参数t设置为时间,建立x方向路径与y方向路径与t的关系,拟合从第11帧
    开始,以前10帧目标位置进行拟合得到x1(t)和y1(t),计算不同时刻处的曲率:
    <mrow> <mi>&kappa;</mi> <mrow> <mo>(</mo> <mi>t</mi> <mo>)</mo> </mrow> <mo>=</mo> <mfrac> <mrow> <msub> <mover> <mi>x</mi> <mo>&CenterDot;</mo> </mover> <mn>1</mn> </msub> <mrow> <mo>(</mo> <mi>t</mi> <mo>)</mo> </mrow> <msub> <mover> <mi>y</mi> <mo>&CenterDot;&CenterDot;</mo> </mover> <mn>1</mn> </msub> <mrow> <mo>(</mo> <mi>t</mi> <mo>)</mo> </mrow> <mo>-</mo> <msub> <mover> <mi>y</mi> <mo>&CenterDot;</mo> </mover> <mn>1</mn> </msub> <mrow> <mo>(</mo> <mi>t</mi> <mo>)</mo> </mrow> <msub> <mover> <mi>x</mi> <mo>&CenterDot;&CenterDot;</mo> </mover> <mn>1</mn> </msub> <mrow> <mo>(</mo> <mi>t</mi> <mo>)</mo> </mrow> </mrow> <msup> <mrow> <mo>(</mo> <msqrt> <mrow> <msubsup> <mover> <mi>x</mi> <mo>&CenterDot;</mo> </mover> <mn>1</mn> <mn>2</mn> </msubsup> <mo>+</mo> <msubsup> <mover> <mi>y</mi> <mo>&CenterDot;</mo> </mover> <mn>1</mn> <mn>2</mn> </msubsup> </mrow> </msqrt> <mo>)</mo> </mrow> <mn>3</mn> </msup> </mfrac> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>34</mn> <mo>)</mo> </mrow> </mrow>
    式中,是相对于路径参数的导
    数;
    曲率半径:
    <mrow> <mi>r</mi> <mrow> <mo>(</mo> <mi>t</mi> <mo>)</mo> </mrow> <mo>=</mo> <mfrac> <mn>1</mn> <mrow> <mi>&kappa;</mi> <mrow> <mo>(</mo> <mi>t</mi> <mo>)</mo> </mrow> </mrow> </mfrac> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>35</mn> <mo>)</mo> </mrow> </mrow>
    角速度为:ω(t)=v×r(t) (36)
    3)基于粒子群优化方法的分层性跟踪
    包括:
    ①计算初始化后的每个粒子的适应度值;
    ②根据适应度值对粒子进行排序,找出种群最优粒子个体;
    ③在前10帧默认运动模型为直线模型,之后根据在线更新的角速度值进行调整,利用
    运动模型分别对第i个粒子的位置进行更新,并令i自增;若运动模型为直线模型,则位置更
    新公式表示为:
    <mrow> <msubsup> <mi>x</mi> <mi>i</mi> <mi>t</mi> </msubsup> <mo>=</mo> <msubsup> <mi>x</mi> <mi>i</mi> <mrow> <mi>t</mi> <mo>-</mo> <mn>1</mn> </mrow> </msubsup> <mo>+</mo> <msubsup> <mi>v</mi> <mi>i</mi> <mi>t</mi> </msubsup> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>37</mn> <mo>)</mo> </mrow> </mrow>
    式中,表示粒子i在t时刻的位置,表示粒子i在t时刻的飞行速度;
    若运动模型为转弯模型,则位置更新公式在x、y轴两个方向表示为:
    <mrow> <msubsup> <mi>x</mi> <mrow> <mi>i</mi> <mn>1</mn> </mrow> <mi>t</mi> </msubsup> <mo>=</mo> <msubsup> <mi>x</mi> <mrow> <mi>i</mi> <mn>1</mn> </mrow> <mrow> <mi>t</mi> <mo>-</mo> <mn>1</mn> </mrow> </msubsup> <mo>+</mo> <mfrac> <mrow> <mi>s</mi> <mi>i</mi> <mi>n</mi> <mrow> <mo>(</mo> <mi>&omega;</mi> <mi>T</mi> <mo>)</mo> </mrow> </mrow> <mi>&omega;</mi> </mfrac> <msubsup> <mi>v</mi> <mrow> <mi>i</mi> <mn>1</mn> </mrow> <mrow> <mi>t</mi> <mo>-</mo> <mn>1</mn> </mrow> </msubsup> <mo>-</mo> <mfrac> <mrow> <mn>1</mn> <mo>-</mo> <mi>c</mi> <mi>o</mi> <mi>s</mi> <mrow> <mo>(</mo> <mi>&omega;</mi> <mi>T</mi> <mo>)</mo> </mrow> </mrow> <mi>&omega;</mi> </mfrac> <msubsup> <mi>v</mi> <mrow> <mi>i</mi> <mn>2</mn> </mrow> <mrow> <mi>t</mi> <mo>-</mo> <mn>1</mn> </mrow> </msubsup> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>38</mn> <mo>)</mo> </mrow> </mrow>
    <mrow> <msubsup> <mi>x</mi> <mrow> <mi>i</mi> <mn>2</mn> </mrow> <mi>t</mi> </msubsup> <mo>=</mo> <msubsup> <mi>x</mi> <mrow> <mi>i</mi> <mn>2</mn> </mrow> <mrow> <mi>t</mi> <mo>-</mo> <mn>1</mn> </mrow> </msubsup> <mo>+</mo> <mfrac> <mrow> <mn>1</mn> <mo>-</mo> <mi>c</mi> <mi>o</mi> <mi>s</mi> <mrow> <mo>(</mo> <mi>&omega;</mi> <mi>T</mi> <mo>)</mo> </mrow> </mrow> <mi>&omega;</mi> </mfrac> <msubsup> <mi>v</mi> <mrow> <mi>i</mi> <mn>1</mn> </mrow> <mrow> <mi>t</mi> <mo>-</mo> <mn>1</mn> </mrow> </msubsup> <mo>+</mo> <mfrac> <mrow> <mi>s</mi> <mi>i</mi> <mi>n</mi> <mrow> <mo>(</mo> <mi>&omega;</mi> <mi>T</mi> <mo>)</mo> </mrow> </mrow> <mi>&omega;</mi> </mfrac> <msubsup> <mi>v</mi> <mrow> <mi>i</mi> <mn>2</mn> </mrow> <mrow> <mi>t</mi> <mo>-</mo> <mn>1</mn> </mrow> </msubsup> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>39</mn> <mo>)</mo> </mrow> </mrow>
    其中:i1代表粒子i在x方向的数据情况,i2代表粒子i在y方向的数据情况;
    ④计算新位置的适应度值;判断新位置是否优于该粒子最佳位置,如果满足条件,用新
    位置替换原来位置,否则,原位置不变;判断新位置是否优于当前最佳位置;如果不满足条
    件,则当前最佳位置不变,否则,用新位置替换当前最佳位置,并调整粒子速度;
    <mrow> <msubsup> <mi>v</mi> <mi>i</mi> <mi>t</mi> </msubsup> <mo>=</mo> <msub> <mi>w</mi> <mn>1</mn> </msub> <msubsup> <mi>v</mi> <mi>i</mi> <mrow> <mi>t</mi> <mo>-</mo> <mn>1</mn> </mrow> </msubsup> <mo>+</mo> <msub> <mi>c</mi> <mn>1</mn> </msub> <mi>r</mi> <mrow> <mo>(</mo> <msubsup> <mi>x</mi> <mi>i</mi> <mrow> <mi>t</mi> <mo>-</mo> <mn>1</mn> </mrow> </msubsup> <mo>-</mo> <msubsup> <mi>x</mi> <mi>i</mi> <mrow> <mi>b</mi> <mi>e</mi> <mi>s</mi> <mi>t</mi> </mrow> </msubsup> <mo>)</mo> </mrow> <mo>+</mo> <msub> <mi>c</mi> <mn>2</mn> </msub> <mi>r</mi> <mrow> <mo>(</mo> <msubsup> <mi>x</mi> <mi>i</mi> <mrow> <mi>t</mi> <mo>-</mo> <mn>1</mn> </mrow> </msubsup> <mo>-</mo> <msup> <mi>x</mi> <mo>*</mo> </msup> <mo>)</mo> </mrow> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>40</mn> <mo>)</mo> </mrow> </mrow>
    其中:x*表示在当前群体中最佳粒子所处位置,表示粒子的最佳位置,w1为惯性参
    数,c1、c2为学习因子,r为[0,1]之间的随机数;
    ⑤判断是否满足迭代终止条件;如果满足,则粒子群优化搜索结束,否则执行③;
    最终,步骤④中得到的群体粒子最佳位置x*为无人机视频中目标跟踪的精确结果。
    5.根据权利要求4所述的一种基于分层模型的无人机目标跟踪方法,所述的1)中:
    种群数目N=10,迭代次数Tpso=30;粒子0时刻的初始飞行速度粒子0时刻的初
    始位置为以上一帧目标为中心进行高斯分布;随机数r1=r2=r;学习因子c1=c2=1。

    说明书

    一种基于分层模型的无人机目标跟踪方法

    技术领域

    本发明属于遥感图像处理技术领域,具体涉及一种基于分层模型的无人机目标跟
    踪方法。

    背景技术

    无人机在民用和军用领域都有着越来越广泛的应用,因此其发展也引起了各国的
    重视,在众多研究方向中视觉能力成为今后无人机发展和竞争的焦点。而无人机视觉跟踪
    作为无人机视觉的主力军,也成为了一项关键技术。据公开资料表明,其中早在2005年的时
    候美国中央佛罗里达大学计算机视觉实验室就开发出了基于MATLAB的COCOA系统,主要用
    于无人机低空航拍视频图像的目标检测与跟踪处理。目前国内的目标跟踪方法大致可以分
    为以下三种:

    一是基于确定性的跟踪方法。该方法以均值漂移为代表,以目标初始位置为起点,
    沿着密度梯度下降最快的方向移动,得到新的位置,再以新的位置作为起点寻找下一个新
    的位置,算法逐次迭代,直到收敛到目标位置。优点是计算简单,速度快;缺点是很容易陷入
    局部最优,在背景干扰、目标遮挡和目标快速运动等情况下目标跟踪的鲁棒性较差。

    二是基于概率性跟踪方法。该方法以粒子滤波为代表,用一组离散的随机粒子来
    近似系统随机变量的概率密度函数,以样本均值代替积分运算,获得状态最小方差估计。优
    点是重采样的引入会使粒子能够将更大的精力投入在最可能成为目标的粒子周围,防止了
    粒子分散;缺点是可能会引起粒子匮乏,导致目标丢失。

    三是基于智能优化算法的跟踪方法。近年来,受均值漂移算法的启发,有学者另辟
    蹊径,采用其它寻优策略的智能优化算法来解决目标跟踪问题,比较典型的有粒子群优化
    算法,通过提高粒子的适应值,不断更新每个粒子自身的状态来控制粒子群的运动,不断迭
    代获得全局最优位置。优点是搜索速度快且搜索效率高,算法简单;缺点是对于离散的优化
    问题处理不佳,容易陷入局部最优。

    其中粒子滤波策略和粒子群优化搜索策略都需要一定数量的粒子,当粒子数目越
    多,准确率就会越高,但是相对应就会增大搜索时间;当粒子数目较少时,搜索时间减小,但
    有些区域可能无法搜索到,由于粒子群优化策略利用迭代的思想,对准确性的影响较小,但
    是粒子滤波会严重影响到跟踪的准确率。

    由于无人机成像环境的特殊性,在平稳飞行的过程中,可得到较为平稳的视频图
    像,但是由于云台运动以及气流等干扰,无人机在飞行过程中可能受到一定的冲击,在这种
    情况下,所成图像帧间目标可能会产生较大运动偏移。

    从上述分析可以看出,现有基于无人机侦察平台的目标跟踪方法在精度和效率之
    间存在矛盾,难以满足现代战场对实时、高精度目标跟踪的需求。

    发明内容

    本发明以无人机视频图像作为研究对象,建立了基于分层模型的无人机目标跟踪
    方法。该方法首先建立初始帧目标的特征表达,然后分两层对目标进行跟踪:第一层先利用
    简单快速的均值漂移方法对目标进行粗跟踪,得到粗跟踪下的目标位置;第二层将粗跟踪
    下的目标位置作为精跟踪的初始位置,利用改进的粒子群优化的搜索方法对目标进行进一
    步精跟踪。通过建立分层模型,得到目标的准确位置。

    本发明的一种基于分层模型的无人机目标跟踪方法,包括以下几个步骤:

    第一步,建立初始帧目标的特征表达。

    为进行无人机视频目标跟踪,本发明要求在无人机侦察平台捕获到地面目标信息
    后,在视频第一帧中进行目标的获取,该任务可由人工标注完成,也可利用自动检测技术。
    获取目标后利用目标区域的基于核函数加权的归一化空间颜色模型作为该目标的特征表
    达。

    第二步,建立以均值漂移为基础的粗跟踪模型。

    本步骤是目标跟踪的第一层,重点建立了以均值漂移为理论基础的粗跟踪模型,
    通过循环迭代计算出目标帧间偏移量,并由此得到粗跟踪中目标的初步估计结果。

    第三步,建立以改进的粒子群优化为基础的精跟踪模型。

    本步骤是目标跟踪的第二层,重点建立了以改进的粒子群优化为理论基础的精跟
    踪模型,通过在线调整角速度进而对粒子的运动规律进行调整,同时根据目标跟踪情况对
    粒子群参数进行自适应性调整,从而得到目标精跟踪结果。

    本发明具有以下优点:

    1)通过建立均值漂移和改进的粒子群优化分层模型的方法对无人机目标进行跟
    踪,该方法优势互补,在不影响实时效率的情况下提高了跟踪的准确性。

    2)引入运动模型的表达,并通过在线更新角速度的方法对当前运动模型进行实时
    调整,从而使粒子搜索更有方向性,提高效率和准确率。

    3)自适应性调整粒子群优化方法中的惯性参数,在一定程度上有利于跳出局部最
    优,防止目标丢失。

    附图说明

    图1为本发明方法流程图;

    图2分层跟踪模型具体流程图;

    图3为目标转弯运动示意图;

    图4本发明实验结果图。

    具体实施方式

    下面结合附图,对本发明的具体实施方法进行详细说明。

    本发明是一种基于分层模型的无人机目标跟踪方法,通过建立初始帧目标的特征
    表达、建立以均值漂移为基础的粗跟踪模型以及建立以改进的粒子群优化为基础的精跟踪
    模型等步骤,实现无人机视频目标的跟踪。图1给出了方法整体流程图,具体实施方法如图2
    所示,包括以下步骤:

    第一步建立初始帧目标的特征表达

    该步骤首先读入视频序列,手动选取要跟踪的目标区域,对所选取的目标区域建
    立基于核函数加权的归一化空间颜色直方图模型,并将该模型的特征向量作为目标区域的
    特征表达。

    基于核函数加权空间颜色模型的核心思想是当目标进行空间位置和形态的变换
    时,不同位置所受的影响程度不同,靠近目标中心的区域受到外界的影响相对于靠近目标
    边缘的区域较小,因此靠近目标中心的区域所具备的权重应大于靠近目标边缘的区域,本
    发明中核函数使用的是Epanechnikov核函数,如式(1):


    其中,K(·)为Epanechnikov核函数,||r||为像素距离目标中心的距离。

    本发明中使用的颜色空间是RGB颜色空间,其中R、G、B三个通道的取值范围均是
    [0,255],每个通道均有256个颜色等级。直方图是一种概率统计方法,统计出现某种颜色组
    合的次数,因此该RGB空间颜色直方图有2563种组合。鉴于该组合数目较大,本发明中采用
    降维的方法来简化,将R、G、B三通道中[0,255]的取值范围转换到[0,16],整个颜色空间就
    可以用163,即4096种组合来表示。最终将三维向量映射为直方图上的一维向量,(R,G,B)=
    256×R+16×G+B,并将其做归一化处理,该归一化后的直方图即可作为初始化阶段中视频
    初始帧目标的特征表达。

    第二步建立以均值漂移为基础的粗跟踪模型

    均值漂移算法是一种基于密度梯度上升的非参数方法,通过迭代运算找到目标位
    置,实现目标跟踪。均值漂移算法计算简单,迭代速度快,因此在初始帧之后的无人机视频
    图像序列中首先使用均值漂移算法对目标进行粗定位。本步骤主要包括两个环节,首先建
    立了均值漂移在跟踪中的数学表达,然后结合具体无人机视频图像对该算法进行分析。

    1)建立均值漂移在跟踪中的数学表达

    ①建立目标模型

    结合核函数加权空间颜色模型,可将目标模型中第u个特征值估计的概率密度表
    示为:


    其中,fu表示目标模型中第u个特征值估计的概率密度,u为直方图中的颜色索引
    值,n为该目标模型的总像素个数,δ[·]为单位冲激响应,k(·)为Epannechnikov核函数,h
    为带宽,x为图像中目标的中心位置,xi为第i个样本的中心位置,b(xi)为位置xi处的颜色量
    化等级,C为系数。

    由概率密度的积分和等于1,可得到系数C为:


    ②建立候选模型

    类似于目标模型,第y个候选模型第u个特征值估计的概率密度可表示为:


    这里的y是当前帧候选区域中心。其他参数和目标模型相同。

    ③建立相似度度量模型

    本发明中相似度的度量采用Bhattacharyya系数:


    其中:p、q为两个样本,m是颜色分块数,pi、qi分别是在p、q中第i部分的成员数。本
    发明中需要判断目标模型和候选模型的相似度,即可表示为:


    其中,相似度ρ的范围在0到1之间。

    ④建立目标粗跟踪的数学表达

    利用均值漂移的思想通过得到偏移量从而到视频中的目标位置,在这个过程中,
    目标位置的确定主要依靠相似度的度量,即判断各个候选区域与目标区域的相似程度,并
    将相似度最高的视为目标。

    当前帧寻找目标的新位置通常是根据上一帧中目标位置的最佳估计值y0得到的,
    因为一般情况下每帧目标位置并不会有太大变化,所以本发明在y0邻域内寻找。因此,必须
    先计算出在当前帧中位于y0处的候选目标的颜色概率密度pu(y0),u=1,2...m,对ρ(y)在pu
    (y0)处进行泰勒展开,可以得到下式:


    将(4)式代入(7)式可得


    其中,

    因为第一项是个常数,所以只有第二项与相似度,即当前目标估计值有关。设第二
    项为f(y),则


    可以看出第二项即是在当前帧中利用核函数和图像像素加权值计算得到的概率
    密度估计,若想得到相似度最高的位置,也就是均值漂移的偏移量,则需求概率密度下降最
    快的位置,此时


    其中:

    其中:k'(·)表示核函数k(·)的导数。

    通过一定的提取变换可得到:


    当时,可得到概率密度梯度方向最大的位置。


    因此,均值漂移偏移量为:


    2)以均值漂移为基础的粗跟踪

    在无人机视频跟踪过程中,通过1)中推导出的均值漂移偏移量中的各个参数在实
    际视频图像中的表示如下:

    y0:上一帧目标的中心位置。

    y:当前帧候选区域中心位置。

    k():Epannechnikov核函数。

    xi:当前帧候选区域的第i个像素点的位置。

    h:h为目标区域各像素点距离中心位置最远的距离。

    n:为样本点的个数,这里代表的是上一帧目标位置处的像素个数。

    本步骤将步骤一中所建立的基于核函数加权的归一化空间颜色直方图模型作为
    均值漂移过程中目标的特征值,通过不断迭代,得到一定精度下的跟踪匹配结果,这里设置
    迭代次数N=50,当候选区域和目标区域的相似度大于一定阈值T时,跳出迭代,这里设置阈
    值T=0.9。

    第三步建立以改进的粒子群优化为基础的精跟踪模型

    由于无人机空中成像的特殊性,摄像机和地面目标同时运动,并且受天气环境影
    响,当经过均值漂移的粗跟踪后,候选目标不一定就是真实目标。此时需判断目标是否丢
    失,若候选目标与实际目标达到较高的相似度,则目标未丢失,利用精跟踪模型对位置进行
    进一步精确定位;若候选目标与真实目标相似度较低,则目标丢失,利用精跟踪模型对目标
    重新进行搜索定位;若精跟踪后仍未搜索到目标,则再次对目标进行手动选取。

    改进的粒子群优化精跟踪模型主要分为以下几个环节:第一,对粒子群优化中的
    参数进行初始化;第二,建立在线更新的运动模型;第三,建立基于粒子群优化算法的自适
    应性参数调整跟踪模型。

    1)参数初始化

    在粒子群优化跟踪模型中,主要存在以下参数:

    粒子群种群数目:N;

    迭代次数:Tpso;

    粒子0时刻的初始飞行速度:

    粒子0时刻的初始位置:

    惯性参数:w1;

    学习因子:c1、c2;

    随机数:r1、r2;

    合适的参数不仅能提高粒子群优化搜索的效率,更能提高其定位精度。通过选取
    不同的参数对50段无人机视频进行测试,得到最优的参数结果。其中,种群数目N=10,通过
    实验可知10个粒子数目已经可以满足搜索需求,不需要过多粒子耗费时间;迭代次数Tpso=
    30;粒子0时刻的初始飞行速度粒子0时刻的初始位置为以上一帧目标为中心进行
    高斯分布;随机数r1=r2=r;学习因子c1=c2=1;惯性参数主要控制搜索能力,如果w1较大,
    则全局搜索能力较强,若w1较小,则局部搜索能力较强。本步骤通过自适应调整惯性参数w1
    的值来控制粒子群优化的搜索能力。

    2)建立运动模型

    无人机视频中目标的运动可大致归结为匀速直线运动和转弯运动。本步骤通过建
    立一种在线调整角速度的方法灵活转化运动模型,进而提高粒子群优化搜索的准确性和效
    率。

    如图3所示,假设目标在(x,y)处转弯,转弯时的速率为v,速度与x轴方向所成夹角
    为θ,本步骤将拐弯处的状态描述为(x,x',y,y'),其中x',y'表示x,y方向上的速度大小。假
    设转弯角速率为ω,则可根据已知条件得到其与时间t的关系。


    x方向上速度的大?。簒'=vx=vcosθ (17)

    y方向上速度的大?。簓'=vy=vsinθ (18)

    x方向上加速度的大?。簒”=v'x=-ω(t)vsinθ (19)

    y方向上加速度的大?。簓”=v'y=ω(t)vcosθ (20)

    建立转弯模型的状态描述方程为:


    式(21)表示的是速度、加速度与位置、速度的关系,若想得到状态转移方程,即式
    (22)所示的形式:

    Xk+1=FkXk+Wk (22)

    其中,Fk为转弯模型的系统状态转移矩阵,Xk+1,Xk表示k+1时刻和k时刻目标的状
    态,即[x,x',y,y']T,Wk为转弯时的系统状态噪声,本发明设置该噪声为随机值,与t无关。要
    知道转弯时的状态,我们的目标就是求解Fk。

    将式(21)代入式(22),化简得到如下形式


    其中A(t)可以表示为


    首先将A(t)离散化,得到


    对式(23)进行拉普拉斯变换

    sL[Fk(s)]-Fk(0)=AFk(s) (26)

    合并同类项为(sI-A)Fk(s)=Fk(0) (27)

    其中I表示单位矩阵。

    已知Fk(0)=(f1(0),f2(0),f3(0),...,fn(0))T=I (28)

    则Fk(s)=(sI-A)-1Fk(0) (29)

    其中


    要求系统状态矩阵Fk(时域下),即Fk=L-1((sI-A)-1) (31)

    通过拉普拉斯逆变换,可得:


    此时就得到转弯下状态矩阵的表达式。当ω=0时,状态矩阵Fk的表达式化简为:


    通过以上两种数学分析,建立了直线运动和转弯运动情况下的状态矩阵,通过上
    一步均值漂移粗跟踪后,再利用粒子群优化的搜索方法进行精跟踪。因在无人机实际飞行
    转弯运动模型中,在图像坐标系中,转弯率并不是一成不变的,为适应实际需求,需要实时
    在线调整角速度,本发明设置在线调整方式如下:

    首先将路径参数t设置为时间,建立x方向路径与y方向路径与t的关系,拟合从第
    11帧开始,以前10帧目标位置进行拟合得到x1(t)和y1(t),计算不同时刻处的曲率:


    式中,是相对于路径参数
    的导数。

    曲率半径:


    角速度为:ω(t)=v×r(t) (36)

    本发明利用该在线更新的方法,得到较为准确的每一时刻的角速度,从而为粒子
    群优化精跟踪中运动模型的选取和更新奠定了基础。

    3)基于粒子群优化方法的分层性跟踪

    在1)中对粒子群优化中所要利用到的参数进行了最优的初始化之后,就需要按照
    粒子群优化中鸟类寻找食物的思想以一定规则寻找无人机视频中的目标。当达到迭代次数
    上限或者超过0.9的相似度并且大于粗跟踪下所得目标位置的相似度时,停止迭代。本步骤
    主要分为以下几个环节:

    ①计算初始化后的每个粒子的适应度值。在粒子群优化精跟踪中,本步骤仍然从
    采用基于核函数加权的归一化空间颜色直方图的特征表达方式,适应度值同样采用
    Bhattacharyya系数的相似度测量方法。

    ②根据适应度值对粒子进行排序,找出种群最优粒子个体。

    ③在前10帧默认运动模型为直线模型,之后根据在线更新的角速度值进行调整,
    本发明利用运动模型分别对第i个粒子的位置进行更新,并令i自增。若运动模型为直线模
    型,则位置更新公式可以表示为:


    式中,表示粒子i在t时刻的位置,表示粒子i在t时刻的飞行速度。

    若运动模型为转弯模型,则位置更新公式在x、y轴两个方向可以表示为:



    其中:i1代表粒子i在x方向的数据情况,i2代表粒子i在y方向的数据情况。

    ④计算新位置的适应度值。判断新位置是否优于该粒子最佳位置,如果满足条件,
    用新位置替换原来位置,否则,原位置不变;判断新位置是否优于当前最佳位置。如果不满
    足条件,则当前最佳位置不变,否则,用新位置替换当前最佳位置,并调整粒子速度。


    其中:x*表示在当前群体中最佳粒子所处位置,表示粒子的最佳位置,w1为惯性
    参数,c1、c2为学习因子,r为[0,1]之间的随机数。

    在不断迭代的过程中,若适应度不断朝着大的方向趋近,那么此时减小w1,加强局
    部区域的搜索,若适应度值一直没有什么变化或者一直较小,则可能陷入局部最优或者丢
    失目标,此时就要调大w1,加强全局搜索能力,跳出局部最优,找到目标。本发明设置每次调
    节的步长为0.5,最大调节到10,最小调节到0。

    ⑤判断是否满足迭代终止条件。如果满足,则粒子群优化搜索结束,否则执行③。

    本步骤主要对均值漂移后得到的粗结果进行进一步精确定位,在运动模型的基础
    上,自适应调整粒子群优化中的惯性参数,调整全局和局部搜索比重,定位到准确的位置。
    通过均值漂移和粒子群优化的分层模型,最终将步骤④中得到的群体粒子最佳位置x*视作
    无人机视频中目标跟踪的精确结果。

    下面主要以两组序列说明本发明的有效性:

    如图4所示,实线代表粗跟踪结果,而虚线代表精跟踪结果,图(a)主要说明在粗跟
    踪将目标丢失的情况下,加大精跟踪的全局搜索能力,从而准确定位到目标,该序列中使用
    分层模型使跟踪精度提高了14.8%;图(b)主要说明在精跟踪已经达到一定的准确率的情
    况下,利用精跟踪的方法进一步提高定位精度,该序列中使用分层模型使跟踪精度提高了
    1.2%。若只单纯使用均值漂移或粒子群优化的方法,难以适应无人机视频的复杂性,尤其
    对于目标丢失的情况,一旦丢失目标就难以再次准确定位,本发明中的分层模型更能满足
    无人机对于跟踪的需求,粗跟踪与精跟踪相辅相成,在有效防止目标丢失的同时增加了目
    标跟踪的准确性。实验表明,本发明对于无人机视频跟踪准确率的提高具有一定的有效性。

    关于本文
    本文标题:一种基于分层模型的无人机目标跟踪方法.pdf
    链接地址://www.4mum.com.cn/p-6079550.html
    关于我们 - 网站声明 - 网站地图 - 资源地图 - 友情链接 - 网站客服 - 联系我们

    [email protected] 2017-2018 www.4mum.com.cn网站版权所有
    经营许可证编号:粤ICP备17046363号-1 
     


    收起
    展开
  • 四川郎酒股份有限公司获第十二届人民企业社会责任奖年度环保奖 2019-05-13
  • 银保监会新规剑指大企业多头融资和过度融资 2019-05-12
  • 韩国再提4国联合申办世界杯 中国网友无视:我们自己来 2019-05-11
  • 中国人为什么一定要买房? 2019-05-11
  • 十九大精神进校园:风正扬帆当有为 勇做时代弄潮儿 2019-05-10
  • 粽叶飘香幸福邻里——廊坊市举办“我们的节日·端午”主题活动 2019-05-09
  • 太原设禁鸣路段 设备在测试中 2019-05-09
  • 拜耳医药保健有限公司获第十二届人民企业社会责任奖年度企业奖 2019-05-08
  • “港独”没出路!“梁天琦们”该醒醒了 2019-05-07
  • 陈卫平:中国文化内涵包含三方面 文化复兴表现在其中 2019-05-06
  • 人民日报客户端辟谣:“合成军装照”产品请放心使用 2019-05-05
  • 【十九大·理论新视野】为什么要“建设现代化经济体系”?   2019-05-04
  • 聚焦2017年乌鲁木齐市老城区改造提升工程 2019-05-04
  • 【专家谈】上合组织——构建区域命运共同体的有力实践者 2019-05-03
  • 【华商侃车NO.192】 亲!楼市火爆,别忘了买车位啊! 2019-05-03
  • 竞彩稳赚秘法 重庆时时分析手机软件 手机赌博赢了几十万 1.992高赔率奖金pk10app 时时安卓计划免费版 内蒙古时时老五星走势图 pk10猜冠军技巧稳赚 赛车稳赚不亏的注码法 dnf铁马有多稳 体育彩票最晚购买时间 北京pk赛车开结果 北京福彩赛车开奖时间 大小单双怎么买稳赚 幸运飞艇倍投 足球在线直播 时时彩稳赚平刷枝术