发明名称 一种基于距离‑多普勒量测的多站雷达人体跟踪方法
摘要 本发明一种基于距离‑多普勒量测的多站雷达人体跟踪方法,包括以下处理步骤:第一步,对各通道接收的回波进行预处理;第二步,计算瞬时距离‑多普勒像;第三步,目标检测;第四步,量测估计;第五步,目标定位;第六步,目标跟踪。本发明方法其有益效果:本发明在瞬时距离‑多普勒域检测人体目标,提高了目标分辨能力;同时运用距离‑多普勒量测进行目标定位和跟踪,达到人体目标精确跟踪目的。
申请公布号 CN106468772A 申请公布日期 2017.03.01
申请号 CN201610846698.8 申请日期 2016.09.23
申请人 南京特艺科技有限公司 发明人 张洪林;李迎雪;王勇
分类号 G01S13/72(2006.01)I;G01S13/58(2006.01)I;G01S13/42(2006.01)I;G01S7/41(2006.01)I 主分类号 G01S13/72(2006.01)I
代理机构 北京慧泉知识产权代理有限公司 11232 代理人 李娜
主权项 一种基于距离‑多普勒量测的多站雷达人体跟踪方法,技术方案包括以下处理步骤:第一步,对各通道接收的回波进行预处理(a)各通道接收回波为解线频调处理后的信号,对信号做一维傅里叶变换,去除残余视频相位,得到目标慢时间‑距离像;(b)运用指数平均背景相消法进行杂波抑制,得到杂波抑制后的慢时间‑距离像;第二步,计算瞬时距离‑多普勒像(a)对杂波抑制后,某一时刻特定时间窗内的慢时间‑距离像进行Keystone变换,消除目标跨距离单元徙动;(b)对Keystone变换后的数据,沿慢时间维做傅里叶变换,得到瞬时距离‑多普勒像I<sup>k</sup>(n,m),其中,I代表瞬时距离‑多普勒像,n=1,2,…,n为距离维采样,n是I的距离维总采样数,m=1,2,…,m为多普勒维采样,m是I的多普勒维总采样数,k=1,2,3为通道标识;第三步,目标检测运用排序统计恒虚警检测器,对第二步得到的瞬时距离‑多普勒像进行目标检测;排序统计恒虚警检测器输出可表示为:<maths num="0001"><math><![CDATA[<mrow><msubsup><mi>I</mi><mi>d</mi><mi>k</mi></msubsup><mrow><mo>(</mo><mi>n</mi><mo>,</mo><mi>m</mi><mo>)</mo></mrow><mo>=</mo><mfenced open = "{" close = ""><mtable><mtr><mtd><mrow><mn>1</mn><mo>,</mo><msup><mrow><mo>&lsqb;</mo><msup><mi>I</mi><mi>k</mi></msup><mrow><mo>(</mo><mi>n</mi><mo>,</mo><mi>m</mi><mo>)</mo></mrow><mo>&rsqb;</mo></mrow><mn>2</mn></msup><mo>&GreaterEqual;</mo><mi>T</mi><mi>h</mi></mrow></mtd></mtr><mtr><mtd><mrow><mn>0</mn><mo>,</mo><msup><mrow><mo>&lsqb;</mo><msup><mi>I</mi><mi>k</mi></msup><mrow><mo>(</mo><mi>n</mi><mo>,</mo><mi>m</mi><mo>)</mo></mrow><mo>&rsqb;</mo></mrow><mn>2</mn></msup><mo>&lt;</mo><mi>T</mi><mi>h</mi></mrow></mtd></mtr></mtable></mfenced><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>1</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0001119625110000011.GIF" wi="781" he="214" /></maths>其中,Th表示排序统计恒虚警检测器估计的自适应阀值;将检测到的目标索引添加入候选目标量测集合M<sup>k</sup>中,<img file="FDA0001119625110000012.GIF" wi="572" he="71" />第四步,量测估计(a)量测范围收缩人体跟踪仅需要人体的躯干量测,人体肢体产生的量测可视为杂波;人体肢体在距离‑多普勒像上扩展面积通常小于人体躯干,定义一个矩形窗:<maths num="0002"><math><![CDATA[<mrow><mi>G</mi><mrow><mo>(</mo><mi>n</mi><mo>,</mo><mi>m</mi><mo>)</mo></mrow><mo>=</mo><munderover><mi>&Sigma;</mi><mrow><mi>p</mi><mo>=</mo><mi>n</mi><mo>-</mo><mfrac><msub><mi>N</mi><mi>r</mi></msub><mn>2</mn></mfrac></mrow><mrow><mi>n</mi><mo>+</mo><mfrac><msub><mi>N</mi><mi>r</mi></msub><mn>2</mn></mfrac></mrow></munderover><munderover><mi>&Sigma;</mi><mrow><mi>q</mi><mo>=</mo><mi>m</mi><mo>-</mo><mfrac><msub><mi>N</mi><mi>f</mi></msub><mn>2</mn></mfrac></mrow><mrow><mi>m</mi><mo>+</mo><mfrac><msub><mi>N</mi><mi>f</mi></msub><mn>2</mn></mfrac></mrow></munderover><msubsup><mi>I</mi><mi>d</mi><mi>k</mi></msubsup><mrow><mo>(</mo><mi>n</mi><mo>+</mo><mi>p</mi><mo>,</mo><mi>m</mi><mo>+</mo><mi>q</mi><mo>)</mo></mrow><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>2</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0001119625110000021.GIF" wi="990" he="214" /></maths>其中,(n,m)∈M<sup>k</sup>,N<sub>r</sub>和N<sub>f</sub>分别为人体躯干在距离‑多普勒像上的扩展宽度;若:G(n,m)<β(N<sub>r</sub>+1)(N<sub>f</sub>+1)  (3)则以(n,m)为中心的目标扩展较小,可将该目标视为杂波,从M<sup>k</sup>删除量测(n,m);其中,β(0<β<1)是一个设定的标量因子;(b)人体距离‑多普勒量测估计将瞬时距离‑多普勒像中人体躯干能量最强点视为量测估计;基于索引集合M<sup>k</sup>,在瞬时距离‑多普勒像I<sup>k</sup>中寻找能量最大值;若I<sup>k</sup>(n,m),(n,m)∈M<sup>k</sup>,对应于能量最强点,将索引(n,m)添加入量测估计集合<img file="FDA0001119625110000022.GIF" wi="74" he="63" />中,并且从M<sup>k</sup>中删除属于矩形窗W(n+p,m+q),<img file="DA00011196251156796.GIF" wi="431" he="127" /><img file="FDA0001119625110000023.GIF" wi="441" he="127" />的量测;然后搜索下一个量测,直至候选量测集合M<sup>k</sup>为空;第五步,目标定位融合各通道中所有距离‑多普勒量测,将所有满足均方误差条件的位置估计作为目标定位点;第六步,目标跟踪将第五步得到的定位点和速度作为目标跟踪系统的输入,运用卡尔曼滤波器,基于最近邻准则,进行多目标跟踪。
地址 211111 江苏省南京市江宁区秣陵街道秣周东路12号