CN116106948A - 网络rtk抗电离层干扰定位方法及相关设备 - Google Patents

网络rtk抗电离层干扰定位方法及相关设备 Download PDF

Info

Publication number
CN116106948A
CN116106948A CN202111322856.7A CN202111322856A CN116106948A CN 116106948 A CN116106948 A CN 116106948A CN 202111322856 A CN202111322856 A CN 202111322856A CN 116106948 A CN116106948 A CN 116106948A
Authority
CN
China
Prior art keywords
terminal
ionosphere
quality factor
rtk
filtering mode
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.)
Pending
Application number
CN202111322856.7A
Other languages
English (en)
Inventor
陈华炎
何锡扬
战兴群
冯绍军
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Qianxun Si Network Zhejiang Co ltd
Original Assignee
Qianxun Si Network Zhejiang Co ltd
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Qianxun Si Network Zhejiang Co ltd filed Critical Qianxun Si Network Zhejiang Co ltd
Priority to CN202111322856.7A priority Critical patent/CN116106948A/zh
Priority to PCT/CN2022/116380 priority patent/WO2023082785A1/zh
Publication of CN116106948A publication Critical patent/CN116106948A/zh
Pending legal-status Critical Current

Links

Images

Classifications

    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
    • G01S19/00—Satellite radio beacon positioning systems; Determining position, velocity or attitude using signals transmitted by such systems
    • G01S19/01—Satellite radio beacon positioning systems transmitting time-stamped messages, e.g. GPS [Global Positioning System], GLONASS [Global Orbiting Navigation Satellite System] or GALILEO
    • G01S19/03—Cooperating elements; Interaction or communication between different cooperating elements or between cooperating elements and receivers
    • G01S19/07—Cooperating elements; Interaction or communication between different cooperating elements or between cooperating elements and receivers providing data for correcting measured positioning data, e.g. DGPS [differential GPS] or ionosphere corrections
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
    • G01S19/00—Satellite radio beacon positioning systems; Determining position, velocity or attitude using signals transmitted by such systems
    • G01S19/01—Satellite radio beacon positioning systems transmitting time-stamped messages, e.g. GPS [Global Positioning System], GLONASS [Global Orbiting Navigation Satellite System] or GALILEO
    • G01S19/03—Cooperating elements; Interaction or communication between different cooperating elements or between cooperating elements and receivers
    • G01S19/07—Cooperating elements; Interaction or communication between different cooperating elements or between cooperating elements and receivers providing data for correcting measured positioning data, e.g. DGPS [differential GPS] or ionosphere corrections
    • G01S19/072—Ionosphere corrections
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
    • G01S19/00—Satellite radio beacon positioning systems; Determining position, velocity or attitude using signals transmitted by such systems
    • G01S19/38—Determining a navigation solution using signals transmitted by a satellite radio beacon positioning system
    • G01S19/39—Determining a navigation solution using signals transmitted by a satellite radio beacon positioning system the satellite radio beacon positioning system transmitting time-stamped messages, e.g. GPS [Global Positioning System], GLONASS [Global Orbiting Navigation Satellite System] or GALILEO
    • G01S19/40—Correcting position, velocity or attitude
    • G01S19/41—Differential correction, e.g. DGPS [differential GPS]
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
    • G01S19/00—Satellite radio beacon positioning systems; Determining position, velocity or attitude using signals transmitted by such systems
    • G01S19/38—Determining a navigation solution using signals transmitted by a satellite radio beacon positioning system
    • G01S19/39—Determining a navigation solution using signals transmitted by a satellite radio beacon positioning system the satellite radio beacon positioning system transmitting time-stamped messages, e.g. GPS [Global Positioning System], GLONASS [Global Orbiting Navigation Satellite System] or GALILEO
    • G01S19/42—Determining position

Abstract

本申请公开了一种网络RTK抗电离层干扰定位方法及相关设备。方法包括：接收服务器播发的差分改正数和质量因子，所述差分改正数和所述质量因子为服务器根据各个基站的观测数据和位置信息以及所述终端上传的概略位置信息生成；根据接收到的质量因子确定RTK滤波模式，所述RTK滤波模式包括无电离层组合滤波模式或非组合滤波模式；采用确定的所述RTK滤波模式对终端的观测数据和终端接收到的差分改正数进行滤波解算，得到所述终端的高精度位置信息。根据本申请实施例，能够通过服务器播发的质量因子进行滤波模式的切换，避免电离层活跃期时电离层误差的干扰，改善电离层活跃期的定位效果。

Description

