CN105486291A - 一种用于无验潮水深测量的动态精密单点定位方法 - Google Patents
一种用于无验潮水深测量的动态精密单点定位方法 Download PDFInfo
- Publication number
- CN105486291A CN105486291A CN201510815230.8A CN201510815230A CN105486291A CN 105486291 A CN105486291 A CN 105486291A CN 201510815230 A CN201510815230 A CN 201510815230A CN 105486291 A CN105486291 A CN 105486291A
- Authority
- CN
- China
- Prior art keywords
- elevation
- water
- depth
- receiver
- water level
- 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
Links
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01C—MEASURING DISTANCES, LEVELS OR BEARINGS; SURVEYING; NAVIGATION; GYROSCOPIC INSTRUMENTS; PHOTOGRAMMETRY OR VIDEOGRAMMETRY
- G01C13/00—Surveying specially adapted to open water, e.g. sea, lake, river or canal
- G01C13/008—Surveying specially adapted to open water, e.g. sea, lake, river or canal measuring depth of open water
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02A—TECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE
- Y02A90/00—Technologies having an indirect contribution to adaptation to climate change
- Y02A90/30—Assessment of water resources
Landscapes
- Life Sciences & Earth Sciences (AREA)
- Hydrology & Water Resources (AREA)
- Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Radar, Positioning & Navigation (AREA)
- Remote Sensing (AREA)
- Position Fixing By Use Of Radio Waves (AREA)
Abstract
本发明公开了一种用于无验潮水深测量的动态精密单点定位方法,涉及动态精密单点定位领域,该方法基于附有水位约束的动态精密单点定位技术,通过开始测试位置的水位值推算得到接收机天线高程,根据内河航道水位变化不大的特点,在观测方程中引入高程约束条件,采用附有高程约束的动态精密单点定位模型得到每个测点的平面和高程信息,进而计算各测点处的水位值,本发明提供的一种用于无验潮水深测量的动态精密单点定位方法,能够加快动态精密单点定位收敛的时间,提高工作效率。
Description
技术领域
本发明涉及动态精密单点定位领域,尤其涉及精密单点定位技术在无验潮水深测量中的应用,加快动态精密单点定位的收敛速度。
技术背景
卫星定位技术在海洋测绘领域已经得到了较为长远的发展,将卫星定位技术与水深测量紧密地结合起来,能够实现自动化的高精度航道水深测量。无验潮水深测量的基本原理是利用高精度的卫星定位技术得到天线中心高程,通过高程转换得到所需坐标系中的标高,根据测深仪数据等信息,最终计算得到水底泥面相对深度基准面的高程。该系统主要由接收机设备、测深系统、水上导航采集软件组成,其中卫星接收机设备的定位精度是测量系统成果的主要影响因素之一。
当前无验潮水深测量主要使用的定位技术为差分定位技术,但差分定位在水深测量中会受到测区离基准站的距离的限制,采用基于实时相位差分的作业距离一般不超过12km,采用事后相位差分的作业距离一般不超过50km,并且存在至少需要两台双频接收机,作业成本较高等缺点。将动态精密单点定位技术应用于无验潮水深测量,能够有效解决离基站较远无法进行定位的问题,它还具有不需要建立基准站、测得的点位结果之间不存在误差积累等优势,但传统的动态精密单点定位技术需要长时间的收敛过程。
发明内容
为了克服现有技术的不足,本发明目的是公开一种附有水位约束的动态精密单点定位新方法,能够较好地应用于内河航道的水深测量中,解决目前差分定位技术受距离限制、作业成本较高;以及传统动态精密单点定位技术收敛时间长的问题。
为达到上述目的,本发明提供如下技术方案:
一种用于无验潮水深测量的动态精密单点定位方法,该方法包括以下步骤:
步骤一:读取验潮站初始位置的水位值d水位,利用初始位置的水位值d水位与高程值H大地高转换公式H大地高=L+d水位+ζ得到接收机天线中心的高程值H大地高;L表示接收机天线到水面的高差;ζ为高程异常值,将理论深度基准面当成似大地水准面,理论深度基准面就是似大地水准面到参考椭球面的高程异常值ζ。
步骤二:从接收机得到的双频伪距和载波相位观测值,列出“消电离层组合”的动态精密单点定位模型:
式中,L1、L2分别是GPS信号调制波L1、L2载波上改正后的相位观测值(单位为m);
P1、P2分别是GPS信号调制波L1、L2载波上改正后的伪距观测值;
f1、f2分别是L1、L2载波的频率;λ1、λ2分别是L1、L2载波的波长;
N1、N2分别是L1、L2载波的整周模糊度;
ρ是卫星到接收机的几何距离;c是光速;
dtr、dts分别是接收机钟差和卫星钟差;
dtrop是对流层延迟项;
分别是“消电离层组合”伪距和相位观测值的残差值;
在动态精密单点定位模型中,每个历元的接收位置参数都是不断变化的,进行Kalman滤波时,设置噪声向量协方差阵的坐标相关参数为(100m)2。将动态精密单点定位模型进行线性化:
V=AX-L
X=[δx,δy,δz,dtr,δtrop,N1,…,Nn]T
式中,接收机r观测的任意一颗卫星i的系数矩阵为:
式中,xi、yi、zi表示卫星i的坐标,Xr、Yr、Zr表示接收机r的坐标,表示卫星到接收机的几何距离;Mi表示天顶方向上的对流层湿延迟投影函数;第一行为载波相位观测方程式系数,第二行为伪距观测方程式系数。
步骤三:在动态精密单点定位模型中加入高程约束条件:
V=RX-W
高程约束条件的系数阵为:
上式中,B0、L0分别表示接收机的概略纬度和经度;表示通过初始位置水位值d水位推算得到的天线中心的高程值H大地高;内河航道水面的短时间波动范围一般在20cm以内,在动态Kalman滤波模型中,将噪声向量协方差阵的高程约束相关参数设置为(0.2m)2;
步骤四:利用动态精密单点定位计算得到接收机天线的三维坐标;通过动态精密单点定位技术可以得到GPS天线中心的高程值H大地高,测量出GPS天线中心到换能器底部的高差L和吃水深度H吃水,结合测深仪得到的水深值H水深,可得到泥面相对于参考椭球面的高程为:
H泥面=H大地高-L-H吃水-H水深
泥面相对于理论深度基准面的高程为:
H=H理论-H泥面
=H理论-(H大地高-L-H吃水-H水深)
=(H理论-H大地高)+L+H吃水+H水深
最终得到深度基准面高程,即水深可以表示为:
H=L+H吃水+H水深-H大地高+ζ
本发明由于采用上述方法,采用精密单点定位技术,达到作业范围不受限制、减少作业设备的优势。采用上述方法后,因内河航道水位值在短时间内变化较小,根据开始测试时已知的水位值反推出接收机天线高程,作为定位高程约束条件,能够有效加快动态精密单点定位的收敛速度,实现了在高程方向的定位结果快速收敛,提高了精密单点定位技术在水深测量领域的应用的效率。
附图说明
图1是方法流程图;
图2是无验潮水深测量中的高程空间结构关系。
具体实施方式
下面将结合附图与具体实施方式对本发明的技术方案作进一步详细说明:
如图1所示为方法流程图。下面使用中海达V30接收机和HD-310测深仪在长江南京段上采集动态GPS/GLONASS双系统、双频伪距和载波观测数据,采样率为1s,卫星截止高度角为10°。采用本发明提出的新方法进行动态解算,并将最终定位结果与传统精密单点定位结果、以及验潮站水位进行比较。
首先从IGS网站下载所需时段的精密星历和精密钟差产品,结合记录的原始观测数据和初始水位值,采用自主编写的附有水位约束的动态精密单点定位软件,对数据进行预处理和各项误差计算后,计算得到每个历元的接收机天线坐标,最终得到水位成果,并与其他方式得到的结果进行对比分析。
如图2所示为无验潮水深测量中的高程空间结构关系。通过卫星定位技术得到高精度的接收机天线的三维坐标,该高程值属于大地高系统,进行一系列高程转换后,可直接得到最终水位结果。反之,可根据水位信息反算得到天线中心的高程值,下面列出该转换方法的计算过程:
H大地高=L+d水位+H理论
在实际工程中,将理论深度基准面当成似大地水准面,所以H理论就是似大地水准面到参考椭球面的高程值,即为高程异常值ζ。则上式可以改写为:
H大地高=L+d水位+ζ
上式中,H大地高表示接收机天线的大地高;L表示接收机天线到水面的高差;d水位表示验潮站水位观测值。
(2)精密单点定位技术无法通过差分对数据进行精确的预处理,再通过相关模型计算出各项误差改正值,并对存在问题的卫星进行标记。数据处理完成后,根据接收机获取的双频伪距和载波相位观测值,经过误差改正后,列出“消电离层组合”的动态精密单点定位模型:
式中,L1、L2分别是GPS信号调制波L1、L2载波上改正后的相位观测值(单位为m);
P1、P2分别是GPS信号调制波L1、L2载波上改正后的伪距观测值;
f1、f2分别是L1、L2载波的频率;λ1、λ2分别是L1、L2载波的波长;
N1、N2分别是L1、L2载波的整周模糊度;
ρ是卫星到接收机的几何距离;c是光速;
dtr、dts分别是接收机钟差和卫星钟差;
dtrop是对流层延迟项;
分别是“消电离层组合”伪距和相位观测值的残差值;
在动态精密单点定位模型中,每个历元的接收位置参数都是不断变化的,进行Kalman滤波时,设置噪声向量协方差阵的坐标相关参数为(100m)2。将动态精密单点定位模型进行线性化:
V=AX-L
X=[δx,δy,δz,dtr,δtrop,N1,…,Nn]T
式中,接收机r观测的任意一颗卫星i的系数矩阵为:
式中,xi、yi、zi表示卫星i的坐标,Xr、Yr、Zr表示接收机r的坐标,表示卫星到接收机的几何距离;Mi表示天顶方向上的对流层湿延迟投影函数;第一行为载波相位观测方程式系数,第二行为伪距观测方程式系数。
添加高程约束条件:
V=RX-W
高程约束条件的系数阵为:
式中,B0、L0分别表示接收机的纬度和经度;表示通过初始位置水位值推算得到的天线中心高程值H大地高。内河的水面波动范围一般在20cm以内,在Kalman滤波的动态模型中,将噪声向量协方差阵的高程约束相关参数设置为(0.2m)2。
利用Kalman滤波即可求得接收机位置、钟差、天顶对流层延迟以及模糊度参数,为水深测量提供高精度的测点坐标。
(3)最终,利用动态精密单点定位计算得到接收机天线的三维坐标,能够进一步得到成果水深。具体过程如下:
通过动态精密单点定位技术可以得到GPS天线中心的高程值H大地高,测量出GPS天线中心到换能器底部的高差L和吃水深度H吃水,结合测深仪得到的水深值H水深,可得到泥面相对于参考椭球面的高程为:
H泥面=H大地高-L-H吃水-H水深
泥面相对于理论深度基准面的高程为:
H=H理论-H泥面
=H理论-(H大地高-L-H吃水-H水深)
=(H理论-H大地高)+L+H吃水+H水深
最终得到深度基准面高程,即成果水深可以表示为:
H=L+H吃水+H水深-H大地高+ζ
Claims (1)
1.一种用于无验潮水深测量的动态精密单点定位方法,其特征在于该方法包括以下步骤:
步骤一:读取验潮站初始位置的水位值d水位,利用初始位置的水位值d水位与高程值H大地高转换公式H大地高=L+d水位+ζ得到接收机天线中心的高程值H大地高;L表示接收机天线到水面的高差;ζ为高程异常值,将理论深度基准面当成似大地水准面,理论深度基准面就是似大地水准面到参考椭球面的高程异常值ζ。
步骤二:从接收机得到的双频伪距和载波相位观测值,列出“消电离层组合”的动态精密单点定位模型:
式中,L1、L2分别是GPS信号调制波L1、L2载波上改正后的相位观测值(单位为m);
P1、P2分别是GPS信号调制波L1、L2载波上改正后的伪距观测值;
f1、f2分别是L1、L2载波的频率;λ1、λ2分别是L1、L2载波的波长;
N1、N2分别是L1、L2载波的整周模糊度;
ρ是卫星到接收机的几何距离;c是光速;
dtr、dts分别是接收机钟差和卫星钟差;
dtrop是对流层延迟项;
分别是“消电离层组合”伪距和相位观测值的残差值;
在动态精密单点定位模型中,每个历元的接收位置参数都是不断变化的,进行Kalman滤波时,设置噪声向量协方差阵的坐标相关参数为(100m)2。将动态精密单点定位模型进行线性化:
V=AX-L
X=[δx,δy,δz,dtr,δtrop,N1,…,Nn]T
式中,接收机r观测的任意一颗卫星i的系数矩阵为:
式中,xi、yi、zi表示卫星i的坐标,Xr、Yr、Zr表示接收机r的坐标,表示卫星到接收机的几何距离;Mi表示天顶方向上的对流层湿延迟投影函数;第一行为载波相位观测方程式系数,第二行为伪距观测方程式系数。
步骤三:在动态精密单点定位模型中加入高程约束条件:
V=RX-W
高程约束条件的系数阵为:
上式中,BO、LO分别表示接收机的概略纬度和经度;表示通过初始位置水位值d水位推算得到的天线中心的高程值H大地高;内河航道水面的短时间波动范围一般在20cm以内,在动态Kalman滤波模型中,将噪声向量协方差阵的高程约束相关参数设置为(0.2m)2;
步骤四:利用动态精密单点定位计算得到接收机天线的三维坐标;通过动态精密单点定位技术可以得到GPS天线中心的高程值H大地高,测量出GPS天线中心到换能器底部的高差L和吃水深度H吃水,结合测深仪得到的水深值H水深,可得到泥面相对于参考椭球面的高程为:
H泥面=H大地高-L-H吃水-H水深
泥面相对于理论深度基准面的高程为:
H=H理论-H泥面
=H理论-(H大地高-L-H吃水-H水深)
=(H理论-H大地高)+L+H吃水+H水深
最终得到深度基准面高程,即水深可以表示为:
H=L+H吃水+H水深-H大地高+ζ。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510815230.8A CN105486291A (zh) | 2015-11-23 | 2015-11-23 | 一种用于无验潮水深测量的动态精密单点定位方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510815230.8A CN105486291A (zh) | 2015-11-23 | 2015-11-23 | 一种用于无验潮水深测量的动态精密单点定位方法 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN105486291A true CN105486291A (zh) | 2016-04-13 |
Family
ID=55673429
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201510815230.8A Pending CN105486291A (zh) | 2015-11-23 | 2015-11-23 | 一种用于无验潮水深测量的动态精密单点定位方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN105486291A (zh) |
Cited By (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106403908A (zh) * | 2016-08-29 | 2017-02-15 | 上海交通大学 | 基于时间序列的水深预测方法和系统 |
CN106802148A (zh) * | 2017-03-01 | 2017-06-06 | 山东科技大学 | 一种水深测量深度基准面的传递方法 |
CN108427133A (zh) * | 2018-03-20 | 2018-08-21 | 交通运输部北海航海保障中心天津海事测绘中心 | 一种基于北斗和声测深技术的验潮仪系统 |
CN108489497A (zh) * | 2018-05-22 | 2018-09-04 | 何竹君 | 一种利用地图防触礁的安全助航方法 |
CN109917494A (zh) * | 2019-02-13 | 2019-06-21 | 上海华测导航技术股份有限公司 | 降雨预报方法、装置、设备和存储介质 |
CN110109163A (zh) * | 2019-04-09 | 2019-08-09 | 天津大学 | 一种附有高程约束的精密单点定位方法 |
CN110109167A (zh) * | 2019-04-18 | 2019-08-09 | 天津大学 | 一种基于高程约束的海上精密定位方法 |
CN111189435A (zh) * | 2020-02-28 | 2020-05-22 | 国家海洋技术中心 | 一种臂挂拖曳式水深测量系统及方法 |
CN111536951A (zh) * | 2020-05-29 | 2020-08-14 | 天津水运工程勘察设计院 | 一种实时动态水深测量系统 |
CN111674522A (zh) * | 2020-05-14 | 2020-09-18 | 中交广州航道局有限公司 | 一种水位测量系统以及一种施工船舶 |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JPH0942960A (ja) * | 1995-08-02 | 1997-02-14 | Taisei Corp | 水中点の高さの測定方法 |
CN104714244A (zh) * | 2015-03-31 | 2015-06-17 | 东南大学 | 一种基于抗差自适应Kalman滤波的多系统动态PPP解算方法 |
CN104729486A (zh) * | 2015-03-16 | 2015-06-24 | 东南大学 | 一种基于似大地水准面精化的无验潮水深测量方法 |
-
2015
- 2015-11-23 CN CN201510815230.8A patent/CN105486291A/zh active Pending
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JPH0942960A (ja) * | 1995-08-02 | 1997-02-14 | Taisei Corp | 水中点の高さの測定方法 |
CN104729486A (zh) * | 2015-03-16 | 2015-06-24 | 东南大学 | 一种基于似大地水准面精化的无验潮水深测量方法 |
CN104714244A (zh) * | 2015-03-31 | 2015-06-17 | 东南大学 | 一种基于抗差自适应Kalman滤波的多系统动态PPP解算方法 |
Non-Patent Citations (1)
Title |
---|
林瑜滢等: "动态精密单点定位传统模型与UofC模型的比较", 《测绘科学技术学报》 * |
Cited By (15)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106403908B (zh) * | 2016-08-29 | 2019-10-08 | 上海交通大学 | 基于时间序列的水深预测方法和系统 |
CN106403908A (zh) * | 2016-08-29 | 2017-02-15 | 上海交通大学 | 基于时间序列的水深预测方法和系统 |
CN106802148A (zh) * | 2017-03-01 | 2017-06-06 | 山东科技大学 | 一种水深测量深度基准面的传递方法 |
CN106802148B (zh) * | 2017-03-01 | 2018-03-06 | 山东科技大学 | 一种水深测量深度基准面的传递方法 |
CN108427133A (zh) * | 2018-03-20 | 2018-08-21 | 交通运输部北海航海保障中心天津海事测绘中心 | 一种基于北斗和声测深技术的验潮仪系统 |
CN108489497A (zh) * | 2018-05-22 | 2018-09-04 | 何竹君 | 一种利用地图防触礁的安全助航方法 |
CN108489497B (zh) * | 2018-05-22 | 2021-07-30 | 何竹君 | 一种利用地图防触礁的安全助航方法 |
CN109917494A (zh) * | 2019-02-13 | 2019-06-21 | 上海华测导航技术股份有限公司 | 降雨预报方法、装置、设备和存储介质 |
CN110109163A (zh) * | 2019-04-09 | 2019-08-09 | 天津大学 | 一种附有高程约束的精密单点定位方法 |
CN110109163B (zh) * | 2019-04-09 | 2022-01-07 | 天津大学 | 一种附有高程约束的精密单点定位方法 |
CN110109167A (zh) * | 2019-04-18 | 2019-08-09 | 天津大学 | 一种基于高程约束的海上精密定位方法 |
CN110109167B (zh) * | 2019-04-18 | 2021-12-07 | 天津大学 | 一种基于高程约束的海上精密定位方法 |
CN111189435A (zh) * | 2020-02-28 | 2020-05-22 | 国家海洋技术中心 | 一种臂挂拖曳式水深测量系统及方法 |
CN111674522A (zh) * | 2020-05-14 | 2020-09-18 | 中交广州航道局有限公司 | 一种水位测量系统以及一种施工船舶 |
CN111536951A (zh) * | 2020-05-29 | 2020-08-14 | 天津水运工程勘察设计院 | 一种实时动态水深测量系统 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN105486291A (zh) | 一种用于无验潮水深测量的动态精密单点定位方法 | |
CN106291639B (zh) | 一种gnss接收机实现定位的方法及装置 | |
CN104102822B (zh) | 一种多频gnss观测值随机特性建模方法 | |
CN108802782B (zh) | 一种惯导辅助的北斗三频载波相位整周模糊度求解方法 | |
CN103868493B (zh) | 一种基于ppp技术的深度基准面大地高测量方法 | |
CN104483691B (zh) | 一种gnss组合精密单点定位方法 | |
CN105182384A (zh) | 一种双模实时伪距差分定位系统和伪距改正数据生成方法 | |
CN104133231B (zh) | 一种基于积分多普勒平滑伪距的导航定位方法 | |
CN104597471A (zh) | 面向时钟同步多天线gnss接收机的定向测姿方法 | |
CN107193023B (zh) | 一种具有闭式解的高精度北斗卫星系统单点定位方法 | |
CN103529482B (zh) | 一种高精度确定载体动态加速度的方法 | |
He | GNSS kinematic position and velocity determination for airborne gravimetry | |
CN105005019A (zh) | 一种基于bds/gps rtk的地下管线测量方法 | |
CN107544075A (zh) | 基于精密单点定位与调和分析估计海潮负荷位移参数方法 | |
CN103852799A (zh) | 一种基于ppp技术的船载重力测量方法 | |
Seepersad | Reduction of initial convergence period in GPS PPP data processing | |
CN116009042A (zh) | 一种单站载波历元间差分实时探测相对形变的方法及系统 | |
CN103529451A (zh) | 一种水面母船校准海底应答器坐标位置的方法 | |
Liu et al. | Real-time precise measurements of ocean surface waves using GNSS variometric approach | |
Chen et al. | Feasibility and performance evaluation of low-cost GNSS devices for sea level measurement based on GNSS-IR | |
CN105158781A (zh) | 一种地球gnss卫星实现外层空间用户定位的方法 | |
Wu et al. | MEMS-based IMU assisted real time difference using raw measurements from smartphone | |
CN109977499B (zh) | 基于位置约束的北斗三频/静力水准仪斜拉桥监测方法 | |
CN102830410B (zh) | 卫星导航中结合多普勒测速的定位方法 | |
CN103389502B (zh) | 基于多个地面基站高精度确定载体加速度的方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
WD01 | Invention patent application deemed withdrawn after publication | ||
WD01 | Invention patent application deemed withdrawn after publication |
Application publication date: 20160413 |