地形起伏对GPS工程控制网高程异常的影响

1、地形起伏对GPS工程控制网高程异常的影响张同刚1基金项目：国家自然科学基金资助项目（40271092） 岑敏仪1 冯义从1 路伯祥1 卢建康2（1.西南交通大学 土木工程学院 地理信息工程中心 四川 成都 610031；2.铁道第二勘测设计院 四川 成都 610031）摘 要：高程异常由短波分量和中长波分量组成。受地形起伏影响显著，特别是在山区和高山区，变化较大，而的变化则相对比较平缓。利用DEM计算GPS控制点的，再结合移除-恢复技术，可以显著提高局部地区的GPS大地高转换为正常高的精度。本文就其核心问题的计算进行了深入的研究，导出了连续型积分计算公式；分别试验了平原、丘陵、和高山地区在不同

2、积分范围和不同参考面高程对计算结果的影响，并分析了这些因素对GPS大地高转化为正常高的影响特点，由此得出一些有益的结论。关键词：高程异常；正常高；全球定位系统；数字高程模型中图分类号: P228文献标识码: AThe Affection of the Rugged Topography on the Height Anomaly of GPS Engineering Control NetworkZHANG Tong-gang1CEN Min-yi1FENG Yi-cong1LU Bai-xiang1LU Jian-kang2（1. College of Civil Engineering,

3、Southwest Jiaotong University, Chengdu Sichuan, 610031; 2. The Second Railways Survey & Design Institute, Chengdu Sichuan, 610031）Abstract: The height anomaly can be separated into two components, , the short wavelength component and , the medium and long wavelength one. varies with topography great

4、ly, especially in the mountain and high mountain area, whileis relatively smoother. The accuracy of the height transformation of GPS height can be improved greatly in local mountain area by employing the remove-and-recover technique after of each GPS point calculated using DEM. We carefully study th

5、e kernel question, which is the calculation of the height anomaly caused by the rugged topography, and educe a continuous integral formula. Then, a series of experiments are designed for studying and analyzing the affection of scope of DEM involved in the integral process and the selection of the re

6、ference height. All of the experiments are performed in the plain, hill, and mountain areas respectively. After analysis of the results, some helpful conclusions are drawn by these tests.Key words: height anomaly；normal height；GPS；DEM 我国西部多为高山丘陵地形，对铁路建设中的高程测量提出了很高的要求，如果采用传统的水准测量方法将遇到很多困难，而解决这个问题的有效途

7、径就是通过将GPS大地高转换为正常高的方法，目前国内外很多学者在这方面进行了大量的研究，并且提出了很多方法。一是拟合法，其本质是用连续的数学曲面来拟合似大地水准面。在平原地区由于似大地水准面比较光滑，选择合适的数学模型就能够满足一般工程建设的要求1，但对高山区等似大地水准面较为粗糙的地区，效果则并不理想。二是重力场模型法2和重力场模型结合内插的方法34，这种方法通过重力场模型求得高程异常值，并结合水准测量和一定的内插方法将GPS测得的大地高转换为正常高。该方法多运用于大范围的GPS高程转换。三是顾及地形改正的拟合法56，这种方法的关键在于首先通过DEM来计算高程异常的短波分量，然后剔除该分量，

8、使得高程异常的剧烈变化变平缓后再进行拟合，最后恢复被剔除的短波分量。试验结果表明这种方法能够适用于起伏较大的高山地区。文献5中给出了高程异常的计算公式，其试验结果显示，利用第三种方法绝大部分点达到了四等水准测量的精度。该文采取了增加待求点与DEM格网点之间的水平距离的办法来处理其积分奇异问题。文献6对这种方法做了进一步的研究和分析，提出了分段积分方法，来解决积分奇异的问题，从而改善了高程异常的计算精度。不足之处是这个公式在计算过程中需要有人工干预，可能会给结果带来一些不确定因素。其他的方法还有神经网络模型，以及神经网络模型和曲面拟合相结合的混合方法7等。目前由于计算机技术和遥感技术的发展，DE

