CN109143289B - 一种gnss单站位移监测方法 - Google Patents
一种gnss单站位移监测方法 Download PDFInfo
- Publication number
- CN109143289B CN109143289B CN201811346260.9A CN201811346260A CN109143289B CN 109143289 B CN109143289 B CN 109143289B CN 201811346260 A CN201811346260 A CN 201811346260A CN 109143289 B CN109143289 B CN 109143289B
- Authority
- CN
- China
- Prior art keywords
- point
- observation
- equation
- satellite
- receiver
- 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.)
- Active
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/13—Receivers
- G01S19/35—Constructional details or hardware or software details of the signal processing chain
- G01S19/37—Hardware or software details of the signal processing chain
Abstract
本发明公开了一种GNSS单站位移监测方法,适用于使用单个GNSS接收机在短时间内测量空间距离较近的两点相对位置。依据对流层延迟等误差在短时间内较强相关性和没有周跳发生时模糊度固定不变的特点,提出GNSS单站位移监测技术,研究对GNSS单站数据时间序列信息进行处理,通过削弱或消除观测量中的时空相关误差,解算精确的单个接收机位移的基线向量,快速测量得到两点的相对位置关系。本发明技术不需要基站和通讯设备,不依赖精密星历,使用单个接收机可达到厘米级定位精度,处理方法简单,节省设备购置成本和提高测量生产效率。
Description
技术领域
本发明涉及测量技术领域,特别涉及一种GNSS(Global Navigation SatelliteSystem,全球导航卫星系统)单站位移监测方法。
背景技术
近年来,随着GNSS接收机性能和解算技术的不断发展进步,GNSS精确定位方法取得了系列进步,其应用领域也得到了极大拓展,GNSS测量技术大致分为静态和动态测量两类。
静态测量需要观测时间长且同步观测接收机较多,数据处理需要专门软件进行解算、平差后才能应用于工作,主要用于高精度控制测量领域。
动态测量观测时间较短,现场可以快速得到测量结果,但因其精度较低主要在导航、测绘领域中放样、测图、变形监测等方面得到广泛应用,动态测量又大致分为以下几种方式:GNSS单点定位(也叫绝对定位)、双差相对定位(RTK)、非差精密单点定位(PPP)。
GNSS单点定位应用只需要一台GNSS接收机即可测量快速测量定位,但只能达到米级定位精度。
双差相对定位(RTK)需要最少一台基站和至少一台流动站接收机,将基站架设在已知点上进行测量,测量时间几秒到几十秒,解算即可得到厘米级精度的结果,目前主要用于测绘领域精度要求较高的地方和工程测量施工放样。
非差精密单点定位(PPP),比静态测量观测时间短,比双差相对定位(RTK)观测时间长,定位精度可以达到厘米级精度,可用于建筑物进行了结构健康监测。
但是,以上三种定位方式都存在局限性:GNSS单点定位因定位精度在米级,限制其一般只能用于导航领域,不能满足测图、放样等工作需要;双差相对定位(RTK)定位精度满足测图、放样等工作需要,但是需要基站作为参考站,所得到的动态定位结果是相对于参考站的位移情况,基站是必需的,要进行精确定位需要至少一台基站接收机和至少一台流动站接收机,要进行一系列基站设置、校正等准备工作后才能进行测量,因此增加了应用成本,作业效率低;非差精密单点定位(PPP)定位不需要基站,定位方式灵活,可用观测数据多,但是PPP定位模式需要接收精密星历,并且受到很多误差因素的影响,数据处理过程变得相对复杂,模糊度的快速固定也没有得到很好解决,与网络RTK相比,PPP定位精度较低,初始化时间更长,以至于至今都没有很好的商业接收机产品问世。
发明内容
本发明的目的在于克服现有技术中所存在的现有GNSS单点定位精度在米级,不能满足测图、放样等工作需要;RTK精度在厘米级,但是需要至少一台基站接收机和至少一台流动站接收机,要进行一系列基站设置、校正等准备工作后才能进行测量,因此增加了应用成本,作业效率低;PPP精度在厘米级,但是需要接收精密星历,并且受到很多误差因素的影响,数据处理过程变得相对复杂,模糊度的快速固定也没有得到很好解决,初始化时间较长的上述不足,提供一种GNSS单站位移监测方法。
为了实现上述发明目的,本发明提供了以下技术方案:
一种GNSS单站位移监测方法,将GNSS接收机由A点移动至B点,A点观测时刻为t1,B点观测时刻为t2,移动过程中保持所述接收机对卫星的跟踪,还
包括在A点处和B点处分别建立该点处所述接收机与卫星相位观测之差的单差观测方程的步骤;
包括将A点处的单差观测方程和B点处的单差观测方程二次求差,得到双差观测方程的步骤;
包括所述接收机利用不同频率信号消除电离层因素对所述双差观测方程影响的步骤;
包括计算改正所述双差观测方程中卫星钟差变化的步骤;
包括将所述双差观测方程线性化,得到误差方程式的步骤;
包括通过最小二乘平差原理,获取精确的基线向量的步骤。
进一步的,所述GNSS包括美国GPS、俄罗斯GLONASS、欧盟GALILEO、中国BeiDou等卫星导航系统。
进一步的,所述卫星包括观测卫星s和参考卫星r,且参考卫星r的高度角大于观测卫星s的高度角。
由于A点到B点的过程中,所述接收机都保持对卫星的跟踪,从A点到B点的距离和时间都较短,很多误差都是时空相关误差,可通过观测值历元间求差的方法得以消除,比如相对论效应、潮汐效应、轨道误差和对流层延迟等,且整周模糊度不会发生变化,均表示为对流层延迟轨道误差简化得到上述双差观测方程,其中,表示双差算子符号,
进一步的,由于电离层延迟变化较快,在数分钟内,即可发生较大变化,故需要使用双频进行消除,根据所述接收机的不同波长λ1和λ2,分别列出两个波长的所述双差观测方程:
其中,IF表示与电离层无关组合。
其中,a1和a2分别是观测卫星s的钟速和钟漂,t0是卫星星历的参考时刻,t1是所述接收机位于A点的时刻,t2是所述接收机位于B点的时刻。
进一步的,将所述双差观测方程线性化,得到误差方程式包括以下步骤:
进一步的,通过观测多颗卫星建立观测误差方程组,利用最小二乘平差原理,解算得到A、B两点基线向量的改正数。
进一步的,还包括利用卫星的不同高度角,对观测值进行加权的步骤,以进一步提高解算的基线向量的改正数的精度。
进一步的,根据A点和B点的初始坐标计算A、B两点的初始基线向量,再利用所述基线向量改正数,得到A点和B点精确的基线向量。
综上所述,由于采用了上述技术方案,本发明的有益效果是:
运用本发明所述的一种GNSS单站位移监测方法,适用于使用单个GNSS接收机在短时间内测量空间距离较近的两点相对位置,依据对流层延迟等误差在短时间内较强相关性和没有周跳发生时模糊度固定不变的特点,提出GNSS单站位移监测技术,研究对GNSS单站数据时间序列信息进行处理,通过削弱或消除观测量中的时空相关误差,解算精确的单个接收机位移的基线向量,快速测量得到两点的相对位置关系;本发明技术不需要基站和通讯设备,不依赖接收精密星历,使用单个接收机即可达到厘米级定位精度,处理方法简单,节省设备购置成本和提高测量生产效率。
附图说明
图1为实施例中GNSS接收机在t1时刻的A点移动至t2时刻的B点的示意图(本发明所述方法外业作业示意图);
图2为实施例中本发明所述方法外业作业流程图;
图3为实施例中本发明所述方法内业作业流程图。
具体实施方式
下面结合试验例及具体实施方式对本发明作进一步的详细描述。但不应将此理解为本发明上述主题的范围仅限于以下的实施例,凡基于本发明内容所实现的技术均属于本发明的范围。
实施例
如图1-3所示,本发明所述的一种基于美国GPS的GNSS单站位移监测方法,包括以下步骤:
步骤一、在A点设置GNSS接收机,启动所述接收机进行数秒短时间的观测,A点与观测卫星s的相位观测可表示为:
同理可得A点与参考卫星r(参考卫星r的高度角大于观测卫星s的高度角)的相位观测方程为:
步骤二、对不同卫星间观测值求差,得到A点与观测卫星s相对于参考卫星r的观测方程为:
步骤三、保持所述接收机对卫星的跟踪,将所述接收机迅速移动到B点(若移动过程中所述接收机与卫星断连,重新执行所述步骤一),进行数秒短时间的观测,则B点与观测卫星s相对于参考卫星r的观测方程为:
步骤四、由于A点到B点的过程中,所述接收机都保持对卫星的跟踪,从A点到B点的距离和时间都较短,很多误差都是时空相关误差,可通过观测值历元间求差的方法得以消除,比如相对论效应、潮汐效应、轨道误差和对流层延迟等,且整周模糊度不会发生变化,均表示为对流层延迟轨道误差因此,对在A点和B点的观测方程二次求差,并可以简化为双差观测方程:
步骤五、由于电离层延迟变化较快,在数分钟内,即可发生米级的较大变化,故需要使用双频进行消除,根据不同波长λ1和λ2,可分别列出双差观测方程:
步骤七、卫星钟差的变化可通过二次多项式模型进行计算改正:
其中,a1和a2分别是观测卫星s的钟速和钟漂,t0是卫星星历的参考时刻,t1是所述接收机位于A点的时刻,t2是所述接收机位于B点的时刻;
步骤八、把已知项放在方程左边,把未知项放在方程右边,得到方程:
将该方程线性化,得到如下方程式:
步骤九、卫星到所述接收机几何距离的变化和卫星钟差的变化体现在相位观测值的变化上,通过观测多颗卫星建立观测误差方程组,利用卫星的不同高度角,对观测值进行加权,利用最小二乘平差原理,解算得到A、B两点基线向量的改正数,根据A点和B点的初始坐标计算A、B两点的初始基线向量,再利用所述基线向量改正数,由此得到A点和B点精确的基线向量:
ΔxAB=xB0-xA0+dx,
ΔyAB=yBO-yAO+dy,
ΔzAB=zB0-zA0+dz。
运用本发明所述的一种GNSS单站位移监测方法,通过编程实现其算法。在实际数据测试中,通过相邻两个点之间的相对位置测试和逐点之间(若干点两两之间)进行相对位置测试,通过本发明所叙述的方法得到的距离值和利用高精度全站仪实测结果进行比较分析,较差差值大部分能达到厘米级;本发明适用于使用单个GNSS接收机在短时间内测量空间距离较近的两点相对位置,本发明技术不需要基站和通讯设备,不依赖接收精密星历,使用单个接收机即可达到厘米级定位精度,处理方法简单,节省设备购置成本和提高测量生产效率。
以上所述仅为本发明的较佳实施例而已,并不用以限制本发明,凡在本发明的精神和原则之内所作的任何修改、等同替换和改进等,均应包含在本发明的保护范围之内。
Claims (8)
1.一种GNSS单站位移监测方法,将GNSS接收机由A点移动至B点,A点观测时刻为t1,B点观测时刻为t2,其特征在于,移动过程中保持所述接收机对卫星的跟踪,还
包括在A点处和B点处分别建立该点处所述接收机与卫星相位观测之差的单差观测方程的步骤;
包括将A点处的单差观测方程和B点处的单差观测方程二次求差,得到双差观测方程的步骤;
包括所述接收机利用不同频率信号消除电离层因素对所述双差观测方程影响的步骤;
包括计算改正所述双差观测方程中卫星钟差变化的步骤;
包括将所述双差观测方程线性化,得到误差方程式的步骤;
包括通过最小二乘平差原理,获取精确的基线向量的步骤:卫星到所述接收机几何距离的变化和卫星钟差的变化体现在相位观测值的变化上,通过观测多颗卫星建立观测误差方程组,利用卫星的不同高度角,对观测值进行加权,利用最小二乘平差原理,解算得到A、B两点基线向量的改正数,根据A点和B点的初始坐标计算A、B两点的初始基线向量,再利用所述基线向量改正数,由此得到A点和B点精确的基线向量。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201811346260.9A CN109143289B (zh) | 2018-11-13 | 2018-11-13 | 一种gnss单站位移监测方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201811346260.9A CN109143289B (zh) | 2018-11-13 | 2018-11-13 | 一种gnss单站位移监测方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN109143289A CN109143289A (zh) | 2019-01-04 |
CN109143289B true CN109143289B (zh) | 2021-09-21 |
Family
ID=64805785
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201811346260.9A Active CN109143289B (zh) | 2018-11-13 | 2018-11-13 | 一种gnss单站位移监测方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN109143289B (zh) |
Families Citing this family (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN112184778B (zh) * | 2020-09-24 | 2022-04-26 | 武汉大学 | 基于gnss动态定位时间序列分割的短期位移探测方法 |
CN113163486B (zh) * | 2021-04-20 | 2022-05-17 | 武汉大学 | 一种基于5g tdoa的相对定位方法及系统 |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US6397147B1 (en) * | 2000-06-06 | 2002-05-28 | Csi Wireless Inc. | Relative GPS positioning using a single GPS receiver with internally generated differential correction terms |
KR20100025973A (ko) * | 2008-08-28 | 2010-03-10 | 한국전자통신연구원 | 단일 지상국에서 위성 궤도 결정 방법 |
CN106125113A (zh) * | 2016-06-20 | 2016-11-16 | 武汉大学 | 一种利用多系统gnss观测值的高精度基线解算方法 |
CN106291639A (zh) * | 2016-08-31 | 2017-01-04 | 和芯星通科技(北京)有限公司 | 一种gnss接收机实现定位的方法及装置 |
CN108363084A (zh) * | 2018-01-18 | 2018-08-03 | 和芯星通科技(北京)有限公司 | 利用卫星定位的方法和装置、卫星导航接收机、存储介质 |
-
2018
- 2018-11-13 CN CN201811346260.9A patent/CN109143289B/zh active Active
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US6397147B1 (en) * | 2000-06-06 | 2002-05-28 | Csi Wireless Inc. | Relative GPS positioning using a single GPS receiver with internally generated differential correction terms |
KR20100025973A (ko) * | 2008-08-28 | 2010-03-10 | 한국전자통신연구원 | 단일 지상국에서 위성 궤도 결정 방법 |
CN106125113A (zh) * | 2016-06-20 | 2016-11-16 | 武汉大学 | 一种利用多系统gnss观测值的高精度基线解算方法 |
CN106291639A (zh) * | 2016-08-31 | 2017-01-04 | 和芯星通科技(北京)有限公司 | 一种gnss接收机实现定位的方法及装置 |
CN108363084A (zh) * | 2018-01-18 | 2018-08-03 | 和芯星通科技(北京)有限公司 | 利用卫星定位的方法和装置、卫星导航接收机、存储介质 |
Non-Patent Citations (2)
Title |
---|
利用高频GPS单站数据解算同震地表位移关键算法研究;李艳辉;《中国优秀硕士学位论文全文数据库基础科学辑》;20170715;正文第4页第2段、第6页至第22页第1段、及图2-1 * |
同震地表位移的单站时间基线模型解算;李艳辉等;《测绘学报》;20180630;第43卷(第6期);全文 * |
Also Published As
Publication number | Publication date |
---|---|
CN109143289A (zh) | 2019-01-04 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN108415049B (zh) | 提高网络rtk双差宽巷模糊度固定正确率的方法 | |
CN109581452B (zh) | 一种gnss参考站载波相位整周模糊度解算方法 | |
US7498979B2 (en) | Fast decimeter-level GNSS positioning | |
CN108802782B (zh) | 一种惯导辅助的北斗三频载波相位整周模糊度求解方法 | |
CN105807300B (zh) | 一种用北斗双频接收机进行高精度动态单点定位的方法 | |
US10422884B2 (en) | Method and system for performing precise point positioning (PPP) ambiguity resolution using GNSS triple frequency signals | |
CN103837879B (zh) | 基于北斗系统民用载波相位组合实现高精度定位的方法 | |
CN108363084A (zh) | 利用卫星定位的方法和装置、卫星导航接收机、存储介质 | |
CN110389364B (zh) | 定位方法、装置、计算机设备和存储介质 | |
CN110687556B (zh) | 一种适用于laas的多径误差模型化方法 | |
CN111596322A (zh) | 一种短基线条件下接收机伪距测量一致性评估方法 | |
CN107966722B (zh) | 一种gnss钟差解算方法 | |
CN110749907A (zh) | 一种基于北斗动定位中接收机的钟差补偿方法及其系统 | |
CN111998849A (zh) | 一种基于惯性导航系统的差分动动定位的方法 | |
CN104316943A (zh) | 一种伪距离和多普勒组合差分定位系统及方法 | |
CN114935770B (zh) | 一种多历元加快精密单点定位收敛速度的方法及装置 | |
CN115421172B (zh) | 一种基于实时与准实时结合的北斗变形监测方法 | |
CN109143289B (zh) | 一种gnss单站位移监测方法 | |
CN112146557A (zh) | 一种基于gnss的实时桥梁变形监测系统及方法 | |
CN110568464B (zh) | 基于bds/gnss多模芯片的精密定位方法及装置 | |
CN110646823A (zh) | 一种基于Helmet后验定权法的GPS\BDS紧组合精密单点定位方法 | |
CN112630811B (zh) | 一种实时ppp-rtk组合定位方法 | |
RU2624268C1 (ru) | Способ определения взаимного положения объектов по сигналам глобальных навигационных спутниковых систем | |
Bisnath | Relative Positioning and Real‐Time Kinematic (RTK) | |
CN115220078A (zh) | 基于载波相位差分的gnss高精度定位方法及导航方法 |
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 | ||
TR01 | Transfer of patent right |
Effective date of registration: 20220819 Address after: No. 16, Jinniu District Road, Jinniu District Road, Chengdu, Sichuan Patentee after: CHINA RAILWAY NO.2 ENGINEERING GROUP Co.,Ltd. Patentee after: Sichuan Chengzheng Engineering Testing Technology Co.,Ltd. Address before: No. 16, Jinniu District Road, Jinniu District Road, Chengdu, Sichuan Patentee before: CHINA RAILWAY NO.2 ENGINEERING GROUP Co.,Ltd. |
|
TR01 | Transfer of patent right |