CN112254719B - 双坐标参数空间目标定位方法 - Google Patents
双坐标参数空间目标定位方法 Download PDFInfo
- Publication number
- CN112254719B CN112254719B CN202011042966.3A CN202011042966A CN112254719B CN 112254719 B CN112254719 B CN 112254719B CN 202011042966 A CN202011042966 A CN 202011042966A CN 112254719 B CN112254719 B CN 112254719B
- Authority
- CN
- China
- Prior art keywords
- target
- coordinate
- station
- latitude
- value
- 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
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01C—MEASURING DISTANCES, LEVELS OR BEARINGS; SURVEYING; NAVIGATION; GYROSCOPIC INSTRUMENTS; PHOTOGRAMMETRY OR VIDEOGRAMMETRY
- G01C21/00—Navigation; Navigational instruments not provided for in groups G01C1/00 - G01C19/00
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01C—MEASURING DISTANCES, LEVELS OR BEARINGS; SURVEYING; NAVIGATION; GYROSCOPIC INSTRUMENTS; PHOTOGRAMMETRY OR VIDEOGRAMMETRY
- G01C21/00—Navigation; Navigational instruments not provided for in groups G01C1/00 - G01C19/00
- G01C21/02—Navigation; Navigational instruments not provided for in groups G01C1/00 - G01C19/00 by astronomical means
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02A—TECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE
- Y02A90/00—Technologies having an indirect contribution to adaptation to climate change
- Y02A90/10—Information and communication technologies [ICT] supporting adaptation to climate change, e.g. for weather forecasting or climate simulation
Landscapes
- Engineering & Computer Science (AREA)
- Radar, Positioning & Navigation (AREA)
- Remote Sensing (AREA)
- Physics & Mathematics (AREA)
- Automation & Control Theory (AREA)
- General Physics & Mathematics (AREA)
- Astronomy & Astrophysics (AREA)
- Position Fixing By Use Of Radio Waves (AREA)
- Radar Systems Or Details Thereof (AREA)
Abstract
本发明公开的一种双坐标多参数的迭代定位方法,属于目标定位技术领域。本发明通过以下技术方案实现:首先,以雷达为站点,获取目标的站点球坐标和WGS‑84坐标组成的双坐标参数:目标在WGS‑84坐标系的经度L、高度H、目标到站点的距离r,站点位置固定为站点在WGS‑84坐标系的经度l、纬度m、高度h,采用初始值计算模块,根据目标双坐标参数和站点位置计算目标纬度初值M01、M02;然后,采用精确值计算模块,根据目标双坐标参数与纬度初值M01和站点位置,通过循环迭代计算目标纬度精确值M11;以同样的方式根据目标双坐标参数与纬度初值M02和站点位置,循环迭代计算目标纬度精确值M12,获取精确的目标定位结果。
Description
技术领域
本发明涉及一种航天器目标监视系统坐标定位方法,尤其涉及一种与地面监控网通信的双坐标多参数的迭代定位方法。
背景技术:
目前日益增多的太空垃圾给在轨运行的航天器带来了严重的安全隐患。为了让航天器安全运行,那就需要完善空间目标监视系统,对运行在航天器附近的目标进行跟踪监视,对可能给航天器运行造成威胁的目标进行处理或者航天器提前做好躲避部署。空间目标是在宇宙空间运行的航天器和空间碎片,空间目标监视系统对空间目标进行探测跟踪、定轨预报、识别编目、侦收分析的情报获取系统。它具有监测非合作目标的能力。其主要功用是:获取空间目标的信息,进行识别、分类、编目,建立目标数据库;空间目标监视系统的任务是对重要空间目标进行精确探测和跟踪,确定可能对航天系统构成威胁的航天器的任务,尺寸、形状和轨道参数等重要目标特性;对目标特性数据进行归类和分发。空间目标监视系统的任务是对重要空间目标进行精确探测和跟踪,确定可能对航天系统构成威胁的航天器的任务,尺寸、形状和轨道参数等重要目标特性;对目标特性数据进行归类和分发。空间目标主要指卫星,包括工作的卫星和不工作的卫星,同时也包括各种空间碎片,如进入空间轨道的助推火箭、保护罩和其他物体,还包括进入地球外层空间的各种宇宙飞行物,如彗星和小行星。目前国际上使用的空间目标监视系统有地面监视系统和天基监视系统,其中地面上典型的有贝克-努恩系统、地基光电深空探测系统(GEODSS)、Maui空间监视站(MSSS)。天基监视系统主要有天基可视计划(SBV)、天基空间监视系统(SBSS)、空间跟踪监视系统(STSS)等。但是这些系统由于距离目标比较远,导致定位精度达到几十米甚至上百米。为了消除太空垃圾带来的安全隐患,在轨航天器需要发展自己的目标监视系统,使其具有跟踪定位目标的能力。现在研究比较多的空间目标定位技术是单星定位和多星交叉定位。但是单星定位收敛时间长、定位精度低,收敛时间需要达到几百秒甚至几千秒才能达到百米的定位精度,多星定位系统中卫星之间的调度比较复杂。空间目标监视系统由数据处理指挥中心(简称系统中心)与若干监测台站(含星载、机载和船载监测系统)组成,它包括探测系统、信息处理系统、通信系统、时间统一系统等4个基本部分。系统中心主要对各监测台站的测量信息进行汇集、处理、分析、存储、发送,提供有关部门使用,并对各台站实施指挥管理,它是全系统的中枢。监测台站主要通过探测系统直接获取空间目标的信息,并进行初步处理。预测空间物体的轨道,对可能发生的碰撞和对己方空间系统的攻击进行告警等。GPS的控制部分由分布在全球的若干跟踪站组成的监控系统所构成,根据其作用的不同,这些跟踪站又被分为主控站、监控站和注入站。主控站的作用是根据各监控站对GPS的观测数据,计算出卫星的星历和卫星钟的改正参数等,并将这些数据通过注入站注入到卫星中去。同时它还对卫星进行控制,向卫星发布指令,当工作卫星出现故障时,调度备用卫星,替代失效的工作卫星工作,另外主控站也具有监控站的功能。监控站的作用是接收卫星信号,监测卫星的工作状态,注入站的作用是将主控站计算出的卫星星历和卫星钟的改正参数等注入到卫星中去。在L1和L2上又分别调制多种信号这些信号主要有C/A码,C/A码又称为粗捕获码它被调制在L1载波上是1MHz的伪随机噪声码PRN码,由于每颗卫星的C/A码都不一样因此经常用它们的PRN信号来区分。C/A码是普通用户用以测定观测站到卫星间距离的一种主要信号。P码又称为精码,它被调制在L1和L2载波上是10MHz的伪随机噪声码。在GPS定位中经常采用下列观测值中的一种或几种进行数据处理以确定出待定点的坐标或待定点之间的基线向量。在进行GPS定位时除了大量地使用上面的观测值进行数据处理以外,还经常使用由上面的观测值通过某些组合而形成的一些特殊观测值。在利用GPS对机动目标进行定位时会受到各种各样因素的影响。影响定位精度的因素可分为以下四大类,与GPS卫星有关的因素包括人为降低,普通用户利用GPS进行导航定位精度。在进行GPS定位时计算在某时刻GPS卫星位置所需的卫星轨道参数是通过各种类型的星历提供的,但不论采用哪种类型的星历所计算出的卫星位置都会与其真实位置有所差异,这就是所谓的卫星星历误差。卫星钟差是GPS卫星上所安装的原子钟的钟面时间与GPS标准时间之间的误差。卫星信号发射天线相位中心偏差是GPS卫星上信号发射天线的标称相位中心与其真实相位中心之间的差异;与传播途径有关的因素,由于地球周围的电离层对电磁波的折射效应使得GPS信号的传播速度发生变化这种变化称为电离层延迟。电磁波所受电离层折射的影响与电磁波的频率以及电磁波传播途径上电子总含量有关。由于地球周围的对流层对电磁波的折射效应使得GPS信号的传播速度发生变化这种变化称为对流层延迟。电磁波所受对流层折射的影响与电磁波传播途径上的温度、湿度和气压有关。由于接收机周围环境的影响使得接收机所接收到的卫星信号中还包含有各种反射和折射信号的影响,这就是所谓的多路径效应。与接收机有关的因素,接收机钟差是GPS接收机所使用的钟的钟面时间与GPS标准时之间的差异,接收机天线相位中心偏差是GPS接收机天线的标称相位中心与其真实的相位中心之间的差异。接收机软件和硬件造成的误差,在进行GPS定位时定位结果还会受到诸如处理与控制软件和硬件等的影响。其他因素,GPS控制部分人为或计算机造成的影响,由于GPS控制部分的问题或用户在进行数据处理时引入的误差,数据处理软件的算法不完善对定位结果的影响等。由于GPS存在100米的误差,在相当多的应用中此误差难以接受。
目标探测的位置主要为目标在站点球坐标系的距离、方位角、俯仰角。站点球坐标系的距离为目标与站点之间的直线距离,站点球坐标系的方位角为目标到站点的射线与站点的正北方向之间的夹角,即站点的正北方向绕站点顺时针旋转的角度,站点球坐标系的俯仰角为目标相对于站点所在直线与站点水平面之间的夹角,向上为正。在目标监视跟踪领域,通常需要将站点球坐标转换为WGS-84坐标(经度、纬度、高度)。WGS-84坐标(WorldGeodeticSystem一1984CoordinateSystem)是一种国际上采用的地心坐标系。坐标原点为地球质心,其地心空间直角坐标系的Z轴指向BIH(国际时间服务机构)1984.O定义的协议地球极(CTP)方向,X轴指向BIH1984.0的零子午面和CTP赤道的交点,Y轴与Z轴、X轴垂直构成右手坐标系,称为1984年世界大地坐标系统。GPS卫星定位测量是用三维地心坐标系——WGS-84坐标系为依据来测定和表示总的空间位置,它既可以用地心空间坐标系(X,Y,Z)表示,也可以用椭球大地坐标系(M,L,H)表示。目前,已有站点球坐标与WGS-84坐标之间的相互转换公式,即利用站点自己的位置,一般为WGS-84坐标,可以将目标在站点球坐标系的位置转为目标在WGS-84坐标系的位置,也可将目标在WGS-84坐标系的位置转为站点球坐标系的位置。这种转换方法的输入参数要么都是站点球坐标或都是WGS-84坐标。但实际中有一种情况:输入的是两个坐标系下的组合坐标,即输入目标在WGS-84坐标系的经度L、高度H、目标到站点的距离r的这种组合坐标,不知道目标在站点球坐标系或目标在WGS-84坐标系的全参数,针对这种包含上述两个坐标系的组合坐标的目标定位,目前在公开文献尚未查到解决方法。
发明内容
为了解决这种组合坐标(目标在WGS-84坐标系的经度L、高度H、目标到站点的距离r)的定位问题,本发明的目的是针对目前空间目标定位技术对目标定位精度低、收敛时间长等问题,提供一种可以提高定位精度,能够更为精确计算空间目标的双坐标参数空间目标定位方法,解决空间目标在站点球坐标系或目标在WGS-84坐标系的全参数的目标定位精度不高的问题。
为了实现上述目的,本发明提出了一种双坐标参数空间目标定位方法,其特征在于包括如下所述的步骤:首先,以雷达侦察站为站点,获取目标的站点球坐标和WGS-84坐标组成的双坐标参数:目标在WGS-84坐标系的经度L、高度H、目标在到站点的距离r,站点位置固定为站点在WGS-84坐标系的经度l、纬度m、高度h;根据上述目标参数和站点位置采用初始值计算模块分别计算出目标的纬度初值M01、M02;然后,根据目标双坐标参数与纬度初值M01、M02和站点位置采用精确值计算模块进行循环迭代计算,分别计算出目标的纬度精确值M11、M12,利用目标的纬度精确值M11、M12获取精确的定位结果。
本发明与现有技术相比具有如下有益效果:
本发明采用初始值计算模块计算目标在WGS-84坐标系的纬度初值M01、M02;能够根据目标的双坐标系的坐标参数(目标在WGS-84坐标系的经度L、高度H、目标到站点的距离r)与站点位置(站点在WGS-84坐标系的经度l、纬度m、高度h),计算目标在WGS-84坐标系的纬度值M01、M02;避免使用误差较大的方位测量值与距离测量值,使用精度极高的距离测量值,使定位精度更高。
本发明采用精确值计算模块,根据目标双坐标参数与纬度初值M01和站点位置,通过循环迭代计算目标纬度精确值M11;根据目标双坐标参数与纬度初值M02和站点位置,通过循环迭代计算目标纬度精确值M12。可以获取精确的定位结果,可用于目标监视系统中已知双坐标系的坐标参数的目标定位,可以解决未知目标在站点球坐标系或在WGS-84坐标系的全参数的目标定位精度不高的问题。
附图说明
为了更清楚地理解本发明,将通过具体实施方案,同时参照附图来描述本发明,其中:
图1双坐标参数空间目标定位。
图2是精确值计算模块的计算流程。
图3是坐标正向转换计算流程。
图4坐标反向转换计算流程。
具体实施方式:
参阅图1。根据本发明,首先,以雷达侦察站为站点,获取目标的站点球坐标和WGS-84坐标组成的双坐标参数:目标在WGS-84坐标系的经度L、高度H、目标在到站点的距离r,站点位置固定为站点在WGS-84坐标系的经度l、纬度m、高度h,采用初始值计算模块,根据目标双坐标参数和站点位置计算目标纬度初值M01、M02;然后,采用精确值计算模块,根据目标双坐标参数与纬度初值M01和站点位置,通过循环迭代计算目标纬度精确值M11;以同样的方式根据目标双坐标参数与纬度初值M02和站点位置,通过循环迭代计算目标纬度精确值M12,获取精确的定位结果。
双坐标参数为目标在WGS-84坐标系的经度L、高度H、目标到站点的距离r,站点位置为站点在WGS-84坐标系的经度l、纬度m、高度h,L的单位为弧度,H的单位为米,r的单位为米,l的单位为弧度,m的单位为弧度,h的单位为米。
初始值计算模块根据目标双坐标参数和站点位置计算目标纬度初值M01、M02,纬度初值M01、纬度初值M02的计算公式为:
其中,a cos是反余弦函数,cos是余弦函数。
精确值计算模块根据目标双坐标参数与纬度初值M01和站点位置,通过循环迭代计算目标纬度精确值M11的计算过程可以参阅图2。
参阅图2。精确值计算模块根据目标双坐标参数与纬度初值M01和站点位置,通过循环迭代计算目标纬度精确值M11的计算步骤为:首先输入目标纬度初值并初始化迭代参数:循环迭代次数i=0,输入目标纬度初值Mi=M01;然后进行坐标正向转换:根据目标的WGS-84坐标(L,Mi,H)和站点位置计算目标的站点球坐标方向参数;再进行坐标反向转换:根据目标的站点球坐标参数和站点位置计算目标的WGS-84坐标(Li、Mi、Hi);最后进行循环迭代是否终止判断:迭代次数i加1,即i=i+1,如果|L-Li|<0.00001且|H-Hi|<0.1或循环迭代次数i>20,则循环迭代终止,输出计算结果Mi作为纬度精确值M11,否则,转入坐标正向转换,开始继续迭代计算。
循环迭代计算目标纬度精确值M11的具体计算步骤为:
1)初始化:精确值计算模块输入目标纬度初值Mi=M01及初始化迭代参数:i=0;
2)坐标正向转换:精确值计算模块根据目标的WGS-84坐标(L、Mi、H)和站点位置(l、m、h),计算目标的站点球坐标方向参数,完成坐标正向转换;计算出的目标的站点球坐标方向参数包括:方位角的正弦值sin ai、方位角的余弦值cos ai、俯仰角的正弦值sin bi、俯仰角的余弦值cos bi;
3)坐标反向转换:精确值计算模块根据目标的站点球坐标参数(目标到站点的距离r和计算的站点球坐标方向参数:方位角的正弦值sin ai、方位角的余弦值cos ai、俯仰角的正弦值sin bi、俯仰角的余弦值cos bi)和站点位置(l、m、h),计算目标的WGS-84坐标(Li、Mi、Hi),完成坐标反向转换;
4)循环迭代是否终止判断:精确值计算模块令迭代次数加1,即i=i+1,如果|L-Li|<0.00001且|H-Hi|<0.1或循环迭代次数i>20,则循环迭代终止,输出计算结果Mi作为纬度精确值M11;否则,转入步骤2)开始继续迭代计算。
同理:精确值计算模块根据目标双坐标参数与纬度初值M02和站点位置,通过循环迭代计算目标纬度精确值M12的计算步骤为:
1)初始化:精确值计算模块输入目标纬度初值Mi=M02及初始化迭代参数:i=0;
2)坐标正向转换:精确值计算模块根据目标的WGS-84坐标(L、Mi、H)和站点位置(l、m、h),计算目标的站点球坐标方向参数,完成坐标正向转换;计算出的目标的站点球坐标方向参数包括:方位角的正弦值sinai、方位角的余弦值cosai、俯仰角的正弦值sinbi、俯仰角的余弦值cosbi;
3)坐标反向转换:精确值计算模块根据目标的站点球坐标参数(目标到站点的距离r和计算的站点球坐标方向参数:方位角的正弦值sinai、方位角的余弦值cosai、俯仰角的正弦值sinbi、俯仰角的余弦值cosbi)和站点位置(l、m、h),计算目标的WGS-84坐标(Li、Mi、Hi),完成坐标反向转换;
4)循环迭代是否终止判断:精确值计算模块令迭代次数加1,即i=i+1,如果|L-Li|<0.00001且|H-Hi|<0.1或循环迭代次数i>20,则循环迭代终止,输出计算结果Mi作为纬度精确值M12;否则,转入步骤2)开始继续迭代计算。
精确值计算模块根据目标的WGS-84坐标(L、Mi、H)和站点位置(l、m、h),计算目标的站点球坐标方向参数的计算过程参阅图3
参阅图3。在坐标正向转换计算中,精确值计算模块首先根据目标的WGS-84坐标(L、Mi、H)计算目标的地心直角坐标(Xi、Yi、Zi);然后,根据目标的地心直角坐标(Xi、Yi、Zi)和站点位置(l、m、h),计算目标的站点直角坐标(xi、yi、zi);最后,根据目标的站点直角坐标(xi、yi、zi),计算目标的站点球坐标方向参数。
精确值计算模块根据目标的WGS-84坐标(L、Mi、H)计算目标的地心直角坐标(Xi、Yi、Zi)的计算公式为:
其中,Ni为目标地心坐标迭代计算的椭球卯酉圈曲率半径。
精确值计算模块根据目标的地心直角坐标(Xi、Yi、Zi)和站点位置(l、m、h),计算目标的站点直角坐标(xi、yi、zi)的计算公式为
其中,Ns为站点地心坐标迭代计算的椭球卯酉圈曲率半径,(Xs、Ys、Zs)为站点的地心直角坐标。
精确值计算模块根据目标的站点直角坐标(xi、yi、zi),计算目标的站点球坐标方向参数的计算公式为:
其中,ri为迭代计算的目标到站点的距离,ai为目标在站点球坐标系的方位角,bi为目标在站点球坐标系的俯仰角。
精确值计算模块根据目标的站点球坐标参数和站点位置(l、m、h),计算目标的WGS-84坐标(Li、Mi、Hi)的计算过程参阅图4。
参阅图4。在坐标反向转换计算中,精确值计算模块根据目标的站点球坐标参数和站点位置(l、m、h),计算目标的WGS-84坐标(Li、Mi、Hi)的计算过程为:首先根据目标的站点球坐标参数,计算目标的站点直角坐标(xi、yi、zi);然后,根据目标的站点直角坐标(xi、yi、zi)和站点位置(l、m、h),计算目标的地心直角坐标(Xi、Yi、Zi);最后,根据目标的地心直角坐标(Xi、Yi、Zi)计算目标的WGS-84坐标(Li、Mi、Hi)。
精确值计算模块根据目标的站点球坐标参数,计算目标的站点直角坐标(xi、yi、zi)的计算公式为:
其中,目标的站点球坐标参数为:目标到站点的距离r和目标的站点球坐标方向参数:方位角的正弦值sinai、方位角的余弦值cosai、俯仰角的正弦值sinbi、俯仰角的余弦值cosbi。
精确值计算模块根据目标的站点直角坐标(xi、yi、zi)和站点位置(l、m、h),计算目标的地心直角坐标(Xi、Yi、Zi)的计算公式为:
其中,Ns为站点地心坐标迭代计算的椭球卯酉圈曲率半径,(Xs、Ys、Zs)为站点的地心直角坐标,sin为正弦函数,cos为余弦函数。
精确值计算模块根据目标的地心直角坐标(Xi、Yi、Zi)计算目标的WGS-84坐标(Li、Mi、Hi)中的经度Li的计算公式为:
其中,atan为反正切函数,为数学常量3.1415926。
精确值计算模块根据目标的地心直角坐标(Xi、Yi、Zi)计算目标的WGS-84坐标(Li、Mi、Hi)中的纬度Mi、高度Hi,需采用循环迭代方式计算。循环迭代计的初始值计算公式和循环迭代计算公式分别为:
和
循环迭代停止条件为|mj-mj-1|<0.0000001且|hj-hj-1|<0.001或迭代次数大于5。迭代停止后,所计算的纬度mj作为目标的WGS-84坐标纬度Mi,即Mi=mj,计算的高度hj作为目标的WGS-84坐标的高度Hi,即Hi=hj,
其中,n0椭球卯酉圈曲率半径循环迭代计算的初值、m0为纬度循环迭代计算的初值、h0为高度循环迭代的初值,nj椭球卯酉圈曲率半径循环迭代计算值、mj-1为上次纬度循环迭代计算值、mj为本次纬度循环迭代计算值、hj-1为上次高度循环迭代计算值、hj为本次高度循环迭代计算值。
Claims (9)
1.一种双坐标参数空间目标定位方法,其特征在于,包括如下步骤:首先,以雷达侦察站为站点,获取目标的站点球坐标和WGS-84坐标组成的双坐标参数:目标在WGS-84坐标系的经度L、高度H、目标在到站点的距离r,站点位置固定为站点在WGS-84坐标系的经度l、纬度m、高度h;根据上述目标参数和站点位置采用初始值计算模块分别计算出目标的纬度初值M01、M02;然后,根据目标双坐标参数与纬度初值M01、M02和站点位置采用精确值计算模块进行循环迭代计算,分别计算出目标的纬度精确值M11、M12,利用目标的纬度精确值M11、M12获取精确的定位结果;
初始值计算模块根据目标双坐标参数和站点位置计算目标纬度初值M01、M02的计算公式为:
其中,acos是反余弦函数,cos是余弦函数。
2.如权利要求1所述的双坐标参数空间目标定位方法,其特征在于:精确值计算模块通过循环迭代计算目标纬度精确值M11、M12的计算步骤为:
1)初始化:精确值计算模块输入目标纬度初值Mi=M01、Mi=M02及初始化迭代参数:i=0;
2)坐标正向转换:精确值计算模块根据目标的WGS-84坐标(L、Mi、H)和站点位置(l、m、h),计算目标的站点球坐标方向参数,完成坐标正向转换,并计算出目标的站点球坐标方向参数的方位角的正弦值sinai、方位角的余弦值cosai、俯仰角的正弦值sinbi、俯仰角的余弦值cosbi;
3)坐标反向转换:精确值计算模块根据目标到站点的距离r和计算的上述方位角的正弦值sinai、方位角的余弦值cosai、俯仰角的正弦值sinbi、俯仰角的余弦值cosbi)和站点位置(l、m、h),计算目标的WGS-84坐标(Li、Mi、Hi),完成坐标反向转换;
4)循环迭代是否终止判断:精确值计算模块令迭代次数加1,即i=i+1,如果|L-Li|<0.00001且|H-Hi|<0.1或循环迭代次数i>20,则循环迭代终止,输出初值Mi=M01的计算结果作为纬度精确值M11,初值Mi=M02的计算结作为纬度精确值M12;否则,转入步骤2)开始继续迭代计算。
3.如权利要求2所述的双坐标参数空间目标定位方法,其特征在于:在坐标正向转换计算中,精确值计算模块首先根据目标的WGS-84坐标(L、Mi、H)计算目标的地心直角坐标(Xi、Yi、Zi);然后,根据目标的地心直角坐标(Xi、Yi、Zi)和站点位置(l、m、h),计算目标的站点直角坐标(xi、yi、zi);最后,根据目标的站点直角坐标(xi、yi、zi),计算目标的站点球坐标方向参数。
4.如权利要求2所述的双坐标参数空间目标定位方法,其特征在于:在坐标反向转换计算中,精确值计算模块首先根据目标的站点球坐标参数,计算目标的站点直角坐标(xi、yi、zi);然后,根据目标的站点直角坐标(xi、yi、zi)和站点位置(l、m、h),计算目标的地心直角坐标(Xi、Yi、Zi);最后,根据目标的地心直角坐标(Xi、Yi、Zi)计算目标的WGS-84坐标(Li、Mi、Hi)。
5.如权利要求4所述的双坐标参数空间目标定位方法,其特征在于:精确值计算模块根据目标的WGS-84坐标(L、Mi、H)计算目标的地心直角坐标(Xi、Yi、Zi)的计算公式为:
其中,Ni为目标地心坐标迭代计算的椭球卯酉圈曲率半径。
6.如权利要求5所述的双坐标参数空间目标定位方法,其特征在于:精确值计算模块根据目标的地心直角坐标(Xi、Yi、Zi)和站点位置(l、m、h),计算目标的站点直角坐标(xi、yi、zi)的计算公式为
其中,Ns为站点地心坐标迭代计算的椭球卯酉圈曲率半径,(Xs、Ys、Zs)为站点的地心直角坐标。
7.如权利要求6所述的双坐标参数空间目标定位方法,其特征在于:精确值计算模块根据目标的站点直角坐标(xi、yi、zi),计算目标的站点球坐标方向参数的计算公式为:
其中,ri为迭代计算的目标到站点的距离,ai为目标在站点球坐标系的方位角,bi为目标在站点球坐标系的俯仰角。
8.如权利要求7所述的双坐标参数空间目标定位方法,其特征在于:精确值计算模块根据目标的站点直角坐标(xi、yi、zi)和站点位置(l、m、h),计算目标的地心直角坐标(Xi、Yi、Zi)的计算公式为:
根据目标的地心直角坐标(Xi、Yi、Zi)计算目标的WGS-84坐标(Li、Mi、Hi)中的经度Li的计算公式为:
其中,Ns为站点地心坐标迭代计算的椭球卯酉圈曲率半径,(Xs、Ys、Zs)为站点的地心直角坐标,sin为正弦函数,cos为余弦函数,atan为反正切函数,π为数学常量3.1415926。
9.如权利要求1所述的双坐标参数空间目标定位方法,其特征在于:精确值计算模块循环迭代计的初始值计算公式和循环迭代计算公式分别为:
和
循环迭代停止条件为|mj-mj-1|<0.0000001且|hj-hj-1|<0.001或迭代次数大于5;迭代停止后,所计算的纬度mj作为目标的WGS-84坐标纬度Mi,即Mi=mj,计算的高度hj作为目标的WGS-84坐标的高度Hii,即Hi=hj,
其中,n0椭球卯酉圈曲率半径循环迭代计算的初值、m0为纬度循环迭代计算的初值、h0为高度循环迭代的初值,nj椭球卯酉圈曲率半径循环迭代计算值、mj-1为上次纬度循环迭代计算值、mj为本次纬度循环迭代计算值、hj-1为上次高度循环迭代计算值、hj为本次高度循环迭代计算值。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202011042966.3A CN112254719B (zh) | 2020-09-28 | 2020-09-28 | 双坐标参数空间目标定位方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202011042966.3A CN112254719B (zh) | 2020-09-28 | 2020-09-28 | 双坐标参数空间目标定位方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN112254719A CN112254719A (zh) | 2021-01-22 |
CN112254719B true CN112254719B (zh) | 2023-09-05 |
Family
ID=74234391
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202011042966.3A Active CN112254719B (zh) | 2020-09-28 | 2020-09-28 | 双坐标参数空间目标定位方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN112254719B (zh) |
Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2008097738A1 (en) * | 2007-02-06 | 2008-08-14 | The Boeing Company | Automatic stereo measurement of a point of interest in a scene |
CN105136128A (zh) * | 2015-08-27 | 2015-12-09 | 中国航天空气动力技术研究院 | 基于两点定位的机体结构测量方法 |
CN105446690A (zh) * | 2015-12-22 | 2016-03-30 | 中国电子科技集团公司第十一研究所 | 具有目标定位功能的信息融合及多信息显示方法 |
CN107328420A (zh) * | 2017-08-18 | 2017-11-07 | 上海木爷机器人技术有限公司 | 定位方法和装置 |
CN107490364A (zh) * | 2017-09-01 | 2017-12-19 | 中国科学院长春光学精密机械与物理研究所 | 一种大角度倾斜成像航空相机对地目标定位方法 |
CN108594275A (zh) * | 2018-04-26 | 2018-09-28 | 桂林电子科技大学 | 一种北斗+gps双模单点定位方法 |
CN110440801A (zh) * | 2019-07-08 | 2019-11-12 | 浙江吉利控股集团有限公司 | 一种定位感知信息获取方法、装置及系统 |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US8004930B2 (en) * | 2008-03-17 | 2011-08-23 | Westerngeco, L.L.C. | Methods and systems for determining coordinates of an underwater seismic component in a reference frame |
-
2020
- 2020-09-28 CN CN202011042966.3A patent/CN112254719B/zh active Active
Patent Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2008097738A1 (en) * | 2007-02-06 | 2008-08-14 | The Boeing Company | Automatic stereo measurement of a point of interest in a scene |
CN105136128A (zh) * | 2015-08-27 | 2015-12-09 | 中国航天空气动力技术研究院 | 基于两点定位的机体结构测量方法 |
CN105446690A (zh) * | 2015-12-22 | 2016-03-30 | 中国电子科技集团公司第十一研究所 | 具有目标定位功能的信息融合及多信息显示方法 |
CN107328420A (zh) * | 2017-08-18 | 2017-11-07 | 上海木爷机器人技术有限公司 | 定位方法和装置 |
CN107490364A (zh) * | 2017-09-01 | 2017-12-19 | 中国科学院长春光学精密机械与物理研究所 | 一种大角度倾斜成像航空相机对地目标定位方法 |
CN108594275A (zh) * | 2018-04-26 | 2018-09-28 | 桂林电子科技大学 | 一种北斗+gps双模单点定位方法 |
CN110440801A (zh) * | 2019-07-08 | 2019-11-12 | 浙江吉利控股集团有限公司 | 一种定位感知信息获取方法、装置及系统 |
Non-Patent Citations (1)
Title |
---|
王前东."双坐标系多参量目标的融合定位".《电讯技术》.2011,第51卷(第6期),第66-70页. * |
Also Published As
Publication number | Publication date |
---|---|
CN112254719A (zh) | 2021-01-22 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
US11860260B2 (en) | Systems and methods for high-integrity satellite positioning | |
Zhu et al. | GNSS position integrity in urban environments: A review of literature | |
US11867823B2 (en) | System and method for determining GNSS positioning corrections | |
Dai et al. | Real-time attitude determination for microsatellite by LAMBDA method combined with Kalman filtering | |
CA2681918A1 (en) | Distance dependant error mitigation in real-time kinematic (rtk) positioning | |
Bahrami et al. | Instantaneous Doppler-aided RTK positioning with single frequency receivers | |
US20230184956A1 (en) | System and method for correcting satellite observations | |
Jokinen | Enhanced ambiguity resolution and integrity monitoring methods for precise point positioning | |
Min et al. | DNN-based approach to mitigate multipath errors of differential GNSS reference stations | |
Martin | GNSS precise point positioning: The enhancement with GLONASS | |
Obst et al. | Probabilistic non-line-of-sight detection in reliable urban GNSS vehicle localization based on an empirical sensor model | |
Ferrão | Positioning with combined gps and glonass observations | |
Ohashi et al. | Multi-layer ionospheric VTEC estimation for GNSS positioning | |
CN112254719B (zh) | 双坐标参数空间目标定位方法 | |
Weisenburger | Effect of constraints and multiple receivers for on-the-fly ambiguity resolution | |
Liu et al. | Actualization Analysis of LEO Opportunistic Doppler Aided GNSS Precise Point Positioning using Moving Horizon Estimation | |
Zocca et al. | A GNSS Multipath and NLoS Mitigation Method for Urban Scenarios Based on Particle Filtering | |
Khojasteh et al. | Introduction to global navigation satellite systems and its errors | |
Noureldin et al. | Global Positioning System | |
Krasuski | Application the GPS code observations in BSSD method for recovery the position of the aircraft | |
Fujita et al. | A Single Frequency PPP Algorithm Based on GR Models | |
Ohashi et al. | Prediction of regional ionospheric models with AR model for GNSS positioning in Japan | |
Iwamoto et al. | Receiver Clock-Based Integrity Monitoring for GNSS Positioning | |
Chang et al. | Vector angle grouping-based solution separation for multipath/NLOS detection and exclusion with the enhancement of doppler test | |
Shah | Using Random Sampling Consensus (RANSAC) to Detect Errors in Global Navigation Satellite Systems (GNSS) Signals and Data |
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 |