GPS高程擬合方法研究及精度對比試驗

張　潘,余代俊,張玉剛,許　馨

(1. 成都理工大學現(xiàn)代測量技術及應用研究所，四川 成都 610059； 2. 西藏自治區(qū)

測繪院，西藏 拉薩 850000)

ZHANG Pan,YU Daijun,ZHANG Yugang,XU Xin

GPS高程擬合方法研究及精度對比試驗

張潘1,2,余代俊1,張玉剛2,許馨2

(1. 成都理工大學現(xiàn)代測量技術及應用研究所，四川 成都 610059； 2. 西藏自治區(qū)

測繪院，西藏 拉薩 850000)

Research on GPS Elevation Fitting Method and Accuracy Comparison Experiment

ZHANG Pan,YU Daijun,ZHANG Yugang,XU Xin

摘要：利用GPS高程擬合技術求得滿足一定精度的正常高是測繪研究的重要問題之一。本文通過某試驗區(qū)GPS控制網(wǎng)數(shù)據(jù)成果，采用加權平均法、二次曲面擬合法求取高程異常值，并將最終成果與三等水準測量成果進行了對比分析得出，當已知點呈區(qū)域面分布時，二次曲面擬合法擬合精度優(yōu)于加權平均擬合法，且擬合精度達毫米級，證明了該模型在實際運用中的優(yōu)越性。

引文格式： 張潘,余代俊,張玉剛,等. GPS高程擬合方法研究及精度對比試驗[J].測繪通報，2015(9)：54-56.DOI:10.13474/j.cnki.11-2246.2015.0278

關鍵詞：高程擬合；加權平均法；二次曲面擬合；高程異常

中圖分類號：P228.4

文獻標識碼：B

文章編號：0494-0911(2015)09-0054-03

收稿日期：2014-07-30

作者簡介：張潘(1987—)，男，碩士，助理工程師，研究方向為攝影測量與遙感。E-mail:980506099@163.com

一、引言

GPS測量可以獲得高精度的三維坐標，它的平面相對定位精度已經(jīng)十分精確，但是它得到的高精度的大地高與我國采用的正常高系統(tǒng)不一致，只具有幾何意義，如何將大地高轉換為正常高，在一定程度上代替?zhèn)鹘y(tǒng)的水準測量作業(yè)，是測繪行業(yè)研究的熱點。

我國采用正常高系統(tǒng)，大地高起算面為參考橢球面，正高起算面為大地水準面，正常高系統(tǒng)起算面為似大地水準面。三者關系如圖1所示。

圖1　大地高、正高和正常高關系

(1)

二、GPS擬合模型

1. 加權平均法

加權平均法，是由內插點周圍部分已知點的高程異常加權平均求得該點的高程異常。設在內插點周圍選n個已知點，高程異常為ξi(i=1，2，…，n)，對應的權為Pi，則內插點j的高程異常為

權Pi可根據(jù)已知點至內插點的距離來計算

式中，di為已知點i至內插點j的水平距離；ε為一小正數(shù)，以防止權函數(shù)的分母趨于0，通常ε取0.01，單位與di單位相同。

2. 二次曲面擬合

曲面擬合法是認為高程異常在一定范圍內變化平緩的前提下，將高程異常(或高程異常差)近似地看作是一定范圍內各點坐標的曲面函數(shù)。設測站點的高程異常與其平面坐標的函數(shù)關系式為

ξi=f(xi,yi)+vi

式中，f(xi,yi)為所選定的曲面函數(shù)，也即ξi的趨勢值；vi為其殘差。當選擇多項式曲面擬合法時

式中，bi(i=0，1，2，…，5)為待定系數(shù)；(xi，yi)為i點的平面坐標

若有n個點，則可以得到下面的矩陣形式

ξ=XB+v

用最小二乘求出bi(i=0，1，2，…，5)6個待定參數(shù)，求出參數(shù)后根據(jù)擬合模型得出高程異常值ξi，并由式(1)得出正常高。

三、試驗分析

1. 控制網(wǎng)概況

試驗區(qū)結合GPS控制網(wǎng)及三等水準進行研究，控制網(wǎng)選取11個點作為已知點，按GPS D級網(wǎng)進行觀測(如圖2所示)，觀測儀器采用Trimble R8靜態(tài)接收機，觀測時段數(shù)大于1.6，每時段時長大于60 min，測量開始與結束時各量取3次儀器高，儀器高前后不超過5 mm。

圖2　GPS觀測網(wǎng)

2. 高程擬合試驗

此次試驗采用加權平均法及二次曲面擬合，其中加權平均法選取方案1為6點(DS05、DS03、DS02、GD07、GD08、GD06)、方案2為7點(DS05、DS03、DS02、GD07、GD08、GD06、GD05)、方案3為6點(DS05、DS03、GD07、GD08、GD06、GD05)，其余點作為檢查點，擬合結果見表1。

二次曲面擬合選取方案1為6點(DS05、DS03、DS02、GD07、GD08、GD06)、方案2為6點(DS05、DS03、GD07、GD08、GD06、GD05)，其余點作為檢查點，擬合結果見表2。

表1　加權平均法檢查點殘差　mm

表2　曲面擬合法檢查點殘差　mm

3. GPS水準精度評定

(1) 內符合精度

式中，n為參與擬合計算的已知GPS水準聯(lián)測點數(shù)。

(2) 外符合精度

式中，n為參與檢核的GPS水準聯(lián)測點數(shù)。

根據(jù)表1、表2計算試驗區(qū)的內符合精度及外符合精度，計算結果見表3。

表3　兩種模型精度對比　mm

四、結論

1) 當測區(qū)已知點呈面狀分布時，加權平均法內符合精度與外符合精度都遠大于二次曲面擬合的符合精度，并通過與三等水準成果的對比，得出二次曲面擬合精度優(yōu)于加權平均法。

2) 二次曲面擬合適用于地形較復雜地區(qū)的高程擬合，擬合時已知點數(shù)量不少于6個，同時在擬合時已知點點位分布對擬合精度影響較大。當已知點包含整個測區(qū)時，擬合精度較穩(wěn)定；當待定點在已知點覆蓋范圍之外時，擬合精度較低。

3) GPS高程擬合精度取決于地形起伏變化均勻程度，用數(shù)學模型求解某一點高程異常時，存在一定的誤差，因此選擇某種模型進行高程擬合時，可根據(jù)測區(qū)實際情況選擇擬合模型。

參考文獻：

[1]吉淵明, 趙水泉. 曲面擬合法求GPS網(wǎng)正常高的幾點認識[J]. 測繪通報, 1998(7): 18-20.

[2]徐紹銓, 張華海, 楊志強, 等. GPS測量原理及應用[M]. 武漢: 武漢測繪科技大學出版社, 1998.

[3]劉基余. 全球定位系統(tǒng)原理及其應用[M]. 北京: 測繪出版社, 1993.

[4]王巖, 劉茂華, 黨永超. 小地區(qū)GPS高程擬合的方法研究與實施[J]. 測繪通報, 2012(S0): 66-68.

[5]沈學標. GPS水準高程擬合精度的分析[J]. 測繪通報, 1998(7): 21-22.

[6]閆偉, 高俊強, 王維. 小地區(qū)GPS高程擬合和水準測量對比試驗[J]. 南京工業(yè)大學學報, 2007, 29(5): 93-96.

[7]何美琳, 文鴻雁, 潘元進, 等. GPS高程擬合的方法比較[J]. 測繪科學, 2013, 38(3): 63-65.

[8]國家測繪局. GB/T 18314—2001 全球定位系統(tǒng)(GPS)測量規(guī)范[S]. 北京: 中國標準出版社, 2001.

[9]張興福,劉成.綜合EGM2008模型和SRTM/DTM2006.0剩余地形模型的GPS高程轉換方法[J]. 測繪學報，2012， 41 (1): 25-32.