本发明公开一种基于测角信息的多目标同一性识别方法，包括：对所有观测站观测到的所有目标进行同一性识别，其中，对两个观测站分别观测到的两个目标进行同一性识别包括：建立目标角度测量模型；计算观测站在地心地固坐标系下的坐标；计算两个目标测角视线的单位向量；计算两个目标测角视线之间公共垂线的长度；计算检验统计量的检测门限；对两个目标进行同一性识别。本发明运算量小，准确度高，能够在复杂环境下实现较高的识别精度。
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所述的基于测角信息的多目标同一性识别方法，其特征在于，第五步中，垂足 的定位误差最大值 的计算公式为：
一种基于测角信息的多目标同一性识别方法
技术领域
[0001] 本发明涉及测量技术领域。更具体地，涉及一种基于测角信息的多目标同一性识别方法。
背景技术
[0002] 通过高精度的测角设备在多个观测站对目标进行测角，利用各个测角视线的交点对目标的位置进行解算，在各类军事和民用领域中得到了广泛应用。基于观测站的测角数据进行正确的目标同一性识别是这种交叉定位方法准确性的前提。该问题的本质就是从来自多信息员多个目标的多个测角视线的交叉组合中剔除虚假交点，找出属于同一个目标的正确组合，即目标同一性识别。针对该问题，已有技术方法主要包括最小距离法、最近领域法、联合概率数据关联法、极大似然法、分层关联算法等。最小距离法计算简单，但在复杂环境下关联度不高。最近领域法和联合概率数据关联法在计算过程中，要求必须已经形成目标航迹，或者使用目标的先验信息以获得目标测角信息预测值，关联结果很大程度上取决于目标测角信息预测精度，使用场景受限。极大似然法和分层关联算法准确率较高，但当观测站和目标数较多时，计算量急剧增大，实时处理难度高，工程化实现较为困难。
[0003] 因此，需要提供一种基于测角信息的多目标同一性识别方法。
发明内容
[0004] 本发明的目的在于针对基于测角信息的目标同一性识别问题，提供一种基于测角信息的多目标同一性识别方法，以避免现有的最近领域法和联合概率数据关联法对目标航迹信息或先验信息的依赖，克服极大似然法和分层关联算法运算量大的问题，同时能够在复杂环境下实现较高的识别精度，实现多个目标的快速同一性识别和交叉定位。
[0005] 为达到上述目的，本发明采用下述技术方案：
[0006] 一种基于测角信息的多目标同一性识别方法，包括：
[0007] 对所有观测站观测到的所有目标进行同一性识别，其中，对两个观测站分别观测到的两个目标进行同一性识别包括：
[0008] 第一步、设两个目标观测分别为观测站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个目标的二维角度信息为
[0009] 第二步、根据观测站Si和观测站Sj的大地坐标位置(Li,λi,Hi)和(Lj,λj,Hj)，计算得到观测站Si和观测站Sj在地心地固坐标系下的坐标 和
[0010] 第三步、分别计算观测站Si观测第mi个目标的二维角度信息 的单位视线向量 和观测站Sj站观测第nj个目标的二维角度信息 的单位视线向量
[0011] 第四步、设观测站Si观测第mi个目标的视线和观测站Sj观测第nj个目标的视线的公共垂线为 公共垂线 在观测站Si观测第mi个目标的视线和观测站Sj观测第nj个目标的视线上的垂足分别为 和 分别计算得到垂足 到观测站Si的距离 和垂足 到观测站Sj的距离 并进一步计算得到垂足 的坐标 和
垂足 的坐标 从而计算得到公共垂线 的长度；
[0012] 第五步、根据观测站Si的方位角误差最大值Δβi和俯仰测角误差最大值Δεi，计算得到垂足 的定位误差最大值 根据观测站Sj的方位角误差最大值Δβj和俯仰测角误差最大值Δεj，计算得到垂足 的定位误差最大值 计算得到检测门限[0013] 第六步、将公共垂线 的长度与检测门限 进行对比以进行同一性识别：如果公共垂线 的长度小于检测门限 则观测站Si观测到的第mi个目标和观测站Sj观测到的第nj个目标为同一个目标；如果公共垂线 长度的大于等于检测门限 则观测站Si观测到的第mi个目标和观测站Sj观测到的第nj个目标为两个不同的目标。
[0014] 优选地，第二步中，计算得到观测站Si和观测站Sj在地心地固坐标系下的坐标和 的计算公式为：
[0015]
[0016]
[0017] 其中，参数Ci＝Eq/(1-e2sin2Li)1/2，参数Cj＝Eq/(1-e2sin2Lj)1/2；e为地球椭球第一偏心率，e2＝0.0066943799013；Eq为赤道半径。
[0018] 优选地，第三步中，分别计算观测站Si观测第mi个目标的二维角度信息 的单位视线向量 和观测站Sj站观测第nj个目标的二维角度信息 的单位视线向量的计算公式为：
[0019]
[0020]
[0021] 其中，
[0022]
[0023]
[0024] 优选地，第四步进一步包括：
[0025] 设观测站Si观测第mi个目标的视线和观测站Sj观测第nj个目标的视线的公共垂线为 该垂线在两条视线上的垂足分别为 和 分别计算得到垂足 到观测站Si的距离 和垂足 到观测站Sj的距离
[0026]
[0027]
[0028] 其中，
[0029]
[0030]
[0031]
[0032]
[0033] 进一步计算得到垂足 的坐标 和垂足 的坐标
[0034]
[0035]
[0036] 计算得到公共垂线 的长度
[0037]
[0038] 优选地，第五步中，
[0039] 垂足 的定位误差最大值 的计算公式为：
[0040]
[0041] 垂足 的定位误差最大值 的计算公式为：
[0042]
[0043] 本发明的有益效果如下：
[0044] 本发明所述技术方案通过构造基于测角视线最小间距的检验统计量，并与根据各个观测站的测角精度构造的检测门限进行比较，对虚假的测角视线交点进行剔除并找出属于同一个目标的正确组合，完成对各个目标的同一性识别，克服了对目标航迹信息或先验信息的依赖，运算量小，准确度高，能够在复杂环境下实现较高的识别精度，实现多个目标的快速同一性识别和交叉定位。
附图说明
[0045] 下面结合附图对本发明的具体实施方式作进一步详细的说明；
[0046] 图1示出基于测角信息的多目标同一性识别方法的流程图。
[0047] 图2示出目标角度观测的示意图。
[0048] 图3示出测角视线公共垂线的示意图。
具体实施方式
[0049] 为了更清楚地说明本发明，下面结合优选实施例和附图对本发明做进一步的说明。附图中相似的部件以相同的附图标记进行表示。本领域技术人员应当理解，下面所具体描述的内容是说明性的而非限制性的，不应以此限制本发明的保护范围。
[0050] 本实施例提供的基于测角信息的多目标同一性识别方法，采用基于观测站测角视线最小间距的检验统计量。两个观测站在某一时刻的目标二维测角数据可以在空间中形成两条测角视线。若两个测角数据来源于同一个目标，则两条视线相交于目标所在位置；否则，两条视线在空间中无法相交，此时可以将两条视线的公共垂线段的长度定义为两条视线的最小间距。若最小间距小于某一个门限值，即可认为两条视线来自于同一个目标，即通过同一性识别；反之则未通过同一性确认，两条视线属于非同源目标。
[0051] 本实施例提供的基于测角信息的多目标同一性识别方法，包括：对所有观测站观测到的所有目标进行同一性识别。
[0052] 其中，如图1所示，对两个观测站分别观测到的两个目标进行同一性识别包括：
[0053] 第一步、建立目标角度测量模型：
[0054] 如图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个目标进行同一性识别。
[0055] 第二步、计算观测站在地心地固坐标系下的坐标：
[0056] 根据观测站Si和观测站Sj的大地坐标位置(Li,λi,Hi)和(Lj,λj,Hj)，计算得到观测站Si和Sj在地心地固坐标系下的坐标 和 公式为：
[0057]
[0058]
[0059] 其中，参数Ci＝Eq/(1-e2sin2Li)1/2，参数Cj＝Eq/(1-e2sin2Lj)1/2；e为地球椭球第一偏心率，e2＝0.0066943799013；Eq为赤道半径。
[0060] 第三步、计算两个目标测角视线的单位向量：
[0061] 分别计算观测站Si观测第mi个目标的二维角度信息 的单位视线向量和观测站Sj站观测第nj个目标的二维角度信息 的单位视线向量 公式为：
[0062]
[0063]
[0064] 其中，
[0065]
[0066]
[0067] 第四步、计算两个目标测角视线之间公共垂线的长度：
[0068] 如图3所示，设观测站Si观测第mi个目标的视线和观测站Sj观测第nj个目标的视线的公共垂线为 公共垂线 在观测站Si观测第mi个目标的视线和观测站Sj观测第nj个目标的视线上的垂足分别为 和 分别计算得到垂足 到观测站Si的距离 和垂足 到观测站Sj的距离 并进一步计算得到垂足 的坐标 和
垂足 的坐标 从而计算得到公共垂线 的长度，具体包括：
[0069] 设观测站Si观测第mi个目标的视线和观测站Sj观测第nj个目标的视线的公共垂线为 该垂线在两条视线上的垂足分别为 和 分别计算得到垂足 到观测站Si的距离 和垂足 到观测站Sj的距离
[0070]
[0071]
[0072] 其中，
[0073]
[0074]
[0075]
[0076]
[0077] 进一步计算得到垂足 的坐标 和垂足 的坐标
[0078]
[0079]
[0080] 计算得到公共垂线 的长度
[0081]
[0082] 第五步、计算检验统计量的检测门限：
[0083] 根据观测站Si的方位角误差最大值Δβi和俯仰测角误差最大值Δεi，计算得到在该 观 测 站 S i 视 线 垂直 方 向 上 的 垂 足 的 定 位 误 差 最 大 值根据观测站Sj的方位角误差最大值Δβj和俯仰测角误差最大值Δεj，计算得到在观测站Sj视线垂直方向上的垂足 的定位误差最大值最终计算得到检测门限
[0084] 第六步、对两个目标进行同一性识别：
[0085] 将公共垂线 的长度作为检验统计量与检测门限 进行对比，对观测站Si观测到的第mi个目标和观测站Sj观测到的第nj个目标进行同一性识别：如果公共垂线 的长度小于检测门限 则认为两条视线属于同一个目标来源，即通过同一性识别，观测站Si观测到的第mi个目标和观测站Sj观测到的第nj个目标为同一个目标；如果公共垂线 的长度大于等于检测门限 则认为两条视线不属于同一个目标来源，即未通过同一性识别，观测站Si观测到的第mi个目标和观测站Sj观测到的第nj个目标为两个不同的目标。
[0086] 显然，本发明的上述实施例仅仅是为清楚地说明本发明所作的举例，而并非是对本发明的实施方式的限定，对于所属领域的普通技术人员来说，在上述说明的基础上还可以做出其它不同形式的变化或变动，这里无法对所有的实施方式予以穷举，凡是属于本发明的技术方案所引伸出的显而易见的变化或变动仍处于本发明的保护范围之列。
