CN102944220B - 重力水准面与gps水准差异分解融合方法 - Google Patents
重力水准面与gps水准差异分解融合方法 Download PDFInfo
- Publication number
- CN102944220B CN102944220B CN201210500297.9A CN201210500297A CN102944220B CN 102944220 B CN102944220 B CN 102944220B CN 201210500297 A CN201210500297 A CN 201210500297A CN 102944220 B CN102944220 B CN 102944220B
- Authority
- CN
- China
- Prior art keywords
- gps
- gravity
- level
- level surface
- residual error
- 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.)
- Expired - Fee Related
Links
Landscapes
- Position Fixing By Use Of Radio Waves (AREA)
Abstract
一种重力大地水准面与GPS水准差异分解融合新方法，属于大地测量学与测量工程技术领域的发明专利。重力大地水准面分辨率高但有系统偏差，GPS水准精度高但分辨率低，实用中必须将两者融合。由于大地水准面形状不规则，非常复杂，然而目前国内国际用多项式、球谐等函数拟合的方法得到的是规则曲面，精度和分辨率提高都很困难。本发明采用趋势函数先把重力水准面与GPS水准的系统偏差移除，残差部分以放大镜原理采用最小曲率Green基函数样条移动插值，形成与重力水准面同等分辨率的残差格网，再把残差格网、系统偏差与重力水准面格网叠加，形成最终的GPS水准重力大地水准面。该发明的优点在于精度高，不受分辨率的约束，能处理任意大小、形状复杂的各级大地水准面，可广泛应用于大地水准面建模中。
Description
技术领域
本发明是一种重力大地水准面与GPS水准差异分解融合新方法，属于“大地测量学与测量工程”学科中的“物理大地测量学”技术领域。
背景技术
随着现代大地测量技术的发展，以及在地球物理学、海洋学等学科中应用的不断扩展，建立高精度高分辨率的大地水准面模型，是解决现代测绘生产、地学科研实验的迫切任务。为了满足地学应用的高要求，需要不断提高大地水准面分辨率和精度。区域大地水准面建模流程是，在全球地球重力场的基础上，结合各种重力观测数据、高分辨率的数字地形模型，经过严密计算形成重力大地水准面，然后再与GPS水准相结合形成GPS水准重力大地水准面。本发明主要解决重力水准面与GPS水准高精度融合问题。
重力大地水准面分辨率高但有系统偏差且精度较低，GPS水准精度高但分辨率低，实用中必须将两者融合。重力大地水准面与GPS水准之间的差异形状不规则，非常复杂，国内国际目前采用多项式函数、球谐函数等拟合方法，然而得到的是规则曲面，因而模型存在规则化的近似。用规则函数去构造一个不规则的曲面存在理论上方法上的缺憾，进一步提高分辨率和模型精度面临巨大困难。本发明采用趋势函数先把重力水准面与GPS水准的系统偏差移除，残差部分采用最小曲率Green基函数样条，在较小的局部范围内利用放大镜的原理，不断对网格结点进行移动观测插值，形成与重力水准面同等分辨率的残差格网，再把残差格网、系统偏差与重力水准面格网叠加，形成最终的GPS水准重力大地水准面。该技术方案与现有采用规则函数对残差进行整体拟合的方法存在本质的区别，可以实现任意不规则曲面的构建而不受到分辨率与精度的限制。本该发明的优点在于精度很高，能处理任意大小、形状复杂的各级大地水准面，可广泛应用于大地水准面建模中。
发明内容
一种重力大地水准面与GPS水准差异分解融合新方法，其基本目标是将高分辨率的重力大地水准面与高精度的GPS水准实现完美融合，从而构造最终的(似)大地水准面格网模型。本发明的两个创新点在于：一是完全抛弃旧的技术路线，不采用函数拟合的方法来拟合重力大地水准面与GPS水准的差异，而是先把趋势项从两者差异中移除，然后再对残差构造格网，最后把趋势项、残差格网、重力大地水准面三者叠加，获得最终的大地水准面模型；二是以结点周围关系密切的GPS水准点对结点进行插值计算，采用放大镜原理仅对局部细微特征进行处理，无需整体建模，从而解决高分辨率与高精度插值问题。
本发明的基本特征在于，针对重力大地水准面与GPS水准之间存在的差异，把该差异分成低频的趋势项和高频的残差两部分分别处理，针对趋势项的构造提出了优化趋势函数；针对残差部分提出了移动格林基样条最小曲率插值方法形成残差格网。
在构建最终GPS重力水准大地水准面时，以高分辨率的重力水准面为基础，恢复趋势项和格网残差，这是全新的技术方案，与差异函数拟合方法存在本质的不同，函数拟合能力是有限的，对高频变化的起伏实现高精度的拟合非常困难，而本发明能对任意复杂的大面积大地水准面进行重力GPS水准融合数据处理，基于最小二乘原理求解参数估值，获得非常高的计算精度。
本发明的关键技术特征和流程包括以下内容：
(A) 由格林基样条插值计算GPS水准点的重力高程异常：
(B) 由GPS水准实测高程异常与插值重力高程异常相减得高程异常差值；
(C) 给出了趋势函数形式，基于最小二乘原理，提取高程异常差异的趋势项，获得趋势函数的系数项；
(D) 把趋势项从高程异常差异中移除，得到高程异常差异的残差；
(E) 针对高程异常差异的残差，采用格林基样条插值方法对格网的结点进行残差插值，得到残差格网模型；
(F) 以重力大地水准面为基础，恢复趋势项，再与残差格网模型进行叠加，获得最终的GPS水准重力大地水准面模型；
附图说明
图1 重力大地水准面与GPS水准融合新方法流程图
具体实施方式
一种重力大地水准面与GPS水准融合的新方法，其特征在于，具体实施中包括以下步骤。
l）重力大地水准面与GPS水准差异计算：
(A) 重力大地水准面模型是规则的格网，GPS水准是任意均衡分布的离散观测点，取得GPS水准点的坐标；
(B) 由格网插值方法，直接插值计算GPS水准点坐标位置的重力水准面高程异常；
(C) 将计算出来的重力水准面高程异常与GPS水准观测的高程异常进行相减，得到重力大地水准面与GPS水准的差异。
2）在重力水准面与GPS水准面的差异的基础上，基于最小二乘原理构造差异的趋势函数，即差异信息的低频部分：
(A) 趋势函数的形式是：
其中，为椭球的扁率；为纬度和经度；为待求参数；
(B) 设矩阵A,X,L为：
；
(C) 解算趋势函数的参数为：。
3）把趋势项部分从重力水准面与GPS水准差异中移除；
(A) 根据GPS水准点的大地经纬度坐标，由步骤3)求得的趋势函数系数，计算趋势差异值；
(B) 把1)计算的实际差异值与趋势差异值相减，即得到移除趋势项的重力大地水准面与GPS水准差异之残差。
4）把移除趋势项的重力大地水准面与GPS水准差异之残差形成与重力大地水准面分辨率相同的规则格网，具体步骤如下：
(A) 设采用最邻近个GPS水准点的残差来构建残差格网；
(B) 设这个GPS水准点的坐标为，；
(C) 设这个GPS水准点的剩余残差为；
(D) 构造Green基函数：，其中，；
(E) 计算权矩阵；
(F) 计算矩阵，其中为结点到个GPS水准点的距离；
(G) 计算结点P应分配的残差值：；
(H) 依次计算其它结点的残差值，形成残差格网。
5）把残差格网、系统偏差与重力水准面格网叠加，形成最终的GPS水准重力大地水准面。
Claims (2)
1.一种重力水准面与GPS水准差异分解融合方法，其特征在于，包括以下步骤：
l)求取重力水准面与GPS水准之间的差异值：
(A)根据GPS水准点的大地坐标，由格林基样条插值计算GPS水准点位置的重力高程异常；
(B)由GPS大地高减去水准正常高程得到GPS水准高程异常；
(C)由重力高程异常减去GPS水准高程异常得到高程异常差异值；
2)从高程异常差异值中提取趋势项：
(A)趋势项具有如下形式：
其中e为椭球的扁率；λ为纬度和经度；ai为待求参数；
(B)设矩阵A,X,L为：
(C)解算趋势函数的参数为：X＝(ATA)-1(ATL)，从而获得趋势项；
3)把趋势项部分从重力水准面与GPS水准差异中移除；
(A)根据GPS水准点的大地经纬度坐标，由步骤2)求得的趋势函数系数，计算趋势差异值；
(B)把1)计算的实际差异值与趋势差异值相减，即得到移除趋势项的重力大地水准面与GPS水准差异之残差；
4)把移除趋势项的重力大地水准面与GPS水准差异之残差形成与重力大地水准面分辨率相同的规则格网：(A)设采用最邻近k个GPS水准点的残差来构建残差格网；
(B)设这k个GPS水准点的坐标为X＝[x1 x2 ... xk]T，Y＝[y1 y2 ... yk]T；
(C)设这k个GPS水准点的剩余残差为Z＝[z1 z2 ... zk]T；
(D)构造Green基函数：其中
(E)计算权矩阵W＝G-1Z；
(F)计算矩阵Gp＝[d01 d02 ... d0k]，其中d0i为结点到k个GPS水准点的距离；
(G)计算结点P应分配的残差值：zp＝GpW；
(H)依次计算其它结点的残差值，形成残差格网；
5)把残差格网、系统偏差与重力水准面格网叠加，形成最终的GPS水准重力大地水准面；
2.根据权利要求1所述的重力水准面与GPS水准差异分解融合方法，应用在GPS水准重力大地水准面模型建立过程中。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201210500297.9A CN102944220B (zh) | 2012-11-30 | 2012-11-30 | 重力水准面与gps水准差异分解融合方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201210500297.9A CN102944220B (zh) | 2012-11-30 | 2012-11-30 | 重力水准面与gps水准差异分解融合方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN102944220A CN102944220A (zh) | 2013-02-27 |
| CN102944220B true CN102944220B (zh) | 2015-01-07 |
Family
ID=47727185
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201210500297.9A Expired - Fee Related CN102944220B (zh) | 2012-11-30 | 2012-11-30 | 重力水准面与gps水准差异分解融合方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN102944220B (zh) |
Families Citing this family (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN104792305B (zh) * | 2015-05-05 | 2017-01-18 | 长沙理工大学 | Gps水准转换模型核函数光滑因子确定方法 |
| CN115098830A (zh) * | 2022-06-28 | 2022-09-23 | 中国人民解放军战略支援部队信息工程大学 | 一种水准面和水准椭球曲率半径差异计算方法及装置 |
Family Cites Families (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US8152366B2 (en) * | 2008-01-23 | 2012-04-10 | University Of Delaware | Estimation of subsurface thermal structure using sea surface height and sea surface temperature |
| CN101713649B (zh) * | 2009-11-05 | 2011-03-23 | 中国测绘科学研究院 | 一种基于扰动重力的似大地水准面陆海无缝拼接方法 |
| CN101839710B (zh) * | 2010-06-12 | 2012-01-04 | 中国测绘科学研究院 | 一种似大地水准面的确定方法 |
-
2012
- 2012-11-30 CN CN201210500297.9A patent/CN102944220B/zh not_active Expired - Fee Related
Also Published As
| Publication number | Publication date |
|---|---|
| CN102944220A (zh) | 2013-02-27 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN102506824B (zh) | 一种城市低空无人机系统生成数字正射影像图的方法 | |
| US7944547B2 (en) | Method and system of generating 3D images with airborne oblique/vertical imagery, GPS/IMU data, and LIDAR elevation data | |
| CN109031301A (zh) | 基于PSInSAR技术的山区地形形变提取方法 | |
| CN109472802B (zh) | 一种基于边缘特征自约束的表面网格模型构建方法 | |
| CN105158760A (zh) | 一种利用InSAR反演地下流体体积变化和三维地表形变的方法 | |
| CN105929398A (zh) | 结合外控点的InSAR高精度高分辨率DEM获取方法 | |
| CN109100719B (zh) | 基于星载sar影像与光学影像的地形图联合测图方法 | |
| CN104535080B (zh) | 大方位失准角下基于误差四元数的传递对准方法 | |
| CN104574512A (zh) | 一种顾及地形语义信息的多尺度dem构建方法 | |
| CN103927418A (zh) | 基于dem的城市道路渠网化排水通道制作方法 | |
| CN105354832A (zh) | 一种山区卫星影像自动配准到地理底图上的方法 | |
| CN110415165A (zh) | 一种极区等间距格网点的平面坐标与大地坐标相互转换方法 | |
| CN114283070B (zh) | 融合无人机影像与激光点云的地形断面制作方法 | |
| CN108427741B (zh) | 一种基于大量高精度控制点的dem相对误差评价方法 | |
| CN106157258B (zh) | 一种星载sar图像几何校正方法 | |
| CN102944220B (zh) | 重力水准面与gps水准差异分解融合方法 | |
| CN111257870A (zh) | 一种利用InSAR监测数据的采煤沉陷积水区水下地形反演方法 | |
| CN113238228B (zh) | 基于水准约束的三维地表形变获取方法、系统及装置 | |
| CN110310370B (zh) | 一种gps与srtm点面融合的方法 | |
| CN111986320A (zh) | 面向智慧城市应用的dem与倾斜摄影模型空间贴合优化算法 | |
| CN111080795A (zh) | 一种基于Delaunay三角网的重力场三维模型构建方法 | |
| CN102930176B (zh) | 大地水准面模型九网格五样条插值方法 | |
| CN116045920A (zh) | 一种基于gf7号立体像对提取dem的方法 | |
| CN115112090A (zh) | 一种基于少量gnss水准点的山区数字高程基准建立方法 | |
| CN114777745A (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 | ||
| CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20150107 Termination date: 20181130 |
|
| CF01 | Termination of patent right due to non-payment of annual fee |