9、M在数据获取、存储和处理速度等方面取得了突破性的进展，国家1：50 000DEM已建立，这就使得顾及地形改正的拟合法的应用前景非常广阔，特别是对于西部地区的铁路建设具有很重要的现实意义，因此有必要对该方法做进一步的研究。由于求解高程异常的精度是影响该方法应用的重要因素，本文就此进行更加深入的试验和研究。本文首先推导高程异常求解的连续积分计算公式，然后详细讨论影响计算高程异常的各种因素，并设计试验方案，分析试验结果，最后得出一些结论和建议。1 高程异常求解的连续积分公式如图 1所示，设测区有一参考面Hr，则高出或者低于Hr的地形对P点的引力位为8： （1）式（1）中：G为万有引力常数（G = 6

10、.67310-8C3S-2g-1），为地球平均质量密度（= 2.67gC-3），H为数字地面模型格网点高程，Hr为参考面高程。由地形起伏引起的高程异常为8： （2）式（2）中 为计算点的正常重力值，即（3）图 1 由地形起伏计算高程异常Fig 1 参考面似大地水准面rr0HrHpzP起伏地面图H式（3）中为测区平均纬度，Hp为计算点高程（单位：km）。文献5对积分公式的积分奇异问题采取了一个简单的近似处理，当（见图 1）小于0.5km时将设为0.5km，来应对积分奇异；文献6给出分段积分计算公式，计算中需人工干预，不利于编程计算，同时带来了计算结果的不确定性。因此有必要进一步完善积分公式。为了

11、便于积分运算，文献6中对式（1）使用泰勒级数展开，最终得到了分段积分公式。本文采用了另一种不同的处理方式，对式（1）直接积分：（4）式（1）（4）是计算高程异常的连续积分公式。注意到式（4）中，当待求点与格网点重合时，此时，使得，导致积分异常。为避免积分异常，在0的DEM格网点的高程值增加10cm，其他DEM格网点并不做这样的改正，这样处理对高程异常积分计算结果的影响为0.7mm。考虑到1：1万DEM格网点高程中误差平原地区为0.51m，山区为2.55m，高山区为510m9，所以这样处理能够满足实际工程的要求，利于编程实现且可避免积分异常。2 DEM范围对高程异常的影响整个测区的范围很大，相应

12、的DEM数据量就非常大，需要多大范围的DEM参与计算，才能使计算的高程异常值满足实际需要的精度？这个问题涉及到要准备的DEM数据量和实际计算效率，它对工程控制网的应用具有重要的现实意义。针对这个问题，本文首先试验与待求点P距离不等的地面点对P的高程异常的影响，然后试验不同地形，不同积分区域对待求点P 的高程异常及其高程异常差的影响。图 2 DEM格网点对P点的高程异常的影响（上半部分Z轴的单位10-5m；下的半部分Z轴的单位为m，DEM格网间距250m）YXZ8 Z4ACP2.1 不同格网点对待求点高程异常值的影响对待求点P，分别计算DEM上各格网对P点的高程异常的影响值，然后以对应的DEM格

13、网为平面坐标，以影响值为Z坐标，在空间三维图上表示，则DEM上每个格网对P点高程异常的影响值分布图如图 2所示，图中上半部分为分布于DEM上各格网点对P点高程异常的影响值，若取整个上半部分区域面上的影响值之和，则可求得由该DEM地形起伏对P点的高程异常。图中下半部分为参与计算P点高程异常的DEM。P点的位置位于DEM的中央。由图 2上半部分可以看出，中央部分（即P点所在位置）有一个明显的凸起，意味着该点附近的点对高程异常的影响较大，而其它每个点的影响要小得多；在X90区域（下称区域C）内单点的影响总体上略微超过X90部分（下称区域A）的点。在参考面一定的情况下，影响计算P点的高程异常因素有二个

14、主要方面，一是格网与P点之间的距离，二是格网高程与参考高程面的高差。在与P点距离相同的情况下，尽管在区域C中点的高程显著大于在区域A中的点，但其对高程异常的影响只是稍大一点。所以距离因素在这二个因素中占据主导地位。试验发现，DEM格网对高程异常的影响并没有随着距离增大而完全为零，随着DEM范围的扩大，这些微小影响的数量也越多，总的影响不能忽视。2.2 不同范围的DEM对高程异常差的影响不同范围大小的DEM计算的试验区格网点高程异常值的变化非常大，而工程控制网主要考虑的是高程点的绝对点位（相对于控制网起算点或平差基准的高差）和相对点位（高程点之间的高差）精度，这种要求在铁路、公路、管线等工程建设

