一种基于iGMAS超快速星历的实时精密单点定位方法与流程

本发明涉及高精度导航技术领域，尤其涉及一种基于igmas超快速星历的实时精密单点定位方法。

背景技术：

为了进一步推动多模gnss(globalnavigationsatellitesystem全球导航卫星系统)兼容和互操作，中国从2012年启动国际gnssigmas(internationalgnssmonitoring&assessmentsystem,监测评估系统)建设。igmas旨在建立我国自主的bds、gps、glonass和galileo导航卫星全弧段、多重覆盖的全球近实时跟踪网，监控gnss导航卫星的运行状况、信号质量和服务性能。

igmas能够向全球用户提供精密星历、钟差、地球定向参数等产品，为卫星导航技术试验提供支持，服务于科学研究以及各类应用，包括ppp(preciouspointpositioning,精密单点定位)技术。

相对于事后处理的精密单点定位，实时精准定位才能满足5g、物联网和无人驾驶等技术发展的需求，其中高质量精密星历、钟差改正产品的实时获取和应用是难点之一。目前国际上开展实时ppp应用主要基于igs(instrumentguidancesystem，仪表引导系统)实时、近实时产品，包括igs实时定位服务、igdg(internet-basedglobaldifferentialgps，基于互联网的全球差分定位系统)实时产品和bnc(bkgntripclient，bkgntrip客户端)软件等，以及西班牙gmv公司、trimble公司等推出的商业实时ppp服务。与国际上比较成熟的igs系统相比，我国主导的igmas建设周期短，跟踪站、分析中心和数据中心数量相对较少，因此igmas产品(特别是超快速产品)的实际应用仍需进一步发展。

目前，国内外应用igmas超快速产品开展的实时精密单点定位研究，主要依靠跟踪站的观测数据模拟实时ppp解算，未能实现igmas超快速数据的实时获取和普遍适用。

技术实现要素：

本发明提供一种基于igmas超快速星历的实时精密单点定位方法，能够利用igmas数据和ppp算法应用，提高双频gnss接收机原始定位的精度。

为达到上述目的，本发明采用如下技术方案：

一种基于igmas超快速星历的实时精密单点定位方法，包括：

采集双频gnss观测数据和igs/igmas提供的精密改正数；

利用所述双频gnss观测数据和实时精密改正数，采用实时ppp算法得到高精度定位。

进一步的，实时精密改正数采用下载软件下载，打通了数据下载与实际应用的通道。

进一步的，双频gnss观测数据包括：双频gnss伪距和载波相位观测值，gnss广播星历、gnss初步定位结果。

进一步的，精密改正数包括：gnss超快速星历、钟差。

进一步的，实时ppp算法包括：

对精密改正数中的gnss超快速星历、钟差进行处理，使得igmas超快速星历、钟差与观测数据采样时刻匹配，进行卫星的位置计算和钟差消除；

对双频gnss观测数据进行粗差和周跳的探测、处理,得到干净可用的观测数据；

将处理后的igmas超快速精密星历、钟差和干净可用的观测数据带入ppp常规模型；

对所述ppp模型采用扩展卡尔曼滤波进行参数估计；

对所述双频gnss观测数据中的各类误差采用各种模型进行处理。

进一步的，参数估计中所涉及的参数包括接收机位置、天顶对流层延迟、接收机钟差、系统间偏差和模糊度。

进一步的，误差包括卫星轨道及钟差、电离层延迟、对流层干延迟、卫星和接收机天线相位中心、相位缠绕、相对论效应、地球自转改正、固体潮、海潮、极潮。

进一步的，误差处理模型包括igmas超快速星历及钟差、双频无电离层组合、hopfield模型+gpt2+vmf、igs14.atx及pco+pcv、wu模型、iers协议、sagnac效应、got4.8模型。

本发明的有益效果是：

现有技术中，应用igmas超快速产品开展的实时精密单点定位研究主要依靠跟踪站的观测数据模拟实时ppp解算，本发明中提供了实时ppp算法，结合观测数据与精密改正数，有效地提升了定位精度，解决了igmas超快速数据的实时获取和普遍适用问题。

附图说明

为了更清楚地说明本发明实施例中的技术方案，下面将对实施例中所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图仅仅是本发明的一些实施例，对于本领域普通技术人员来讲，在不付出创造性劳动的前提下，还可以根据这些附图获得其它的附图。

图1是igmas数据下载应用用户界面；

图2是igmas数据实时下载程序设计流程图；

图3是实时ppp技术方案流程图；

图4是接收机自身定位和应用igmas超快速星历的ppp定位误差对比图；