网络RTK抗电离层干扰定位方法及相关设备
技术领域
本申请属于卫星定位技术领域，尤其涉及一种网络RTK抗电离层干扰定位方法、装置、系统、设备及存储介质。
背景技术
目前，在网络RTK(Real-time kinematic，实时动态载波相位差分)定位技术中，服务器可以对电离层误差和对流层误差进行建模，并根据终端的概略位置向终端播发相应的差分改正数。终端可以通过差分改正数消除电离层误差和对流层误差的影响，实现高精度定位。
然而，在电离层活跃期，服务器在电离层误差建模时存在残余误差，从而导致终端在使用服务器所播发的差分改正数后，仍存在较大的电离层误差，使得终端无法固定整周模糊度或者整周模糊度固定错误，无法进行高精度定位。
发明内容
本申请实施例提供了一种网络RTK抗电离层干扰定位方法、装置、系统、设备及存储介质，能够解决网络RTK定位在电离层活跃期时无法实现高精度定位的技术问题。
第一方面，本申请实施例提供一种网络RTK抗电离层干扰定位方法，应用于终端，方法包括：
接收服务器播发的差分改正数和质量因子，差分改正数和质量因子为服务器根据各个基站的观测数据和位置信息以及终端上传的概略位置信息生成，质量因子用于表征生成差分改正数所用的大气误差模型精度；
根据接收到的质量因子确定RTK滤波模式，RTK滤波模式包括无电离层组合滤波模式或非组合滤波模式；
采用确定的RTK滤波模式对终端的观测数据和终端接收到的差分改正数进行滤波解算，得到终端的高精度位置信息。
通过服务器播发表征大气误差模型精度的质量因子，能够使得终端根据该质量因子从无电离层组合滤波模式或非组合滤波模式中选择RTK滤波模式，在电离层活跃期时能够选择消除电离层误差的无电离层组合滤波模式进行滤波解算，以消除电离层误差的影响，实现电离层活跃期的高精度定位。
在一些实施例中，根据接收到的质量因子确定RTK滤波模式，包括：
根据接收到的质量因子判断当前历元的电离层是否活跃；
在当前历元的电离层活跃时，获取第一时间段内电离层活跃的历元数量；
在电离层活跃的历元数量达到预设第一历元阈值时，将RTK滤波模式确定为无电离层组合滤波模式。
终端可以根据接收到的质量因子来判断当前历元的电离层是否活跃，并确定采用的RTK滤波模式为无电离层组合模式，以消除电离层活跃时电离层误差的影响，从而得到终端的高精度位置信息。
在一些实施例中，将RTK滤波模式确定为无电离层组合滤波模式之后，还包括：
在当前历元的电离层不活跃时，获取第二时间段内电离层不活跃的历元数量；
在电离层不活跃的历元数量达到预设第二历元阈值时，RTK滤波模式切换为非组合滤波模式。
终端在根据质量因子确定当前历元的电离层不活跃时，采用的RTK滤波模式为非组合滤波模式能够消除各种误差对整周模糊度的影响，从而快速固定整周模糊度，计算出终端的高精度位置信息。
在一些实施例中，质量因子包括对应多个卫星的电离层质量因子和对应多个卫星的对流层质量因子；根据接收到的质量因子判断当前历元的电离层是否活跃，包括：
将接收到的多个卫星的电离层质量因子和多个卫星的对流层质量因子分别进行从小到大排序；
判断第一分位的电离层质量因子是否大于第二分位的对流层质量因子，以及第一分位的电离层质量因子是否大于预设最小活跃阈值，其中第一分位小于或等于第二分位；
根据接收到的质量因子判断当前历元的电离层是否活跃之后，还包括：
在第一分位的电离层质量因子大于第二分位的对流层质量因子，且第一分位的电离层质量因子大于预设最小活跃阈值时，确定当前历元的电离层活跃；或者
将接收到的多个卫星的电离层质量因子和多个卫星的对流层质量因子分别进行从大到小排序；
判断第三分位的电离层质量因子是否大于第四分位的对流层质量因子，以及第三分位的电离层质量因子是否大于预设最小活跃阈值，其中第三分位大于或等于第四分位；
根据接收到的质量因子判断当前历元的电离层是否活跃之后，还包括：
在第三分位的电离层质量因子大于第四分位的对流层质量因子，且第三分位的电离层质量因子大于预设最小活跃阈值时，确定当前历元的电离层活跃。
通过将电离层质量因子与对流层质量因子进行对比，能够排除基站观测、基站网络拓扑结构等其他因素造成的公共误差的影响，从而准确地判断当前电离层的活跃状态。
第二方面，本申请实施例提供一种网络RTK抗电离层干扰定位方法，应用于服务器，方法包括：
根据基站网络中各个基站的观测数据和位置信息，生成多条基线的电离层误差和对流层误差，基线为两个基站相连形成；
基于多条基线的电离层误差和对流层误差以及各个基站的位置信息进行建模，得到大气误差模型；
根据大气误差模型的模型参数计算基站网络中多个格网点位置的差分改正数以及多个格网点位置的质量因子；
根据终端发送的概略位置信息从多个格网点位置中确定距离终端概略位置最近的目标格网点；
向终端播发目标格网点的差分改正数和质量因子，以使终端根据目标格网点的质量因子确定RTK滤波模式，并根据RTK滤波模式采用差分改正数进行滤波解算以得到终端的高精度位置信息，RTK滤波模式包括无电离层组合滤波模式或非组合滤波模式。
服务器能够根据终端发送的概略位置向终端播发目标格网点的差分改正数和质量因子，以使终端根据质量因子确定RTK滤波模式后，根据差分改正数和终端的观测数据进行滤波解算，得到终端的高精度位置信息。
在一些实施例中，根据基站网络中各个基站的观测数据和位置信息，生成多条基线的电离层误差和对流层误差，包括：
根据基站网络中各个基站的观测数据和位置信息，解算多条基线的双差模糊度；
根据多条基线的双差模糊度计算多条基线的电离层误差值和对流层误差值。
服务器可以根据各个基站的观测数据得到各条基线的双差模糊度，并反算出每条基线的差分改正数，从而根据基线的差分改正数进行大气误差模型的建模。
在一些实施例中，根据大气误差模型的模型参数计算基站网络中多个格网点位置的差分改正数以及多个格网点位置的质量因子，包括：
根据大气误差模型的模型参数、各条基线对应基站的位置信息以及各条基线的电离层误差值和对流层误差值计算得到大气误差模型中各条基线的观测值残差；
根据各条基线的观测值残差计算模型参数的精度因子矩阵；
根据精度因子矩阵和基站网络中多个格网点的位置信息，分别计算得到每个格网点位置的质量因子。
服务器可以确定各个格网点的差分改正数，并根据各条基线实际观测得到的差分改正数和大气误差模型所计算得出的差分改正数计算得出精度因子矩阵，从而根据精度因子矩阵确定每个格网点的质量因子。
第三方面，本申请实施例提供一种网络RTK抗电离层干扰的定位装置，装置包括：
接收模块，用于接收服务器播发的差分改正数和质量因子，差分改正数和质量因子为服务器根据各个基站的观测数据和位置信息以及终端上传的概略位置信息生成，质量因子用于表征生成差分改正数所用的大气误差模型精度；
判定模块，用于根据接收到的质量因子确定RTK滤波模式，RTK滤波模式包括无电离层组合滤波模式或非组合滤波模式；
定位模块，用于采用确定的RTK滤波模式对终端的观测数据和终端接收到的差分改正数进行滤波解算，得到终端的高精度位置信息；
或者，装置包括：
误差模块，用于根据基站网络中各个基站的观测数据和位置信息，生成多条基线的电离层误差和对流层误差，基线为两个基站相连形成；
建模模块，用于基于多条基线的电离层误差和对流层误差以及各个基站的位置信息进行建模，得到大气误差模型；
计算模块，用于根据大气误差模型的模型参数计算基站网络中多个格网点位置的差分改正数以及多个格网点位置的质量因子；
匹配模块，用于根据终端发送的概略位置信息从多个格网点位置中确定距离终端概略位置最近的目标格网点；
播发模块，用于向终端播发目标格网点的差分改正数和质量因子，以使终端根据目标格网点的质量因子确定RTK滤波模式，并根据RTK滤波模式采用差分改正数进行滤波解算以得到终端的高精度位置信息，RTK滤波模式包括无电离层组合滤波模式或非组合滤波模式。
第四方面，本申请实施例提供了一种网络RTK抗电离层干扰定位系统，包括服务器和终端，系统包括：
服务器根据基站网络中各个基站的观测数据和位置信息，生成多条基线的电离层误差和对流层误差，基线为两个基站相连形成；
服务器基于多条基线的电离层误差和对流层误差以及各个基站的位置信息进行建模，得到大气误差模型；
服务器根据大气误差模型的模型参数计算基站网络中多个格网点位置的差分改正数以及多个格网点位置的质量因子；
服务器根据终端发送的概略位置信息从多个格网点位置中确定距离终端概略位置最近的目标格网点；
服务器向终端播发目标格网点的差分改正数和质量因子；
终端接收服务器播发的差分改正数和质量因子，差分改正数和质量因子为服务器根据各个基站的观测数据和位置信息以及终端上传的概略位置信息生成，质量因子用于表征生成差分改正数所用的大气误差模型精度；
终端根据接收到的质量因子确定RTK滤波模式，RTK滤波模式包括无电离层组合滤波模式或非组合滤波模式；
终端采用确定的RTK滤波模式对终端的观测数据和终端接收到的差分改正数进行滤波解算，得到终端的高精度位置信息。
第五方面，本申请实施例提供了一种网络RTK抗电离层干扰定位设备，网络RTK抗电离层干扰定位设备包括：处理器以及存储有计算机程序指令的存储器；处理器执行计算机程序指令时实现如上的网络RTK抗电离层干扰定位方法。
第六方面，本申请实施例提供了一种计算机存储介质，计算机存储介质上存储有计算机程序指令，计算机程序指令被处理器执行时实现如上的网络RTK抗电离层干扰定位方法。
与现有技术相比，本申请实施例提供的网络RTK抗电离层干扰定位方法，在服务器对电离层误差和对流层误差进行建模时，生成表征大气误差模型精度的质量因子，并将质量因子与差分改正数一同播发给终端。终端可以根据该质量因子判断当前的差分改正数是否受到电离层活跃的影响，并根据判断结果从无电离层组合滤波模式或非组合滤波模式中选择合适的RTK滤波模式，对终端的观测数据和终端接收的差分改正数进行滤波解算，以在电离层活跃时消除电离层误差的影响，实现高精度定位。通过服务器所播发的质量因子，可以根据电离层的活跃状态进行滤波模式的切换，避免电离层活跃期时电离层误差的干扰，改善电离层活跃期的定位效果。
附图说明
为了更清楚地说明本申请实施例的技术方案，下面将对本申请实施例中所需要使用的附图作简单地介绍，显而易见地，下面所描述的附图仅仅是本申请的一些实施例，对于本领域普通技术人员来讲，在不付出创造性劳动的前提下，还可以根据这些附图获得其他的附图。
图1是本申请一实施例提供的网络RTK抗电离层干扰定位方法的流程示意图；
图2是本申请另一实施例提供的网络RTK抗电离层干扰定位方法的流程示意图；
图3是本申请又一实施例提供的网络RTK抗电离层干扰定位方法的流程示意图；
图4是图2实施例中S210的细化流程示意图；
图5是本申请一实施例提供的网络RTK抗电离层干扰定位方法的流程示意图；
图6是本申请另一实施例提供的网络RTK抗电离层干扰定位方法的流程示意图；
图7是本申请又一实施例提供的网络RTK抗电离层干扰定位方法的流程示意图；
图8是本申请一实施例提供的网络RTK抗电离层干扰定位装置的结构示意图；
图9是本申请另一实施例提供的网络RTK抗电离层干扰定位装置的结构示意图；
图10是本申请一实施例提供的网络RTK抗电离层干扰定位装置设备的硬件结构示意图。
具体实施方式
下面将详细描述本申请的各个方面的特征和示例性实施例，为了使本申请的目的、技术方案及优点更加清楚明白，以下结合附图及具体实施例，对本申请进行进一步详细描述。应理解，此处所描述的具体实施例仅意在解释本申请，而不是限定本申请。对于本领域技术人员来说，本申请可以在不需要这些具体细节中的一些细节的情况下实施。下面对实施例的描述仅仅是为了通过示出本申请的示例来提供对本申请的更好的理解。
需要说明的是，在本文中，诸如第一和第二等之类的关系术语仅仅用来将一个实体或者操作与另一个实体或操作区分开来，而不一定要求或者暗示这些实体或操作之间存在任何这种实际的关系或者顺序。而且，术语“包括”、“包含”或者其任何其他变体意在涵盖非排他性的包含，从而使得包括一系列要素的过程、方法、物品或者设备不仅包括那些要素，而且还包括没有明确列出的其他要素，或者是还包括为这种过程、方法、物品或者设备所固有的要素。在没有更多限制的情况下，由语句“包括……”限定的要素，并不排除在包括所述要素的过程、方法、物品或者设备中还存在另外的相同要素。
需要说明的是，在不冲突的情况下，本申请中的实施例及实施例中的特征可以相互组合。下面将结合附图对实施例进行详细描述。
在现代GNSS(Global navigation satellite system，全球卫星导航系统)高精度实时定位中，网络RTK定位技术已经成为目前的主流应用。在网络RTK定位技术中，各个基站可以分别获取电离层和对流层的大气延迟信息，并计算生成该基站所对应的电离层误差和对流层误差。根据终端所发送的概略位置，能够根据各个基站所组成的基站网络进行大气误差模型的建模，并根据建立的大气误差模型计算得到概略位置的差分改正数后发送至终端，以使终端对观测数据和差分改正数进行滤波解算，得到终端的高精度位置信息。
然而，在电离层活跃期，对电离层的误差进行建模时，存在较大的观测值误差。当终端接收到差分改正数后，由于仍存在较大的电离层误差，使得终端无法固定整周模糊度或者错误固定整周模糊度，从而导致无法得到高精度的定位结果。
为了解决上述技术问题，本申请实施例提供了一种网络RTK抗电离层干扰定位方法、装置、系统、设备及存储介质。下面首先对本申请实施例所提供的网络RTK抗电离层干扰定位方法进行介绍。
图1示出了本申请一个实施例提供的网络RTK抗电离层干扰定位方法的流程示意图。网络RTK抗电离层干扰定位方法应用于终端，终端包括：
S110，接收服务器播发的差分改正数和质量因子，所述差分改正数和所述质量因子为服务器根据各个基站的观测数据和位置信息以及所述终端上传的概略位置信息生成，所述质量因子用于表征生成所述差分改正数所用的大气误差模型精度；
S120，根据接收到的质量因子确定RTK滤波模式，所述RTK滤波模式包括无电离层组合滤波模式或非组合滤波模式；
S130，采用确定的所述RTK滤波模式对终端的观测数据和终端接收到的差分改正数进行滤波解算，得到所述终端的高精度位置信息。
服务器可以和各个基站以及终端进行数据通信，以实现数据传输。其中，基站的观测数据可以包括观测到的载波、伪距观测值等，基站的位置信息可以包括基站坐标、天线信息编码等。终端的观测数据可以包括与基站相同历元观测到的载波、伪距观测值等。可以理解的是，基站所获取到的数据还可以包括星历数据，终端所获取到的数据同样可以包括星历数据。
在本实施例中，各个基站能够将其观测数据和位置信息发送至服务器，终端则可以将其概略位置信息发送至服务器。服务器根据基站观测数据和位置信息以及终端上传的概略位置信息，能够生成该终端对应的差分改正数和质量因子，终端在接收到差分改正数和质量因子后，能够根据质量因子选择相应的RTK滤波模式，从而在电离层活跃时选择能够有效消除电离层残余误差的RTK滤波模式，还能够在确定的RTK滤波模式下，根据差分改正数与自身的观测数据进行滤波解算，以固定整周模糊度，并得到终端的高精度位置信息。通过服务器播发表征大气误差模型精度的质量因子，能够使得终端根据该质量因子从无电离层组合滤波模式或非组合滤波模式中选择RTK滤波模式，在电离层活跃期时能够选择消除电离层误差的无电离层组合滤波模式进行滤波解算，以消除电离层误差的影响，实现电离层活跃期的高精度定位。
在S110中，各个基站可以上传观测到的观测数据和位置信息，终端可以上传其概略位置信息。服务器在接收到各个基站的观测数据和位置信息以及终端的概略位置信息后，可以生成差分改正数和质量因子，并向终端播发差分改正数和质量因子。
在S120中，终端在接收到差分改正数和质量因子后，可以根据接收到的质量因子确定进行滤波解算的RTK滤波模式。RTK滤波模式可以为无电离层组合滤波模式或非组合滤波模式。其中，若RTK滤波模式采用无电离层组合滤波模式，能够在电离层的残余误差较大时消除电离层误差的影响，实现电离层活跃期的终端高精度定位；若RTK滤波模式采用非组合滤波模式，则能够在电离层不活跃时实现终端的高精度定位。
作为一个可选实施例，请参照图2，为了在电离层活跃时消除电离层误差的影响，上述S120，可以包括：
S210，根据接收到的质量因子判断当前历元的电离层是否活跃；
S220，在当前历元的电离层活跃时，获取第一时间段内电离层活跃的历元数量；
S230，在电离层活跃的历元数量达到预设第一历元阈值时，将所述RTK滤波模式确定为无电离层组合滤波模式。
在本实施例中，终端可以根据接收到的质量因子来判断当前历元的电离层是否活跃，若确定当前历元的电离层活跃，则获取第一时间段内电离层活跃的历元数量，并在达到预设值时确定采用的RTK滤波模式为无电离层组合模式，从而根据无电离层组合模式对终端的观测数据和终端接收到的差分改正数进行滤波解算，以消除电离层活跃时电离层误差的影响，从而得到终端的高精度位置信息。
在S210中，终端可以在接收到质量因子后，根据预设的判断规则对质量因子进行判断，以确定当前历元的电离层是否活跃。在满足预设的判断规则时，即确定当前历元的电离层活跃。
在S220中，在终端确定当前历元的电离层活跃时，可以对当前历元为止的第一时间段内，电离层活跃的历元数量进行计数，以获取第一时间段内电离层活跃的历元数量。预设的判断规则可以是第一时间段内电离层活跃的历元数量是否达到预设第一历元阈值。
例如，在终端确定当前历元电离层活跃时，可以将当前历元对应的电离层活跃值设置为“活跃”或者“1”，并存储在数据库中。同样地，在终端确定当前历元电离层为不活跃时，可以将当前历元对应的电离层活跃值设置为“不活跃”或者“0”，并存储在数据库中。在终端确定当前历元电离层活跃时，还可以从数据库中获取过去第一时间段内电离层活跃值为“活跃”或者“1”的数量，从而确定第一时间段内电离层活跃的历元数量。
在S230中，终端在获取第一时间段内电离层活跃的历元数量后，可以将电离层活跃的历元数量与预设第一历元阈值进行比较，若电离层活跃的历元数量达到预设第一历元阈值时，则终端可以将采用的RTK滤波模式确定为无电离层组合滤波模式。
无电离层组合滤波模式中，终端可以根据自身的观测数据以及无电离层组合公式计算得到无电离层组合观测值。无电离层组合公式如下：
Figure BDA0003345977080000111
Figure BDA0003345977080000112
其中，P1,P2分别为第一频、第二频的伪距观测值(单位：m)，φ1,φ2分别为第一频、第二频的载波观测值(单位：周)，PIF,φIF分别为伪距和载波的无电离层组合观测值。通过获取无电离层组合观测值，并用无电离层组合观测值代替伪距观测值和载波观测值进行滤波解算，能够在电离层活跃时消除电离层的残余误差，得到高精度的终端位置信息。
可以理解的是，在上述实施例中，预设的判断规则还可以是判断第一时间段内电离层连续活跃的历元数量是否达到预设第一历元阈值。在该判断规则中，每当终端确定当前历元的电离层活跃时，可以从当前历元开始，获取数量为预设第一历元阈值的历史历元电离层活跃值。例如，在预设第一历元阈值为10时，可以从当前历元开始往前，获取连续10个历元的电离层活跃值，若该10个历元的电离层活跃值均为“活跃”或者“1”，则终端可以将RTK滤波模式确定为无电离层组合滤波模式。若存在至少一个历元的电离层活跃值为“不活跃”或者“0”，则不对RTK滤波模式进行调整。
作为一个可选实施例，请参照图3，为了在电离层不活跃时实现高精度定位，上述S230之后，可以包括：
S310，在当前历元的电离层不活跃时，获取第二时间段内电离层不活跃的历元数量；
S320，在电离层不活跃的历元数量达到预设第二历元阈值时，所述RTK滤波模式切换为非组合滤波模式。
在本实施例中，终端在确定当前历元的电离层不活跃时，获取第二时间段内电离层不活跃的历元数量，并在第二时间段内电离层不活跃的历元数量达到相应的预设值时确定采用的RTK滤波模式为非组合滤波模式。终端通过非组合滤波模式对终端的观测数据和终端接收到的差分改正数进行滤波解算，能够消除各种误差对整周模糊度的影响，从而快速固定整周模糊度，计算出终端的高精度位置信息。
在S310中，终端确定当前历元的电离层不活跃时，可以对当前历元为止的第二时间段内，电离层活跃的历元数量进行计数，以获取第二时间段内电离层不活跃的历元数量。例如，在终端确定当前历元电离层不活跃时，可以将当前历元对应的电离层活跃值设置为“不活跃”或者“0”，并存储在数据库中。在终端确定当前历元电离层不活跃时，还可以从数据库中获取过去第二时间段内电离层活跃值为“不活跃”或者“0”的数量，从而确定第二时间段内电离层活跃的历元数量。
在S320中，终端在获取第二时间段内电离层不活跃的历元数量后，可以将电离层不活跃的历元数量与预设第二历元阈值进行比较，若电离层不活跃的历元数量达到预设第二历元阈值时，则终端可以将采用的RTK滤波模式确定为非组合滤波模式。
同样地，上述预设的判断规则还可以是电离层连续不活跃的历元数量是否达到预设第二历元阈值。即，在终端确定当前历元的电离层不活跃时，可以获取当前历元之前的连续多个历元的电离层活跃值，若存在预设第二历元阈值个连续的历元，其电离层活跃值为“不活跃”或者“0”，则终端可以将RTK滤波模式切换为非组合滤波模式。若在连续多个历元的电离层活跃值中，存在至少一个电离层活跃值为“活跃”或者“1”，则不对RTK滤波模式进行调整。
作为一个可选实施例，请参照图4，为了在根据质量因子判断电离层是否活跃时消除公共误差的影响，质量因子可以包括对应多个卫星的电离层质量因子和对应多个卫星的对流层质量因子，上述S210，可以包括：
S410，将接收到的多个卫星的电离层质量因子和多个卫星的对流层质量因子分别进行从小到大排序；
S420，判断所述第一分位的电离层质量因子是否大于所述第二分位的对流层质量因子，以及所述第一分位的电离层质量因子是否大于预设最小活跃阈值，其中第一分位小于或等于第二分位；
S210之后，还包括：
S430，在所述第一分位的电离层质量因子大于所述第二分位的对流层质量因子，且所述第一分位的电离层质量因子大于预设最小活跃阈值时，确定当前历元的电离层活跃；
作为另一个可选实施例，上述S210，可以包括：
S440，将接收到的多个卫星的电离层质量因子和多个卫星的对流层质量因子分别进行从大到小排序；
S450，判断所述第三分位的电离层质量因子是否大于所述第四分位的对流层质量因子，以及所述第三分位的电离层质量因子是否大于预设最小活跃阈值，其中第三分位大于或等于第四分位；
S210之后，还包括：
S460，在所述第三分位的电离层质量因子大于所述第四分位的对流层质量因子，且所述第三分位的电离层质量因子大于预设最小活跃阈值时，确定当前历元的电离层活跃。
在本实施例中，终端所接收到的质量因子，包括多个卫星的电离层质量因子和多个卫星的对流层质量因子。终端在接收到多个卫星的电离层质量因子和多个卫星的对流层质量因子后，可以将电离层质量因子进行排序，将对流层质量因子也进行排序。在排序后可以获取特定分位的电离层质量因子以及特定分位的对流层质量因子，并根据获取到的电离层质量因子、对流层质量因子以及预设最小活跃阈值的大小来确定当前历元的电离层是否活跃。在计算质量因子时，还会存在其他影响因素，导致得到的质量因子存在一定误差。通过将电离层质量因子与对流层质量因子进行对比，能够排除基站观测、基站网络拓扑结构等其他因素造成的公共误差的影响，从而准确地判断当前电离层的活跃状态。
在S410中，终端可以将接收到的多个卫星的电离层质量因子按照从小到大进行排序，并将接收到的多个卫星的对流层质量因子按照从小到大进行排序。
在S420中，终端在排序后的电离层质量因子中，可以确定第一分位的电离层质量因子；在排序后的对流层质量因子中，则可以确定第二分位的对流层质量因子。其中第一分位的分位值可以为小于或等于第二分位的分位值。例如，第一分位可以为75％，第二分位可以为90％。其中，第一分位为75％的电离层质量因子σI,75，是指由小到大的排序中，小于该电离层质量因子的电离层质量因子数量占总数量的75％。同样地，第二分位为90％的对流层质量因子σT,90，是指由小到大的排序中，小于该对流层质量因子的对流层质量因子数量占总数量的90％。
在获取到第一分位的电离层质量因子和第二分位的对流层质量因子后，终端可以对第一分位的电离层质量因子和第二分位的对流层质量因子进行比较，判断第一分位的电离层质量因子是否大于第二分位的对流层质量因子。终端还可以对第一分位的电离层质量因子是否大于预设最小活跃阈值进行判断。
在S430中，终端可以在第一分位的电离层质量因子大于第二分位的对流层质量因子，并且第一分位的电离层质量因子大于预设最小活跃阈值时，确定当前历元的电离层活跃。
例如，在获取到第一分位为75％的电离层质量因子σI,75以及第二分位为90％的对流层质量因子σT,90后，可以确定预先设置的最小活跃阈值Ilow，并判断是否σI,75＞σT,90且σI,75＞Ilow。在σI,75＞σT,90且σI,75＞Ilow时，可以确定当前历元的电离层活跃，否则可以确定当前历元的电离层平静或不活跃。
可以理解的是，影响质量因子的因素不仅仅是电离层建模误差或者对流层建模误差，还包括基站网络的几何拓扑结构以及基站的观测噪声。若仅仅只是将电离层质量因子与预设的最小活跃阈值进行大小比较，难以判断此时电离层是否为活跃状态。因此，还可以将电离层质量因子与对流层质量因子进行对比，以排除基站观测、基站网络拓扑结构等其他因素造成的公共误差的影响。
在S440中，终端可以将接收到的多个卫星的电离层质量因子按照从大到小进行排序，并将接收到的多个卫星的对流层质量因子按照从大到小进行排序。
在S450中，终端在排序后的电离层质量因子中，可以确定第三分位的电离层质量因子；在排序后的对流层质量因子中，则可以确定第四分位的对流层质量因子。其中第三分位的分位值大于或等于第四分位的分位值。
在S460中，终端可以在第三分位的电离层质量因子大于第四分位的对流层质量因子，并且第三分位的电离层质量因子大于预设最小活跃阈值时，确定当前历元的电离层活跃。
在S130中，终端可以根据确定的RTK滤波模式对终端的观测数据和终端接收到的差分改正数进行滤波解算，以得到终端的高精度位置信息。例如，终端根据接收到的质量因子，可以在确定电离层为活跃期时，确定RTK滤波模式为无电离层组合滤波模式；在确定电离层不为活跃期时，确定RTK滤波模式为非组合滤波模式。
在电离层不为活跃期时，终端能够通过非组合滤波模式对差分改正数和终端的观测数据进行滤波解算，从而在在电离层不活跃时期得到终端的高精度位置信息。
在电离层为活跃期时，终端能够通过无电离层组合滤波模式，根据两个频率的载波观测值和伪距观测值计算出无电离层组合观测值，从而使得终端在滤波解算时能够消除电离层误差的影响，在电离层活跃期也能够根据差分改正数和终端的观测数据得到终端的高精度位置信息。
图5示出了本申请另一个实施例提供的网络RTK抗电离层干扰定位方法的流程示意图。网络RTK抗电离层干扰定位方法应用于服务器，网络RTK抗电离层干扰定位方法包括：
S510，根据基站网络中各个基站的观测数据和位置信息，生成多条基线的电离层误差和对流层误差，所述基线为两个所述基站相连形成；
S520，基于多条基线的电离层误差和对流层误差以及各个基站的位置信息进行建模，得到大气误差模型；
S530，根据所述大气误差模型的模型参数计算所述基站网络中多个格网点位置的差分改正数以及多个格网点位置的质量因子；
S540，根据终端发送的概略位置信息从多个格网点位置中确定距离终端概略位置最近的目标格网点；
S550，向所述终端播发所述目标格网点的差分改正数和质量因子，以使所述终端根据所述目标格网点的质量因子确定RTK滤波模式，并根据所述RTK滤波模式采用差分改正数进行滤波解算以得到所述终端的高精度位置信息，所述RTK滤波模式包括无电离层组合滤波模式或非组合滤波模式。
在本实施例中，服务器能够根据基站网络中各个基站发送的数据生成多条基线的电离层误差和对流层误差，并根据基线的误差进行建模，得到大气误差模型。通过大气误差模型可以计算出基站网络覆盖的区域内各个格网点位置的差分改正数和质量因子。在接收到终端发送的概略位置后，可以根据概略位置从多个格网点中确定距离最近的格网点，并将该格网点对应的差分改正数和质量因子播发至终端，以使终端根据质量因子确定RTK滤波模式后，根据差分改正数和终端的观测数据进行滤波解算，得到终端的高精度位置信息。
在S510中，服务器可以接收基站网络中各个基站分别发送的观测数据和位置信息。在基站网络的各个基站中，将任意两个基站相连即可形成一条基线，根据该两个基站的观测数据和位置信息即可生成该条基线的电离层误差和对流层误差。在将各个基站进行两两相连后，即可生成基站网络中多条基线分别对应的电离层误差和对流层误差。
作为一个可选实施例，请参照图6，为了生成基线的电离层误差和对流层误差，上述S510，可以包括：
S610，根据基站网络中各个基站的观测数据和位置信息，解算多条基线的双差模糊度；
S620，根据多条基线的双差模糊度计算多条基线的电离层误差值和对流层误差值。
在本实施例中，服务器可以获取基站网络中各个基站的观测数据和位置信息。在选择两个基站相连形成一条基线时，根据这两个基站的观测数据和位置信息即可计算得到该条基线的双差模糊度。在将基线的双差模糊度进行固定后，可以反算出该基线的电离层误差值和对流层误差值。通过对每条基线依次进行反算，即可得到多条基线的电离层误差值和对流层误差值。可以理解的是，基线的电离层误差值和对流层误差值即为基线对应的差分改正数。
在S610中，通过多个基站两两相连，即可形成多条基线，多条基线以及基站可以构成网形，从而形成基站网络。根据各个基站的观测数据和位置信息，即可解算出对应的基线的双差模糊度。
在S620中，根据基线固定的双差模糊度，可以通过相应的公式计算得到基线的电离层误差值和对流层误差值。反算公式如下：
Figure BDA0003345977080000171
Figure BDA0003345977080000172
其中，
Figure BDA0003345977080000173
为电离层误差值，
Figure BDA0003345977080000174
为对流层误差值。f1和f2分别为基站观测信号的第一和第二频率(单位：MHz)，
Figure BDA0003345977080000175
和
Figure BDA0003345977080000176
分别为第一频率和第二频率恢复了模糊度的载波观测值(单位：m)，
Figure BDA0003345977080000177
为双差真空几何距离值(单位：m)，
Figure BDA0003345977080000178
可以根据卫星坐标和基站精确坐标计算得到。
在S520中，服务器在分别计算得到多条基线的电离层误差和对流层误差后，可以根据多条基线的电离层误差和对流层误差以及各个基站的位置信息，对该基站网络所覆盖的区域进行大气误差模型的建模，并解出大气误差模型的模型参数。在建立大气误差模型后，对于该区域内的特定位置，可以根据特定位置的坐标以及大气误差模型的模型参数计算得出该特定位置的电离层误差和对流层误差。
例如，大气误差模型的建模方程可以为：
L＝AX；
其中，L为大气误差模型的观测值，对于电离层误差和对流层误差，可以分别进行建模。即，电离层误差的观测值为
Figure BDA0003345977080000181
对流层误差的观测值为
Figure BDA0003345977080000182
m为用于建立大气误差模型的基线数量。A为与基线的坐标位置相关的系数矩阵，X为待求解的模型参数。
在确定系数矩阵A以及观测值L后，可以对上述建模方程进行求解，以确定大气误差模型的模型参数。例如，可以采用最小二乘法对上述建模方程进行求解，得到大气误差模型的模型参数为：
Figure BDA0003345977080000183
在S530中，终端可以根据基站网络的范围确定范围内的多个格网点位置，并根据上述大气误差模型的模型参数以及每个格网点的位置坐标确定每个格网点位置对应的差分改正数和质量因子。其中，格网点位置的差分改正数即该格网点位置的电离层误差和对流层误差，质量因子为该格网点位置下各个卫星的电离层建模误差质量因子和各个卫星的对流层建模误差质量因子。
在根据上述建模方程得到大气误差模型的模型参数后，若想要对某个特定位置的电离层误差和对流层误差进行计算，可以通过差分改正数计算公式得到，差分改正数计算公式如下：
Figure BDA0003345977080000184
其中，Bcorrection为格网点位置的坐标对应的系数矩阵，Lcorrection则为解算后得到的格网点位置的电离层误差值或对流层误差值。
可以理解的是，在基站网络的范围内，可以通过内插值算法生成多个格网点位置，并通过上述差分改正数计算公式计算每个格网点位置对应的电离层误差值或对流层误差值，即计算得到每个格网点位置对应的差分改正数。内插值算法可以是反距离加权法、线性内插法、最小二乘配置法等。在计算出各个格网点位置的差分改正数后，若服务器接收到终端发送的概略位置，即可确定与概略位置最接近的格网点位置，并向终端播发该最接近的格网点位置的差分改正数和质量因子。
作为一个可选实施例，请参照图7，为了生成格网点位置的质量因子，上述S530，可以包括：
S710，根据所述大气误差模型的模型参数、各条基线对应基站的位置信息以及各条基线的电离层误差值和对流层误差值计算得到所述大气误差模型中各条基线的观测值残差；
S720，根据所述各条基线的观测值残差计算所述模型参数的精度因子矩阵；
S730，根据所述精度因子矩阵和所述基站网络中多个格网点的位置信息，分别计算得到每个格网点位置的质量因子。
在本实施例中，在确定基站网络范围内的各个格网点位置，并计算得到每个格网点位置分别对应的电离层误差和对流层误差，即每个格网点位置分别对应的差分改正数后，还可以计算得到每个格网点位置下，各个卫星的电离层建模误差质量因子和各个卫星的对流层建模误差质量因子。
在S710中，服务器在解算出大气误差模型的模型参数后，还可以将各个基线的坐标位置代入模型参数中，得到大气误差模型下基线位置对应的电离层误差值和对流层误差值。服务器根据大气误差模型下基线的电离层误差值和对流层误差值以及上述根据基线的双差模糊度反算得到的电离层误差值和对流层误差值，即可确定在大气误差模型下，实际的误差值与模型计算得到的误差值之间的观测值残差。观测值残差的计算公式如下：
Figure BDA0003345977080000191
上述公式中，A为与基线的坐标位置相关的系数矩阵，L为该基线根据双差模糊度反算得到的电离层误差值或者对流层误差值，V为大气误差模型的观测值残差。
在S720中，根据上述各条基线的观测值残差，可以计算得到相应的单位权标准差，单位权标准差的计算公式为：
Figure BDA0003345977080000201
其中，σ0为单位权标准差，n为基线个数，k为大气误差模型中建模参数的个数。
在计算得到单位权标准差后，可以根据单位权标准差计算大气误差模型建模参数的精度因子矩阵。精度因子矩阵的计算公式如下：
Figure BDA0003345977080000202
在S730中，在计算得到精度因子矩阵后，服务器可以根据基站网络中格网点的位置信息，计算得到该格网点位置的质量因子。对每个格网点的位置信息进行计算后即可得到每个格网点位置的质量因子。计算公式如下：
Figure BDA0003345977080000203
其中，Bcorrection为格网点位置位置的坐标对应的系数矩阵，D为上述精度因子矩阵。
可以理解的是，每个格网点位置的质量因子可以包括电离层建模误差质量因子和对流层建模误差质量因子。在计算大气误差模型的观测值残差时，若输入的观测值为各条基线的电离层误差值，则最终计算得到的质量因子为各个卫星的电离层建模误差质量因子。同样地，若输入的观测值为各条基线的对流层误差值，则最终计算得到的质量因子为各个卫星的对流层建模误差质量因子。
在S540中，服务器在确定基站网络中每个格网点位置的差分改正数和质量因子后，可以根据终端发送的概略位置信息确定终端对应的概略位置，并从多个格网点位置中确定与终端的概略位置最接近的格网点，作为目标格网点。
在S550中，服务器在根据终端的概略位置确定目标格网点后，可以获取目标格网点的差分改正数和质量因子，并播发给终端。终端在结合搜到目标格网点的差分改正数和质量因子后，可以根据质量因子确定RTK滤波模式，并根据确定的RTK滤波模式采用差分改正数以及终端的观测数据进行滤波解算，从而得到终端的高精度位置信息。
基于上述实施例提供的网络RTK抗电离层干扰定位方法，相应地，本申请还提供了网络RTK抗电离层干扰定位装置的具体实现方式。请参见以下实施例。
首先参见图8，本申请实施例提供的网络RTK抗电离层干扰定位装置800包括以下模块：
接收模块801，用于接收服务器播发的差分改正数和质量因子，所述差分改正数和所述质量因子为服务器根据各个基站的观测数据和位置信息以及所述终端上传的概略位置信息生成，所述质量因子用于表征生成所述差分改正数所用的大气误差模型精度；
判定模块802，用于根据接收到的质量因子确定RTK滤波模式，所述RTK滤波模式包括无电离层组合滤波模式或非组合滤波模式；
定位模块803，用于采用确定的所述RTK滤波模式对终端的观测数据和终端接收到的差分改正数进行滤波解算，得到所述终端的高精度位置信息；
或者，请参见图9，所述装置包括：
误差模块901，用于根据基站网络中各个基站的观测数据和位置信息，生成多条基线的电离层误差和对流层误差，所述基线为两个所述基站相连形成；
建模模块902，用于基于多条基线的电离层误差和对流层误差以及各个基站的位置信息进行建模，得到大气误差模型；
计算模块903，用于根据所述大气误差模型的模型参数计算所述基站网络中多个格网点位置的差分改正数以及多个格网点位置的质量因子；
匹配模块904，用于根据终端发送的概略位置信息从多个格网点位置中确定距离终端概略位置最近的目标格网点；
播发模块905，用于向所述终端播发所述目标格网点的差分改正数和质量因子，以使所述终端根据所述目标格网点的质量因子确定RTK滤波模式，并根据所述RTK滤波模式采用差分改正数进行滤波解算以得到所述终端的高精度位置信息，所述RTK滤波模式包括无电离层组合滤波模式或非组合滤波模式。
在本实施例中，各个基站能够将其观测数据和位置信息发送至服务器，终端则可以将其概略位置信息发送至服务器。服务器根据基站观测数据和位置信息以及终端上传的概略位置信息，能够生成该终端对应的差分改正数和质量因子，终端在接收到差分改正数和质量因子后，能够根据质量因子选择相应的RTK滤波模式，从而在电离层活跃时选择能够有效消除电离层残余误差的RTK滤波模式，还能够在确定的RTK滤波模式下，根据差分改正数与自身的观测数据进行滤波解算，以固定整周模糊度，并得到终端的高精度位置信息。通过服务器播发表征大气误差模型精度的质量因子，能够使得终端根据该质量因子从无电离层组合滤波模式或非组合滤波模式中选择RTK滤波模式，在电离层活跃期时能够选择消除电离层误差的无电离层组合滤波模式进行滤波解算，以消除电离层误差的影响，实现电离层活跃期的高精度定位。
作为本申请的一种实现方式，为了在电离层活跃时消除电离层误差的影响，上述判定模块802还可以包括：
判定单元，用于根据接收到的质量因子判断当前历元的电离层是否活跃；
第一获取单元，用于在当前历元的电离层活跃时，获取第一时间段内电离层活跃的历元数量；
第一确定单元，用于在电离层活跃的历元数量达到预设第一历元阈值时，将所述RTK滤波模式确定为无电离层组合滤波模式。
作为本申请的一种实现方式，为了在电离层不活跃时实现高精度定位，上述判定模块802还可以包括：
第二获取单元，用于在当前历元的电离层不活跃时，获取第二时间段内电离层不活跃的历元数量；
第二确定单元，用于在电离层不活跃的历元数量达到预设第二历元阈值时，所述RTK滤波模式切换为非组合滤波模式。
作为本申请的一种实现方式，为了根据质量因子判断电离层是否活跃，质量因子包括对应多个卫星的电离层质量因子和对应多个卫星的对流层质量因子，上述判定单元还可以包括：
第一排序子单元，用于将接收到的多个卫星的电离层质量因子和多个卫星的对流层质量因子分别进行从小到大排序；
第一比较子单元，用于判断所述第一分位的电离层质量因子是否大于所述第二分位的对流层质量因子，以及所述第一分位的电离层质量因子是否大于预设最小活跃阈值，其中第一分位小于或等于第二分位；
所述判定模块802还包括：
第一活跃子单元，用于在所述第一分位的电离层质量因子大于所述第二分位的对流层质量因子，且所述第一分位的电离层质量因子大于预设最小活跃阈值时，确定当前历元的电离层活跃；
或者，上述判定单元还可以包括：
第二排序子单元，用于将接收到的多个卫星的电离层质量因子和多个卫星的对流层质量因子分别进行从大到小排序；
第二排序子单元，用于判断所述第三分位的电离层质量因子是否大于所述第四分位的对流层质量因子，以及所述第三分位的电离层质量因子是否大于预设最小活跃阈值，其中第三分位大于或等于第四分位；
所述判定模块802还包括：
第二排序子单元，用于在所述第三分位的电离层质量因子大于所述第四分位的对流层质量因子，且所述第三分位的电离层质量因子大于预设最小活跃阈值时，确定当前历元的电离层活跃。
作为本申请的一种实现方式，为了生成基线的电离层误差和对流层误差，上述误差模块901还可以包括：
固定单元，用于根据基站网络中各个基站的观测数据和位置信息，解算多条基线的双差模糊度；
误差单元，用于根据多条基线的双差模糊度计算多条基线的电离层误差值和对流层误差值。
作为本申请的一种实现方式，为了生成格网点位置的质量因子，上述计算模块903还可以包括：
残差单元，用于根据所述大气误差模型的模型参数、各条基线对应基站的位置信息以及各条基线的电离层误差值和对流层误差值计算得到所述大气误差模型中各条基线的观测值残差；
精度单元，用于根据所述各条基线的观测值残差计算所述模型参数的精度因子矩阵；
计算单元，用于根据所述精度因子矩阵和所述基站网络中多个格网点的位置信息，分别计算得到每个格网点位置的质量因子。
本申请实施例提供的网络RTK抗电离层干扰的定位装置能够实现图1至图7的方法实施例中的各个过程，为避免重复，这里不再赘述。
本申请一个实施例提供的网络RTK抗电离层干扰定位系统，网络RTK抗电离层干扰定位系统包括服务器和终端，系统包括：
所述服务器根据基站网络中各个基站的观测数据和位置信息，生成多条基线的电离层误差和对流层误差，所述基线为两个所述基站相连形成；
所述服务器基于多条基线的电离层误差和对流层误差以及各个基站的位置信息进行建模，得到大气误差模型；
所述服务器根据所述大气误差模型的模型参数计算所述基站网络中多个格网点位置的差分改正数以及多个格网点位置的质量因子；
所述服务器根据终端发送的概略位置信息从多个格网点位置中确定距离终端概略位置最近的目标格网点；
所述服务器向所述终端播发所述目标格网点的差分改正数和质量因子；
所述终端接收服务器播发的差分改正数和质量因子，所述差分改正数和所述质量因子为服务器根据各个基站的观测数据和位置信息以及所述终端上传的概略位置信息生成，所述质量因子用于表征生成所述差分改正数所用的大气误差模型精度；
所述终端根据接收到的质量因子确定RTK滤波模式，所述RTK滤波模式包括无电离层组合滤波模式或非组合滤波模式；
所述终端采用确定的所述RTK滤波模式对终端的观测数据和终端接收到的差分改正数进行滤波解算，得到所述终端的高精度位置信息。
本申请实施例提供的网络RTK抗电离层干扰的定位系统能够实现图1至图7的方法实施例中的各个过程，为避免重复，这里不再赘述。
图10示出了本申请实施例提供的网络RTK抗电离层干扰定位设备的硬件结构示意图。
在网络RTK抗电离层干扰定位设备可以包括处理器1001以及存储有计算机程序指令的存储器1002。
具体地，上述处理器1001可以包括中央处理器(CPU)，或者特定集成电路(Application Specific Integrated Circuit，ASIC)，或者可以被配置成实施本申请实施例的一个或多个集成电路。
存储器1002可以包括用于数据或指令的大容量存储器。举例来说而非限制，存储器1002可包括硬盘驱动器(Hard Disk Drive，HDD)、软盘驱动器、闪存、光盘、磁光盘、磁带或通用串行总线(Universal Serial Bus，USB)驱动器或者两个或更多个以上这些的组合。在合适的情况下，存储器1002可包括可移除或不可移除(或固定)的介质。在合适的情况下，存储器1002可在综合网关容灾设备的内部或外部。在特定实施例中，存储器1002是非易失性固态存储器。
存储器可包括只读存储器(ROM)，随机存取存储器(RAM)，磁盘存储介质设备，光存储介质设备，闪存设备，电气、光学或其他物理/有形的存储器存储设备。因此，通常，存储器包括一个或多个编码有包括计算机可执行指令的软件的有形(非暂态)计算机可读存储介质(例如，存储器设备)，并且当该软件被执行(例如，由一个或多个处理器)时，其可操作来执行参考根据本公开的一方面的方法所描述的操作。
处理器1001通过读取并执行存储器1002中存储的计算机程序指令，以实现上述实施例中的任意一种网络RTK抗电离层干扰定位方法。
在一个示例中，网络RTK抗电离层干扰定位设备还可包括通信接口1003和总线1010。其中，如图10所示，处理器1001、存储器1002、通信接口1003通过总线1010连接并完成相互间的通信。
通信接口1003，主要用于实现本申请实施例中各模块、装置、单元和/或设备之间的通信。
总线1010包括硬件、软件或两者，将网络RTK抗电离层干扰定位设备的部件彼此耦接在一起。举例来说而非限制，总线可包括加速图形端口(AGP)或其他图形总线、增强工业标准架构(EISA)总线、前端总线(FSB)、超传输(HT)互连、工业标准架构(ISA)总线、无限带宽互连、低引脚数(LPC)总线、存储器总线、微信道架构(MCA)总线、外围组件互连(PCI)总线、PCI-Express(PCI-X)总线、串行高级技术附件(SATA)总线、视频电子标准协会局部(VLB)总线或其他合适的总线或者两个或更多个以上这些的组合。在合适的情况下，总线1010可包括一个或多个总线。尽管本申请实施例描述和示出了特定的总线，但本申请考虑任何合适的总线或互连。
该网络RTK抗电离层干扰定位设备可以基于上述实施例，从而实现结合图1至图9描述的网络RTK抗电离层干扰定位方法和装置。
另外，结合上述实施例中的网络RTK抗电离层干扰定位方法，本申请实施例可提供一种计算机存储介质来实现。该计算机存储介质上存储有计算机程序指令；该计算机程序指令被处理器执行时实现上述实施例中的任意一种网络RTK抗电离层干扰定位方法。
需要明确的是，本申请并不局限于上文所描述并在图中示出的特定配置和处理。为了简明起见，这里省略了对已知方法的详细描述。在上述实施例中，描述和示出了若干具体的步骤作为示例。但是，本申请的方法过程并不限于所描述和示出的具体步骤，本领域的技术人员可以在领会本申请的精神后，作出各种改变、修改和添加，或者改变步骤之间的顺序。
以上的结构框图中所示的功能块可以实现为硬件、软件、固件或者它们的组合。当以硬件方式实现时，其可以例如是电子电路、专用集成电路(ASIC)、适当的固件、插件、功能卡等等。当以软件方式实现时，本申请的元素是被用于执行所需任务的程序或者代码段。程序或者代码段可以存储在机器可读介质中，或者通过载波中携带的数据信号在传输介质或者通信链路上传送。“机器可读介质”可以包括能够存储或传输信息的任何介质。机器可读介质的例子包括电子电路、半导体存储器设备、ROM、闪存、可擦除ROM(EROM)、软盘、CD-ROM、光盘、硬盘、光纤介质、射频(RF)链路，等等。代码段可以经由诸如因特网、内联网等的计算机网络被下载。
还需要说明的是，本申请中提及的示例性实施例，基于一系列的步骤或者装置描述一些方法或系统。但是，本申请不局限于上述步骤的顺序，也就是说，可以按照实施例中提及的顺序执行步骤，也可以不同于实施例中的顺序，或者若干步骤同时执行。
上面参考根据本公开的实施例的方法、装置和计算机程序产品的流程图和/或框图描述了本公开的各方面。应当理解，流程图和/或框图中的每个方框以及流程图和/或框图中各方框的组合可以由计算机程序指令实现。这些计算机程序指令可被提供给通用计算机、专用计算机、或其它可编程数据处理装置的处理器，以产生一种机器，使得经由计算机或其它可编程数据处理装置的处理器执行的这些指令使能对流程图和/或框图的一个或多个方框中指定的功能/动作的实现。这种处理器可以是但不限于是通用处理器、专用处理器、特殊应用处理器或者现场可编程逻辑电路。还可理解，框图和/或流程图中的每个方框以及框图和/或流程图中的方框的组合，也可以由执行指定的功能或动作的专用硬件来实现，或可由专用硬件和计算机指令的组合来实现。
以上，仅为本申请的具体实施方式，所属领域的技术人员可以清楚地了解到，为了描述的方便和简洁，上述描述的系统、模块和单元的具体工作过程，可以参考前述方法实施例中的对应过程，在此不再赘述。应理解，本申请的保护范围并不局限于此，任何熟悉本技术领域的技术人员在本申请揭露的技术范围内，可轻易想到各种等效的修改或替换，这些修改或替换都应涵盖在本申请的保护范围之内。

