基于NTRIP测定任意坐标系成果的RTD方法技术

一种基于NTRIP测定任意坐标系成果的RTD方法，对源节点进行扩展，用于更详尽地传达流动站需求，不需要修改NTRIP协议，只需在NtripCaster软件中加入判断及算法，或在NtripCaster和流动站之间增加NtripProxy中继模块；需要实时得到参心坐标系下坐标，只需在登录服务器时选择相应的源节点并在手薄上填入非涉密七参数；需要实时得到正常高，只需在登录服务器时选择相应的源节点；需要实时得到ECEF坐标，只需在登录服务器时选择相应的源节点。作业方式基本不变，在符合保密政策的前提下，使RTD实时测量放样参心坐标系坐标和正常高得以实现，显著提高了测绘生产效率。

【技术实现步骤摘要】

【技术保护点】
一种基于NTRIP测定任意坐标系成果的RTD方法，其特征在于：在NtripCaster服务器和NtripClient之间插入一个NtripProxy中继模块，NtripProxy中继模块所在服务器记为NtripProxy服务器，NtripClient表示流动站；测定过程包括以下步骤，步骤1，定义源节点和辅助七参数，计算辅助坐标系到目标坐标系的转换参数作为非涉密七参数；所述辅助七参数，包括3个平移量ΔX1、ΔY1、ΔZ1，3个旋转参数εX1、εY1、εZ1，1个尺度缩放因子m1，得到下述布尔莎模型，XYZFZ=(1+m1)1ϵZ1-ϵY1-ϵZ11ϵX1ϵY1-ϵX11XYZECEF+ΔX1ΔY1ΔZ1,其中，XYZFZ表示辅助坐标系中的坐标，XYZECEF表示ECEF坐标系中的坐标，ECEF坐标系表示地心地固坐标系；所述计算辅助坐标系到目标坐标系的转换参数作为非涉密七参数，包括设ECEF坐标系到目标坐标系的转化参数为ΔX0、ΔY0、ΔZ0、εX0、εY0、εZ0、m0，定义辅助坐标系到目标坐标系的转换参数为ΔX、ΔY、ΔZ、εX、εY、εZ、m，通过下式计算，ΔXΔYΔZ=ΔX0ΔY0ΔZ0+(1+m)1ϵZ-ϵY-ϵZ1ϵXϵY-ϵX1ΔX1ΔY1ΔZ1m=m0-m11+m1ϵX=ϵX0-ϵX1ϵY=ϵY0-ϵY1ϵZ=ϵZ0-ϵZ1步骤2，生成网络RTK差分电文，包含以下子步骤，步骤2.1，NtripClient登录NtripProxy服务器，连接时选取源节点，并将初步定位得到的 ECEF坐标系下的流动站概略位置XrovECEF=XrovECEFYrovECEFZrovECEF发往NtripProxy服务器，NtripProxy服务器将XrovECEF转发给NtripCaster服务器，NtripCaster服务器计算给出参考站坐标XrefECEF＝(XrefECEF,YrefECEF,ZrefECEF)，并基于XrefECEF生成伪距改正数Λ及变率如下，Λ＝[dP1dP2…dPk]TΛ·=(Λ(t)-Λ(t-1))/dt,其中，k代表卫星颗数，为各颗卫星的初始伪距改正数；代表第i颗星到参考站的几何距离，Pi代表某时刻对第i颗星的伪距观测值，dt代表GNSS信号采样间隔；NtripCaster服务器将伪距改正数Λ及变率编制为电文并发送给NtripProxy服务器；步骤2.2，NtripProxy服务器执行以下操作，步骤2.2.1，将电文解码，还原出伪距改正数Λ，步骤2.2.2，将NtripClient发送的概略位置XrovECEF按下述公式转换为经纬度大地高，LrovECEF=arctan(YrovECEFXrovECEF)BrovECEF=arctan(ZrovECEF+e′2bsin3θ(XrovECEF2+YrovECEF2)-e2acos3θ)HrovECEF=XrovECEF2+YrovECEF2cosBrovECEF-N其中N=a1-e2sin2BrovECEF,θ=arctan(aZrovECEFbXrovECEF2+YrovECEF2),e2=a2-b2a2,a为ECEF坐标系的椭球长半轴长，b为椭球短半轴长；步骤2.2.3，利用LrovECEF、BrovECEF，结合似大地水准面精化的成果，内插出NtripClient处的高程异常值ξrovECEF，令NtripClient的伪大地高Hr′ovECEF＝HrovECEF?ξrovECEF，将NtripClient 的伪大地坐标ΦrovECEF′=LrovECEFBrovECEFHrovECEF′重新转换为空间直角坐标如下，XrovECEF′=XrovECEF′YrovECEF′ZrovECEF′=(N+HrovECEF′)cosBrovECEFcosLrovEC...

【技术特征摘要】

【专利技术属性】
技术研发人员：章迪，郭际明，巢佰崇，罗年学，
申请(专利权)人：武汉大学，
类型：发明
国别省市：