发明名称 一种基于空域分布特征的飞鸟群目标跟踪方法
摘要 本发明公开了一种基于空域分布特征的飞鸟群目标跟踪方法,本发明提出的飞鸟群目标跟踪方法基于“探鸟雷达”实时获取的鸟情信息,首先确定当前时刻飞鸟群目标集合的预估中心,并利用隶属度关系建立多个备选的群目标量测集合;然后,提取出每个备选群目标量测集合的空域分布特征,进而估计出每个备选群目标集合与前一时刻群目标集合的空域分布相似度;最后,以相似度最大的群目标量测集合修正群预估中心。本发明基于空域分布特征的群目标跟踪方法适用于密集飞行的飞鸟群目标,不会受到群目标量测缺失或空间位置偏移的影响。
申请公布号 CN103913740B 申请公布日期 2015.11.18
申请号 CN201410152711.0 申请日期 2014.04.16
申请人 中国民航科学技术研究院 发明人 陈唯实;李敬
分类号 G01S13/66(2006.01)I;G01S7/41(2006.01)I 主分类号 G01S13/66(2006.01)I
代理机构 北京永创新实专利事务所 11121 代理人 赵文颖
主权项 一种基于空域分布特征的飞鸟群目标跟踪方法,包括如下步骤:步骤1,群目标中心预估;设飞鸟群目标在前一个扫描周期中的群目标集合为G<sub>k‑1</sub>,中心为O<sub>k‑1</sub>,前一个扫描周期中包含q个量测信息,飞鸟群运动速度为v,则飞鸟群目标在当前扫描周期中的群预估中心<img file="FDA00004919068000000110.GIF" wi="65" he="72" />为:<maths num="0001" id="cmaths0001"><math><![CDATA[<mrow><msubsup><mi>O</mi><mi>k</mi><mo>*</mo></msubsup><mo>=</mo><msub><mi>O</mi><mrow><mi>k</mi><mo>-</mo><mn>1</mn></mrow></msub><mo>+</mo><mi>v</mi><mo>&CenterDot;</mo><mi>&Delta;t</mi><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>1</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000491906800000011.GIF" wi="1178" he="80" /></maths>其中,k为当前扫描周期,k‑1为前一个扫描周期,Δt为雷达的扫描周期步长;O<sub>k‑1</sub>、<img file="FDA0000491906800000012.GIF" wi="97" he="81" />v为均为二维向量;步骤2,建立备选群目标集合;基于步骤1中获得的群目标预估中心<img file="FDA0000491906800000013.GIF" wi="66" he="81" />和当前扫描周期中的量测信息M<sub>k</sub>={m<sub>i</sub>|i=1,2,...,n}        (2)其中:M<sub>k</sub>表示当前扫描周期中的量测信息集合,m<sub>i</sub>表示当前扫描周期中的量测信息,n表示量测信息的数量;建立N个备选群目标集合,每个备选群目标集合<img file="FDA00004919068000000111.GIF" wi="65" he="68" />的中心<img file="FDA00004919068000000112.GIF" wi="66" he="67" />由下式计算<maths num="0002" id="cmaths0002"><math><![CDATA[<mrow><msubsup><mi>O</mi><mi>k</mi><mi>j</mi></msubsup><mo>=</mo><msubsup><mi>O</mi><mi>k</mi><mo>*</mo></msubsup><mo>+</mo><msubsup><mi>&eta;</mi><mi>k</mi><mi>j</mi></msubsup><mo>,</mo><mi>j</mi><mo>=</mo><mn>1,2</mn><mo>,</mo><mo>.</mo><mo>.</mo><mo>.</mo><mo>,</mo><mi>N</mi><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>3</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000491906800000014.GIF" wi="1303" he="80" /></maths>其中:<img file="FDA0000491906800000015.GIF" wi="63" he="79" />表示群目标中心随机偏移向量;每个备选群目标集合中包含的元素需满足以下条件<maths num="0003" id="cmaths0003"><math><![CDATA[<mrow><msubsup><mi>G</mi><mi>k</mi><mi>j</mi></msubsup><mo>=</mo><mo>{</mo><msub><mi>m</mi><mi>i</mi></msub><mo>|</mo><mo>|</mo><msub><mi>m</mi><mi>i</mi></msub><mo>-</mo><msubsup><mi>O</mi><mi>k</mi><mi>j</mi></msubsup><mo>|</mo><mo>&lt;</mo><mi>&theta;</mi><mo>}</mo><mo>,</mo><mi>i</mi><mo>=</mo><mn>1,2</mn><mo>,</mo><mo>.</mo><mo>.</mo><mo>.</mo><mo>,</mo><mi>n</mi><mo>;</mo><mi>j</mi><mo>=</mo><mn>1,2</mn><mo>,</mo><mo>.</mo><mo>.</mo><mo>.</mo><mo>,</mo><mi>N</mi><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>4</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000491906800000016.GIF" wi="1556" he="102" /></maths>式中,θ为距离阈值;获得N个备选群目标集合,每个备选群目标集合<img file="FDA00004919068000000113.GIF" wi="67" he="68" />中包含p<sub>j</sub>个量测信息,且满足p<sub>j</sub>≤n;步骤3,空域分布特征提取;获取<img file="FDA00004919068000000114.GIF" wi="68" he="73" />中任意两个量测信息之间的差矢量模组成的距离集合,记为<maths num="0004" id="cmaths0004"><math><![CDATA[<mrow><mi>L</mi><mo>{</mo><msubsup><mi>G</mi><mi>k</mi><mi>j</mi></msubsup><mo>}</mo><mo>=</mo><msubsup><mrow><mo>{</mo><msubsup><mi>l</mi><mi>s</mi><mrow><mi>k</mi><mo>,</mo><mi>j</mi></mrow></msubsup><mo>}</mo></mrow><mrow><mi>s</mi><mo>=</mo><mn>1</mn></mrow><msubsup><mi>C</mi><msub><mi>p</mi><mi>j</mi></msub><mn>2</mn></msubsup></msubsup><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>5</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000491906800000017.GIF" wi="1207" he="113" /></maths>式中,<img file="FDA0000491906800000018.GIF" wi="632" he="84" />且a<b,s表示空域分布特征集合<img file="FDA0000491906800000019.GIF" wi="68" he="78" />中每个元素的角标,<maths num="0005" id="cmaths0005"><math><![CDATA[<mrow><msubsup><mi>C</mi><msub><mi>p</mi><mi>j</mi></msub><mn>2</mn></msubsup><mo>=</mo><msub><mi>p</mi><mi>j</mi></msub><mo>&CenterDot;</mo><mrow><mo>(</mo><msub><mi>p</mi><mi>j</mi></msub><mo>-</mo><mn>1</mn><mo>)</mo></mrow><mo>/</mo><mn>2</mn></mrow>]]></math><img file="FDA0000491906800000021.GIF" wi="386" he="82" /></maths>表示元素的总数;步骤4,相似度计算;按照公式(6)提取前一个扫描周期中的群目标集合G<sub>k‑1</sub>的空域分布特征集合为:<maths num="0006" id="cmaths0006"><math><![CDATA[<mrow><mi>L</mi><mo>{</mo><msub><mi>G</mi><mrow><mi>k</mi><mo>-</mo><mn>1</mn></mrow></msub><mo>}</mo><mo>=</mo><msubsup><mrow><mo>{</mo><msubsup><mi>l</mi><mi>t</mi><mrow><mi>k</mi><mo>-</mo><mn>1</mn></mrow></msubsup><mo>}</mo></mrow><mrow><mi>t</mi><mo>=</mo><mn>1</mn></mrow><msubsup><mi>C</mi><mi>q</mi><mn>2</mn></msubsup></msubsup><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>6</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA00004919068000000217.GIF" wi="1284" he="114" /></maths>式中:t表示空域分布特征集合G<sub>k‑1</sub>中每个元素的角标,<img file="FDA0000491906800000023.GIF" wi="318" he="79" />表示元素的总数,q为群目标集合G<sub>k‑1</sub>中的量测个数,计算其与当前扫描周期中N个备选群目标集合空域分布特征集合<img file="FDA00004919068000000218.GIF" wi="414" he="87" />内元素两两之间的距离,得到距离矩阵<maths num="0007" id="cmaths0007"><math><![CDATA[<mrow><msubsup><mi>R</mi><mrow><mi>k</mi><mo>-</mo><mn>1</mn><mo>,</mo><mi>k</mi></mrow><mi>j</mi></msubsup><mo>=</mo><mo>[</mo><msub><mi>r</mi><mi>st</mi></msub><mo>]</mo><mo>,</mo><mn>1</mn><mo>&le;</mo><mi>s</mi><mo>&le;</mo><msubsup><mi>C</mi><msub><mi>p</mi><mi>j</mi></msub><mn>2</mn></msubsup><mo>,</mo><mn>1</mn><mo>&le;</mo><mi>t</mi><mo>&le;</mo><msubsup><mi>C</mi><mi>q</mi><mn>2</mn></msubsup><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>7</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000491906800000024.GIF" wi="1424" he="92" /></maths>其中,rst表示<img file="FDA0000491906800000025.GIF" wi="163" he="113" />和<img file="FDA0000491906800000026.GIF" wi="158" he="112" />内元素之间的绝对差;设集合<img file="FDA0000491906800000027.GIF" wi="172" he="113" />和<img file="FDA0000491906800000028.GIF" wi="162" he="113" />之间的距离为<maths num="0008" id="cmaths0008"><math><![CDATA[<mrow><msubsup><mi>&mu;</mi><mrow><mi>k</mi><mo>-</mo><mn>1</mn><mo>,</mo><mi>k</mi></mrow><mi>j</mi></msubsup><mo>=</mo><mi>max</mi><mrow><mo>(</mo><mi>max</mi><mrow><mo>(</mo><msub><mi>R</mi><mi>min</mi></msub><mo>)</mo></mrow><mo>,</mo><mi>max</mi><mrow><mo>(</mo><msub><mi>C</mi><mi>min</mi></msub><mo>)</mo></mrow><mo>)</mo></mrow><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>8</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000491906800000029.GIF" wi="1386" he="83" /></maths>其中,R<sub>min</sub>和C<sub>min</sub>分别为距离矩阵<img file="FDA00004919068000000210.GIF" wi="123" he="74" />中每行和每列中最小元素组成的集合;设相似度为<maths num="0009" id="cmaths0009"><math><![CDATA[<mrow><msubsup><mi>&lambda;</mi><mrow><mi>k</mi><mo>-</mo><mn>1</mn><mo>,</mo><mi>k</mi></mrow><mi>j</mi></msubsup><mo>=</mo><mi>exp</mi><mo>{</mo><mo>-</mo><mi>&xi;</mi><mo>&CenterDot;</mo><msubsup><mi>&mu;</mi><mrow><mi>k</mi><mo>-</mo><mn>1</mn><mo>,</mo><mi>k</mi></mrow><mi>j</mi></msubsup><mo>}</mo><mo>,</mo><mi>j</mi><mo>=</mo><mn>1,2</mn><mo>,</mo><mo>.</mo><mo>.</mo><mo>.</mo><mo>,</mo><mi>N</mi><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>9</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA00004919068000000211.GIF" wi="1406" he="94" /></maths>其中,ξ为调整系数;设λ<sub>max</sub>表示相似度<img file="FDA00004919068000000212.GIF" wi="110" he="82" />的最大值,λ<sub>max</sub>对应的备选群目标集合<img file="FDA00004919068000000213.GIF" wi="104" he="81" />作为当前扫描周期的群目标量测集合,其中心坐标为<img file="FDA00004919068000000214.GIF" wi="138" he="79" />步骤5,群目标中心修正;利用步骤4中获得的相似度系数λ<sub>max</sub>和中心坐标<img file="FDA00004919068000000215.GIF" wi="112" he="81" />修正群预估中心,如下式<maths num="0010" id="cmaths0010"><math><![CDATA[<mrow><msub><mi>O</mi><mi>k</mi></msub><mo>=</mo><msubsup><mi>O</mi><mi>k</mi><mo>*</mo></msubsup><mo>+</mo><msub><mi>&lambda;</mi><mi>max</mi></msub><mo>&CenterDot;</mo><mrow><mo>(</mo><msubsup><mover><mi>G</mi><mo>&OverBar;</mo></mover><mi>k</mi><msub><mi>&lambda;</mi><mi>max</mi></msub></msubsup><mo>-</mo><msubsup><mi>O</mi><mi>k</mi><mo>*</mo></msubsup><mo>)</mo></mrow><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>10</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA00004919068000000216.GIF" wi="1299" he="95" /></maths>最后得到修正后k时刻飞鸟群目标估计中心O<sub>k</sub>,完成k时刻的飞鸟群目标跟踪,返回步骤1,进入k+1时刻,依次循环直至完成飞鸟群目标跟踪。
地址 100028 北京市朝阳区西坝河北里甲24号