Claims (11)

1.一种网络RTK抗电离层干扰定位方法，应用于终端，其特征在于，所述方法包括：
接收服务器播发的差分改正数和质量因子，所述差分改正数和所述质量因子为服务器根据各个基站的观测数据和位置信息以及所述终端上传的概略位置信息生成，所述质量因子用于表征生成所述差分改正数所用的大气误差模型精度；
根据接收到的质量因子确定RTK滤波模式，所述RTK滤波模式包括无电离层组合滤波模式或非组合滤波模式；
采用确定的所述RTK滤波模式对终端的观测数据和终端接收到的差分改正数进行滤波解算，得到所述终端的高精度位置信息。
2.根据权利要求1所述的网络RTK抗电离层干扰定位方法，其特征在于，所述根据接收到的质量因子确定RTK滤波模式，包括：
根据接收到的质量因子判断当前历元的电离层是否活跃；
在当前历元的电离层活跃时，获取第一时间段内电离层活跃的历元数量；
在电离层活跃的历元数量达到预设第一历元阈值时，将所述RTK滤波模式确定为无电离层组合滤波模式。
3.根据权利要求2所述的网络RTK抗电离层干扰定位方法，其特征在于，所述将所述RTK滤波模式确定为无电离层组合滤波模式之后，还包括：
在当前历元的电离层不活跃时，获取第二时间段内电离层不活跃的历元数量；
在电离层不活跃的历元数量达到预设第二历元阈值时，所述RTK滤波模式切换为非组合滤波模式。
4.根据权利要求2所述的网络RTK抗电离层干扰定位方法，其特征在于，所述质量因子包括对应多个卫星的电离层质量因子和对应多个卫星的对流层质量因子；所述根据接收到的质量因子判断当前历元的电离层是否活跃，包括：
将接收到的多个卫星的电离层质量因子和多个卫星的对流层质量因子分别进行从小到大排序；
判断第一分位的电离层质量因子是否大于第二分位的对流层质量因子，以及所述第一分位的电离层质量因子是否大于预设最小活跃阈值，其中第一分位小于或等于第二分位；
所述根据接收到的质量因子判断当前历元的电离层是否活跃之后，还包括：
在所述第一分位的电离层质量因子大于所述第二分位的对流层质量因子，且所述第一分位的电离层质量因子大于预设最小活跃阈值时，确定当前历元的电离层活跃；或者
将接收到的多个卫星的电离层质量因子和多个卫星的对流层质量因子分别进行从大到小排序；
判断第三分位的电离层质量因子是否大于第四分位的对流层质量因子，以及所述第三分位的电离层质量因子是否大于预设最小活跃阈值，其中第三分位大于或等于第四分位；
所述根据接收到的质量因子判断当前历元的电离层是否活跃之后，还包括：
在所述第三分位的电离层质量因子大于所述第四分位的对流层质量因子，且所述第三分位的电离层质量因子大于预设最小活跃阈值时，确定当前历元的电离层活跃。
5.一种网络RTK抗电离层干扰定位方法，应用于服务器，其特征在于，所述方法包括：
根据基站网络中各个基站的观测数据和位置信息，生成多条基线的电离层误差和对流层误差，所述基线为两个所述基站相连形成；
基于多条基线的电离层误差和对流层误差以及各个基站的位置信息进行建模，得到大气误差模型；
根据所述大气误差模型的模型参数计算所述基站网络中多个格网点位置的差分改正数以及多个格网点位置的质量因子；
根据终端发送的概略位置信息从多个格网点位置中确定距离终端概略位置最近的目标格网点；
向所述终端播发所述目标格网点的差分改正数和质量因子，以使所述终端根据所述目标格网点的质量因子确定RTK滤波模式，并根据所述RTK滤波模式采用差分改正数进行滤波解算以得到所述终端的高精度位置信息，所述RTK滤波模式包括无电离层组合滤波模式或非组合滤波模式。
6.根据权利要求5所述的网络RTK抗电离层干扰定位方法，其特征在于，所述根据基站网络中各个基站的观测数据和位置信息，生成多条基线的电离层误差和对流层误差，包括：
根据基站网络中各个基站的观测数据和位置信息，解算多条基线的双差模糊度；
根据多条基线的双差模糊度计算多条基线的电离层误差值和对流层误差值。
7.根据权利要求5所述的网络RTK抗电离层干扰定位方法，其特征在于，所述根据所述大气误差模型的模型参数计算所述基站网络中多个格网点位置的差分改正数以及多个格网点位置的质量因子，包括：
根据所述大气误差模型的模型参数、各条基线对应基站的位置信息以及各条基线的电离层误差值和对流层误差值计算得到所述大气误差模型中各条基线的观测值残差；
根据所述各条基线的观测值残差计算所述模型参数的精度因子矩阵；
根据所述精度因子矩阵和所述基站网络中多个格网点的位置信息，分别计算得到每个格网点位置的质量因子。
8.一种网络RTK抗电离层干扰的定位装置，其特征在于，所述装置包括：
接收模块，用于接收服务器播发的差分改正数和质量因子，所述差分改正数和所述质量因子为服务器根据各个基站的观测数据和位置信息以及所述终端上传的概略位置信息生成，所述质量因子用于表征生成所述差分改正数所用的大气误差模型精度；
判定模块，用于根据接收到的质量因子确定RTK滤波模式，所述RTK滤波模式包括无电离层组合滤波模式或非组合滤波模式；
定位模块，用于采用确定的所述RTK滤波模式对终端的观测数据和终端接收到的差分改正数进行滤波解算，得到所述终端的高精度位置信息；
或者，所述装置包括：
误差模块，用于根据基站网络中各个基站的观测数据和位置信息，生成多条基线的电离层误差和对流层误差，所述基线为两个所述基站相连形成；
建模模块，用于基于多条基线的电离层误差和对流层误差以及各个基站的位置信息进行建模，得到大气误差模型；
计算模块，用于根据所述大气误差模型的模型参数计算所述基站网络中多个格网点位置的差分改正数以及多个格网点位置的质量因子；
匹配模块，用于根据终端发送的概略位置信息从多个格网点位置中确定距离终端概略位置最近的目标格网点；
播发模块，用于向所述终端播发所述目标格网点的差分改正数和质量因子，以使所述终端根据所述目标格网点的质量因子确定RTK滤波模式，并根据所述RTK滤波模式采用差分改正数进行滤波解算以得到所述终端的高精度位置信息，所述RTK滤波模式包括无电离层组合滤波模式或非组合滤波模式。
9.一种网络RTK抗电离层干扰定位系统，包括服务器和终端，其特征在于，所述系统包括：
所述服务器根据基站网络中各个基站的观测数据和位置信息，生成多条基线的电离层误差和对流层误差，所述基线为两个所述基站相连形成；
所述服务器基于多条基线的电离层误差和对流层误差以及各个基站的位置信息进行建模，得到大气误差模型；
所述服务器根据所述大气误差模型的模型参数计算所述基站网络中多个格网点位置的差分改正数以及多个格网点位置的质量因子；
所述服务器根据终端发送的概略位置信息从多个格网点位置中确定距离终端概略位置最近的目标格网点；
所述服务器向所述终端播发所述目标格网点的差分改正数和质量因子；
所述终端接收服务器播发的差分改正数和质量因子，所述差分改正数和所述质量因子为服务器根据各个基站的观测数据和位置信息以及所述终端上传的概略位置信息生成，所述质量因子用于表征生成所述差分改正数所用的大气误差模型精度；
所述终端根据接收到的质量因子确定RTK滤波模式，所述RTK滤波模式包括无电离层组合滤波模式或非组合滤波模式；
所述终端采用确定的所述RTK滤波模式对终端的观测数据和终端接收到的差分改正数进行滤波解算，得到所述终端的高精度位置信息。
10.一种网络RTK抗电离层干扰定位设备，其特征在于，所述网络RTK抗电离层干扰定位设备包括：处理器以及存储有计算机程序指令的存储器；
所述处理器执行所述计算机程序指令时实现如权利要求1-7中任一项所述的网络RTK抗电离层干扰定位方法。
11.一种计算机存储介质，其特征在于，所述计算机存储介质上存储有计算机程序指令，所述计算机程序指令被处理器执行时实现如权利要求1-7中任一项所述的网络RTK抗电离层干扰定位方法。
CN202111322856.7A 2021-11-09 2021-11-09 网络rtk抗电离层干扰定位方法及相关设备 Pending CN116106948A (zh)