图5是应用igmas超快速星历的ppp算法在enu方向的误差收敛图。

具体实施方式

为使本领域技术人员更好地理解本发明的技术方案，下面结合具体实施方式对本发明作进一步详细描述。

本发明实施例提供一种基于igmas超快速星历的实时精密单点定位方法，流程如图3所示，包括：

s1、利用下载软件igmasdownload实时下载igmas超快速星历，下载文件的命名格式为isuwwwwd_hh.sp3/clk.z(wwww为北斗周；d为星期：0为星期日，1-6表示星期一至六；hh为小时，分为00，06，12，18；sp3/clk分别表示星历和钟差文件)。igmasdownload软件可以实现对igmas提供的精密改正数的实时或分时下载，下载方式可以选择实时或者事后批量下载。

igmasdownload程序实时刷新igmas数据中心产品文件，下载到指定文件夹并自动解压，以备实时ppp算法读取应用。同时，事后模式可以根据用户需求批量下载其他产品，包括跟踪站地心坐标、地球自转参数、大气环境参数、频间偏差信息、电离层闪烁指数、民用监测评估结果、完好性产品。下载应用用户界面如图1所示，下载流程如图2所示。

s2、双频gnss接收机获取双频伪距、载波相位观测数据、广播星历和原始定位结果。

s3、对igmas精密星历数据与gnss观测数据进行预处理，包括粗差探测消除，接收机钟跳和模糊度周跳探测。

s4、利用处理后的精密星历与gnss观测数据，采用ppp常规模型对定位结果进行解算，采用双频伪距、载波相位观测值的无电离层组合作为函数模型，其基础模型表达式为

其中，分别为伪距、载波相位无电离层组合观测值，为接收机和卫星间几何距离，s为卫星编号，r为接收机编号，为对流层延迟，δtr、δt^s分别表示接收机和卫星钟差，c是光速，δm为多路径延迟，为包含卫星和接收机端初始相位和硬件延迟的无电离层组合模糊度，和分别表示两种组合观测值，即伪距和载波相的观测噪声及未被模型化的误差。

s5、采用卡尔曼滤波进行参数估计，待估计参数包括接收机位置、天顶对流层延迟、接收机钟差、系统间偏差和模糊度。再利用估计后的参数对双频gnss接收机得到的原始定位结果的三维位置参数(x,y,z)进行静态处理。观测值参数估计方程和非线性展开为：

yk＝hkx+vk(7)

式中，yk是新获得的量测值；hk是量测矩阵；x是待估计参数；vk是量测噪声；和表示去电离层组合伪距观和载波相位观测值；(xr,0,yr,0,zr,0)为测站位置；(x^s,y^s,z^s)为经过精密星历修正后的gnss卫星位置；为测站与gnss卫星的空间距离；f1表示gps卫星信号频率；c为光速；εp,if和εφ,if分别表示去电离层组合伪距观测方程误差和载波相位观测方程误差。

s6、对双频gnss观测数据中的各类误差采用各种模型进行处理，处理后得到单点高精度定位参数。

如图3所示，本实施例中引入北斗伴侣m2检验算法有效性，北斗伴侣m2自身测量精度高，通过对比，能检测本发明的测算结果是否满足高精度条件。

具体的，各项误差处理方法如表1所示：

表1各项误差处理方法

本发明使用novateloem617双频接收机(单点定位精度：水平5m，高程10m)实时接收gps观测数据，截止高度角设为10°，使用北斗伴侣m2rtk固定解(定位精度：水平0.02m+1ppm，高程0.04m+1ppm)作为定位结果参考标准。在定位参考点的设计上，通过将novatel天线和北斗伴侣m2架设到两台三脚架并调整抵消其高度差，预先测量水平和高度误差并在算法中予以修正。使用igmas超快速星历和钟差产品，对gps实测数据进行实时精密单点定位实验，并与基于igs超快速星历的实时ppp定位进行对比，其中定位结果对比如表2和图4、图5所示。

表2接收机和实时ppp定位误差均方根统计

本发明的有益效果是：

本发明中提供了实时ppp算法，结合观测数据与精密改正数，有效地提升了定位精度，解决i了gmas超快速数据的实时获取和普遍适用问题。

本发明将双频gnss接收机的原始定位精度提高了一个量级，得到了分米级的高精度定位。

以上所述，仅为本发明的具体实施方式，但本发明的保护范围并不局限于此，任何熟悉本技术领域的技术人员在本发明揭露的技术范围内，可轻易想到的变化或替换，都应涵盖在本发明的保护范围之内。因此，本发明的保护范围应该以权利要求的保护范围为准。