15、的带状控制网中更明显。因此，在工程控制网的GPS高程转换中，主要是使用高程异常差，因而需要计算和分析不同范围大小的DEM对试验区格网点高程异常差的影响。为比较分析的需要，本文以平原、丘陵和高山的DEM为代表性地区，试验方案如下： 1DEM为240240，间距250m的规则格网点；2为保证每个待求格网点均不在DEM的边缘部分，取DEM中央部分6060个格网作为试验区域；3每次实验在试验区域的计算格网点上，以上、下、左、右各分别增加10、20、90个格网的DEM区域（称计算区域）来计算其高程异常影响值；4相邻两次不同面积的计算区域计算的高程异常值相减，求出不同范围大小的DEM对试验区格网点的高程异

16、常差的影响值。实验发现，当平原和丘陵地区的DEM计算区域大于（40+40）（40+40）格网面积时，高程异常差的变化很小，DEM计算区域从（40+40）（40+40）增加到（50+50）（50+50）格网面积时，平原和丘陵二个DEM试验区域（6060格网面积，区域面积1515km2）的高程异常差变化的最大值分别为0.4mm和1.3mm；而在高山地区，当DEM计算区域从（70+70）（70+70）增加到（80+80）（80+80）格网面积时，试验区域的高程异常差变化最大值为55.6mm，参与计算的DEM范围继续增大，给高程异常差带来的影响会越来越小。图 3 高山地区DEM半径为17.5km-20

17、km之间的高程异常较差（格网间距250m,等值线间隔0.005m）图3表示DEM区域从（70+70）（70+70）增加到（80+80）（80+80）格网面积时，计算的高山地区试验区域高程异常差的等值线图，从图 3可以明显看出高程异常差在X方向的变化要大于Y方向上的变化，这一点与实际地形坡度相对应；区域A、C高程异常差总体上的变化在0.160.4cm/km之间，区域B高程异常差的变化较大，在0.61.2cm/km之间，该地区的高差很大，在600m左右。由此可以看出，GPS水准点如果选取在高程相近的地方，有利于减弱地形起伏给高程异常差带来的影响，进而有利于最后的高程转换精度。相同的格网面积，如果D

18、EM的格网间距不同，那么参与积分计算的实际地面面积就不一样。这样就带来了一个新的问题，以实际地面面积还是以DEM的格网面积为准？为此，在上述平原、丘陵和高山的DEM上重采样，格网间距分别为125米和62.5米，在新的DEM上重新进行上述试验。从高山试验区的试验结果发现，如果保持实际地面面积不变，不同分辨率的DEM之间的计算结果相差非常小，差值范围在00.08cm/km之间。如果保持格网面积不变，相应的实际地面面积分别减少为原来的1/4和1/16，不同实际地面面积的DEM，计算的高程异常差结果差别显著，其范围在0.52cm/km之间。在平原和丘陵试验区均有相似的结论。通过这些试验结果可以看出影响

19、高程异常差的是实际地面面积，而不是DEM的格网面积，DEM的分辨率对计算高程异常差的影像可以忽略。综合分析以上实验结果，可知在平原和丘陵地区参加计算高程异常的DEM区域范围，一般比测区范围上下左右大于10km即可，在高山地区应大于18km。3 参考面高程的选择（4）式中，选择不同的参考面对高程异常和高程异常差有多大的影响？如何选择一个合适的参考面高程？这是使用者普遍关心的问题。为此本文设计了两组方案进行试验，并对结果数据进行了分析。3.1 参考面高程对高程异常差的影响本文选取特大山区QL为试验区，DEM为432480格网面积, 格网间隔250m，山区最低处在400m以下，最高处在2 500m以

