CN103727920B - 基于大地水准面模型测量水准高差的方法 - Google Patents
基于大地水准面模型测量水准高差的方法 Download PDFInfo
- Publication number
- CN103727920B CN103727920B CN201310738247.9A CN201310738247A CN103727920B CN 103727920 B CN103727920 B CN 103727920B CN 201310738247 A CN201310738247 A CN 201310738247A CN 103727920 B CN103727920 B CN 103727920B
- Authority
- CN
- China
- Prior art keywords
- control points
- elevation
- geoid
- formula
- point
- 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
- 238000000034 method Methods 0.000 title claims abstract description 32
- 238000005259 measurement Methods 0.000 claims abstract description 15
- 238000005516 engineering process Methods 0.000 claims description 7
- 230000003068 static effect Effects 0.000 claims description 3
- 230000005484 gravity Effects 0.000 description 22
- 230000009897 systematic effect Effects 0.000 description 8
- 230000007812 deficiency Effects 0.000 description 3
- 238000004364 calculation method Methods 0.000 description 2
- 238000010586 diagram Methods 0.000 description 2
- 238000013507 mapping Methods 0.000 description 2
- 230000005856 abnormality Effects 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 238000005290 field theory Methods 0.000 description 1
- 238000004519 manufacturing process Methods 0.000 description 1
- 238000003672 processing method Methods 0.000 description 1
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 description 1
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01C—MEASURING DISTANCES, LEVELS OR BEARINGS; SURVEYING; NAVIGATION; GYROSCOPIC INSTRUMENTS; PHOTOGRAMMETRY OR VIDEOGRAMMETRY
- G01C5/00—Measuring height; Measuring distances transverse to line of sight; Levelling between separated points; Surveyors' levels
Landscapes
- 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
本发明公开了一种基于大地水准面模型测量水准高差的方法,它采用大地水准面模型计算出相邻控制点的高程异常差,通过GNSS测量得到相邻控制点的大地高差,二者相减,即可转换为水准高差,进而计算各个控制点的高程。本发明可应用于多种地形,地形包括平原、丘陵、山区等,其各个控制点的间距可达5~8km,工作效率可提高6至10倍,可达到三、四等水准测量的精度;且只需首末端两个控制点的高程,不需要进行几何水准测量,具有工作强度小、效率高的特点。
Description
技术领域
本发明涉及测绘技术领域,具体地指一种基于大地水准面模型测量水准高差的方法。
背景技术
现有技术中,水准高差的测量一般采用几何水准的方法进行;几何水准是指用水准仪和水准尺测定地面上两点间高差的方法。在地面两点间安置水准仪,观测竖立在两点上的水准标尺,按尺上读数推算两点间的高差。通常由水准原点或任一已知高程点出发,沿选定的水准路线逐站测定各点的高程,一站的距离为20~200m。
采用几何水准的方法进行水准高差的测量,其工作强度大、效率低,严重制约了测绘保障。
也可以采用重力似大地水准面精化方法,获得各个点的高程异常,但其有以下缺点:
1)重力似大地水准面及其存在的不足:
重力似大地水准面是利用重力测量数据、数字地形模型和参考重力位模型,应用移去-恢复技术,仅采用重力方法计算获得的似大地水准面模型,因为重力和高程数据分辨率较高,因此它具有高分辨率的特点;
但是,重力似大地水准面也因为各种因素的影响,使其计算精度较低:一是部分地区地面重力数据稀少甚至是空白;二是重力资料分布不规则且我国地形变化剧烈,在进行重力数据格网化时必定会带来一定插值误差;三是参考重力位模型与我国重力场存在系统差,位模型与重力数据的长波误差也不尽一致,这些因素都将对计算得到的重力似大地水准面产生系统性影响;四是重力似大地水准面与GPS/水准似大地水准面之间存在较大的系统偏差,不利于GPS测量的应用;
2)重力似大地水准面与GPS/水准的拟合及其存在的不足:
用重力场理论确定的高分辨率的重力似大地水准面受地面重力数据误差及其分布密度等影响而含有系统偏差,GPS/水准的精度虽然很高,但受水准测量方法的限制,难以实现在时间和空间上达到一定的密度要求,因此,如何将这两类似大地水准面综合得到一个消除系统误差的高精度的区域似大地水准面有不同的处理方法,通过采用多项式拟合多面函数拟合等方法,将高分辨率但精度较低的重力似大地水准面拟合到高精度但数量较少的GPS(似)大地水准面上,从而消除重力似大地水准面的系统偏差,重力似大地水准面与GPS/水准似大地水准面的拟合,能够较好的消除两者之间的系统偏差,提高似大地水准面的精度,但是该方法仍然存在一些不足:
(1)重力似大地水准面与GPS/水准似大地水准面之间的系统偏差非常复杂,是不可能用某种数学曲面代替的,这种曲面拟合显然误差较大;
(2)因为GPS/水准点一般比较稀疏,用于曲面拟合的数据较少,不能精细表现复杂的似大地水准面,难以达到很好的拟合效果,而要保证拟合精度,就必须布测较密集的GPS/水准点,这样势必造成过多地增加水准测量工作量,增大了外业劳动强度,加大了项目的周期和生产成本;
(3)拟合曲面是一种平滑的数学曲面,会导致在拟合过程中损失部分高频信息。
发明内容
本发明的目的就是要解决上述背景技术的不足,提供一种可提高测量工作效率且精度高的基于大地水准面模型测量水准高差的方法。
本发明的技术方案为:一种基于大地水准面模型测量水准高差的方法,其特征在于,它包括以下步骤:
1)、在待测量区域选取多个控制点,多个控制点连成多个闭合环,其中首端闭合环的首端控制点与末端闭合环的末端控制点为已知高程的控制点;
2)、利用EGM2008全球大地水准面模型或省域似大地水准面模型,计算各个控制点的高程异常;
3)、将相邻控制点的高程异常相减,即可得到相邻控制点的高程异常差;
4)、通过GNSS测量技术,得到各个控制点的大地高,将相邻控制点的大地高相减,即可得到相邻控制点的大地高差;
5)、将上述得到的相邻控制点的大地高差减去高程异常差,即得到相邻控制点的水准高差。
步骤2)中,采用邻域9点法计算各个控制点的高程异常。
采用邻域9点法计算各个控制点的高程异常的步骤为:
a)、选取任一控制点,以该控制点为中心,选取该控制点邻域方格的9个已知高程异常的数据点,构成核函数:
Qj=[(X-Xj)2+(Y-Yj)2+0.8]0.5(式1)
其中,(Xj,Yj)为第j个已知高程异常的数据点在EGM2008全球大地水准面模型或省域似大地水准面模型的坐标,(X,Y)为9个已知高程异常的数据点在EGM2008全球大地水准面模型或省域似大地水准面模型对应的坐标;
b)、当j为1时,将(X,Y)代入到式1中,可得到Q11、Q12…Q19;
当j为2时,将(X,Y)代入到式1中,可得到Q21、Q22…Q29;
……
当j为9时,将(X,Y)代入到式1中,可得到Q91、Q92…Q99;
c)、利用式2得到未知参数K值,其中式2中z为已知高程异常的数据点对应的高程异常;
d)、解得式2中的K值,并利用式3得到ZP的值,ZP即为待求控制点的高程异常;类似a)、b)两个步骤,计算QP1、QP2…QP9;
进一步地,所述控制点的选取满足GNSS静态定位测量要求,控制点之间的间距为5~8km。
进一步地,所述构成闭合环的控制点数量为3至6个。
本发明提出了采用大地水准面模型计算出相邻控制点的高程异常差,通过GNSS测量得到相邻控制点的大地高差,二者相减,即可转换为水准高差,进而计算各个控制点的高程。
本发明的各个控制点的间距可达5~8km,工作效率可提高6至10倍,可达到三、四等水准测量的精度;且只需首末端两个控制点的高程,不需要进行几何水准测量,具有工作强度小、效率高的特点。
附图说明
图1为本发明中控制点的选取示意图;
图2为本发明中控制点邻域方格的已知高程异常的数据点选取示意图。
具体实施方式
下面结合附图和具体实施例对本发明作进一步的详细说明。
本发明所描述的一种基于大地水准面模型测量水准高差的方法,其特征在于,它包括以下步骤:
1)、在待测量区域选取多个控制点,多个控制点连成多个闭合环,其中首端闭合环的首端控制点与末端闭合环的末端控制点为已知高程的控制点;
如图1所示,例如在此待测区域选取10个控制点,其中BM01\BM02为已知高程的控制点,P01\P02……P07\P08为未知高程的控制点。
另外,BM01-P01-P02构成闭合环,其余类似。构成闭合环的控制点数量可以不少于3个,但应不多于6个;
2)、利用EGM2008全球大地水准面模型或省域似大地水准面模型,采用邻域9点法计算各个控制点的高程异常;
3)、将相邻控制点的高程异常相减,即可得到相邻控制点的高程异常差;
此步骤中:如果直接采用高程异常加上大地高的方法来计算高程,不同高程系统的偏差,会导致得到的高程误差较大,无法满足等级水准的要求。而高程异常差既不包含不同高程系统的偏差(一般偏差可达到0.2m~0.3m),且由于大地水准面模型的相关性,模型的相关性误差经差分也已经消除,故得到的“高程异常差”的精度较高;
4)、通过GNSS测量技术,得到各个控制点的大地高,将相邻控制点的大地高相减,即可得到相邻控制点的大地高差;
5)、将上述得到的相邻控制点的大地高差减去高程异常差,即得到相邻控制点的水准高差。
本实施例中:采用邻域9点法计算各个控制点的高程异常的步骤为:
a)、选取任一控制点,以该控制点为中心,选取该控制点邻域方格的9个已知高程异常的数据点,构成核函数:
Qj=[(X-Xj)2+(Y-Yj)2+0.8]0.5(式1)
其中,(Xj,Yj)为第j个已知高程异常的数据点在EGM2008全球大地水准面模型或省域似大地水准面模型的坐标,(X,Y)为9个已知高程异常的数据点在EGM2008全球大地水准面模型或省域似大地水准面模型对应的坐标;
b)、当j为1时,将(X,Y)代入到式1中,可得到Q11、Q12…Q19;
当j为2时,将(X,Y)代入到式1中,可得到Q21、Q22…Q29;
……
当j为9时,将(X,Y)代入到式1中,可得到Q91、Q92…Q99;
c)、利用式2得到未知参数K值,其中式2中z为已知高程异常的数据点对应的高程异常;类似a)、b)两个步骤,计算QP1、QP2…QP9;
d)、解得式2中的K值,并利用式3得到ZP的值,ZP即为待求控制
点的高程异常,类似a)、b)两个步骤,计算QP1、QP2…QP9;
其中ZP为GNSS控制点的高程异常。
另外,在利用GNSS测量技术中,各个控制点的观测时间不小于60分钟,并构成GNSS水准线路;解算GNSS基线后,进行GNSS控制网平差,得到各个控制点的大地高,将相邻控制点的大地高相减,计算相邻控制点的“大地高差”。
利用上述方法计算得到的水准高差,可采用条件平差方法,根据高等级水准点的高程、水准高差、间距,可计算高程闭合差,并得到各个控制点的平差后高程。
本发明提出了采用大地水准面模型计算出相邻控制点的高程异常差,通过GNSS测量得到相邻控制点的大地高差,二者相减,即可转换为水准高差,进而计算各个控制点的高程。
本发明可应用于多种地形,地形包括平原、丘陵、山区等,其各个控制点的间距可达5~8km,工作效率可提高6至10倍,可达到三、四等水准测量的精度;且只需两个控制点的高程,不需要进行几何水准测量,具有工作强度小、效率高的特点。
实施例:某项水利工程,水准线路长度为75.8km,设置有13个控制点,其中HBJT02、HBJT24为已知高程点,采用GNSS静态定位测量联测各个控制点,并采用EGM2008全球大地水准面模型按照上述方式计算高程异常,按上述方法计算转换的水准高差,为了验证本方法的精度,以三等水准联测各个控制点,转换的水准高差与三等水准高差比较,从下表可见,高差较差均在限差之内。
表1
表2
高程闭合差/m | 0.0405 |
四等限差/m | 0.174 |
表3
经统计,在一定距离(不大于10km)内,本实施例得到的高程异常差的精度较高(可达到0.01m~0.02m),而大地高差的精度可达到0.02m~0.03m,根据误差传播律,则获得的水准高差精度为0.022m~0.036m。
采用本发明实施方式比较简便易行,不需要进行大量的重力测量、水准测量,以及重力似大地水准面确定与GPS/水准的拟合等复杂计算,可代替传统几何水准测量的技术,具有重要的应用价值。
Claims (5)
1.一种基于大地水准面模型测量水准高差的方法,其特征在于,它包括以下步骤:
1)、在待测量区域选取多个控制点,多个控制点连成多个闭合环,其中首端闭合环的首端控制点与末端闭合环的末端控制点为已知高程的控制点;
2)、利用EGM2008全球大地水准面模型或省域似大地水准面模型,计算各个控制点的高程异常;
3)、将相邻控制点的高程异常相减,即可得到相邻控制点的高程异常差;
4)、通过GNSS测量技术,得到各个控制点的大地高,将相邻控制点的大地高相减,即可得到相邻控制点的大地高差;
5)、将上述得到的相邻控制点的大地高差减去高程异常差,即得到相邻控制点的水准高差。
2.根据权利要求1所述的基于大地水准面模型测量水准高差的方法,其特征在于:步骤2)中,采用邻域9点法计算各个控制点的高程异常。
3.根据权利要求2所述的基于大地水准面模型测量水准高差的方法,其特征在于,采用邻域9点法计算各个控制点的高程异常的步骤为:
a)、选取任一控制点,以该控制点为中心,选取该控制点邻域方格的9个已知高程异常的数据点,构成核函数:
Qj=[(X-Xj)2+(Y-Yj)2+0.8]0.5(式1)
其中,(Xj,Yj)为第j个已知高程异常的数据点在EGM2008全球大地水准面模型或省域似大地水准面模型的坐标,(X,Y)为9个已知高程异常的数据点在EGM2008全球大地水准面模型或省域似大地水准面模型对应的坐标;
b)、当j为1时,将(X,Y)代入到式1中,可得到Q11、Q12…Q19;
当j为2时,将(X,Y)代入到式1中,可得到Q21、Q22…Q29;
……
当j为9时,将(X,Y)代入到式1中,可得到Q91、Q92…Q99;
c)、利用式2得到未知参数K值,其中式2中z为已知高程异常的数据点对应的高程异常,类似a)、b)两个步骤,计算QP1、QP2…QP9;
d)、解得式2中的K值,并利用式3得到ZP的值,ZP即为待求控制点的高程异常;
4.根据权利要求1所述的基于大地水准面模型测量水准高差的方法,其特征在于:所述控制点的选取满足GNSS静态定位测量要求,控制点之间的间距为5~8km。
5.根据权利要求1所述的基于大地水准面模型测量水准高差的方法,其特征在于:构成所述闭合环的控制点数量为3至6个。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201310738247.9A CN103727920B (zh) | 2013-12-27 | 2013-12-27 | 基于大地水准面模型测量水准高差的方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201310738247.9A CN103727920B (zh) | 2013-12-27 | 2013-12-27 | 基于大地水准面模型测量水准高差的方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN103727920A CN103727920A (zh) | 2014-04-16 |
CN103727920B true CN103727920B (zh) | 2016-03-16 |
Family
ID=50452100
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201310738247.9A Active CN103727920B (zh) | 2013-12-27 | 2013-12-27 | 基于大地水准面模型测量水准高差的方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN103727920B (zh) |
Families Citing this family (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104573396A (zh) * | 2015-01-27 | 2015-04-29 | 国家测绘地理信息局大地测量数据处理中心 | 水准网独立最小闭合环、最短附合线生成算法 |
CN105157672B (zh) * | 2015-08-24 | 2017-08-01 | 国家电网公司 | 一种用于杆塔间的高差测量仪 |
CN105953772A (zh) * | 2016-04-22 | 2016-09-21 | 深圳市华信天线技术有限公司 | 一种高程异常的测定方法及测定装置 |
CN108761510B (zh) * | 2018-04-09 | 2020-07-24 | 中铁二局集团有限公司 | 利用基于地形改正的重力场模型进行水准高差测量的方法 |
CN108917706B (zh) * | 2018-06-26 | 2020-10-16 | 中铁大桥勘测设计院集团有限公司 | 一种gnss水准测量方法 |
CN111895980A (zh) * | 2020-06-23 | 2020-11-06 | 中铁第一勘察设计院集团有限公司 | 基于北斗gnss的铁路控制测量控制点装置及控制方法 |
CN117367367A (zh) * | 2023-12-04 | 2024-01-09 | 中国铁路设计集团有限公司 | 基于北斗/gnss的超长距离跨河、跨海高程传递测量方法 |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101839710A (zh) * | 2010-06-12 | 2010-09-22 | 中国测绘科学研究院 | 一种似大地水准面计算的优化方法 |
CN102175213A (zh) * | 2011-01-30 | 2011-09-07 | 中国测绘科学研究院 | Gps水准高程传递方法 |
CN102589517A (zh) * | 2012-01-12 | 2012-07-18 | 东南大学 | 基于egm2008的区域似大地水准面精化方法 |
-
2013
- 2013-12-27 CN CN201310738247.9A patent/CN103727920B/zh active Active
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101839710A (zh) * | 2010-06-12 | 2010-09-22 | 中国测绘科学研究院 | 一种似大地水准面计算的优化方法 |
CN102175213A (zh) * | 2011-01-30 | 2011-09-07 | 中国测绘科学研究院 | Gps水准高程传递方法 |
CN102589517A (zh) * | 2012-01-12 | 2012-07-18 | 东南大学 | 基于egm2008的区域似大地水准面精化方法 |
Non-Patent Citations (1)
Title |
---|
"GPS水准及其在测绘工程中的应用";邸国辉等;《地理空间信息》;20060228;第4卷(第1期);第6-8页 * |
Also Published As
Publication number | Publication date |
---|---|
CN103727920A (zh) | 2014-04-16 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN103727920B (zh) | 基于大地水准面模型测量水准高差的方法 | |
CN105629263B (zh) | 一种对流层大气延迟误差估计改正方法和改正系统 | |
CN104091064B (zh) | 基于优化解空间搜索法的PS‑DInSAR地表形变测量参数估计方法 | |
CN104330073B (zh) | 一种烟囱倾斜度的测量方法 | |
CN103217688B (zh) | 一种基于不规则三角网机载激光雷达点云平差计算方法 | |
CN105182366A (zh) | 一种基于实测气象参数的对流层天顶延迟改正方法 | |
CN101839710B (zh) | 一种似大地水准面的确定方法 | |
CN102175213B (zh) | Gps水准高程传递方法 | |
CN103196425B (zh) | 特长隧道横向贯通误差的估测方法 | |
CN101825662B (zh) | 短距测量接地网接地电阻的测量方法及其测量结构 | |
CN104035138A (zh) | 一种全球及局部海洋扰动重力的精确快速计算方法 | |
CN104613932A (zh) | 利用垂线偏差与重力异常确定似大地水准面模型的方法 | |
CN102589517A (zh) | 基于egm2008的区域似大地水准面精化方法 | |
CN104776827B (zh) | Gps高程异常数据的粗差探测方法 | |
CN106885586A (zh) | 无需专门定标场的卫星雷达高度计绝对定标方法 | |
CN105180811A (zh) | 基于同名特征地物的移动测量系统激光扫描仪标定方法 | |
CN103293554A (zh) | 一种基于地震数据的静校正方法 | |
CN104111061A (zh) | 一种基础资料缺乏地区测量控制点成果获取方法 | |
CN104360396B (zh) | 一种海上井间tti介质三种初至波走时层析成像方法 | |
CN101957193B (zh) | 一种海岛礁高程传递的优化方法 | |
CN103837113A (zh) | 一种相对坐标测量方法 | |
CN105573963B (zh) | 一种电离层水平不均匀结构重构方法 | |
CN104181571A (zh) | 实现弱或无cors信号区地面点的精密坐标与高程快速测量的方法 | |
CN105866386A (zh) | 一种利用军事地图测算土壤流失量的方法 | |
CN106815607B (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 | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant |