CN101957193A - 一种海岛礁高程传递的优化方法 - Google Patents
一种海岛礁高程传递的优化方法 Download PDFInfo
- Publication number
- CN101957193A CN101957193A CN 201010198624 CN201010198624A CN101957193A CN 101957193 A CN101957193 A CN 101957193A CN 201010198624 CN201010198624 CN 201010198624 CN 201010198624 A CN201010198624 A CN 201010198624A CN 101957193 A CN101957193 A CN 101957193A
- Authority
- CN
- China
- Prior art keywords
- gravity
- data
- height
- land
- difference
- 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
Images
Landscapes
- Position Fixing By Use Of Radio Waves (AREA)
- Testing Or Calibration Of Command Recording Devices (AREA)
Abstract
一种海岛礁高程传递的优化方法,属于大地测量学与测量工程技术领域的发明专利。水准测量是目前高程传递最直接有效精度最高并且使用最为广泛的方法,然后在一些特殊区域却影响了其使用范围,例如在海岛礁高程传递中,水准测量因在海面上无法稳定使用水准测量设备而使得利用水准仪直接测量高程的方法不可用。基于此,本发明提出了一种海岛礁高程传递的优化方法,该方法通过陆地、海岸带、海岛重力数据的高度集成,通过高程基准系统差改正的空间重力异常数据计算陆地与所传递海岛礁一致的高精度整体重力似大地水准面数值模型,最后利用计算区域内的已知高程点,通过测量已知高程的GPS大地高和待测海岛礁点的GPS大地高,完成海岛礁的高程传递。
Description
技术领域
本发明是一种海岛礁高程传递的优化方法,属于“大地测量学与测量工程”学科中的“物理大地测量学”技术领域。
背景技术
目前,使用最多精度最高,并且使用最为广泛的高程测量方法为水准测量,通常规定一、二等水准测量需要进行精密水准测量的同时需在相应的水准路线上测量重力数据,而三、四等以及等外水准测量则无需重力测量。我国的高程基准网就是基于水准测量建立起来的。事实上,由于水准测量固有的局限性使得在跨海高程传递中不能发挥作用。目前,在近距离跨海高程传递中(一般小于10km)使用三角高程传递也能达到较高的精度效果,而在更远的距离的跨海高程传递中由于大气折光等多种因素影响,使得结果精度无法保证。基于(似)大地水准面的高程传递原理是利用陆地和海洋重力异常数据,通过物理大地测量原理计算陆地与所传递海岛礁一致的高精度整体重力似大地水准面数值模型,然后利用计算区域内的已知高程点,通过测量已知高程的GPS大地高和待测海岛礁点的GPS大地高,完成海岛礁的高程传递。但目前该方法存在两个问题。第一,卫星测高数据可以认为是全球基准一致的数据,而重力异常数据由于使用了局部高程基准而使得其与全球数据存在系统偏差;第二,在海岸带区域,由于卫星测高数据精度较低,并且可用数据较少,因此需要利用所有一切可利用资源才能满足精度要求。
本发明中,充分利用方法(利用GPS水准数据和重力场模型确定高程基准垂直偏差的方法,具体参见M.,et.al.Determination ofgeopotential differences between local vertical datums and realizationof a world height system.Stud.Geophys.Geod 2001(45):127-132)的优点计算位差,利用位差改正局部重力异常数据,从而使得陆海重力异常数据基准统一;在充分研究海岸带数据特点的基础上,提出了更为严密的海岸带重力场数据集成方法,该方法利用高精度数据控制低精度数据的思想,通过大地测量专业熟知的平差理论实现几种数据的高度集成应用,同时本方法可以避免高精度数据被低精度数据污染的可能性。本发明将在厘米级精度内实现海岛礁高程传递。
发明内容
1、一种海岛礁高程传递的优化方法,其特征在于,包括以下步骤:
1)各类重力数据的基准统一;
2)各类重力数据的集成;
4)利用上述步骤中得到的最终陆地空间重力异常与海洋空间重力异常简单几何拼接处理,形成陆海一致的海岛礁空间重力异常数据;
5)按现有方法对陆海一致的海岛礁空间重力异常数据进行格网化处理,形成陆海一致的空间重力异常格网数值模型;
6)利用空间重力异常格网数值模型利用现有Molodensky公式直接计算重力似大地水准面;
7)以重力似大地水准面为高程基准面,通过陆地高程已知点和海岛高程未知点的GPS大地高观测,得到两点的高程差,计算得到未知点的高程。
2、本发明内容1所述步骤2)中,关于各类重力数据的集成,其特征在于,提出了一种基于经典平差理论实现海岸带各类重力数据集成的方法,其步骤为:
1)本发明所指的各类重力数据,包括卫星测高数据、航空重力数据、船测重力数据、地面实测重力数据,将各类重力数据转化为空间重力异常,并统一到同一边界面;
2)将各类重力数据移去地形改正项;
3)计算卫星测高数据和航空重力数据的重力异常水平梯度矢量;
4)以地面重力异常计算得到水平梯度矢量作为真值,航空重力水平梯度矢量作为观测量,建立函数模型,按经典平差方法计算航空重力梯度矢量改正数;
5)利用水平梯度矢量改正数完成对航空重力异常数据的校正,并利用统计分析方法完成与地面重力异常的系统差计算,并进行系统差改正;
6)用同样的经典平差方法利用船测重力数据完成对卫星测高数据的校正:
7)将重力异常集成数据恢复地形改正则分别得到陆地和海洋空间重力异常数据。
3、一种海岛礁高程传递的优化方法,其特征在于,基于方法的重力异常高程基准系统差迭代改正方法,其方法步骤为:
2)陆地空间重力异常的系统差改正,得到新的陆地空间重力异常;
3)利用新的陆地空间异常计算陆地重力似大地水准面;
4)基于GPS水准数据与陆地重力似大地水准面,计算新系统差;
5)计算得到的新系统差与1)步骤计算的系统差相减得到系统差改正数;
6)利用系统差改正数校正2)步骤中计算得到的新的陆地空间重力异常数据,然后利用此数据继续计算重力似大地水准面,再次与GPS水准点比较,得到新的系统差改正数,依次循环,直到系统差改正数小于设计要求为止,此时,可以得到最终改正过的陆地空间重力异常。
附图说明
图1海岸带重力数据集成计算流程
具体实施方式
基于似大地水准面的高程传递原理是以重力似大地水准面为高程基准面,以通过陆地高程已知点和海岛高程未知点的GPS大地高观测,得到两点的高程差,计算得到未知点的高程,从而实现高程传递。如果区域似大地水准面精度较高,则可以借助GPS技术测量得到的高精度大地高求解得到测点高精度高程值。因此基于似大地水准面的高程传递方法的关键在于确定高精度的似大地水准面。在陆地区域,重力数据一般具有较好的精度和分辨率并且基准一致而使数据集成问题变得简单。当完全处于空旷海域,由于卫星测高数据具有较高精度亦使问题单一化。而在海岸带区域,由于卫星测高数据精度较低,地面实测重力数据与船测重力数据相对较少,有些区域会有航空重力数据,使得海岸带数据集成处理变得非常复杂。本发明的两点创新在于:1、基于方法的重力异高程基准系统差迭代改正算法(图1),可以消除局部高程基准与全球高程基准之间的系统差;2、海岸带重力数据集成处理,将经典平差方法引入到重力数据的集成计算(图2)。具体实施方法和步骤为:
一、各类重力数据的基准统一;
本发明所指的各类重力数据的基准统一,是将全部参与大地水准面精化的陆海重力数据转换到基准一致的地心坐标系中,并采用的统一潮汐基准,我国一般采用无潮汐基准。
二、各类重力数据的集成;
1)本发明所指的各类重力数据,包括卫星测高数据、航空重力数据、船测重力数据、地面实测重力数据,将各类重力数据转化为空间重力异常,并统一到同一边界面,本发明所指的边界面为地球表面;另外本集成数据为离散点形式存在。
2)将各类重力数据移去地形改正项;
对地形改正项的计算既可采用平面地形改正公式,又可采用球面地形改正公式。然后直接在各类重力数据上减去地形改正项,这样做会使各类重力数据变得较为平滑。
3)计算卫星测高数据和航空重力数据的重力异常水平梯度矢量;
本步骤所指的卫星测高数据和航空重力数据已经移去了地形改正项,即是上两步骤基础上的结果。卫星测高数据和航空重力数据的重力异常水平梯度矢量算法为相邻两点数据的差分即为梯度矢量。
4)以地面重力异常计算得到水平梯度矢量作为真值,航空重力水平梯度矢量作为观测量,建立函数模型,按经典平差方法计算航空重力梯度矢量改正数;
地面重力异常数据一般精度要远高于航空重力数据,地面重力异常的水平梯度矢量算法为地面两点间的重力异常,理当精度也高于航空重力梯度矢量数据,因此将其作为控制数据(真值)航空重力水平梯度矢量数据作为观测值,完全参照水准网平差方法建立函数模型和求解方式,以距离为定权,求得航空重力梯度矢量的改正数。平差方法及过程和原理完全与水准网平差一致。
5)利用水平梯度矢量改正数完成对航空重力异常数据的校正,并利用统计分析方法完成与地面重力异常的系统差计算,并进行系统差改正;
利用上述计算得到改正数完成对航空重力水平梯度矢量的改正,并按水平梯度矢量校正航空重力异常;利用统计方法计算地面重力异常与经校正的航空重力异常系统偏差,并进行系统偏差改正。
6)用同样的经典平差方法利用船测重力数据完成对卫星测高数据的校正;
按上述4)-5)步骤中相同方法完成对卫星测高数据的校正,即将船测重力数据的水平梯度矢量作为真值,对卫星测高数据进行相应处理。
7)将重力异常集成数据恢复地形改正则分别得到陆地和海洋空间重力异常数据。
在上述3)-7)步骤中,各类重力数据已经移去地形改正,因此这步骤就要恢复地形改正项,即在重力数据上加上地形改正项,则分别得到陆地和海洋空间重力异常。
2)陆地空间重力异常的系统差改正,得到新的陆地空间重力异常;
3)利用新的陆地空间异常计算陆地重力似大地水准面;
4)基于GPS水准数据与陆地重力似大地水准面,计算新系统差;
5)计算得到的新系统差与1)步骤计算的系统差相减得到系统差改正数;
6)利用系统差改正数校正2)步骤中计算得到的新的陆地空间重力异常数据,然后利用此数据继续计算重力似大地水准面,再次与GPS水准点比较,得到新的系统差改正数,依次循环,直到系统差改正数小于设计要求为止,此时,可以得到最终改正过的陆地空间重力异常(离散点形式)。
四、利用上述步骤中得到的最终陆地空间重力异常与海洋空间重力异常简单几何拼接处理,形成陆海一致的海岛礁空间重力异常数据;
本发明所指的简单拼接是指将陆地空间重力异常与海洋空间重力异常,将这些以离散点形式存在的文件,经统一格式后(例如格式可以为:点号,经度,纬度,空间重力异常值。只要格式统一即可),将两个格式统一的数据拷贝到一个文件,则完成简单拼接。
五、按现有方法对陆海一致的海岛礁空间重力异常数据进行格网化处理,形成陆海一致的空间重力异常格网数值模型;
六、利用空间重力异常格网数值模型利用现有Molodensky公式直接计算重力似大地水准面;
七、以重力似大地水准面为高程基准面,通过陆地高程已知点和海岛高程未知点的GPS大地高观测,得到两点的高程差,计算得到未知点的高程。
未知点高程HB计算公式为:HB=HA+[(hB-ζB)-(hA-ζA)]
其中,B代表未知高程点,A代表已知高程点,hA、hB分别代表已知点和未知点的大地高,ζA、ζH分别代表已知点和未知点的高程异常,HA、HB分别代表已知点和未知点的高程。
Claims (4)
1.一种海岛礁高程传递的优化方法,其特征在于,包括以下步骤:
1)各类重力数据的基准统一;
2)各类重力数据的集成;
4)利用上述步骤中得到的最终陆地空间重力异常与海洋空间重力异常简单几何拼接处理,形成陆海一致的海岛礁空间重力异常数据;
5)按现有方法对陆海一致的海岛礁空间重力异常数据进行格网化处理,形成陆海一致的空间重力异常格网数值模型;
6)利用空间重力异常格网数值模型利用现有Molodensky公式直接计算重力似大地水准面;
7)以重力似大地水准面为高程基准面,通过陆地高程已知点和海岛高程未知点的GPS大地高观测,得到两点的高程差,计算得到未知点的高程。
2.根据权利要求1所述步骤2)中,关于各类重力数据的集成,其特征在于,提出了一种基于经典平差理论实现海岸带各类重力数据集成的方法,其步骤为:
1)本发明所指的各类重力数据,包括卫星测高数据、航空重力数据、船测重力数据、地面实测重力数据,将各类重力数据转化为空间重力异常,并统一到同一边界面;
2)将各类重力数据移去地形改正项;
3)计算卫星测高数据和航空重力数据的重力异常水平梯度矢量;
4)以地面重力异常计算得到水平梯度矢量作为真值,航空重力水平梯度矢量作为观测量,建立函数模型,按经典平差方法计算航空重力梯度矢量改正数;
5)利用水平梯度矢量改正数完成对航空重力异常数据的校正,并利用统计分析方法完成与地面重力异常的系统差计算,并进行系统差改正;
6)利用上述3)-5)步骤中的方法,利用船测重力数据完成对卫星测高数据的校正;
7)将重力异常集成数据恢复地形改正则分别得到陆地和海洋空间重力异常数据。
2)陆地空间重力异常的系统差改正,得到新的陆地空间重力异常;
3)利用新的陆地空间异常计算陆地重力似大地水准面;
4)基于GPS水准数据与陆地重力似大地水准面,计算新系统差;
5)计算得到的新系统差与1)步骤计算的系统差相减得到系统差改正数;
6)利用系统差改正数校正2)步骤中计算得到的新的陆地空间重力异常数据,然后利用此数据继续计算重力似大地水准面,再次与GPS水准点比较,得到新的系统差改正数,依次循环,直到系统差改正数小于设计要求为止,此时,可以得到最终改正过的陆地空间重力异常。
4.一种海岛礁高程传递的优化方法,其应用范围不仅仅限于海岛礁高程传递,在单一陆地和单一海洋上,本方法能被独立使用。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN2010101986240A CN101957193B (zh) | 2010-06-12 | 2010-06-12 | 一种海岛礁高程传递的优化方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN2010101986240A CN101957193B (zh) | 2010-06-12 | 2010-06-12 | 一种海岛礁高程传递的优化方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN101957193A true CN101957193A (zh) | 2011-01-26 |
CN101957193B CN101957193B (zh) | 2012-02-22 |
Family
ID=43484665
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN2010101986240A Expired - Fee Related CN101957193B (zh) | 2010-06-12 | 2010-06-12 | 一种海岛礁高程传递的优化方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN101957193B (zh) |
Cited By (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102230795A (zh) * | 2011-03-24 | 2011-11-02 | 鲍李峰 | 利用重力位差实现跨海高程基准传递 |
CN103245324A (zh) * | 2012-02-06 | 2013-08-14 | 中国测绘科学研究院 | 海岛遥感测图高程精度控制与修正方法及系统 |
CN104567802A (zh) * | 2015-01-06 | 2015-04-29 | 山东科技大学 | 集成船载重力和gnss的测线式陆海高程传递方法 |
CN112985348A (zh) * | 2021-02-08 | 2021-06-18 | 雷军珍 | 利用重力、gnss-rtk和重力场模型实现高程测量方法 |
CN113819882A (zh) * | 2021-09-09 | 2021-12-21 | 江苏海洋大学 | 一种跨海高程点间重力位差计算方法 |
Families Citing this family (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102841386A (zh) * | 2012-09-25 | 2012-12-26 | 中国石油化工股份有限公司 | 一种海洋重力调查提取重力异常的方法 |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101493324A (zh) * | 2009-03-10 | 2009-07-29 | 长安大学 | 基于cqg2000的区域似大地水准面精化方法 |
CN101713649A (zh) * | 2009-11-05 | 2010-05-26 | 中国测绘科学研究院 | 一种基于扰动重力的似大地水准面陆海无缝拼接方法 |
-
2010
- 2010-06-12 CN CN2010101986240A patent/CN101957193B/zh not_active Expired - Fee Related
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101493324A (zh) * | 2009-03-10 | 2009-07-29 | 长安大学 | 基于cqg2000的区域似大地水准面精化方法 |
CN101713649A (zh) * | 2009-11-05 | 2010-05-26 | 中国测绘科学研究院 | 一种基于扰动重力的似大地水准面陆海无缝拼接方法 |
Non-Patent Citations (6)
Title |
---|
《北京测绘》 20050331 王岩等 三种高程传递方法的精度分析 , 第1期 2 * |
《武汉大学学报(信息科学版)》 20011231 李建成等 长距离跨海高程基准传递方法的研究 第26卷, 第6期 2 * |
《海洋测绘》 20100930 郑童瑜等 无控制岛礁高程的确定 第30卷, 第5期 2 * |
《海洋测绘》 20110331 王鹏等 跨海高程传递方法与实践 第31卷, 第2期 2 * |
《港工技术》 20050630 姜晨光等 国家高程系统向岛礁精确传递问题的研究 , 第2期 2 * |
《现代测绘》 20070331 李凤斌等 长距离跨海高程基准传递方法及精度 第30卷, 第2期 2 * |
Cited By (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102230795A (zh) * | 2011-03-24 | 2011-11-02 | 鲍李峰 | 利用重力位差实现跨海高程基准传递 |
CN102230795B (zh) * | 2011-03-24 | 2013-03-13 | 鲍李峰 | 利用重力位差实现跨海高程基准传递 |
CN103245324A (zh) * | 2012-02-06 | 2013-08-14 | 中国测绘科学研究院 | 海岛遥感测图高程精度控制与修正方法及系统 |
CN103245324B (zh) * | 2012-02-06 | 2015-04-08 | 中国测绘科学研究院 | 海岛遥感测图高程精度控制与修正方法及系统 |
CN104567802A (zh) * | 2015-01-06 | 2015-04-29 | 山东科技大学 | 集成船载重力和gnss的测线式陆海高程传递方法 |
CN112985348A (zh) * | 2021-02-08 | 2021-06-18 | 雷军珍 | 利用重力、gnss-rtk和重力场模型实现高程测量方法 |
CN112985348B (zh) * | 2021-02-08 | 2022-08-16 | 雷军珍 | 利用重力、gnss-rtk和重力场模型实现高程测量方法 |
CN113819882A (zh) * | 2021-09-09 | 2021-12-21 | 江苏海洋大学 | 一种跨海高程点间重力位差计算方法 |
CN113819882B (zh) * | 2021-09-09 | 2023-06-16 | 江苏海洋大学 | 一种跨海高程点间重力位差计算方法 |
Also Published As
Publication number | Publication date |
---|---|
CN101957193B (zh) | 2012-02-22 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN101957193B (zh) | 一种海岛礁高程传递的优化方法 | |
CN102506824B (zh) | 一种城市低空无人机系统生成数字正射影像图的方法 | |
CN102175213B (zh) | Gps水准高程传递方法 | |
CN108896040B (zh) | 天空海一体化水下潜器惯性/重力组合导航方法和系统 | |
CN101839710A (zh) | 一种似大地水准面计算的优化方法 | |
CN101493324A (zh) | 基于cqg2000的区域似大地水准面精化方法 | |
Xu | Application of GPS-RTK technology in the land change survey | |
CN110440762B (zh) | 一种多旋翼无人机山区航测成图的网格化像控点布设方法 | |
CN102236108A (zh) | 一种磁性地表三维地形改正方法 | |
CN108919304A (zh) | 一种基于参考平面的移动测量系统中pos误差补偿方法 | |
CN103727920A (zh) | 基于大地水准面模型测量水准高差的方法 | |
Oluyori et al. | Computations of geoid undulation from comparison of GNSS/levelling with EGM 2008 for geodetic applications | |
CN102589528A (zh) | 一种多时相影像的海岛岸线量测方法 | |
CN102230795A (zh) | 利用重力位差实现跨海高程基准传递 | |
CN104181571A (zh) | 实现弱或无cors信号区地面点的精密坐标与高程快速测量的方法 | |
CN115423955B (zh) | 一种基于多源数据的最优深度基准面大地高模型构建方法 | |
CN115346128A (zh) | 一种光学立体卫星dem高程改正和融合方法 | |
Xu et al. | Approach to minimizing the influence of changeable observational environment on long-distance sea-crossing trigonometric leveling | |
CN106125149A (zh) | 点质量模型中浅层高分辨率点质量最佳埋藏深度确定方法 | |
Mo et al. | Evaluation of AW3D30 Elevation Accuracy in China. | |
CN105651310B (zh) | 一种新型测绘仪倾斜量修正方法 | |
CN113358092B (zh) | 一种确定我国高程基准的垂直偏差的大数据数值算法 | |
Sulaiman et al. | Geoid model estimation without additive correction using KTH approach for Peninsular Malaysia | |
CN102944220A (zh) | 重力水准面与gps水准差异分解融合方法 | |
CN113450348B (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 | ||
C17 | Cessation of patent right | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20120222 Termination date: 20120612 |