20、上，高差超过2 100m。ABCD图 4 同一区域选用不同高程参考面引起的高程异常的等值线图（ABCD分别为参考面高程为500，1 000，1 500，2 000m时的高程异常等值线，等值线间距A:0.05m,B:0.02m,C:0.1m,D:0.2m，坐标单位km）试验分别选择了500m、1 000m、1 500m、2 000m为参考面高程，以计算点四周约18km（90个格网）DEM区域为计算区域，使用前述连续积分计算公式分别计算试验区QL中央252300格网面积的高程异常值，然后计算各点在不同高程面之间（分别为500-1 000，1 000-1 500，1 500-2 000m）的高程异常

21、差，试验结果如图4、图5所示。图 5 同高程参考面之间的高程异常差的等值线图（I，II，III分别是参考面高程500m-1 000m，1 000m-1 500m，1 500m-2 000m之间的高程异常较差的等值线图，等值线间距0.1m，坐标单位km）IIIIII一方面，图5中，四组试验结果差别非常明显，相同点处的高程异常在不同的参考高程面下的高程异常值数值相差超过1m（见 图5），因而不同的参考面高程对高程异常值的影响十分显著。另一方面，从图5.I，II，III中的任意一幅图都可以看出，不同参考高程面对高程异常差的影响是非常大的，最大已经超过了1cm/km；比较这三幅图，可以看出这三个等值线

22、图几乎相同，这是由于每幅图中选取的两个参考高程面之差均相同（等于500m）造成的。经过更多的试验发现每个点的高程异常值随参考高程面的变化呈近似的线性变化，但不同位置的点，其变化率不同。而正是这个变化率不同导致了采用不同参考面高程引起了计算的高程异常差的变化。通过以上分析可知，参考高程面是不能任意选取的。3.2 参考面高程对拟合高程的影响不同的参考高程面计算的任意两个点间的高程异常差并不相同，人们自然会想到需要进一步探讨究竟如何采用参考高程面才是合适和满足工程测量要求的。因此，需要试验采用不同参考面分别计算GPS控制点的高程异常，在所有GPS控制点的大地高中除去地形起伏产生的，并选择其中若干个G

23、PS点作为检查点，用剩下的GPS点作为控制点，采用多项式拟合，以求取试验区高程异常的拟合函数，再恢复地形起伏对GPS控制点高程异常的影响值，即采用移除恢复技术求取GPS检查点的高程，以公式表示：然后，再求GPS检查点的实测正常高与拟合高程的不符值，最后通过统计这些不符值来计算其高程的拟合精度，从中便可以发现采用哪一种形式的参考高程面对拟合结果更有利一些。试验选取了两个高山区，其中：（1） 试验区QL，（DEM参数见3.1节），平均高程为1 500m，布设了8个GPS公共控制点；（2） 试验区JL，DEM为360372格网面积，格网间距250m，平均高程为2 500m。该试验区海拔较高，最低处为

24、1 765m，最高处为3 100m，高差1 335m，布设了12个GPS公共控制点。在两个试验区试验，首先选择不同的参考面高程，如表1第2行所列，分别计算GPS点的高程异常，然后以其中4个GPS控制点作为拟合计算点，用二次多项式进行最小二乘高程拟合，其余的GPS点作为检查点。由移除恢复技术获得的检查点高程与实测高程的不符值，可得到每一个参考面的统计精度，结果列于表1，实验数据显示当参考面高程为平均高程面时，拟合效果最好。表 1 检查点高程不符值的统计精度试验区试验区QL试验区JL参考面高程 /m01000150030000150025003500每公里高差中误差/cmkm-110.67.15.89.27.97.77.58.64 结束语通过实验验算和数据分析，可以得出以下几点结论：1) 利用地形起伏计算高程异常的连续积分公式在计算过程中，无需人工干预，避免了人为误差，利于编程自动计算；2) DEM的计算区域对计算高程异常的影响，会随着范围的增大而逐渐减小。在本次研究的算例中，平原和丘陵地区参与计算高程异常的DEM的范围，在测区周边各增10km左右即可；高山地区则需增加18km左右。因此，实用中可根据实际需要准备适当大小的DEM数据，以免由于盲目增加大面积DEM数据而使得测量成本增加；3) 计算各点高程异