Priority Applications (2)

Application Number Priority Date Filing Date Title
CN202111322856.7A CN116106948A (zh) 2021-11-09 2021-11-09 网络rtk抗电离层干扰定位方法及相关设备
PCT/CN2022/116380 WO2023082785A1 (zh) 2021-11-09 2022-08-31 网络rtk抗电离层干扰定位方法、装置、系统、设备及存储介质

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202111322856.7A CN116106948A (zh) 2021-11-09 2021-11-09 网络rtk抗电离层干扰定位方法及相关设备

Publications (1)

Publication Number Publication Date
CN116106948A true CN116106948A (zh) 2023-05-12

Family

ID=86266050

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202111322856.7A Pending CN116106948A (zh) 2021-11-09 2021-11-09 网络rtk抗电离层干扰定位方法及相关设备

Country Status (2)

Country Link
CN (1) CN116106948A (zh)
WO (1) WO2023082785A1 (zh)

Families Citing this family (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN115276769A (zh) * 2022-07-22 2022-11-01 福建鼎旸信息科技股份有限公司 一种卫星数据精准播发与服务覆盖方法及系统
CN116931007B (zh) * 2023-08-31 2023-12-08 腾讯科技（深圳）有限公司 电离层延迟处理方法、装置、设备及存储介质

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN107797126A (zh) * 2017-09-26 2018-03-13 东南大学 基于星型网络的bds/gps广播式网络rtk算法
CN108415049A (zh) * 2018-04-19 2018-08-17 千寻位置网络有限公司 提高网络rtk双差宽巷模糊度固定正确率的方法
CN109828288A (zh) * 2019-01-23 2019-05-31 东南大学 一种基于区域cors的实时电离层建模与监测方法
CN112748455A (zh) * 2019-10-29 2021-05-04 千寻位置网络有限公司 考虑电离层活跃度的网络rtk解算方法及系统

Family Cites Families (16)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
RU2237257C2 (ru) * 2002-01-25 2004-09-27 Закрытое акционерное общество "Конструкторское бюро навигационных систем" Способ устранения влияния тропосферных и ионосферных ошибок измерения в одночастотных приёмниках спутниковой навигации
US7868820B2 (en) * 2005-09-09 2011-01-11 Trimble Navigation Limited Ionosphere modeling apparatus and methods
CN101509970A (zh) * 2009-03-20 2009-08-19 东南大学 一种网络差分改正信息的生成与发布方法
CN102323572B (zh) * 2011-06-14 2015-07-01 北京航空航天大学 一种卫星导航信号电离层差分改正数估计方法
US9557419B2 (en) * 2012-12-18 2017-01-31 Trimble Inc. Methods for generating accuracy information on an ionosphere model for satellite navigation applications
CN103605145B (zh) * 2013-12-04 2017-02-08 北极星云空间技术股份有限公司 基于gnss多频数据和cors实现网络实时动态定位的方法
CN103760572B (zh) * 2014-01-13 2016-09-07 东南大学 一种基于区域cors的单频ppp电离层加权方法
CN106093992A (zh) * 2016-08-24 2016-11-09 南京莱斯信息技术股份有限公司 一种基于cors的亚米级组合定位导航系统及导航方法
US10578748B2 (en) * 2016-10-14 2020-03-03 Eung Gi Paek High-performance time transfer using time reversal (T3R)
CN107229061B (zh) * 2017-07-18 2019-09-03 武汉大学 一种基于低轨卫星的星地差分实时精密定位方法
CN108169776B (zh) * 2017-11-23 2022-01-21 中国科学院光电研究院 基于背景模型和实测数据的电离层延迟误差修正方法
CN111142133B (zh) * 2018-11-02 2022-06-24 千寻位置网络有限公司 基于多个连续运行参考站的后处理定位方法和系统
CN110531392A (zh) * 2019-04-03 2019-12-03 南京航空航天大学 一种基于ppp算法的高精度定位方法和系统
CN111929705A (zh) * 2020-08-24 2020-11-13 中国人民解放军海军工程大学 基于动态基准站的定位方法、装置、电子设备及存储介质
CN111983654B (zh) * 2020-08-24 2021-03-19 中国矿业大学 一种基于gnss的北极区域电离层相位闪烁因子构建方法
CN112835082A (zh) * 2021-01-05 2021-05-25 广州星际互联科技有限公司 一种gnss区域增强电离层和对流层大气产品质量指标计算方法、电子设备、存储介质

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN107797126A (zh) * 2017-09-26 2018-03-13 东南大学 基于星型网络的bds/gps广播式网络rtk算法
CN108415049A (zh) * 2018-04-19 2018-08-17 千寻位置网络有限公司 提高网络rtk双差宽巷模糊度固定正确率的方法
CN109828288A (zh) * 2019-01-23 2019-05-31 东南大学 一种基于区域cors的实时电离层建模与监测方法
CN112748455A (zh) * 2019-10-29 2021-05-04 千寻位置网络有限公司 考虑电离层活跃度的网络rtk解算方法及系统

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
黄丁发，张勤，张小红，周乐韬: "卫星导航定位原理", 31 January 2015, 武汉大学出版社, pages: 157 - 158 *

Also Published As

Publication number Publication date
WO2023082785A1 (zh) 2023-05-19

Similar Documents

Publication Publication Date Title
CN108931915B (zh) 利用导航卫星的授时方法和装置、计算机可读存储介质
EP3430438B1 (en) Estimation of inter-frequency bias for ambiguity resolution in global navigation satellite system receivers
CN111045034B (zh) 基于广播星历的gnss多系统实时精密时间传递方法及系统
US10114123B2 (en) Accuracy and performance of the hybrid positioning system
CN116106948A (zh) 网络rtk抗电离层干扰定位方法及相关设备
US9291716B2 (en) Method and apparatus for determining physical characteristics of a receiving device in a navigation system
EP3106897A1 (en) Gnss receiver with an on-board capability to implement an optimal error correction mode
CN114791613A (zh) 星历预报方法和装置
CN103852774A (zh) Gnss接收器及其用于伪距和位置确定的计算方法
CN110456397B (zh) 一种多天线超短基线定位监测方法、装置及存储介质
CN112327340B (zh) 终端定位精度评估方法、装置、设备以及介质
CN114791614A (zh) 钟差预报方法和装置
CN111352137B (zh) 一种顾及广播星历误差的多模gnss异步rtk定位方法
CN114721012B (zh) 一种电离层延迟效应改正方法及装置
CN110068848B (zh) 一种高性能rtk处理技术方法
US20110181462A1 (en) System and Method for Positioning with GNSS Using Multiple Integer Candidates
Rahman et al. ECEF position accuracy and reliability in the presence of differential correction latency
CN108254762B (zh) 伪距差分定位方法及系统
CN113534206B (zh) 基于北斗地基增强系统的接入虚拟参考站快速选择方法
CN110045634B (zh) 一种gnss基准站的非差误差建模方法
KR101799876B1 (ko) Ｄｇｎｓｓ 측위 방법 및 이를 이용하는 ｇｎｓｓ 수신기
CN114152961A (zh) 一种导航系统的周跳处理方法和装置
CN111948682A (zh) 一种基于bds三频组合的伪距和载波相位随机模型构建方法
CN111158021A (zh) 电离层行扰的估计方法及系统、预警终端
CN116990840A (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