[发明专利]一种根据GGA定位信息快速计算最近格网点的方法

权利要求书

1.一种根据GGA定位信息快速计算最近格网点的方法，其特征在于包括：

S100、计算基准站的精确坐标、概略坐标和基准站时钟相对于卫星时钟的钟差改正数，计算基准站的精确坐标与概略坐标之间的差值作为基准站的差分数据；计算流动站的概略坐标；

S200、基准站通过网络将基准站的精确坐标和差分数据发送到控制中心，流动站通过网络将流动站的概略坐标发送给控制中心；

S300、控制中心根据流动站的概略坐标和各基准站的精确坐标，为流动站选取参考基准站组，并由参考基准站组的各个基准站的差分数据计算虚拟基准站的虚拟差分数据；

S400、控制中心将虚拟基准站的虚拟差分数据通过网络发送给对应的流动站，流动站根据得到的虚拟基准站的虚拟差分数据计算自身的精确坐标；

S500、基准站和流动站分别根据各自的精确坐标生成各自的卫星定位GGA数据并分别通过网络将各自的卫星定位GGA数据发送给控制中心；

S600、控制中心对基准站和流动站的各自的卫星定位GGA数据进行加偏获得基准站和流动站各自的GCJ02坐标系的经度坐标和纬度坐标；

S700、根据流动站的GCJ02坐标系的经度坐标和纬度坐标计算精确坐标定位位置的多个最近格网点坐标。

2.根据权利要求1所述的方法，其特征在于：

步骤S100中，在基准站上安置卫星接收机进行长时间观测得到基准站的精确坐标。

3.根据权利要求2所述的方法，其特征在于：

步骤S100中，求解以下公式获得基准站的概略坐标和基准站时钟相对于卫星时钟的钟差改正数：

其中，c为电磁波在真空中的速度；Δt_i表示第i个卫星至基准站的测距信号时延，i＝1,2,…,n；τ表示基准站时钟相对于卫星时钟的钟差改正数；(x_i,y_i,z_i)表示第i个卫星的坐标；(x₀,y₀,z₀)为基准站的概略坐标；n为卫星个数，n大于等于4。

4.根据权利要求3所述的方法，其特征在于：

步骤S100中，计算基准站的精确坐标与概略坐标之间的差值作为差分数据。

5.根据权利要求4所述的方法，其特征在于：

步骤S300包括以下步骤：

S310、控制中心为流动站选择参考基准站组；

S320、根据以下公式计算虚拟基准站的虚拟差分数据：

式中，ΔX为虚拟基准站的虚拟差分数据；ΔX_j表示参考基准站组内第j个基准站的差分数据，j＝1,2,…,m，m为参考基准站组内基准站的总数；W_j表示为参考基准站组内第j个基准站的权重；d_j表示参考基准站组内第j个基准站至流动站的距离，(x_j,y_j,z_j)为第j个基准站的精确坐标；p为调整参数。

6.根据权利要求5所述的方法，其特征在于：

步骤S400包括以下步骤：控制中心将虚拟基准站的虚拟差分数据通过网络发送给对应的流动站，流动站将接收到的虚拟基准站的虚拟差分数据加入到流动站的概略坐标即得到流动站的精确坐标；

步骤S600中以下述公式对卫星定位GGA数据进行加偏：

其中，B和L分别为卫星定位GGA数据中WGS-84椭球的大地纬度坐标和大地经度坐标；B_GCJ02和L_GCJ02分别为GCJ02坐标系的经度坐标和纬度坐标；ΔB为卫星定位GGA数据中的大地纬度坐标至GCJ02坐标系的纬度坐标的偏移量，ΔL为卫星定位GGA数据中的大地经度坐标至GCJ02坐标系的经度坐标的偏移量；a和e分别为卫星定位GGA数据中WGS-84椭球的长半轴和第一偏心率；B′＝B-35，L′＝L-105。