CN107607904B - 一种基于测角信息的多目标同一性识别方法 - Google Patents
一种基于测角信息的多目标同一性识别方法 Download PDFInfo
- Publication number
- CN107607904B CN107607904B CN201710794563.6A CN201710794563A CN107607904B CN 107607904 B CN107607904 B CN 107607904B CN 201710794563 A CN201710794563 A CN 201710794563A CN 107607904 B CN107607904 B CN 107607904B
- Authority
- CN
- China
- Prior art keywords
- observation station
- target
- observation
- intersection point
- station
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
Landscapes
- Pharmaceuticals Containing Other Organic And Inorganic Compounds (AREA)
- Position Fixing By Use Of Radio Waves (AREA)
Abstract
本发明公开一种基于测角信息的多目标同一性识别方法,包括:对所有观测站观测到的所有目标进行同一性识别,其中,对两个观测站分别观测到的两个目标进行同一性识别包括:建立目标角度测量模型;计算观测站在地心地固坐标系下的坐标;计算两个目标测角视线的单位向量;计算两个目标测角视线之间公共垂线的长度;计算检验统计量的检测门限;对两个目标进行同一性识别。本发明运算量小,准确度高,能够在复杂环境下实现较高的识别精度。
Description
技术领域
本发明涉及测量技术领域。更具体地,涉及一种基于测角信息的多目标同一性识别方法。
背景技术
通过高精度的测角设备在多个观测站对目标进行测角,利用各个测角视线的交点对目标的位置进行解算,在各类军事和民用领域中得到了广泛应用。基于观测站的测角数据进行正确的目标同一性识别是这种交叉定位方法准确性的前提。该问题的本质就是从来自多信息员多个目标的多个测角视线的交叉组合中剔除虚假交点,找出属于同一个目标的正确组合,即目标同一性识别。针对该问题,已有技术方法主要包括最小距离法、最近领域法、联合概率数据关联法、极大似然法、分层关联算法等。最小距离法计算简单,但在复杂环境下关联度不高。最近领域法和联合概率数据关联法在计算过程中,要求必须已经形成目标航迹,或者使用目标的先验信息以获得目标测角信息预测值,关联结果很大程度上取决于目标测角信息预测精度,使用场景受限。极大似然法和分层关联算法准确率较高,但当观测站和目标数较多时,计算量急剧增大,实时处理难度高,工程化实现较为困难。
因此,需要提供一种基于测角信息的多目标同一性识别方法。
发明内容
本发明的目的在于针对基于测角信息的目标同一性识别问题,提供一种基于测角信息的多目标同一性识别方法,以避免现有的最近领域法和联合概率数据关联法对目标航迹信息或先验信息的依赖,克服极大似然法和分层关联算法运算量大的问题,同时能够在复杂环境下实现较高的识别精度,实现多个目标的快速同一性识别和交叉定位。
为达到上述目的,本发明采用下述技术方案:
一种基于测角信息的多目标同一性识别方法,包括:
对所有观测站观测到的所有目标进行同一性识别,其中,对两个观测站分别观测到的两个目标进行同一性识别包括:
第一步、设两个目标观测分别为观测站Si观测到的第mi个目标和观测站Sj观测到的第nj个目标,i、j=1,2,…,K,K为观测站个数,mi=1,2,…,Mi,Mi为观测站Si观测到的目标个数,nj=1,2,…,Nj,Nj为观测站Sj观测到的目标个数;观测站Si和观测站Sj的大地坐标位置分别为(Li,λi,Hi)和(Lj,λj,Hj);设β和ε分别表示在直角坐标系下观测站观测目标的方位角和俯仰角,观测站Si观测第mi个目标的二维角度信息为观测站Sj观测第nj个目标的二维角度信息为
第二步、根据观测站Si和观测站Sj的大地坐标位置(Li,λi,Hi)和(Lj,λj,Hj),计算得到观测站Si和观测站Sj在地心地固坐标系下的坐标和
第三步、分别计算观测站Si观测第mi个目标的二维角度信息的单位视线向量和观测站Sj站观测第nj个目标的二维角度信息的单位视线向量
第四步、设观测站Si观测第mi个目标的视线和观测站Sj观测第nj个目标的视线的公共垂线为公共垂线在观测站Si观测第mi个目标的视线和观测站Sj观测第nj个目标的视线上的垂足分别为和分别计算得到垂足到观测站Si的距离和垂足到观测站Sj的距离并进一步计算得到垂足的坐标和垂足的坐标从而计算得到公共垂线的长度;
第五步、根据观测站Si的方位角误差最大值Δβi和俯仰测角误差最大值Δεi,计算得到垂足的定位误差最大值根据观测站Sj的方位角误差最大值Δβj和俯仰测角误差最大值Δεj,计算得到垂足的定位误差最大值计算得到检测门限
第六步、将公共垂线的长度与检测门限进行对比以进行同一性识别:如果公共垂线的长度小于检测门限则观测站Si观测到的第mi个目标和观测站Sj观测到的第nj个目标为同一个目标;如果公共垂线长度的大于等于检测门限则观测站Si观测到的第mi个目标和观测站Sj观测到的第nj个目标为两个不同的目标。
优选地,第二步中,计算得到观测站Si和观测站Sj在地心地固坐标系下的坐标和的计算公式为:
其中,参数Ci=Eq/(1-e2sin2Li)1/2,参数Cj=Eq/(1-e2sin2Lj)1/2;e为地球椭球第一偏心率,e2=0.0066943799013;Eq为赤道半径。
优选地,第三步中,分别计算观测站Si观测第mi个目标的二维角度信息的单位视线向量和观测站Sj站观测第nj个目标的二维角度信息的单位视线向量的计算公式为:
其中,
优选地,第四步进一步包括:
设观测站Si观测第mi个目标的视线和观测站Sj观测第nj个目标的视线的公共垂线为该垂线在两条视线上的垂足分别为和分别计算得到垂足到观测站Si的距离和垂足到观测站Sj的距离
其中,
进一步计算得到垂足的坐标和垂足的坐标
计算得到公共垂线的长度
优选地,第五步中,
垂足的定位误差最大值的计算公式为:
垂足的定位误差最大值的计算公式为:
本发明的有益效果如下:
本发明所述技术方案通过构造基于测角视线最小间距的检验统计量,并与根据各个观测站的测角精度构造的检测门限进行比较,对虚假的测角视线交点进行剔除并找出属于同一个目标的正确组合,完成对各个目标的同一性识别,克服了对目标航迹信息或先验信息的依赖,运算量小,准确度高,能够在复杂环境下实现较高的识别精度,实现多个目标的快速同一性识别和交叉定位。
附图说明
下面结合附图对本发明的具体实施方式作进一步详细的说明;
图1示出基于测角信息的多目标同一性识别方法的流程图。
图2示出目标角度观测的示意图。
图3示出测角视线公共垂线的示意图。
具体实施方式
为了更清楚地说明本发明,下面结合优选实施例和附图对本发明做进一步的说明。附图中相似的部件以相同的附图标记进行表示。本领域技术人员应当理解,下面所具体描述的内容是说明性的而非限制性的,不应以此限制本发明的保护范围。
本实施例提供的基于测角信息的多目标同一性识别方法,采用基于观测站测角视线最小间距的检验统计量。两个观测站在某一时刻的目标二维测角数据可以在空间中形成两条测角视线。若两个测角数据来源于同一个目标,则两条视线相交于目标所在位置;否则,两条视线在空间中无法相交,此时可以将两条视线的公共垂线段的长度定义为两条视线的最小间距。若最小间距小于某一个门限值,即可认为两条视线来自于同一个目标,即通过同一性识别;反之则未通过同一性确认,两条视线属于非同源目标。
本实施例提供的基于测角信息的多目标同一性识别方法,包括:对所有观测站观测到的所有目标进行同一性识别。
其中,如图1所示,对两个观测站分别观测到的两个目标进行同一性识别包括:
第一步、建立目标角度测量模型:
如图2所示,设基于测角信息进行同一性识别的两个目标观测分别来自观测站Si和观测站Sj(i、j=1,2,…,K,K为观测站个数),观测站Si和观测站Sj的大地坐标位置分别为(Li,λi,Hi)和(Lj,λj,Hj),大地坐标系中的三个坐标变量依次代表大地纬度、大地经度和大地高度。设β和ε分别表示在直角坐标系下观测站观测目标的方位角和俯仰角:观测站Si观测第mi个目标的二维角度信息为(其中和分别表示在直角坐标系下观测站Si观测第mi个目标的方位角和俯仰角),mi=1,2,…,Mi,Mi为观测站Si观测到的目标个数,观测站Sj观测第nj个目标的二维角度信息为nj=1,2,…,Nj,Nj为观测站Sj观测到的目标个数,基于测角信息对观测站Si观测到的第mi个目标和观测站Sj观测到的第nj个目标进行同一性识别。
第二步、计算观测站在地心地固坐标系下的坐标:
根据观测站Si和观测站Sj的大地坐标位置(Li,λi,Hi)和(Lj,λj,Hj),计算得到观测站Si和Sj在地心地固坐标系下的坐标和公式为:
其中,参数Ci=Eq/(1-e2sin2Li)1/2,参数Cj=Eq/(1-e2sin2Lj)1/2;e为地球椭球第一偏心率,e2=0.0066943799013;Eq为赤道半径。
第三步、计算两个目标测角视线的单位向量:
分别计算观测站Si观测第mi个目标的二维角度信息的单位视线向量和观测站Sj站观测第nj个目标的二维角度信息的单位视线向量公式为:
其中,
第四步、计算两个目标测角视线之间公共垂线的长度:
如图3所示,设观测站Si观测第mi个目标的视线和观测站Sj观测第nj个目标的视线的公共垂线为公共垂线在观测站Si观测第mi个目标的视线和观测站Sj观测第nj个目标的视线上的垂足分别为和分别计算得到垂足到观测站Si的距离和垂足到观测站Sj的距离并进一步计算得到垂足的坐标和垂足的坐标从而计算得到公共垂线的长度,具体包括:
设观测站Si观测第mi个目标的视线和观测站Sj观测第nj个目标的视线的公共垂线为该垂线在两条视线上的垂足分别为和分别计算得到垂足到观测站Si的距离和垂足到观测站Sj的距离
其中,
进一步计算得到垂足的坐标和垂足的坐标
计算得到公共垂线的长度
第五步、计算检验统计量的检测门限:
根据观测站Si的方位角误差最大值Δβi和俯仰测角误差最大值Δεi,计算得到在该观测站Si视线垂直方向上的垂足的定位误差最大值 根据观测站Sj的方位角误差最大值Δβj和俯仰测角误差最大值Δεj,计算得到在观测站Sj视线垂直方向上的垂足的定位误差最大值 最终计算得到检测门限
第六步、对两个目标进行同一性识别:
将公共垂线的长度作为检验统计量与检测门限进行对比,对观测站Si观测到的第mi个目标和观测站Sj观测到的第nj个目标进行同一性识别:如果公共垂线的长度小于检测门限则认为两条视线属于同一个目标来源,即通过同一性识别,观测站Si观测到的第mi个目标和观测站Sj观测到的第nj个目标为同一个目标;如果公共垂线的长度大于等于检测门限则认为两条视线不属于同一个目标来源,即未通过同一性识别,观测站Si观测到的第mi个目标和观测站Sj观测到的第nj个目标为两个不同的目标。
显然,本发明的上述实施例仅仅是为清楚地说明本发明所作的举例,而并非是对本发明的实施方式的限定,对于所属领域的普通技术人员来说,在上述说明的基础上还可以做出其它不同形式的变化或变动,这里无法对所有的实施方式予以穷举,凡是属于本发明的技术方案所引伸出的显而易见的变化或变动仍处于本发明的保护范围之列。
Claims (5)
1.一种基于测角信息的多目标同一性识别方法,其特征在于,包括:
对所有观测站观测到的所有目标进行同一性识别,其中,对两个观测站分别观测到的两个目标进行同一性识别包括:
第一步、设两个目标观测分别为观测站Si观测到的第mi个目标和观测站Sj观测到的第nj个目标,i、j=1,2,…,K,K为观测站个数,mi=1,2,…,Mi,Mi为观测站Si观测到的目标个数,nj=1,2,…,Nj,Nj为观测站Sj观测到的目标个数;观测站Si和观测站Sj的大地坐标位置分别为(Li,λi,Hi)和(Lj,λj,Hj);设β和ε分别表示在直角坐标系下观测站观测目标的方位角和俯仰角,观测站Si观测第mi个目标的二维角度信息为观测站Sj观测第nj个目标的二维角度信息为
第二步、根据观测站Si和观测站Sj的大地坐标位置(Li,λi,Hi)和(Lj,λj,Hj),计算得到观测站Si和观测站Sj在地心地固坐标系下的坐标和
第三步、分别计算观测站Si观测第mi个目标的二维角度信息的单位视线向量和观测站Sj站观测第nj个目标的二维角度信息的单位视线向量
第四步、设观测站Si观测第mi个目标的视线和观测站Sj观测第nj个目标的视线的公共垂线为公共垂线在观测站Si观测第mi个目标的视线和观测站Sj观测第nj个目标的视线上的垂足分别为和分别计算得到垂足到观测站Si的距离和垂足到观测站Sj的距离并进一步计算得到垂足的坐标和垂足的坐标从而计算得到公共垂线的长度;
第五步、根据观测站Si的方位角误差最大值Δβi和俯仰测角误差最大值Δεi,计算得到垂足的定位误差最大值根据观测站Sj的方位角误差最大值Δβj和俯仰测角误差最大值Δεj,计算得到垂足的定位误差最大值计算得到检测门限
第六步、将公共垂线的长度与检测门限进行对比以进行同一性识别:如果公共垂线的长度小于检测门限则观测站Si观测到的第mi个目标和观测站Sj观测到的第nj个目标为同一个目标;如果公共垂线长度的大于等于检测门限则观测站Si观测到的第mi个目标和观测站Sj观测到的第nj个目标为两个不同的目标。
2.根据权利要求1所述的基于测角信息的多目标同一性识别方法,其特征在于,第二步中,计算得到观测站Si和观测站Sj在地心地固坐标系下的坐标和的计算公式为:
其中,参数Ci=Eq/(1-e2sin2Li)1/2,参数Cj=Eq/(1-e2sin2Lj)1/2;e为地球椭球第一偏心率,e2=0.0066943799013;Eq为赤道半径。
3.根据权利要求2所述的基于测角信息的多目标同一性识别方法,其特征在于,第三步中,分别计算观测站Si观测第mi个目标的二维角度信息的单位视线向量和观测站Sj站观测第nj个目标的二维角度信息的单位视线向量的计算公式为:
其中,
4.根据权利要求3所述的基于测角信息的多目标同一性识别方法,其特征在于,第四步进一步包括:
设观测站Si观测第mi个目标的视线和观测站Sj观测第nj个目标的视线的公共垂线为该垂线在两条视线上的垂足分别为和分别计算得到垂足到观测站Si的距离和垂足到观测站Sj的距离
其中,
进一步计算得到垂足的坐标和垂足的坐标
计算得到公共垂线的长度
5.根据权利要求4所述的基于测角信息的多目标同一性识别方法,其特征在于,第五步中,
垂足的定位误差最大值的计算公式为:
垂足的定位误差最大值的计算公式为:
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710794563.6A CN107607904B (zh) | 2017-09-06 | 2017-09-06 | 一种基于测角信息的多目标同一性识别方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710794563.6A CN107607904B (zh) | 2017-09-06 | 2017-09-06 | 一种基于测角信息的多目标同一性识别方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN107607904A CN107607904A (zh) | 2018-01-19 |
CN107607904B true CN107607904B (zh) | 2019-07-19 |
Family
ID=61057379
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201710794563.6A Active CN107607904B (zh) | 2017-09-06 | 2017-09-06 | 一种基于测角信息的多目标同一性识别方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN107607904B (zh) |
Families Citing this family (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110232717A (zh) * | 2019-06-10 | 2019-09-13 | 北京壹氢科技有限公司 | 一种适用于多对多目标识别的目标同一性识别方法 |
CN110501685B (zh) * | 2019-08-23 | 2021-06-04 | 北京电子工程总体研究所 | 一种基于雷达信号相位补偿的多帧相参积累方法 |
CN112989865B (zh) * | 2019-12-02 | 2023-05-30 | 山东浪潮科学研究院有限公司 | 一种基于人头姿态判断的人群注意力焦点判断方法 |
CN113791407B (zh) * | 2021-07-28 | 2023-10-13 | 中国人民解放军国防科技大学 | 一种基于椭球模型的双站测向交叉定位方法 |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
KR101282718B1 (ko) * | 2012-12-28 | 2013-07-05 | 한국항공우주연구원 | 선형배열 영상센서와 자세제어 센서 간의 절대 오정렬 보정방법 |
CN103901458A (zh) * | 2014-03-14 | 2014-07-02 | 中国测绘科学研究院 | 一种动态定位系统及动态定位方法 |
CN104391272A (zh) * | 2014-10-27 | 2015-03-04 | 深圳大学 | 利用测向数据进行目标定位的方法与系统 |
CN106772237A (zh) * | 2016-11-18 | 2017-05-31 | 中国船舶重工集团公司第七二四研究所 | 长基线时差定位干涉仪测向辅助解模糊方法 |
CN107289918A (zh) * | 2017-07-24 | 2017-10-24 | 北京电子工程总体研究所 | 一种便携式寻北测距组合测量设备 |
-
2017
- 2017-09-06 CN CN201710794563.6A patent/CN107607904B/zh active Active
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
KR101282718B1 (ko) * | 2012-12-28 | 2013-07-05 | 한국항공우주연구원 | 선형배열 영상센서와 자세제어 센서 간의 절대 오정렬 보정방법 |
CN103901458A (zh) * | 2014-03-14 | 2014-07-02 | 中国测绘科学研究院 | 一种动态定位系统及动态定位方法 |
CN104391272A (zh) * | 2014-10-27 | 2015-03-04 | 深圳大学 | 利用测向数据进行目标定位的方法与系统 |
CN106772237A (zh) * | 2016-11-18 | 2017-05-31 | 中国船舶重工集团公司第七二四研究所 | 长基线时差定位干涉仪测向辅助解模糊方法 |
CN107289918A (zh) * | 2017-07-24 | 2017-10-24 | 北京电子工程总体研究所 | 一种便携式寻北测距组合测量设备 |
Also Published As
Publication number | Publication date |
---|---|
CN107607904A (zh) | 2018-01-19 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN107607904B (zh) | 一种基于测角信息的多目标同一性识别方法 | |
CN104729485B (zh) | 一种基于车载全景影像与街景地图匹配的视觉定位方法 | |
CN104299244B (zh) | 基于单目相机的障碍物检测方法及装置 | |
Liu et al. | Automated matching of multiple terrestrial laser scans for stem mapping without the use of artificial references | |
CN108061889A (zh) | Ais与雷达角度系统偏差的关联方法 | |
CN104077809B (zh) | 基于结构性线条的视觉slam方法 | |
CN106295512B (zh) | 基于标识的多纠正线室内视觉数据库构建方法以及室内定位方法 | |
CN103697855B (zh) | 一种基于海天线检测的船体水平姿态测量方法 | |
US11113896B2 (en) | Geophysical sensor positioning system | |
CN101914881A (zh) | 一种高速铁路基桩控制网cpiii的快速测量方法 | |
CN105571636B (zh) | 一种用于定位目标的方法及测量设备 | |
CN106092059A (zh) | 一种基于多点拟合的结构物平面位移监测方法 | |
CN109345617B (zh) | 一种基于长条带多站点云的链式高精度拼接与平差方法 | |
CN112085003A (zh) | 公共场所异常行为自动识别方法及装置、摄像机设备 | |
Chu et al. | GPS refinement and camera orientation estimation from a single image and a 2D map | |
CN106370160A (zh) | 一种机器人室内定位系统和方法 | |
CN108153979A (zh) | 基于InSAR的形变信息提取方法、终端及存储介质 | |
Yao et al. | Integrity monitoring for Bluetooth low energy beacons RSSI based indoor positioning | |
CN108445505A (zh) | 线环境下基于激光雷达的特征显著性检测方法 | |
CN109459030A (zh) | 一种基于地标的行人定位修正方法及系统 | |
CN112800938A (zh) | 无人驾驶车辆检测侧面落石发生的方法及装置 | |
CN103954220A (zh) | 撞桥试验中船体运动状态数字图像测量方法 | |
CN106651957A (zh) | 基于模板的单目视觉目标空间定位方法 | |
CN104239678A (zh) | 一种实现干涉仪测向定位的方法和装置 | |
CN105547282A (zh) | 一种用于移动定位目标的方法及测量设备 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |