CN112327340A - 终端定位精度评估方法、装置、设备以及介质 - Google Patents
终端定位精度评估方法、装置、设备以及介质 Download PDFInfo
- Publication number
- CN112327340A CN112327340A CN202110011939.8A CN202110011939A CN112327340A CN 112327340 A CN112327340 A CN 112327340A CN 202110011939 A CN202110011939 A CN 202110011939A CN 112327340 A CN112327340 A CN 112327340A
- Authority
- CN
- China
- Prior art keywords
- satellite
- value
- pseudo
- terminal
- observation
- 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.)
- Granted
Links
- 238000011156 evaluation Methods 0.000 title claims abstract description 125
- 239000011159 matrix material Substances 0.000 claims abstract description 157
- 238000000034 method Methods 0.000 claims abstract description 75
- 238000012216 screening Methods 0.000 claims abstract description 19
- 238000003491 array Methods 0.000 claims description 28
- 238000004590 computer program Methods 0.000 claims description 17
- 238000010276 construction Methods 0.000 claims description 12
- 230000006870 function Effects 0.000 claims description 12
- 230000007704 transition Effects 0.000 claims description 9
- 238000001914 filtration Methods 0.000 claims 1
- 238000005516 engineering process Methods 0.000 abstract description 9
- 230000005540 biological transmission Effects 0.000 description 24
- 238000004891 communication Methods 0.000 description 16
- 238000012937 correction Methods 0.000 description 14
- 238000004364 calculation method Methods 0.000 description 13
- 238000010586 diagram Methods 0.000 description 11
- 238000001514 detection method Methods 0.000 description 10
- 238000005259 measurement Methods 0.000 description 8
- 238000012545 processing Methods 0.000 description 8
- 238000013178 mathematical model Methods 0.000 description 4
- 230000009286 beneficial effect Effects 0.000 description 3
- 238000012854 evaluation process Methods 0.000 description 3
- 230000004927 fusion Effects 0.000 description 3
- 230000010365 information processing Effects 0.000 description 3
- 238000010295 mobile communication Methods 0.000 description 2
- 238000012544 monitoring process Methods 0.000 description 2
- 241000287196 Asthenes Species 0.000 description 1
- 238000013473 artificial intelligence Methods 0.000 description 1
- 238000002425 crystallisation Methods 0.000 description 1
- 230000008025 crystallization Effects 0.000 description 1
- 239000005433 ionosphere Substances 0.000 description 1
- 238000013507 mapping Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 230000003287 optical effect Effects 0.000 description 1
- 230000008054 signal transmission Effects 0.000 description 1
- 238000012360 testing method Methods 0.000 description 1
- 239000005436 troposphere Substances 0.000 description 1
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/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
-
- 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
-
- 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/13—Receivers
- G01S19/14—Receivers specially adapted for specific applications
Landscapes
- Engineering & Computer Science (AREA)
- Radar, Positioning & Navigation (AREA)
- Remote Sensing (AREA)
- Computer Networks & Wireless Communication (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Position Fixing By Use Of Radio Waves (AREA)
Abstract
本申请实施例提供了一种终端定位精度评估方法、装置、设备以及介质,该方法涉及导航定位技术,可以应用在电子地图中,该方法包括:获取卫星导航星历,获取移动终端与卫星之间的初始伪距观测值,根据初始伪距观测值和卫星导航星历,确定卫星对应的卫星位置和卫星钟差;根据卫星位置、卫星钟差、移动终端对应的终端位置和终端钟差,对初始伪距观测值进行筛选,得到目标伪距观测值;根据目标伪距观测值和终端位置,构建移动终端与卫星之间的伪距观测矩阵;根据伪距观测矩阵对终端位置进行评估,确定终端位置对应的位置精度指标;位置精度指标用于表征移动终端的定位精度。采用本申请实施例,可以提高移动终端的位置评估准确性。
Description
技术领域
本申请涉及互联网技术领域,尤其涉及一种终端定位精度评估方法、装置、设备以及介质。
背景技术
全球卫星导航系统(Global Navigation Satellite System)已经形成了全球空间定位的强大应用。利用全球卫星导航系统可以对智能手机进行定位,基于智能手机的定位信息可以为使用该智能手机的用户提供相应的服务(例如,打车服务、为用户推荐附近的餐馆、酒店等)。
然而,由于现代城市的特殊环境,如高楼大厦形成城市峡谷,立交桥形成定位盲区,使得卫星导航定位的全面应用受到阻碍,在利用全球卫星导航系统对智能手机进行定位的过程中,造成智能手机的定位位置与智能手机的实际位置出现偏差,从而导致智能手机的定位准确性过低。
发明内容
本申请实施例提供一种终端定位精度评估方法、装置、设备以及介质,可以提高移动终端的位置评估准确性。
本申请实施例一方面提供了一种终端定位精度评估方法,包括:
获取卫星导航星历,获取移动终端与卫星之间的初始伪距观测值,根据初始伪距观测值和卫星导航星历,确定卫星对应的卫星位置和卫星钟差;
获取移动终端对应的终端位置和终端钟差,根据卫星位置、卫星钟差、终端位置以及终端钟差,对初始伪距观测值进行筛选,得到目标伪距观测值;
根据目标伪距观测值和终端位置,构建移动终端与卫星之间的伪距观测矩阵;
根据伪距观测矩阵对终端位置进行评估,确定终端位置对应的位置精度指标;该位置精度指标用于表征移动终端的定位精度。
本申请实施例一方面提供了一种终端定位精度评估装置,包括:
位置确定模块,用于获取卫星导航星历,获取移动终端与卫星之间的初始伪距观测值,根据初始伪距观测值和卫星导航星历,确定卫星对应的卫星位置和卫星钟差;
观测值获取模块,用于获取移动终端对应的终端位置和终端钟差,根据卫星位置、卫星钟差、终端位置以及终端钟差,对初始伪距观测值进行筛选,得到目标伪距观测值;
观测矩阵构建模块,用于根据目标伪距观测值和终端位置,构建移动终端与卫星之间的伪距观测矩阵;
精度评估模块,用于根据伪距观测矩阵对终端位置进行评估,确定终端位置对应的位置精度指标;该位置精度指标用于表征移动终端的定位精度。
其中,位置确定模块包括:
发射时间确定单元,用于获取移动终端接收到卫星发射信号的接收时间,根据初始伪距观测值、接收时间以及光速值,确定卫星发射信号的发射时间;
卫星位置确定单元,用于根据发射时间以及卫星导航星历中所包含的参数列表,确定卫星对应的卫星位置和卫星钟差。
其中,观测值获取模块包括:
距离估计值确定单元,用于根据卫星位置、卫星钟差、终端位置以及终端钟差,确定移动终端与卫星之间的距离估计值;
目标观测值确定单元,用于根据距离估计值和初始伪距观测值之间的残差值,对初始伪距观测值进行筛选,得到目标伪距观测值。
其中,卫星的数量为N个,N个卫星属于G个卫星系统,每个卫星均对应一个距离估计值和一个初始伪距观测值,N和G均为正整数;
目标观测值确定单元包括:
卫星分组子单元,用于获取N个卫星分别对应的残差值,根据G个卫星系统,对N个卫星分别对应的残差值进行分组,得到G个第一序列;一个卫星对应的残差值是由该卫星对应的距离估计值和初始伪距观测值所确定的;
排序子单元,用于获取G个第一序列中的第l个第一序列,对第l个第一序列中所包含的残差值进行排序,将排序后的第l个第一序列确定为第二序列,获取第二序列对应的初始掩码数组;l为小于G的正整数;
序列更新子单元,用于根据第二序列对应的中位数与第二序列中所包含的残差值,对第二序列进行更新,得到第三序列;
序列更新子单元,还用于根据第三序列对应的下四分位数,以及第三序列对应的上四分位数,对第三序列进行更新,得到第四序列;
掩码数组更新子单元,用于根据第二序列对应的下四分位数、第二序列对应的上四分位数以及第四序列对应的中位数,对初始掩码数组进行更新,得到目标掩码数组;
筛选子单元,用于根据目标掩码数组中的非零元素,对第二序列中的残差值所对应的初始伪距观测值进行筛选,得到第二序列对应的目标伪距观测值。
其中,序列更新子单元具体用于:
根据第三序列对应的下四分位数和上四分位数,确定第三序列中所包含的元素数值的有效取值范围;
将第三序列中处于有效取值范围的元素数值,确定为第四序列。
其中,目标伪距观测值的数量为u个,u个目标伪距观测值是指移动终端分别与G个卫星系统中的u个卫星之间的距离观测值,u和G均为正整数;
观测矩阵构建模块包括:
观测值分组单元,用于根据G个卫星系统,对u个目标伪距观测值进行分组,得到G个观测值数组;
参考卫星确定单元,用于获取u个卫星分别对应的载噪比,根据载噪比确定G个观测值数组分别对应的参考卫星,将参考卫星对应的目标伪距观测值确定为参考伪距观测值;参考卫星是指每个观测值数组所对应的卫星中具有最大载噪比的卫星;
观测差值确定单元,用于将u个卫星中除参考卫星之外的卫星确定为非参考卫星,根据参考伪距观测值与非参考卫星对应的目标伪距观测值,确定G个观测值数组关联的距离观测差值;一个距离观测差值用于表征一个观测值数组中的一个非参考卫星,与该观测值数组对应的参考卫星之间的观测距离;
估计差值确定单元,用于根据终端位置、光速值、非参考卫星对应的卫星位置和卫星钟差,以及参考卫星对应的卫星位置和卫星钟差,确定G个观测值数组关联的距离估计差值;一个距离估计差值用于表征一个观测值数组中的一个非参考卫星,与该观测值数组对应的参考卫星之间的估计距离;
观测方程构建单元,用于根据距离观测差值和距离估计差值,构建G个观测值数组分别对应的伪距观测方程,将G个观测值数组分别对应的伪距观测方程组合为伪距观测矩阵。
其中,评估模块包括:
观测值误差获取单元,用于若卫星数量u大于数量阈值,则获取u个卫星分别对应的伪距观测值误差;
观测权阵构建单元,用于根据G个观测值数组关联的非参考卫星的伪距观测值误差,以及G个观测值数组关联的参考卫星的伪距观测值误差,生成观测值误差矩阵,获取观测值误差矩阵中的对角元素构建观测权阵;
协方差矩阵确定单元,用于获取伪距观测矩阵针对终端位置的微分矩阵,根据观测权阵和微分矩阵,确定移动终端对应的位置协方差矩阵;
单位权中误差确定单元,用于根据观测权阵、距离观测差值所组成的矩阵以及卫星数量u,确定u个卫星对应的单位权中误差;
精度指标确定单元,用于根据位置协方差矩阵和单位权中误差,确定终端位置对应的位置精度指标。
其中,位置精度指标包括终端水平定位精度和终端高程定位精度;
精度指标确定单元包括:
评估数值获取子单元,用于根据位置协方差矩阵,确定移动终端对应的过渡评估矩阵,获取过渡评估矩阵对应的第一水平评估数值、第二水平评估数值以及高程评估数值;
水平定位精度确定子单元,用于根据第一水平评估数值与第二水平评估数值之和的平方根,以及单位权中误差,确定终端位置对应的终端水平定位精度;
高程定位精度确定子单元,用于根据高程评估数值的平方根与单位权中误差,确定终端位置对应的终端高程定位精度。
其中,该装置还包括:
卫星数量判断模块,用于若卫星数量u小于数量阈值,则根据距离观测差值对终端位置进行评估;
正确结果确定模块,用于当距离观测差值小于或等于误差阈值时,确定终端位置对应的位置精度评估结果为定位正确结果;
错误结果确定模块,用于当距离观测值大于误差阈值时,确定终端位置对应的位置精度评估结果为定位错误结果。
其中,该装置还包括:
差值序列构建模块,用于若卫星数量u等于数量阈值,且伪距观测矩阵中的伪距观测方程所对应的卫星属于相同卫星系统,则根据距离观测差值构建差值序列;
评估结果获取模块,用于获取差值序列对应的方均根,根据方均根对终端位置进行评估,得到终端位置对应的位置精度评估结果。
其中,该装置还包括:
评估结果确定模块,还用于根据位置精度指标确定终端位置对应的位置精度评估结果;
功能服务提供模块,用于若位置精度评估结果为定位正确结果,则根据终端位置为移动终端提供应用功能服务;
终端位置修正模块,用于若位置精度评估结果为定位错误结果,则根据无线局域网络和移动终端中的传感器,对终端位置进行修正。
其中,上述获取卫星导航星历,包括:
请求发送单元,用于向卫星定位基站服务器发送星历获取请求,以使卫星定位基站服务器根据星历获取请求,在卫星星历数据库中获取卫星导航星历;
导航星历接收单元,用于接收卫星定位基站服务器返回的卫星导航星历;卫星导航星历以二进制流的形式进行传输。
本申请实施例一方面提供了一种计算机设备,包括存储器和处理器,存储器与处理器相连,存储器用于存储计算机程序,处理器用于调用计算机程序,以使得该计算机设备执行本申请实施例中上述一方面提供的方法。
本申请实施例一方面提供了一种计算机可读存储介质,计算机可读存储介质中存储有计算机程序,计算机程序适于由处理器加载并执行,以使得具有处理器的计算机设备执行本申请实施例中上述一方面提供的方法。
根据本申请的一个方面,提供了一种计算机程序产品或计算机程序,该计算机程序产品或计算机程序包括计算机指令,该计算机指令存储在计算机可读存储介质中。计算机设备的处理器从计算机可读存储介质读取该计算机指令,处理器执行该计算机指令,使得该计算机设备执行上述一方面提供的方法。
本申请实施例可以通过获取卫星导航星历,以及移动终端与卫星之间的初始伪距观测值,该卫星导航星历包括用于确定卫星位置的参数列表,进而根据卫星导航星历和初始伪距观测值可以确定卫星对应的卫星位置和卫星钟差;根据卫星位置、卫星钟差、移动终端的终端位置和终端钟差,可以对初始伪距观测值进行筛选,得到目标伪距观测值,进而可以根据目标伪距观测值和终端位置,构建移动终端与卫星之间的伪距观测矩阵,根据伪距观测矩阵对终端位置进行评估,可以确定终端位置对应的位置精度指标,该位置精度指标用于表征移动终端的定位精度。可见,通过卫星导航星历和初始伪距观测值计算得到卫星位置和卫星钟差,剔除初始伪距观测值中的观测值粗差得到目标伪距观测值,进而可以根据目标伪距观测值和终端位置构建伪距观测矩阵,基于该伪距观测矩阵可以得到终端位置对应的位置精度指标,该位置精度指标的数值大小可以用于准确评估移动终端的定位质量,进而可以提高移动终端的位置评估准确性。
附图说明
为了更清楚地说明本申请实施例或现有技术中的技术方案,下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍,显而易见地,下面描述中的附图仅仅是本申请的一些实施例,对于本领域普通技术人员来讲,在不付出创造性劳动的前提下,还可以根据这些附图获得其他的附图。
图1是本申请实施例提供的一种网络架构的结构示意图;
图2是本申请实施例提供的一种终端位置精度评估场景示意图;
图3是本申请实施例提供的一种终端定位精度评估方法的流程示意图;
图4是本申请实施例提供的一种移动终端与卫星定位基站服务器之间的通信示意图;
图5是本申请实施例提供的一种伪距观测值的确定示意图;
图6是本申请实施例提供的一种根据卫星导航星历确定卫星位置和卫星钟差的示意图;
图7是本申请实施例提供的一种终端定位精度评估方法的流程示意图;
图8是本申请实施例提供的一种构建伪距观测矩阵的流程示意图;
图9是本申请实施例提供的一种评估移动终端定位精度的流程示意图;
图10是本申请实施例提供的一种终端定位精度评估装置的结构示意图;
图11是本申请实施例提供的一种计算机设备的结构示意图。
具体实施方式
下面将结合本申请实施例中的附图,对本申请实施例中的技术方案进行清楚、完整地描述,显然,所描述的实施例仅仅是本申请一部分实施例,而不是全部的实施例。基于本申请中的实施例,本领域普通技术人员在没有作出创造性劳动前提下所获得的所有其他实施例,都属于本申请保护的范围。
本申请实施例涉及以下几个概念:
位置服务(LBS,Location Based Services):LBS服务是无线运营公司为用户提供的一种与位置有关的服务;LBS服务是利用各类型的定位技术来获取定位设备当前的所在位置,通过移动互联网向定位设备提供信息资源和基础服务。LBS服务中融合了移动通讯、互联网络、空间定位、位置信息、大数据等多种信息技术,利用移动互联网络服务平台进行数据更新和交互,使用户可以通过空间定位来获取相应的服务。
全球卫星导航系统(Global Navigation Satellite System):全球卫星导航系统也称为全球导航卫星系统,可以在地球表面或近地空间的任何地点为用户提供全天候的3维坐标和速度以及时间信息的空基无线电导航定位系统。常见的全球卫星导航系统包括美国全球定位系统(Global Positioning System,GPS)、中国北斗卫星导航系统(BeiDouNavigation Satellite System,BDS)、俄罗斯格洛纳斯卫星导航系统(GLONASS)和欧盟伽利略定位系统(GALILEO)四大卫星导航系统。最早出现的是美国的GPS(GlobalPositioning System),现阶段技术最完善的也是GPS系统。随着近年来BDS、GLONASS系统在亚太地区的全面服务开启,尤其是BDS系统在民用领域发展越来越快。卫星导航系统已经在航空、航海、通信、人员跟踪、消费娱乐、测绘、授时、车辆监控管理和汽车导航与信息服务等方面广泛使用,而且总的发展趋势是为实时应用提供高精度服务。
CORS(Continuously Operating Reference Stations,连续运行参考站,也可以称为卫星定位基站):CORS系统是卫星定位技术、计算机网络技术、数字通讯技术等高新科技多方位、深度结晶的产物。CORS系统由基准站网、数据处理中心、数据传输系统、定位导航数据播发系统、用户应用系统五个部分组成,各基准站与监控分析中心间通过数据传输系统连接成一体,形成专用网络。
请参见图1,图1是本申请实施例提供的一种网络架构的结构示意图。如图1所示,该网络架构可以包括服务器10d和用户终端集群,该用户终端集群可以包括一个或者多个用户终端,这里不对用户终端的数量进行限制。如图1所示,该用户终端集群可以具体包括用户终端10a、用户终端10b以及用户终端10c等。其中,服务器10d可以是独立的物理服务器,也可以是多个物理服务器构成的服务器集群或者分布式系统,还可以是提供云服务、云数据库、云计算、云函数、云存储、网络服务、云通信、中间件服务、域名服务、安全服务、CDN、以及大数据和人工智能平台等基础云计算服务的云服务器。用户终端10a、用户终端10b以及用户终端10c等均可以包括:智能手机、平板电脑、笔记本电脑、掌上电脑、移动互联网设备(mobile internet device,MID)、可穿戴设备(例如智能手表、智能手环等)以及自动驾驶系统中的车载电脑等具有定位功能的移动终端。如图1所示,用户终端10a、用户终端10b以及用户终端10c等可以分别与服务器10d进行网络连接,以便于每个用户终端可以通过该网络连接与服务器10d之间进行数据交互。
以用户终端10a为例,该用户终端10a可以称为移动终端,即可以在移动中使用的计算机设备;用户终端10a可以理解为一个综合信息处理平台,具有非常丰富的通信方式。例如,可以通过全球移动通信系统(Global System for Mobile Communication,GSM)、码分多址(Code Division Multiple Access,CDMA)、宽带码分多址(Wideband CodeDivision Multiple Access,WCDMA)、增强型数据速率GSM演进技术(Enhanced Data Ratefor GSM Evolution,EDGE)、4G(第四代通讯技术)、5G(第五代通讯技术)等无线运营网通讯,也可以通过无线局域网络(Wireless-Fidelity,WiFi),蓝牙和红外线进行通信;另外,用户终端10a集成有全球卫星导航系统定位芯片,该全球卫星导航系统定位芯片可以用于处理卫星信号,并对用户终端10a的使用用户进行精准定位,可以用于位置服务。服务器10d可以是指CORS服务器,用户终端10a根据集成的全球卫星导航系统定位芯片,获取到该用户终端10a的终端位置后,可以对该用户终端10a的终端位置进行定位精度评估。在对终端位置进行定位精度评估的过程中,用户终端10a可以向CORS服务器发送星历获取请求,CORS服务器在接收到用户终端10a发送的星历获取请求后,可以向该用户终端10a播发实时导航星历(也可以称为两行轨道数据);其中,该卫星导航星历是用于描述太空飞行体位置和速度的表达式,即卫星导航星历可以是指卫星运行随时间变化的精确位置或轨迹表,卫星导航星历可以理解为时间的函数。
用户终端10a可以接收CORS服务器播发的卫星导航星历,从各卫星系统的原始数据中获取该用户终端10a与各卫星之间的几何距离测量值(可以称为伪距观测值);用户终端10a可以根据伪距观测值和卫星导航星历计算各个卫星分别对应的卫星位置和卫星钟差,还可以对伪距观测值进行处理,剔除观测值粗差;进而可以根据处理后的伪距观测值和用户终端10a的终端位置构建伪距观测矩阵,基于该伪距观测矩阵对终端位置进行评估,输出终端位置对应的位置精度指标,该位置精度指标可以用于表示用户终端10a所对应的终端位置的定位精度,位置精度指标的数值越大,表示终端位置越不准确;该位置精度指标可以用于辅助无线局域网络、传感器等,对用户终端10a进行融合定位,以提高用户终端10a的定位准确性。
请一并参见图2,图2是本申请实施例提供的一种终端位置精度评估场景示意图。以电子地图应用场景为例,对终端位置精度评估过程进行描述。如图2所示,移动终端20a可以为上述图1的用户终端集群中的任意一个用户终端,该移动终端20a中安装有电子地图应用,当移动终端20a的使用者想要查询某地点时,可以启动移动终端20a中的电子地图应用,移动终端20a可以响应针对电子地图应用的启动操作,在移动终端20a的终端屏幕中显示该电子地图应用的首页,该电子地图应用的首页可以显示移动终端20a所在的终端位置:位置S1,以及位置S1所在的附近地理区域。
移动终端20a在获取到位置S1后,可以向CORS服务器20b(可以是指上述图1所对应实施例中的服务器10d)发送星历获取请求,CORS服务器20b在接收到移动终端20a发送的星历获取请求后,可以从卫星星历数据库20c中读取各卫星系统对应的实时导航星历(也可以称为卫星导航星历),并向移动终端20a播发实时导航星历,移动终端20a可以接收CORS服务器20b播发的实时导航星历。与此同时,移动终端20a还可以向各卫星系统请求卫星发射信号的原始数据,从原始数据中获取移动终端20a分别与各卫星系统中的每个卫星之间的几何距离观测值,本申请实施例中可以将移动终端20a与卫星之间的几何距离观测值称为伪距观测值,此时的伪距观测值是未进行筛选处理的初始值,因此也可以称为初始伪距观测值。如图2所示,假设存在n个卫星(n为正整数),移动终端20a与卫星1之间的初始伪距观测值为:观测值1,移动终端20a与卫星2之间的初始伪距观测值为:观测值2,移动终端20a与卫星3之间的初始伪距观测值为:观测值3,……,移动终端20a与卫星n之间的初始伪距观测值为:观测值n,因此,移动终端20a可以获取n个卫星的初始伪距观测值如序列20d所示。
移动终端20a可以根据序列20d中所包含的初始伪距观测值,以及自身接收到卫星发射信号的接收时间t1,计算出每个卫星发射信号的发射时间;由于星历是指卫星运行随时间变化的精确位置或轨迹表,因此将卫星发射信号的发射时间输入至实时导航星历,可以计算得到n个卫星分别对应的卫星位置,以得到序列20e。其中,接收时间t1可以是指移动终端20a中的电子地图应用的启动时间,卫星1对应的卫星位置可以表示为:位置1,卫星2对应的卫星位置可以表示为:位置2,卫星3对应的卫星位置可以表示为:位置3,……,卫星n对应的卫星位置可以表示为:位置n。
由于序列20d中所包含的初始伪距观测值可能存在粗差,因此移动终端20a可以对序列20d中所包含的初始伪距观测值进行筛选,得到包含目标伪距观测值的序列20f。其中,粗差可以是指粗大误差,比在平常观测条件下所可能出现的最大误差还要大的误差,粗差可以是由测试人员粗心大意或仪器故障所造成的差错;移动终端20a可以采用粗差探测法对初始伪距观测值进行筛选,粗差探测法可以包括但不限于:数据探测法、拟准检定法、抗差估计法。如图2所示,若卫星4对应的初始伪距观测值为观测值粗差,则可以剔除卫星4对应的初始伪距观测值:位置4。
移动终端20a可以根据终端位置S1,以及包含目标伪距观测值的序列20f,构建卫星伪距观测矩阵20g。移动终端20a可以根据位置S1和卫星位置,确定卫星对应的距离估计差值,进而可以根据距离估计差值构建伪距观测矩阵20g。需要说明的是,初始伪距观测值的筛选过程以及伪距观测矩阵的构建过程均需要考虑卫星钟差、终端钟差、真空中的光速值等信息;其中,钟差可以是指时钟读数与真实系统时间之间的差异,如卫星钟差可以是指卫星钟与真实系统时间之间的差异,即各卫星分别对应的钟差;终端钟差可以是指终端钟与真实系统时间之间的差异,即移动终端20a对应的钟差。
进一步地,移动终端20a可以根据伪距观测矩阵20g对位置S1进行精度评估,输出位置S1对应的位置精度指标,该位置精度指标可以包括水平定位精度A1和高程定位精度A2。当水平定位精度A1和高程定位精度A2大于精度阈值(例如,精度阈值为B1)时,可以确定位置S1与移动终端20a的真实位置之间存在较大的差异,即移动终端20a的定位质量较差,可以结合移动终端20a对应的无线局域网络和传感器等,对位置S1进行修正,得到修正后的位置S2。
可选的,当水平定位精度A1和高程定位精度A2小于或等于精度阈值时,可以确定位置S1与移动终端20a的真实位置之间存在非常小的差异,将位置S1作为移动终端20a的真实位置。移动终端20a的使用者在电子地图应用中搜索地点1时,可以将位置S1作为该使用者的当前位置,查询位置S1与所输入地点1之间的路线。
本申请实施例中,根据各卫星所对应的目标伪距观测值,以及移动终端20a对应的位置S1,构建各卫星对应的伪距观测矩阵,根据伪距观测矩阵对位置S1进行精度评估,以确定移动终端的定位质量,可以有效提高移动终端的位置定位精度。
请参见图3,图3是本申请实施例提供的一种终端定位精度评估方法的流程示意图。该终端定位精度评估方法可以由移动终端执行,该移动终端可以为上述图1的用户终端集群中的任一用户终端。如图3所示,该终端定位精度评估方法可以包括以下步骤:
步骤S101,获取卫星导航星历,获取移动终端与卫星之间的初始伪距观测值,根据初始伪距观测值和卫星导航星历,确定卫星对应的卫星位置和卫星钟差。
具体的,移动终端根据所集成的全球卫星定位导航系统定位芯片,获取到自身所在的定位位置(终端位置)后,可以对移动终端对应的终端位置进行定位精度评估,以确定终端位置的准确性。移动终端可以向卫星定位基站服务器发送星历获取请求,卫星定位基站服务器接收到移动终端发送的星历获取请求后,可以从卫星星历数据库中获取实时的卫星导航星历,并将实时的卫星导航星历以二进制流的形式传输至移动终端,此时的移动终端可以接收卫星定位基站服务器播发的卫星导航星历;其中,该卫星定位基站服务器也可以称为CORS服务器,该卫星导航星历可以包括一组用于确定卫星位置的参数列表,卫星星历数据库中可以包括所有卫星的导航星历。
请一并参见图4,图4是本申请实施例提供的一种移动终端与卫星定位基站服务器之间的通信示意图。如图4所示,移动终端30a可以包括终端定位精度评估单元30b和无线局域网络连接模块30e,该终端定位精度评估单元30b可以用于对移动终端的终端位置进行精度评估,该无线局域网络连接模块30e可以用于为移动终端30a连接WiFi/4G/5G网络。移动终端30a可以与CORS服务器30c(即卫星定位基站服务器)进行网络通信连接,移动终端30a中的终端定位精度评估单元30b可以通过无线局域网络连接模块30e所提供的WiFi/4G/5G网络,向CORS服务器30c发送星历获取请求;CORS服务器30c接收到移动终端30a发送的星历获取请求后,可以从卫星星历数据库30d中获取实时导航星历(即上述卫星导航星历);该实时导航星历可以是指与星历获取请求的发送时间相关联的导航星历,如星历获取请求的发送时间为T1,该实时导航星历可以为T0至T2时间范围内的卫星导航星历,T1处于T0至T2时间范围。CORS服务器30c可以通过WiFi/4G/5G网络,将实时导航星历以二进制流的形式传输至终端定位精度评估单元30b,此时的移动终端30b可以接收CORS服务器30c所播发的实时导航星历,即从CORS服务器30c获取卫星导航星历。
与此同时,移动终端还可以获取卫星对应的原始数据,在原始数据中可以提取当前移动终端与卫星之间的几何距离测量值,此时的几何距离测量值可以称为初始伪距观测值,其中,初始伪距观测值可以是指含有时钟误差和大气层折射延迟的距离测量值,而非移动终端与卫星之间的“真实距离”,因此将其成为伪距。换言之,移动终端可以从卫星对应的原始数据中直接提取初始伪距观测值。
请一并参见图5,图5是本申请实施例提供的一种伪距观测值的确定示意图。如图5所示的卫星可以包括北斗卫星、全球定位系统卫星(GPS卫星)以及格洛纳斯卫星(GLONASS卫星);移动终端与卫星之间的初始伪距观测值可以包括:北斗卫星与移动终端之间的初始伪距观测值,GPS卫星与移动终端之间的初始伪距观测值,GLONASS卫星与移动终端之间的初始伪距观测值。其中,北斗卫星与移动终端之间的初始伪距观测值可以表示为北斗卫星与移动终端之间的几何距离测量值,GPS卫星与移动终端之间的初始伪距观测值可以表示为GPS卫星与移动终端之间的几何距离测量值,GLONASS卫星与移动终端之间的初始伪距观测值可以表示为GLONASS卫星与移动终端之间的几何距离测量值。
进一步地,移动终端可以根据初始伪距观测值和卫星导航星历,计算卫星所对应的卫星位置和卫星钟差,即移动终端可以获取自身接收到卫星发射信号的接收时间,根据初始伪距观测值、接收时间以及光速值,确定卫星发射信号的发射时间,进而可以根据发射时间和卫星导航星历中所包含的参数列表,确定卫星对应的卫星位置和卫星钟差。其中,卫星的数量可以为n个(n为正整数),一个卫星可以对应一个初始伪距观测值、一个卫星位置以及一个卫星钟差,每个卫星所对应的卫星位置和卫星钟差可以是不同的。
其中,计算卫星位置和卫星钟差的具体过程可以包括:若移动终端在时刻tRx接收到n个卫星发射的信号,n颗卫星信号发射时刻可以为:
则n个卫星对应的初始伪距观测值可以表示为:
由上述公式(1)可以得出n个卫星的信号发射时刻为:
其中,上述公式(1)和公式(2)中的c可以表示为真空中的光速值,tRx可以表示为移动终端接收到n个卫星发射信号的接收时间,可以表示为卫星i与移动终端之间的初始伪距观测值,可以表示为卫星i的发射信号的发射时间。
其中,ri可以表示为卫星i对应的卫星位置,该卫星位置以三维坐标的形式进行表示,dti可以表示为卫星i对应的卫星钟差。
需要说明的是,本申请实施例中的n个卫星可以属于不同的卫星系统,如n个卫星中可以包括p颗GPS卫星,q颗GLONASS卫星,j颗GALILEO卫星,k颗北斗卫星。请一并参见图6,图6是本申请实施例提供的一种根据卫星导航星历确定卫星位置和卫星钟差的示意图。如图6所示,移动终端40a通过WiFi/4G/5G网络向CORS服务器发送星历获取请求,CORS服务器在接收到移动终端40a所发送的星历获取请求后,可以从卫星星历数据库中读取各个卫星系统分别对应的卫星导航星历。如卫星系统可以包括全球定位系统、北斗卫星定位系统、格洛纳斯卫星导航系统以及伽利略定位系统,则CORS服务器可以在卫星星历数据库中获取北斗星历参数表、全球定位系统星历参数表、格洛纳斯星历参数表以及伽利略星历参数表,并将上述四个卫星系统分别对应的星历参数表,通过WiFi/4G/5G网络以二进制流的形式传输至移动终端40a。移动终端40a可以接收CORS服务器播发的北斗星历参数表、全球定位系统星历参数表、格洛纳斯星历参数表以及伽利略星历参数表,即获取各卫星系统所对应的卫星导航系统。其中,移动终端40a可以包括卫星信息处理单元40b,该卫星信息处理单元40b可以根据CORS服务器所播发的星历参数表,以及卫星发射信号的发射时间,计算各卫星分别对应的卫星位置、卫星钟差等信息。
步骤S102,获取移动终端对应的终端位置和终端钟差,根据卫星位置、卫星钟差、终端位置以及终端钟差,对初始伪距观测值进行筛选,得到目标伪距观测值。
具体的,由于移动终端获取到的初始伪距观测值可能存在粗差,即初始伪距观测值与真实几何距离之间存在很大差异,因此移动终端可以对各卫星分别对应的初始伪距观测值进行筛选,将初始伪距观测值中所包含的观测值粗差进行剔除,得到目标伪距观测值。可以理解的是,目标伪距观测值是指对初始伪距观测值进行筛选后所得到的伪距观测值,即目标伪距观测值的数量小于或等于初始伪距观测值的数量;当目标伪距观测值的数量与初始伪距观测值的数量相等时,表示初始伪距观测值中不存在观测值粗差;目标伪距观测值的数量小于初始伪距观测值的数量时,表示初始伪距观测值中存在观测值粗差。例如,当卫星i与移动终端之间的初始伪距观测值为观测值粗差时,可以将卫星i与移动终端之间的初始伪距观测值进行剔除;当卫星i与移动终端之间的初始伪距观测值不是观测值粗差时,可以将卫星i与移动终端之间的初始伪距观测值确定为目标伪距观测值。
其中,移动终端可以采用粗差探测法对初始伪距观测值进行筛选,即移动终端可以获取自身对应的终端位置和终端钟差,并根据卫星位置、终端位置、卫星钟差、终端钟差,确定移动终端与各卫星之间的距离估计值,进而可以根据初始伪距观测值和距离估计值之间的残差值,对初始伪距观测值进行筛选,得到剔除观测值粗差后所得到的目标伪距观测值。其中,初始伪距观测值可以是指移动终端与卫星之间的测量值,距离估计值可以是指根据终端位置、终端钟差、卫星位置以及卫星终端等信息建立数学模型,所计算得到的移动终端与卫星之间的估计值。移动终端通过对n个卫星分别对应的残差值进行处理,剔除初始伪距观测值中的观测值粗差,得到粗差探测处理后的目标伪距观测值,该目标伪距观测值的数量可以为u个,即经过粗差探测处理后保留了u个卫星对应的目标伪距观测值,u为小于或等于n的正整数。
步骤S103,根据目标伪距观测值和终端位置,构建移动终端与卫星之间的伪距观测矩阵。
具体的,移动终端可以根据u个卫星所属的卫星系统,对u个卫星所对应的目标伪距观测值进行分组,得到G个观测值数组,本申请实施例中的G可以为数值4,每个卫星系统均可以对应一个观测值数组。移动终端可以从u个卫星中选择每个卫星系统所对应的参考卫星,即参考卫星的数量等于观测值数组的数量,参考卫星对应的目标伪距观测值可以称为参考伪距观测值,将u个卫星中除参考卫星之外的其余卫星确定为非参考卫星;移动终端可以根据参考卫星对应的参考伪距观测值、非参考卫星对应的目标伪距观测值、终端位置、卫星位置以及卫星钟差构建每个卫星系统所对应的伪距观测方程,进而可以将所有卫星系统所对应的伪距观测方程进行组合,构成伪距观测矩阵。其中,参考卫星可以是指每个观测值数组所对应的卫星中具有最大载噪比的卫星,载噪比也可以称为载波-噪声功率谱密度比(Carrier-to-noise-density ratio),载噪比可以是指解调前的射频信号载波(Carrier)功率与噪声功率的比值。
举例来说,u取值为12时,此时的u个卫星分别表示为卫星1至卫星12,卫星1对应的目标伪距观测值可以表示为伪距观测值1,卫星2对应的目标伪距观测值可以表示为伪距观测值2,……,卫星12对应的目标伪距观测值可以表示为伪距观测值12。若卫星1、卫星2以及卫星3属于GPS卫星,则可以将伪距观测值1、伪距观测值2以及伪距观测值3构成观测值数组1;若卫星4、卫星5以及卫星6属于GLONASS卫星,则可以将伪距观测值4、伪距观测值5以及伪距观测值6构成观测值数组2;若卫星7、卫星8以及卫星9属于GALILEO卫星,则可以将伪距观测值7、伪距观测值8以及伪距观测值9构成观测值数组3;若卫星10、卫星11以及卫星12属于北斗卫星,则可以将伪距观测值10、伪距观测值11以及伪距观测值12构成观测值数组4。当卫星1为GPS卫星中具有最大载噪比的卫星时,可以将卫星1作为GPS卫星的参考卫星;当卫星4为GLONASS卫星中具有最大载噪比的卫星时,可以将卫星4作为GLONASS卫星的参考卫星;当卫星7为GALILEO卫星中具有最大载噪比的卫星时,可以将卫星7作为GALILEO卫星的参考卫星;当卫星10为北斗卫星中具有最大载噪比的卫星时,可以将卫星10作为北斗卫星的参考卫星。
移动终端可以根据观测值数组1中伪距观测值2与伪距观测值1之间的差值、移动终端、卫星1的卫星位置和卫星钟差、卫星2的卫星位置和卫星钟差,构建卫星2的伪距观测方程;同理,可以构建卫星3的伪距观测方程,进而得到GPS卫星对应的伪距观测方程。以此类推,可以构建GLONASS卫星对应的伪距观测方程、GALILEO卫星对应的伪距观测方程以及北斗卫星对应的伪距观测方程,进而可以将所有伪距观测方程组成伪距观测矩阵。
步骤S104,根据伪距观测矩阵对终端位置进行评估,确定终端位置对应的位置精度指标;该位置精度指标用于表征移动终端的定位精度。
具体的,移动终端可以根据伪距观测矩阵中所包含的伪距观测方程的数量,确定终端位置的精度评估方式,根据精度评估方式对终端位置进行精度评估,可以确定终端位置对应的位置精度指标,该位置精度指标可以用于表征移动终端的位置定位精度;位置精度指标的数值越大,表示移动终端的定位质量越差,终端位置越不准确。
可选的,当伪距观测方程的数量小于数量阈值(数量阈值的具体数值可以进行认为设定,本申请实施例中,可以将数量阈值设置为3)时,可以直接根据参考卫星所对应的参考伪距观测值与非参考卫星所对应的目标伪距观测值之间的差值,对终端位置进行精度评估,当差值大于预先设定的误差阈值时,可以确定终端位置对应的位置精度评估结果为定位正确结果,即终端位置定位正确;当差值大于预先设定的误差阈值时,可以确定终端位置对应的位置精度评估结果为定位错误结果,即终端位置定位错误。
可选的,当伪距观测方程的数量等于数量阈值时,可以根据参考卫星所对应的参考伪距观测值与非参考卫星所对应的目标伪距观测值之间的差值,构建差值序列,并计算差值序列的方均根(Root Meam Rquare,RMS),根据方均根的数值对终端位置进行精度评估,当方均根小于或等于预先设定的阈值时,可以确定终端位置对应的位置精度评估结果为定位正确结果,即终端位置定位正确;当方均根大于预先设定的误差阈值时,可以确定终端位置对应的位置精度评估结果为定位错误结果,即终端位置定位错误。需要说明的是,不同精度评估方式中所设置的误差阈值可以存在差异,该误差阈值均可以根据具体需求来进行人为设定。
可选的,当伪距观测方程的数量大于数量阈值时,此时的位置精度指标可以包括终端水平定位精度和终端高程定位精度,即移动终端可以计算单位权中误差和位置协方差矩阵,并根据单位权中误差和位置协方差矩阵,确定终端位置对应的终端水平定位精度和终端高程定位精度。其中,单位权中误差和位置协方差矩阵可以是基于伪距观测矩阵所确定的。
可选的,移动终端可以根据位置精度指标确定终端位置对应的位置精度评估结果;若位置精度评估结果为定位正确结果,则根据终端位置为移动终端提供应用功能服务;若位置精度评估结果为定位错误结果,则根据无线局域网络和移动终端中的传感器,对终端位置进行修正。例如,在智能手机电子地图应用场景中,当移动终端对应的位置精度评估结果为定位正确结果时,可以基于终端位置为移动终端提供打车、查询地点、查询公交、查询地铁、路线导航等功能服务;当移动终端对应的位置精度评估结果为定位错误结果时,可以基于无线局域网络(WiFi)和传感器进行融合定位,对终端位置进行修正,进而可以提高移动终端的位置评估准确性。当然,对于修正后的终端位置,同样需要对其进行精度评估,在修正后的终端位置所对应的位置精度评估结果为定位正确结果时,基于修正后的终端位置重新为移动终端提供打车、查询地点、查询公交、查询地铁、路线导航等功能服务。
可选的,在自动驾驶场景中,此时的移动终端可以为自动驾驶系统中的车载电脑,当移动终端对应的位置精度评估结果为定位正确结果时,可以基于终端位置实时记录自动驾驶车辆的行驶轨迹,并基于终端位置对自动驾驶车辆进行位置跟踪;当移动终端对应的位置精度评估结果为定位错误结果时,可以基于无线局域网络和传感器进行融合定位,对终端位置进行修正,得到修正后的终端位置,根据修正后的终端位置实时记录自动驾驶车辆的行驶轨迹,并基于终端位置对自动驾驶车辆进行位置跟踪。可选的,在车机导航系统中,可以根据位置精度评估结果,以获取移动终端较为准确的位置信息,基于该位置信息可以为用户提供更为准确的路线导航,进而可以提高导航准确率。
本申请实施例中,移动终端可以通过卫星导航星历和初始伪距观测值计算得到卫星位置和卫星钟差,剔除初始伪距观测值中的观测值粗差得到目标伪距观测值,进而可以根据目标伪距观测值和终端位置构建伪距观测矩阵,基于该伪距观测矩阵可以得到终端位置对应的位置精度指标,该位置精度指标的数值大小可以用于准确评估移动终端的定位质量,进而可以提高移动终端的位置评估准确性;根据位置精度指标,结合无线局域网络和传感器对移动终端进行融合定位,并基于定位位置为移动终端提供应用功能服务,可以增强应用功能服务的推荐准确性;本申请实施例提供的方法可以适用于智能手机、车机导航系统等场景中,适用性高。
请参见图7,图7是本申请实施例提供的一种终端定位精度评估方法的流程示意图。该终端定位精度评估方法可以由移动终端执行,该移动终端可以为上述图1的用户终端集群中的任一用户终端。如图7所示,该终端定位精度评估方法可以包括以下步骤:
步骤S201,发送星历请求和获取实时卫星导航星历,获取初始伪距观测值。
步骤S202,根据初始伪距观测值计算卫星信号发射时刻。
步骤S203,计算卫星位置和钟差。
其中,步骤S201-步骤S203的具体实现方式可以参见上述图3所对应实施例中对步骤S101的描述,这里不再进行赘述。
步骤S204,处理卫星初始伪距观测值。
具体的,移动终端可以采用粗差探测法对卫星的初始伪距观测值进行粗差探测,剔除初始伪距观测值中所包含的观测值粗差,即移动终端可以根据卫星位置、卫星钟差、终端位置以及终端钟差,确定移动终端与卫星之间的距离估计值,进而可以根据距离估计值和初始伪距观测值之间的残差值,对初始伪距观测值进行筛选,得到目标伪距观测值。
其中,卫星的数量可以为N个,N个卫星属于G个卫星系统,每个卫星均对应一个距离估计值和一个初始伪距观测值,N和G均为正整数,本申请实施例中,G个卫星系统可以包括:GPS卫星系统、GLONASS卫星系统、GALILEO卫星系统、北斗卫星系统,即此时的G可以取值为4。移动终端可以获取N个卫星分别对应的残差值,根据G个卫星系统对N个卫星分别对应的残差值进行分组,得到G个第一序列,其中,一个卫星对应的残差值可以是由该卫星对应的距离估计值和初始伪距观测值所确定的。其中,初始伪距观测值的数学模型可以表示为:
其中,上述公式(4)中的可以表示为卫星i的初始伪距观测值;ru可以表示为移动终端的真实位置,该真实位置ru可以表示为三维坐标(x u ,y u ,z u );ri可以表示为卫星i的卫星位置,该卫星位置ri可以表示为三维坐标(x i ,y i ,z i );dtr可以表示为移动终端的终端钟差;dti可以表示为卫星i的卫星钟差,c为真空中的光速值,ζ为误差改正数(该误差改正数可以包括电离层、对流层以及地球自转改正,ζ的取值可以由经验模型计算得到);||ru-ri||可以用于表示真实位置ru与卫星位置ri之间的欧几里得距离。需要说明的是,移动终端所获得的终端位置可以是指还未进行精度评估的位置信息,该公式(4)中的真实位置ru对于移动终端本身而言是未知的。
其中,可以表示为通过上述公式(4)所示的数学模型,计算得到移动终端与卫星i之间的几何距离;可以表示为移动终端与卫星i之间的残差值。移动终端可以根据上述公式(5),计算N个卫星分别与移动终端之间的残差值,即可以得到N个残差值。
若N个卫星包括p颗GPS卫星,q颗GLONASS卫星,j颗GALILEO卫星,k颗北斗卫星,则可以将N个卫星分成四组,即将N个残差值划分为GPS、GLONASS、GALILEO、北斗四个序列,此时的四个序列可以称为G个第一序列。其中,每个第一序列均可以进行观测值粗差探测处理,下面以GPS序列为例,对观测值粗差探测处理过程进行具体描述。
移动终端可以从G个第一序列中获取第l个第一序列,此时的第l个第一序列可以为GPS序列,并对第l个第一序列中所包含的残差值进行排序,将排序后的第l个第一序列确定为第二序列,获取第二序列对应的初始掩码数组,l为小于G的正整数。其中,第二序列可以理解为对GPS序列中所包含的残差值进行排序后的序列,此时的第二序列可以表示为:
mask[p]={0,0,……,0}(7)
进一步地,移动终端可以根据第二序列对应的中位数与第二序列中所包含的残差值,对第二序列进行更新,得到第三序列;进而可以根据第三序列对应的下四分位数和第三序列对应的上四分位数,对第三序列进行更新,得到第四序列;移动终端可以根据第二序列对应的下四分位数、第二序列对应的上四分位数以及第四序列对应的中位数,对初始掩码数组进行更新,得到目标掩码数组;根据目标掩码数组中的非零元素,对第二序列中的残差值所对应的初始伪距观测值进行筛选,得到第二序列对应的目标伪距观测值。其中,移动终端可以分别计算第二序列的下四分位数、中位数以及下四分位数,下四分位数可以表示为第二序列中处于25%位置处的残差值,中位数可以表示为第二序列中处于50%位置处的残差值,上四分位数可以表示为第二序列中处于75%位置处的残差值。
其中,上述公式(8)中的Int表示为取整运算,Dec表示为取小数部分。
移动终端可以根据第二序列的上四分位数和下四分位数,对第二序列进行粗差探测,即对第二序列中的每个元素(残差值)均做下述处理:当或时,mask[i]=mask[i]+1,否则不作处理,其中,i为小于q的正整数,此次操作是对初始掩码数组的更新处理过程。
进一步地,移动终端还可以计算第三序列分别对应的下四分位数和下四分位数,其计算过程与上述第二序列的下四分位数和下四分位数的计算过程类似,这里不再进行赘述。第三序列的下四分位数可以表示为,第三序列的上四分位数可以表示为。
移动终端在得到第四序列后,可以计算第四序列的中位数,该中位数的计算过程与上述第二序列的中位数M2/4的计算过程类似;随后可以作如下处理:当时,mask[i]=mask[i]+1,否则不作处理,此次操作是对上述更新后的掩码数组再次进行更新处理的过程,经过此次操作后,可以得到目标掩码数组。
移动终端可以根据目标掩码数组mask[p]剔除伪距观测值粗差,若mask[i]≥1,则剔除初始伪距观测值;若mask[i]=0,则保留初始伪距观测值,可以得到GPS序列对应的目标伪距观测值。根据上述一系列操作过程,可以分别获取GLONASS序列、GALILEO序列以及北斗序列所对应的目标伪距观测值。
步骤S205,利用卫星目标伪距观测值和终端位置构建伪距观测矩阵。
具体的,经过筛选后的目标伪距观测值可以表示为:,其中,1 ≤ai≤N,i=1,2,……,u,即目标伪距观测值的数量为u个,u为小于N的正整数,在上述筛选过程中,共剔除了N-u个初始伪距观测值。移动终端可以根据G个卫星系统(可以包括GPS卫星系统、GLONASS卫星系统、GALILEO卫星系统以及北斗卫星系统),对u个目标伪距观测值进行分组,得到G个观测值数组,该G个观测值数组可以分别表示为:
其中,上述公式(12)中,可以表示为GPS卫星系统对应的观测值数组,可以表示为GLONASS卫星系统对应的观测值数组,可以表示为GALILEO卫星系统对应的观测值数组,可以表示为北斗卫星系统对应的观测值数组,p+q+j+k=N,b+z+x+f=u。
进一步地,移动终端可以获取u个卫星分别对应的载噪比,根据载噪比确定G个观测值数组分别对应的参考卫星,将参考卫星对应的目标伪距观测值确定为参考伪距观测值,该参考卫星是指每个观测值数组所对应的卫星中具有最大载噪比的卫星。换言之,移动终端可以在各个卫星系统中分别选择出载噪比最大的卫星作为参考卫星,可以将GPS卫星系统中的参考卫星确定为p1,将GLONASS卫星系统中的参考卫星确定为q1,将GALILEO卫星系统中的参考卫星确定为j1,将北斗卫星系统中的参考卫星确定为k1;参考卫星p1对应的目标伪距观测值为,参考卫星q1对应的目标伪距观测值为,参考卫星j1对应的目标伪距观测值为,参考卫星k1对应的目标伪距观测值为。
进一步地,移动终端可以将u个卫星中除参考卫星之外的卫星确定为非参考卫星,根据参考伪距观测值与非参考卫星对应的目标伪距观测值,确定G个观测值数组关联的距离观测差值,其中,一个距离观测差值用于表征一个观测值数组中的一个非参考卫星,与该观测值数组对应的参考卫星之间的观测距离;移动终端可以根据终端位置、光速值、非参考卫星对应的卫星位置和卫星钟差,以及参考卫星对应的卫星位置和卫星钟差,确定G个观测值数组关联的距离估计差值,其中,一个距离估计差值用于表征一个观测值数组中的一个非参考卫星,与该观测值数组对应的参考卫星之间的估计距离;进而可以根据距离观测差值和距离估计差值,构建G个观测值数组分别对应的伪距观测方程,将G个观测值数组分别对应的伪距观测方程组合为伪距观测矩阵。
其中,GPS卫星系统对应的伪距观测方程可以表示为:
其中,上述公式(13)中,可以表示为GPS卫星系统中的卫星pi与参考卫星p1之间的距离观测差值,可以表示为GPS卫星系统中的卫星pi与参考卫星p1之间的距离估计差值;表示终端位置,rpi表示卫星pi对应的卫星位置,rp1表示参考卫星p1对应的卫星位置;dtp1表示参考卫星p1对应的卫星钟差,dtpi表示卫星pi对应的卫星钟差;ζpi表示为卫星pi的误差改正数,ζp1表示为参考卫星p1的误差改正数,其中,c·(dtp1-dtpi)和ζpi-ζp1均为已知值。GPS卫星系统对应的伪距观测方程的数量可以表示为b-1。
同理,GLONASS卫星系统对应的伪距观测方程可以表示为:
其中,上述公式(14)中,可以表示为GLONASS卫星系统中的卫星qi与参考卫星q1之间的距离观测差值,可以表示为GLONASS卫星系统中的卫星qi与参考卫星q1之间的距离估计差值;rqi表示卫星qi对应的卫星位置,rq1表示参考卫星q1对应的卫星位置;dtq1表示参考卫星q1对应的卫星钟差,dtqi表示卫星qi对应的卫星钟差;ζqi表示为卫星qi对应的误差改正数,ζq1表示为参考卫星q1的误差改正数,其中,c·(dtq1-dtqi)和ζqi-ζq1均为已知值。GLONASS卫星系统对应的伪距观测方程的数量可以表示为z-1。
同理,GALILEO卫星系统对应的伪距观测方程可以表示为:
其中,上述公式(15)中,可以表示为GALILEO卫星系统中的卫星ji与参考卫星j1之间的距离观测差值,可以表示为GALILEO卫星系统中的卫星ji与参考卫星j1之间的距离估计差值;rji表示卫星ji对应的卫星位置,rj1表示参考卫星j1对应的卫星位置;dtj1表示参考卫星j1对应的卫星钟差,dtji表示卫星ji对应的卫星钟差;ζji表示为卫星ji的误差改正数,ζj1表示为参考卫星j1的误差改正数,其中,c·(dtj1-dtji)和ζji-ζj1均为已知值。GALILEO卫星系统对应的伪距观测方程的数量可以表示为x-1。
同理,北斗卫星系统对应的伪距观测方程可以表示为:
其中,上述公式(15)中,可以表示为北斗卫星系统中的卫星ki与参考卫星k1之间的距离观测差值,可以表示为北斗卫星系统中的卫星ki与参考卫星k1之间的距离估计差值;rki表示卫星ki对应的卫星位置,rk1表示参考卫星k1对应的卫星位置;dtk1表示参考卫星k1对应的卫星钟差,dtki表示卫星ki对应的卫星钟差;ζki表示为卫星ki的误差改正数,ζki表示为参考卫星k1的误差改正数,其中,c·(dtk1-dtki)和ζki-ζk1均为已知值。北斗卫星系统对应的伪距观测方程的数量可以表示为f-1。
进一步地,将上述公式(13)-公式(16)所表示的伪距观测方程进行组合,可以得到伪距观测矩阵O,该伪距观测矩阵O可以表示为:
步骤S206,评估移动终端位置精度。
具体的,移动终端在构建了伪距观测矩阵后,可以根据该伪距观测矩阵对移动终端对应的终端位置进行精度评估,得到终端位置对应的位置精度评估结果。移动终端可以卫星数量u确定终端位置所采用的精度评估方式。
其中,若卫星数量u小于数量阈值(该数量阈值可以预先进行人为设定,例如,该数量阈值可以设置为3),则移动终端可以直接根据距离观测差值对终端位置进行评估;当距离观测差值小于或等于误差阈值(可以根据实际需求预先进行人为设定)时,确定终端位置对应的位置精度评估结果为定位正确结果;当距离观测值大于误差阈值时,确定终端位置对应的位置精度评估结果为定位错误结果。
下面均以数量阈值取值为3进行说明,当卫星数量u小于3时,此时的卫星数量可以为1个或者2个,当卫星数量u为1时,移动终端无法评估终端位置的定位精度;当卫星数量u为2,且2个卫星属于同一个卫星系统时,移动终端可以直接根据、、中的其中一个值(即距离观测差值)判断移动终端的定位质量,即根据距离观测阈值来评估移动终端的定位准确性;当卫星数量为2,但2个卫星不属于同一个卫星系统时,移动终端无法评估终端位置的定位精度。若距离观测差值大于误差阈值,则表示移动终端的定位质量较差,该移动终端的终端位置所对应的位置精度评估结果为定位错误结果;若距离观测差值小于或等于误差阈值,则表示移动终端的定位质量良好,该移动终端的终端位置所对应的位置精度评估结果为定位正确结果。
可选的,若卫星数量u等于数量阈值,且伪距观测矩阵中的伪距观测方程所对应的卫星属于相同卫星系统,则移动终端可以根据距离观测差值构建差值序列,并获取差值序列对应的方均根,根据方均根对终端位置进行评估,得到终端位置对应的位置精度评估结果。当卫星数量u等于3,且u个卫星均属于同一卫星系统时,假设u个卫星均属于GPS卫星系统,则可以构建差值序列V,该差值序列V可以表示为:
移动终端可以计算差值序列V的方均根(RMS),RMS的计算方式如下述公式(19)所示:
根据RMS与预先设定的误差阈值之间的大小关系,确定终端位置对应的位置精度评估结果。若卫星数量u等于3,但是这3个卫星不属于同一个卫星系统,则移动终端无法评估移动终端的定位精度。
可选的,若卫星数量u大于数量阈值,则移动终端可以获取u个卫星分别对应的伪距观测值误差(也可以称为伪距观测值误差模型),该伪距观测值误差可以表示为:
其中,上述公式(20)中的CN0可以表示为卫星ai的载噪比。
移动终端可以根据G个观测值数组关联的非参考卫星的伪距观测值误差,以及G个观测值数组关联的参考卫星的伪距观测值误差,生成观测值误差矩阵S,获取观测值误差矩阵S中的对角元素构建观测权阵W。其中,观测值误差矩阵S如下述公式(21)所示:
其中,diag(S)可以表示为将观测值误差矩阵S转换为对角矩阵,此时的对角矩阵即为观测权阵W。
进一步地,移动终端可以获取伪距观测矩阵针对终端位置的微分矩阵H,并根据观测权阵W和微分矩阵H,确定移动终端对应的位置协方差矩阵P。其中,微分矩阵H是指伪距观测矩阵O关于终端位置的微分矩阵,如下述公式(22)所示:
P=(HTW-1H)-1(23)
其中,HT可以是指微分矩阵H的转置矩阵,W-1可以是指观测权阵W的逆矩阵。
移动终端可以根据观测权阵W、距离观测差值所组成的矩阵以及卫星数量u,确定u个卫星对应的单位权中误差σ0,进而可以根据位置协方差矩阵P和单位权中误差σ0,确定终端位置对应的位置精度指标。其中,该单位权中误差σ0如下述公式(24)所示:
其中,上述位置精度指标可以包括终端水平定位精度和终端高程定位精度;移动终端可以根据位置协方差矩阵P,确定移动终端对应的过渡评估矩阵Pu,进而可以获取过渡评估矩阵Pu对应的第一水平评估数值Pu(0,0)、第二水平评估数值Pu(1,1)以及高程评估数值Pu(2,2);根据第一水平评估数值Pu(0,0)与第二水平评估数值Pu(1,1)之和的平方根,以及单位权中误差σ0,可以确定终端位置对应的终端水平定位精度;根据高程评估数值的平方根与单位权中误差σ0,确定终端位置对应的终端高程定位精度。其中,过渡评估矩阵Pu可以如下述公式(25)所示:
步骤S207,输出终端位置精度指标。
具体的,在计算得到上述终端水平定位精度和终端高程定位精度后,可以在移动终端的终端屏幕中输出终端水平定位精度和终端高程定位精度,根据终端水平定位精度和终端高程定位精度来确定终端位置的定位质量,可以提高终端位置的精度评估准确性。
请一并参见图8,图8是本申请实施例提供的一种构建伪距观测矩阵的流程示意图。如图8所示,移动终端可以通过下述步骤S301-步骤S306,实现伪距观测矩阵的构建过程。
步骤S301,移动终端在筛选出目标伪距观测值之后,目标伪距观测值可以表示为:,其中,1≤ai≤N,i=1,2,……,u,进而可以通过上述公式(20)为每个目标伪距观测值设置伪距观测值误差模型(也可以称为伪距观测值误差)。
步骤S302,移动终端可以将u个目标伪距观测值(一个卫星对应一个目标伪距观测值)分为GPS卫星系统、GLONASS卫星系统、GALILEO卫星系统以及北斗卫星系统分别对应的观测值数组,如上述公式(12)所示;进而可以获取u个卫星分别对应的载噪比,从上述四个卫星系统中分别选择载噪比最大的卫星作为参考卫星,如GPS卫星系统、GLONASS卫星系统、GALILEO卫星系统以及北斗卫星系统中的参考卫星可以依次表示为p1、q1、j1、k1。
步骤S303,移动终端可以根据参考卫星,分别构建GPS伪距观测方程、GLONASS伪距观测方程、GALILEO伪距观测方程以及北斗伪距观测方程,具体构建方式可以参见上述公式(13)-公式(16),这里不再进行赘述。
步骤S304,移动终端可以将公式(13)-公式(16)所表示的四个卫星系统对应的伪距观测方程进行组合,构成如公式(17)所示的伪距观测方程。
步骤S305,移动终端可以根据每个目标伪距观测值分别对应的伪距观测值误差模型,生成观测值误差矩阵,进而将观测值误差矩阵转换为对角矩阵,即获取观测值误差矩阵中的对角元素,构成观测权阵,其中,观测权阵的计算过程如上述公式(21)所示。
步骤S306,移动终端可以根据上述步骤S304所构建的伪距观测矩阵,计算该伪距观测矩阵关于终端位置的微分矩阵,微分矩阵的计算过程如上述公式(22)所示。
需要说明的是,在构建伪距观测矩阵的过程中无需使用到观测权阵和微分矩阵,该观测矩阵和微分矩阵是用于评估移动终端的终端位置的。
请一并参见图9,图9是本申请实施例提供的一种评估移动终端定位精度的流程示意图。如图9所示,移动终端定位精度的评估过程可以通过步骤S401-步骤S406来实现。
步骤S401,目标伪距观测值的数量为u时,表示u个目标伪距观测值可以对应u个卫星,移动终端可以对卫星数量u进行判断,若卫星数量u小于3(此时的3可以是指预先设定的数量阈值),则执行步骤S402,移动终端无法评估终端位置的定位精度。若卫星数量u等于3,且3个卫星均属于同一个卫星系统,则执行步骤S403,通过计算方均根(RMS)值,并根据RMS值与预先设定的误差阈值,确定终端位置的准确性;若卫星数量u等于3,但3个卫星不属于同一个卫星系统,则执行步骤S402。若卫星数量大于3,则执行步骤S404-步骤S406,移动终端可以根据上述图8所对应实施例中确定的观测权阵、距离观测差值所组成的矩阵以及卫星数量u,计算单位权中误差,如上述公式(24)所示;进而可以根据上述图8所对应实施例中确定的观测矩阵和微分矩阵,计算终端位置对应的位置协方差矩阵,如上述(23)所示;根据单位权中误差和位置协方差矩阵,可以计算得到终端位置对应的终端水平定位精度和终端高程定位精度。其中,终端位置的精度评估过程可以参见上述步骤S206中的描述,这里不再进行赘述。可选的,若卫星数量u小于3,则此时的卫星数量u可以为1或2,当卫星数量u为2,且2个卫星属于同一个卫星系统时,移动终端可以根据非参考卫星的目标伪距观测值和参考卫星的参考伪距观测值之间的距离观测差值,确定移动终端的定位质量;当然,当卫星数量为1,或者卫星数量为2但是2个卫星不属于同一个卫星系统时,移动终端是无法评估终端定位精度的。
本申请实施例中,移动终端可以通过卫星导航星历和初始伪距观测值计算得到卫星位置和卫星钟差,剔除初始伪距观测值中的观测值粗差得到目标伪距观测值,进而可以根据目标伪距观测值和终端位置构建伪距观测矩阵,基于该伪距观测矩阵可以得到终端位置对应的位置精度指标,该位置精度指标的数值大小可以用于准确评估移动终端的定位质量,进而可以提高移动终端的位置评估准确性;根据位置精度指标,结合无线局域网络和传感器对移动终端进行融合定位,并基于定位位置为移动终端提供应用功能服务,可以增强应用功能服务的推荐准确性;基于卫星伪距观测距离,实现对移动终端的定位精度评估过程,计算过程简单,可以提高移动终端的定位精度评估效率。
请参见图10,图10是本申请实施例提供的一种终端定位精度评估装置的结构示意图。如图10所示,该终端定位精度评估装置1可以包括:位置确定模块11,观测值获取模块12,观测矩阵构建模块13,精度评估模块14;
位置确定模块11,用于获取卫星导航星历,获取移动终端与卫星之间的初始伪距观测值,根据初始伪距观测值和卫星导航星历,确定卫星对应的卫星位置和卫星钟差;
观测值获取模块12,用于获取移动终端对应的终端位置和终端钟差,根据卫星位置、卫星钟差、终端位置以及终端钟差,对初始伪距观测值进行筛选,得到目标伪距观测值;
观测矩阵构建模块13,用于根据目标伪距观测值和终端位置,构建移动终端与卫星之间的伪距观测矩阵;
精度评估模块14,用于根据伪距观测矩阵对终端位置进行评估,确定终端位置对应的位置精度指标;该位置精度指标用于表征移动终端的定位精度。
其中,位置确定模块11,观测值获取模块12,观测矩阵构建模块13,精度评估模块14的具体功能实现方式可以参见上述图3所对应实施例中的步骤S101-步骤S104,这里不再进行赘述。
在一些可行的实施方式中,位置确定模块11可以包括:请求发送单元111,导航星历接收单元112,发射时间确定单元113,卫星位置确定单元114;
请求发送单元111,用于向卫星定位基站服务器发送星历获取请求,以使卫星定位基站服务器根据星历获取请求,在卫星星历数据库中获取卫星导航星历;
导航星历接收单元112,用于接收卫星定位基站服务器返回的卫星导航星历;卫星导航星历以二进制流的形式进行传输。
发射时间确定单元113,用于获取移动终端接收到卫星发射信号的接收时间,根据初始伪距观测值、接收时间以及光速值,确定卫星发射信号的发射时间;
卫星位置确定单元114,用于根据发射时间以及卫星导航星历中所包含的参数列表,确定卫星对应的卫星位置和卫星钟差。
其中,请求发送单元111,导航星历接收单元112,发射时间确定单元113,卫星位置确定单元114的具体功能实现方式可以参见上述图3所对应实施例中的步骤S101,这里不再进行赘述。
在一些可行的实施方式中,观测值获取模块12可以包括:距离估计值确定单元121,目标观测值确定单元122;
距离估计值确定单元121,用于根据卫星位置、卫星钟差、终端位置以及终端钟差,确定移动终端与卫星之间的距离估计值;
目标观测值确定单元122,用于根据距离估计值和初始伪距观测值之间的残差值,对初始伪距观测值进行筛选,得到目标伪距观测值。
其中,距离估计值确定单元121,目标观测值确定单元122的具体功能实现方式可以参见上述图3所对应实施例中的步骤S102,这里不再进行赘述。
在一些可行的实施方式中,卫星的数量为N个,N个卫星属于G个卫星系统,每个卫星均对应一个距离估计值和一个初始伪距观测值,N和G均为正整数;
目标观测值确定单元122可以包括:卫星分组子单元1221,排序子单元1222,序列更新子单元1223,掩码数组更新子单元1224,筛选子单元1225;
卫星分组子单元1221,用于获取N个卫星分别对应的残差值,根据G个卫星系统,对N个卫星分别对应的残差值进行分组,得到G个第一序列;一个卫星对应的残差值是由该卫星对应的距离估计值和初始伪距观测值所确定的;
排序子单元1222,用于获取G个第一序列中的第l个第一序列,对第l个第一序列中所包含的残差值进行排序,将排序后的第l个第一序列确定为第二序列,获取第二序列对应的初始掩码数组;l为小于G的正整数;
序列更新子单元1223,用于根据第二序列对应的中位数与第二序列中所包含的残差值,对第二序列进行更新,得到第三序列;
序列更新子单元1223,还用于根据第三序列对应的下四分位数,以及第三序列对应的上四分位数,对第三序列进行更新,得到第四序列;
掩码数组更新子单元1224,用于根据第二序列对应的下四分位数、第二序列对应的上四分位数以及第四序列对应的中位数,对初始掩码数组进行更新,得到目标掩码数组;
筛选子单元1225,用于根据目标掩码数组中的非零元素,对第二序列中的残差值所对应的初始伪距观测值进行筛选,得到第二序列对应的目标伪距观测值。
其中,序列更新子单元1223具体用于:
根据第三序列对应的下四分位数和上四分位数,确定第三序列中所包含的元素数值的有效取值范围;
将第三序列中处于有效取值范围的元素数值,确定为第四序列。
其中,卫星分组子单元1221,排序子单元1222,序列更新子单元1223,掩码数组更新子单元1224,筛选子单元1225的具体功能实现方式可以参见上述图7所对应实施例中的步骤S204,这里不再进行赘述。
在一些可行的实施方式中,目标伪距观测值的数量为u个,u个目标伪距观测值是指移动终端分别与G个卫星系统中的u个卫星之间的距离观测值,u和G均为正整数;
观测矩阵构建模块13可以包括:观测值分组单元131,参考卫星确定单元132,观测差值确定单元133,估计差值确定单元134,观测方程构建单元135;
观测值分组单元131,用于根据G个卫星系统,对u个目标伪距观测值进行分组,得到G个观测值数组;
参考卫星确定单元132,用于获取u个卫星分别对应的载噪比,根据载噪比确定G个观测值数组分别对应的参考卫星,将参考卫星对应的目标伪距观测值确定为参考伪距观测值;参考卫星是指每个观测值数组所对应的卫星中具有最大载噪比的卫星;
观测差值确定单元133,用于将u个卫星中除参考卫星之外的卫星确定为非参考卫星,根据参考伪距观测值与非参考卫星对应的目标伪距观测值,确定G个观测值数组关联的距离观测差值;一个距离观测差值用于表征一个观测值数组中的一个非参考卫星,与该观测值数组对应的参考卫星之间的观测距离;
估计差值确定单元134,用于根据终端位置、光速值、非参考卫星对应的卫星位置和卫星钟差,以及参考卫星对应的卫星位置和卫星钟差,确定G个观测值数组关联的距离估计差值;一个距离估计差值用于表征一个观测值数组中的一个非参考卫星,与该观测值数组对应的参考卫星之间的估计距离;
观测方程构建单元135,用于根据距离观测差值和距离估计差值,构建G个观测值数组分别对应的伪距观测方程,将G个观测值数组分别对应的伪距观测方程组合为伪距观测矩阵。
其中,观测值分组单元131,参考卫星确定单元132,观测差值确定单元133,估计差值确定单元134,观测方程构建单元135的具体功能实现方式可以参见上述图7所对应实施例中的步骤S205,这里不再进行赘述。
在一些可行的实施方式中,评估模块14可以包括:观测值误差获取单元142,观测权阵构建单元143,协方差矩阵确定单元144,单位权中误差确定单元145,精度指标确定单元146;
观测值误差获取单元142,用于若卫星数量u大于数量阈值,则获取u个卫星分别对应的伪距观测值误差;
观测权阵构建单元143,用于根据G个观测值数组关联的非参考卫星的伪距观测值误差,以及G个观测值数组关联的参考卫星的伪距观测值误差,生成观测值误差矩阵,获取观测值误差矩阵中的对角元素构建观测权阵;
协方差矩阵确定单元144,用于获取伪距观测矩阵针对终端位置的微分矩阵,根据观测权阵和微分矩阵,确定移动终端对应的位置协方差矩阵;
单位权中误差确定单元145,用于根据观测权阵、距离观测差值所组成的矩阵以及卫星数量u,确定u个卫星对应的单位权中误差;
精度指标确定单元146,用于根据位置协方差矩阵和单位权中误差,确定终端位置对应的位置精度指标。
其中,观测值误差获取单元142,观测权阵构建单元143,协方差矩阵确定单元144,单位权中误差确定单元145,精度指标确定单元146的具体功能实现方式可以上述图7所对应实施例中的步骤S206,这里不再进行赘述。
在一些可行的实施方式中,位置精度指标包括终端水平定位精度和终端高程定位精度;
精度指标确定单元146可以包括:评估数值获取子单元1461,水平定位精度确定子单元1462,高程定位精度确定子单元1463;
评估数值获取子单元1461,用于根据位置协方差矩阵,确定移动终端对应的过渡评估矩阵,获取过渡评估矩阵对应的第一水平评估数值、第二水平评估数值以及高程评估数值;
水平定位精度确定子单元1462,用于根据第一水平评估数值与第二水平评估数值之和的平方根,以及单位权中误差,确定终端位置对应的终端水平定位精度;
高程定位精度确定子单元1463,用于根据高程评估数值的平方根与单位权中误差,确定终端位置对应的终端高程定位精度。
其中,评估数值获取子单元1461,水平定位精度确定子单元1462,高程定位精度确定子单元1463的具体功能实现方式可以参见上述图7所对应实施例中的步骤S405-步骤S406,这里不再进行赘述。
在一些可行的实施方式中,该终端定位精度评估装置1还可以包括:卫星数量u判断模块15,正确结果确定模块16,错误结果确定模块17;
卫星数量u判断模块15,用于若卫星数量u小于数量阈值,则根据距离观测差值对终端位置进行评估;
正确结果确定模块16,用于当距离观测差值小于或等于误差阈值时,确定终端位置对应的位置精度评估结果为定位正确结果;
错误结果确定模块17,用于当距离观测值大于误差阈值时,确定终端位置对应的位置精度评估结果为定位错误结果。
其中,卫星数量u判断模块15,正确结果确定模块16,错误结果确定模块17的具体功能实现方式可以参见上述图7所对应实施例中的步骤S206,这里不再进行赘述。
在一些可行的实施方式中,终端定位精度评估装置1还可以包括:差值序列构建模块18,评估结果获取模块19;
差值序列构建模块18,用于若卫星数量u等于数量阈值,且伪距观测矩阵中的伪距观测方程所对应的卫星属于相同卫星系统,则根据距离观测差值构建差值序列;
评估结果获取模块19,用于获取差值序列对应的方均根,根据方均根对终端位置进行评估,得到终端位置对应的位置精度评估结果。
其中,差值序列构建模块18,评估结果获取模块19的具体功能实现方式可以参见上述图7所对应实施例中的步骤S206,这里不再进行赘述。
在一些可行的实施方式中,终端定位精度评估装置1还可以包括:评估结果确定模块19,功能服务提供模块20,终端位置修正模块21;
评估结果确定模块19,还用于根据位置精度指标确定终端位置对应的位置精度评估结果;
功能服务提供模块20,用于若位置精度评估结果为定位正确结果,则根据终端位置为移动终端提供应用功能服务;
终端位置修正模块21,用于若位置精度评估结果为定位错误结果,则根据无线局域网络和移动终端中的传感器,对终端位置进行修正。
其中,评估结果确定模块19,功能服务提供模块20,终端位置修正模块21的具体功能实现方式可以参见上述图3所对应实施例中的步骤S104,这里不再进行赘述。
本申请实施例中,移动终端可以通过卫星导航星历和初始伪距观测值计算得到卫星位置和卫星钟差,剔除初始伪距观测值中的观测值粗差得到目标伪距观测值,进而可以根据目标伪距观测值和终端位置构建伪距观测矩阵,基于该伪距观测矩阵可以得到终端位置对应的位置精度指标,该位置精度指标的数值大小可以用于准确评估移动终端的定位质量,进而可以提高移动终端的位置评估准确性;根据位置精度指标,结合无线局域网络和传感器对移动终端进行融合定位,并基于定位位置为移动终端提供应用功能服务,可以增强应用功能服务的推荐准确性;基于卫星伪距观测距离,实现对移动终端的定位精度评估过程,计算过程简单,可以提高移动终端的定位精度评估效率;本申请实施例提供的方法可以适用于智能手机、车机导航系统等场景中,适用性高。
请参见图11,图11是本申请实施例提供的一种计算机设备的结构示意图。如图11所示,该计算机设备1000可以包括:处理器1001,网络接口1004和存储器1005,此外,上述计算机设备1000还可以包括:用户接口1003,和至少一个通信总线1002。其中,通信总线1002用于实现这些组件之间的连接通信。其中,用户接口1003可以包括显示屏(Display)、键盘(Keyboard),可选用户接口1003还可以包括标准的有线接口、无线接口。可选的,网络接口1004可以包括标准的有线接口、无线接口(如WI-FI接口)。存储器1005可以是高速RAM存储器,也可以是非不稳定的存储器(non-volatile memory),例如至少一个磁盘存储器。可选的,存储器1005还可以是至少一个位于远离前述处理器1001的存储装置。如图11所示,作为一种计算机可读存储介质的存储器1005中可以包括操作系统、网络通信模块、用户接口模块以及设备控制应用程序。
在如图11所示的计算机设备1000中,网络接口1004可提供网络通讯功能;而用户接口1003主要用于为用户提供输入的接口;而处理器1001可以用于调用存储器1005中存储的设备控制应用程序,以实现:
获取卫星导航星历,获取移动终端与卫星之间的初始伪距观测值,根据初始伪距观测值和卫星导航星历,确定卫星对应的卫星位置和卫星钟差;该卫星导航星历包括用于确定卫星位置的参数列表;
获取移动终端对应的终端位置和终端钟差,根据卫星位置、卫星钟差、终端位置以及终端钟差,对初始伪距观测值进行筛选,得到目标伪距观测值;
根据目标伪距观测值和终端位置,构建移动终端与卫星之间的伪距观测矩阵;
根据伪距观测矩阵对终端位置进行评估,确定终端位置对应的位置精度指标;该位置精度指标用于表征移动终端的定位精度。
应当理解,本申请实施例中所描述的计算机设备1000可执行前文图3、图7-图9任一个所对应实施例中对终端定位精度评估方法的描述,也可执行前文图10所对应实施例中对终端定位精度评估装置1的描述,在此不再赘述。另外,对采用相同方法的有益效果描述,也不再进行赘述。
此外,这里需要指出的是:本申请实施例还提供了一种计算机可读存储介质,且计算机可读存储介质中存储有前文提及的终端定位精度评估装置1所执行的计算机程序,且计算机程序包括程序指令,当处理器执行程序指令时,能够执行前文图3、图7-图9任一个所对应实施例中对终端定位精度评估方法的描述,因此,这里将不再进行赘述。另外,对采用相同方法的有益效果描述,也不再进行赘述。对于本申请所涉及的计算机可读存储介质实施例中未披露的技术细节,请参照本申请方法实施例的描述。作为示例,程序指令可被部署在一个计算设备上执行,或者在位于一个地点的多个计算设备上执行,又或者,在分布在多个地点且通过通信网络互连的多个计算设备上执行,分布在多个地点且通过通信网络互连的多个计算设备可以组成区块链系统。
此外,需要说明的是:本申请实施例还提供了一种计算机程序产品或计算机程序,该计算机程序产品或者计算机程序可以包括计算机指令,该计算机指令可以存储在计算机可读存储介质中。计算机设备的处理器从计算机可读存储介质读取该计算机指令,处理器可以执行该计算机指令,使得该计算机设备执行前文图3、图7-图9任一个所对应实施例中对终端定位精度评估方法的描述,因此,这里将不再进行赘述。另外,对采用相同方法的有益效果描述,也不再进行赘述。对于本申请所涉及的计算机程序产品或者计算机程序实施例中未披露的技术细节,请参照本申请方法实施例的描述。
需要说明的是,对于前述的各个方法实施例,为了简单描述,故将其都表述为一系列的动作组合,但是本领域技术人员应该知悉,本申请并不受所描述的动作顺序的限制,因为依据本申请,某一些步骤可以采用其他顺序或者同时进行。其次,本领域技术人员也应该知悉,说明书中所描述的实施例均属于优选实施例,所涉及的动作和模块并不一定是本申请所必须的。
本申请实施例方法中的步骤可以根据实际需要进行顺序调整、合并和删减。
本申请实施例装置中的模块可以根据实际需要进行合并、划分和删减。
本领域普通技术人员可以理解实现上述实施例方法中的全部或部分流程,是可以通过计算机程序来指令相关的硬件来完成,计算机程序可存储于一计算机可读取存储介质中,该程序在执行时,可包括如上述各方法的实施例的流程。其中,存储介质可为磁碟、光盘、只读存储器(Read-Only Memory,ROM)或随机存储器(Random Access Memory,RAM)等。
以上所揭露的仅为本申请较佳实施例而已,当然不能以此来限定本申请之权利范围,因此依本申请权利要求所作的等同变化,仍属本申请所涵盖的范围。
Claims (15)
1.一种终端定位精度评估方法,其特征在于,包括:
获取卫星导航星历,获取移动终端与卫星之间的初始伪距观测值,根据所述初始伪距观测值和所述卫星导航星历,确定所述卫星对应的卫星位置和卫星钟差;
获取所述移动终端对应的终端位置和终端钟差,根据所述卫星位置、所述卫星钟差、所述终端位置以及所述终端钟差,对所述初始伪距观测值进行筛选,得到目标伪距观测值;
根据所述目标伪距观测值和所述终端位置,构建所述移动终端与所述卫星之间的伪距观测矩阵;
根据所述伪距观测矩阵对所述终端位置进行评估,确定所述终端位置对应的位置精度指标;所述位置精度指标用于表征所述移动终端的定位精度。
2.根据权利要求1所述的方法,其特征在于,所述根据所述初始伪距观测值和所述卫星导航星历,确定所述卫星对应的卫星位置和卫星钟差,包括:
获取所述移动终端接收到卫星发射信号的接收时间,根据所述初始伪距观测值、所述接收时间以及光速值,确定所述卫星发射信号的发射时间;
根据所述发射时间以及所述卫星导航星历中所包含的参数列表,确定所述卫星对应的卫星位置和卫星钟差。
3.根据权利要求1-2任一项所述的方法,其特征在于,所述根据所述卫星位置、所述卫星钟差、所述终端位置以及所述终端钟差,对所述初始伪距观测值进行筛选,得到目标伪距观测值,包括:
根据所述卫星位置、所述卫星钟差、所述终端位置以及所述终端钟差,确定所述移动终端与所述卫星之间的距离估计值;
根据所述距离估计值和所述初始伪距观测值之间的残差值,对所述初始伪距观测值进行筛选,得到目标伪距观测值。
4.根据权利要求3所述的方法,其特征在于,所述卫星的数量为N个,N个卫星属于G个卫星系统,每个卫星均对应一个距离估计值和一个初始伪距观测值,N和G均为正整数;
所述根据所述距离估计值和所述初始伪距观测值之间的残差值,对所述初始伪距观测值进行筛选,得到目标伪距观测值,包括:
获取N个卫星分别对应的残差值,根据所述G个卫星系统,对所述N个卫星分别对应的残差值进行分组,得到G个第一序列;一个卫星对应的残差值是由该卫星对应的距离估计值和初始伪距观测值所确定的;
获取所述G个第一序列中的第l个第一序列,对所述第l个第一序列中所包含的残差值进行排序,将排序后的第l个第一序列确定为第二序列,获取所述第二序列对应的初始掩码数组;l为小于G的正整数;
根据所述第二序列对应的中位数与所述第二序列中所包含的残差值,对所述第二序列进行更新,得到第三序列;
根据所述第三序列对应的下四分位数,以及所述第三序列对应的上四分位数,对所述第三序列进行更新,得到第四序列;
根据所述第二序列对应的下四分位数、所述第二序列对应的上四分位数以及所述第四序列对应的中位数,对所述初始掩码数组进行更新,得到目标掩码数组;
根据所述目标掩码数组中的非零元素,对所述第二序列中的残差值所对应的初始伪距观测值进行筛选,得到所述第二序列对应的目标伪距观测值。
5.根据权利要求4所述的方法,其特征在于,所述根据所述第三序列对应的下四分位数,以及所述第三序列对应的上四分位数,对所述第三序列进行更新,得到第四序列,包括:
根据所述第三序列对应的下四分位数,以及所述第三序列对应的上四分位数,确定所述第三序列中所包含的元素数值的有效取值范围;
将所述第三序列中处于所述有效取值范围的元素数值,确定为所述第四序列。
6.根据权利要求1所述的方法,其特征在于,所述目标伪距观测值的数量为u个,u个目标伪距观测值是指所述移动终端分别与G个卫星系统中的u个卫星之间的距离观测值,u和G均为正整数;
所述根据所述目标伪距观测值和所述终端位置,构建所述移动终端与所述卫星之间的伪距观测矩阵,包括:
根据所述G个卫星系统,对u个目标伪距观测值进行分组,得到G个观测值数组;
获取所述u个卫星分别对应的载噪比,根据所述载噪比确定所述G个观测值数组分别对应的参考卫星,将所述参考卫星对应的目标伪距观测值确定为参考伪距观测值;所述参考卫星是指每个观测值数组所对应的卫星中具有最大载噪比的卫星;
将所述u个卫星中除所述参考卫星之外的卫星确定为非参考卫星,根据所述参考伪距观测值与所述非参考卫星对应的目标伪距观测值,确定所述G个观测值数组关联的距离观测差值;一个距离观测差值用于表征一个观测值数组中的一个非参考卫星,与该观测值数组对应的参考卫星之间的观测距离;
根据所述终端位置、光速值、所述非参考卫星对应的卫星位置和卫星钟差,以及所述参考卫星对应的卫星位置和卫星钟差,确定所述G个观测值数组关联的距离估计差值;一个距离估计差值用于表征一个观测值数组中的一个非参考卫星,与该观测值数组对应的参考卫星之间的估计距离;
根据所述距离观测差值和所述距离估计差值,构建所述G个观测值数组分别对应的伪距观测方程,将所述G个观测值数组分别对应的伪距观测方程组合为伪距观测矩阵。
7.根据权利要求6所述的方法,其特征在于,所述根据所述伪距观测矩阵对所述终端位置进行评估,确定所述终端位置对应的位置精度指标,包括:
若卫星数量u大于数量阈值,则获取所述u个卫星分别对应的伪距观测值误差;
根据所述G个观测值数组关联的非参考卫星的伪距观测值误差,以及所述G个观测值数组关联的参考卫星的伪距观测值误差,生成观测值误差矩阵,获取所述观测值误差矩阵中的对角元素构建观测权阵;
获取所述伪距观测矩阵针对所述终端位置的微分矩阵,根据所述观测权阵和所述微分矩阵,确定所述移动终端对应的位置协方差矩阵;
根据所述观测权阵、所述距离观测差值所组成的矩阵以及所述卫星数量u,确定所述u个卫星对应的单位权中误差;
根据所述位置协方差矩阵和所述单位权中误差,确定所述终端位置对应的位置精度指标。
8.根据权利要求7所述的方法,其特征在于,所述位置精度指标包括终端水平定位精度和终端高程定位精度;
所述根据所述位置协方差矩阵和所述单位权中误差,确定所述终端位置对应的位置精度指标,包括:
根据所述位置协方差矩阵,确定所述移动终端对应的过渡评估矩阵,获取所述过渡评估矩阵对应的第一水平评估数值、第二水平评估数值以及高程评估数值;
根据所述第一水平评估数值与所述第二水平评估数值之和的平方根,以及所述单位权中误差,确定所述终端位置对应的终端水平定位精度;
根据所述高程评估数值的平方根与所述单位权中误差,确定所述终端位置对应的终端高程定位精度。
9.根据权利要求7所述的方法,其特征在于,还包括:
若所述卫星数量u小于所述数量阈值,则根据所述距离观测差值对所述终端位置进行评估;
当所述距离观测差值小于或等于误差阈值时,确定所述终端位置对应的位置精度评估结果为定位正确结果;
当所述距离观测值大于所述误差阈值时,确定所述终端位置对应的位置精度评估结果为定位错误结果。
10.根据权利要求7所述的方法,其特征在于,还包括:
若所述卫星数量u等于所述数量阈值,且所述伪距观测矩阵中的伪距观测方程所对应的卫星属于相同卫星系统,则根据所述距离观测差值构建差值序列;
获取所述差值序列对应的方均根,根据所述方均根对所述终端位置进行评估,得到所述终端位置对应的位置精度评估结果。
11.根据权利要求1所述的方法,其特征在于,还包括:
根据所述位置精度指标确定所述终端位置对应的位置精度评估结果;
若所述位置精度评估结果为定位正确结果,则根据所述终端位置为所述移动终端提供应用功能服务;
若所述位置精度评估结果为定位错误结果,则根据无线局域网络和所述移动终端中的传感器,对所述终端位置进行修正。
12.根据权利要求1所述的方法,其特征在于,所述获取卫星导航星历,包括:
向卫星定位基站服务器发送星历获取请求,以使所述卫星定位基站服务器根据所述星历获取请求,在卫星星历数据库中获取所述卫星导航星历;
接收所述卫星定位基站服务器返回的所述卫星导航星历;所述卫星导航星历以二进制流的形式进行传输。
13.一种终端定位精度评估装置,其特征在于,包括:
位置确定模块,用于获取卫星导航星历,获取移动终端与卫星之间的初始伪距观测值,根据所述初始伪距观测值和所述卫星导航星历,确定所述卫星对应的卫星位置和卫星钟差;
观测值获取模块,用于获取所述移动终端对应的终端位置和终端钟差,根据所述卫星位置、所述卫星钟差、所述终端位置以及所述终端钟差,对所述初始伪距观测值进行筛选,得到目标伪距观测值;
观测矩阵构建模块,用于根据所述目标伪距观测值和所述终端位置,构建所述移动终端与所述卫星之间的伪距观测矩阵;
精度评估模块,用于根据所述伪距观测矩阵对所述终端位置进行评估,确定所述终端位置对应的位置精度指标;所述位置精度指标用于表征所述移动终端的定位精度。
14.一种计算机设备,其特征在于,包括存储器和处理器;
所述存储器与所述处理器相连,所述存储器用于存储计算机程序,所述处理器用于调用所述计算机程序,以使得所述计算机设备执行权利要求1至12任一项所述的方法。
15.一种计算机可读存储介质,其特征在于,所述计算机可读存储介质中存储有计算机程序,所述计算机程序适于由处理器加载并执行,以使得具有所述处理器的计算机设备执行权利要求1至12任一项所述的方法。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110011939.8A CN112327340B (zh) | 2021-01-06 | 2021-01-06 | 终端定位精度评估方法、装置、设备以及介质 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110011939.8A CN112327340B (zh) | 2021-01-06 | 2021-01-06 | 终端定位精度评估方法、装置、设备以及介质 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN112327340A true CN112327340A (zh) | 2021-02-05 |
CN112327340B CN112327340B (zh) | 2021-04-13 |
Family
ID=74301510
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202110011939.8A Active CN112327340B (zh) | 2021-01-06 | 2021-01-06 | 终端定位精度评估方法、装置、设备以及介质 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN112327340B (zh) |
Cited By (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN113093250A (zh) * | 2021-06-04 | 2021-07-09 | 腾讯科技(深圳)有限公司 | 伪距观测数据修复方法、定位信息确定方法及装置 |
CN113281796A (zh) * | 2021-07-23 | 2021-08-20 | 腾讯科技(深圳)有限公司 | 位置确定方法、速度确定方法、装置、设备和存储介质 |
CN113608248A (zh) * | 2021-06-25 | 2021-11-05 | 北京建筑大学 | 北斗5g融合的高精度巡检人员定位方法及相关设备 |
CN113821816A (zh) * | 2021-11-18 | 2021-12-21 | 杭州格物智安科技有限公司 | 一种基于位置的区块链共识方法、系统及装置 |
CN115268252A (zh) * | 2022-08-05 | 2022-11-01 | 腾讯科技(深圳)有限公司 | 时间管理方法、装置、计算机、可读存储介质及程序产品 |
CN116520371A (zh) * | 2023-04-28 | 2023-08-01 | 银河航天(成都)通信有限公司 | 提高移动终端的差分定位精度的方法、装置及存储介质 |
Citations (18)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101799524A (zh) * | 2009-07-10 | 2010-08-11 | 中国测绘科学研究院 | 全球导航卫星系统接收机自主完备性监测的方法 |
CN101950025A (zh) * | 2010-08-11 | 2011-01-19 | 中国电子科技集团公司第二十研究所 | 用于局域增强系统的数据质量监测方法 |
CN102981174A (zh) * | 2012-09-12 | 2013-03-20 | 航天东方红卫星有限公司 | 一种gps天线相位中心变化对相对定位精度的修正方法 |
CN103076618A (zh) * | 2012-12-26 | 2013-05-01 | 北京空间飞行器总体设计部 | 一种星载测量型gnss接收机地面验证和性能评测方法 |
CN103941270A (zh) * | 2014-02-28 | 2014-07-23 | 北京邮电大学 | 一种多系统融合定位的方法及装置 |
CN104215977A (zh) * | 2014-09-09 | 2014-12-17 | 泰斗微电子科技有限公司 | 一种基于卫星导航系统的精度评估方法及系统 |
CN107064972A (zh) * | 2017-05-23 | 2017-08-18 | 中国测绘科学研究院 | 一种gnss原始观测数据变换方法 |
CN107229061A (zh) * | 2017-07-18 | 2017-10-03 | 武汉大学 | 一种基于低轨卫星的星地差分实时精密定位方法 |
CN108919321A (zh) * | 2018-05-18 | 2018-11-30 | 长安大学 | 一种基于尝试法的gnss定位粗差探测方法 |
CN109358348A (zh) * | 2018-11-08 | 2019-02-19 | 电子科技大学 | 一种多模卫星导航接收机直接定位解算方法 |
CN110286395A (zh) * | 2019-07-11 | 2019-09-27 | 国网上海市电力公司 | 一种北斗系统定位精度计算方法 |
CN110515097A (zh) * | 2019-09-02 | 2019-11-29 | 江苏省测绘工程院 | 应用于基准站的gnss卫星观测粗差剔除方法和装置 |
CN110927665A (zh) * | 2019-11-09 | 2020-03-27 | 武汉中海庭数据技术有限公司 | 一种千寻云迹虚拟参考站定位精度的验证方法及系统 |
CN111025341A (zh) * | 2019-11-22 | 2020-04-17 | 中国科学院上海天文台 | 卫星轨道的误差精化方法 |
US20200166652A1 (en) * | 2016-06-06 | 2020-05-28 | Brian G. Agee | Multi-Subband Methods for Reduced Complexity, Wideband Blind Resilient Detection and Geo-Observable Estimation of Global Navigation Satellite Signals |
CN111443364A (zh) * | 2020-03-27 | 2020-07-24 | 惠州市德赛西威智能交通技术研究院有限公司 | 一种车载定位接收机精度测试方法及其系统 |
CN111626406A (zh) * | 2020-05-19 | 2020-09-04 | 辽宁工程技术大学 | 一种基于一次差分的bds钟差预测方法 |
CN111736188A (zh) * | 2020-03-09 | 2020-10-02 | 腾讯科技(深圳)有限公司 | 一种卫星定位方法、装置、电子设备和存储介质 |
-
2021
- 2021-01-06 CN CN202110011939.8A patent/CN112327340B/zh active Active
Patent Citations (18)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101799524A (zh) * | 2009-07-10 | 2010-08-11 | 中国测绘科学研究院 | 全球导航卫星系统接收机自主完备性监测的方法 |
CN101950025A (zh) * | 2010-08-11 | 2011-01-19 | 中国电子科技集团公司第二十研究所 | 用于局域增强系统的数据质量监测方法 |
CN102981174A (zh) * | 2012-09-12 | 2013-03-20 | 航天东方红卫星有限公司 | 一种gps天线相位中心变化对相对定位精度的修正方法 |
CN103076618A (zh) * | 2012-12-26 | 2013-05-01 | 北京空间飞行器总体设计部 | 一种星载测量型gnss接收机地面验证和性能评测方法 |
CN103941270A (zh) * | 2014-02-28 | 2014-07-23 | 北京邮电大学 | 一种多系统融合定位的方法及装置 |
CN104215977A (zh) * | 2014-09-09 | 2014-12-17 | 泰斗微电子科技有限公司 | 一种基于卫星导航系统的精度评估方法及系统 |
US20200166652A1 (en) * | 2016-06-06 | 2020-05-28 | Brian G. Agee | Multi-Subband Methods for Reduced Complexity, Wideband Blind Resilient Detection and Geo-Observable Estimation of Global Navigation Satellite Signals |
CN107064972A (zh) * | 2017-05-23 | 2017-08-18 | 中国测绘科学研究院 | 一种gnss原始观测数据变换方法 |
CN107229061A (zh) * | 2017-07-18 | 2017-10-03 | 武汉大学 | 一种基于低轨卫星的星地差分实时精密定位方法 |
CN108919321A (zh) * | 2018-05-18 | 2018-11-30 | 长安大学 | 一种基于尝试法的gnss定位粗差探测方法 |
CN109358348A (zh) * | 2018-11-08 | 2019-02-19 | 电子科技大学 | 一种多模卫星导航接收机直接定位解算方法 |
CN110286395A (zh) * | 2019-07-11 | 2019-09-27 | 国网上海市电力公司 | 一种北斗系统定位精度计算方法 |
CN110515097A (zh) * | 2019-09-02 | 2019-11-29 | 江苏省测绘工程院 | 应用于基准站的gnss卫星观测粗差剔除方法和装置 |
CN110927665A (zh) * | 2019-11-09 | 2020-03-27 | 武汉中海庭数据技术有限公司 | 一种千寻云迹虚拟参考站定位精度的验证方法及系统 |
CN111025341A (zh) * | 2019-11-22 | 2020-04-17 | 中国科学院上海天文台 | 卫星轨道的误差精化方法 |
CN111736188A (zh) * | 2020-03-09 | 2020-10-02 | 腾讯科技(深圳)有限公司 | 一种卫星定位方法、装置、电子设备和存储介质 |
CN111443364A (zh) * | 2020-03-27 | 2020-07-24 | 惠州市德赛西威智能交通技术研究院有限公司 | 一种车载定位接收机精度测试方法及其系统 |
CN111626406A (zh) * | 2020-05-19 | 2020-09-04 | 辽宁工程技术大学 | 一种基于一次差分的bds钟差预测方法 |
Non-Patent Citations (1)
Title |
---|
龚栎澎等: "星间频间伪距定位模型的精度分析及程序实现", 《矿山测量》 * |
Cited By (11)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN113093250A (zh) * | 2021-06-04 | 2021-07-09 | 腾讯科技(深圳)有限公司 | 伪距观测数据修复方法、定位信息确定方法及装置 |
CN113093250B (zh) * | 2021-06-04 | 2021-08-27 | 腾讯科技(深圳)有限公司 | 伪距观测数据修复方法、定位信息确定方法及装置 |
CN113608248A (zh) * | 2021-06-25 | 2021-11-05 | 北京建筑大学 | 北斗5g融合的高精度巡检人员定位方法及相关设备 |
CN113608248B (zh) * | 2021-06-25 | 2023-06-13 | 北京建筑大学 | 北斗5g融合的高精度巡检人员定位方法及相关设备 |
CN113281796A (zh) * | 2021-07-23 | 2021-08-20 | 腾讯科技(深圳)有限公司 | 位置确定方法、速度确定方法、装置、设备和存储介质 |
CN113281796B (zh) * | 2021-07-23 | 2021-10-15 | 腾讯科技(深圳)有限公司 | 位置确定方法、速度确定方法、装置、设备和存储介质 |
CN113821816A (zh) * | 2021-11-18 | 2021-12-21 | 杭州格物智安科技有限公司 | 一种基于位置的区块链共识方法、系统及装置 |
CN113821816B (zh) * | 2021-11-18 | 2022-03-25 | 杭州格物智安科技有限公司 | 一种基于位置的区块链共识方法、系统及装置 |
CN115268252A (zh) * | 2022-08-05 | 2022-11-01 | 腾讯科技(深圳)有限公司 | 时间管理方法、装置、计算机、可读存储介质及程序产品 |
CN116520371A (zh) * | 2023-04-28 | 2023-08-01 | 银河航天(成都)通信有限公司 | 提高移动终端的差分定位精度的方法、装置及存储介质 |
CN116520371B (zh) * | 2023-04-28 | 2024-05-28 | 银河航天(成都)通信有限公司 | 提高移动终端的差分定位精度的方法、装置及存储介质 |
Also Published As
Publication number | Publication date |
---|---|
CN112327340B (zh) | 2021-04-13 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN112327340B (zh) | 终端定位精度评估方法、装置、设备以及介质 | |
US8255160B2 (en) | Integrated mobile terminal navigation | |
US7948433B2 (en) | Calculation method for network-specific factors in a network of reference stations for a satellite-based positioning system | |
CN111736185B (zh) | 终端定位方法、装置及计算机可读存储介质和终端设备 | |
RU2004134213A (ru) | Определение местоположения оконечного устройства беспроводной связи в смешанной системе определения местоположения | |
US11221415B2 (en) | Positioning device and positioning method | |
EP3106897A1 (en) | Gnss receiver with an on-board capability to implement an optimal error correction mode | |
US11275180B2 (en) | Positioning device and positioning method | |
CN113050142B (zh) | 终端设备的定位方法、装置、电子设备及可读存储介质 | |
CN114384569A (zh) | 终端定位方法、装置、设备以及介质 | |
CN112556696B (zh) | 一种对象定位方法、装置、计算机设备以及存储介质 | |
CN113835109B (zh) | 终端的定位方法、装置、电子设备、存储介质及程序产品 | |
CN112924997B (zh) | 目标终端定位方法、装置和存储介质及电子设备 | |
US20110187592A1 (en) | Satellite based position of a cellular terminal | |
KR20130036145A (ko) | 이동 정보 결정 장치, 수신기 및 그에 의한 방법 | |
CN112904390A (zh) | 定位方法、装置、计算机设备和存储介质 | |
CN113281796B (zh) | 位置确定方法、速度确定方法、装置、设备和存储介质 | |
CN113917509B (zh) | 一种双差模糊度固定方法、设备以及可读存储介质 | |
CN115616637A (zh) | 一种基于三维格网多径建模的城市复杂环境导航定位方法 | |
CN112731268B (zh) | 一种差分数据的处理方法和定位跟踪系统 | |
CN117008166A (zh) | 定位质量评估方法、装置、设备、存储介质和程序产品 | |
CN114167469A (zh) | 一种基于5g/gnss组合的车辆导航信息监测方法和装置 | |
CN115242297B (zh) | 移动终端的运动参数确定方法、装置、设备和存储介质 | |
CN111596330B (zh) | 定位方法、装置、终端和存储介质 | |
CN116893433B (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 | ||
GR01 | Patent grant | ||
GR01 | Patent grant | ||
REG | Reference to a national code |
Ref country code: HK Ref legal event code: DE Ref document number: 40038747 Country of ref document: HK |