陳福金

（珠海市測繪院，廣東 珠海 519000）

1 引言

近年來，GPS 技術(shù)已廣泛應(yīng)用于各行各業(yè)。在測繪行業(yè)，采用GPS 技術(shù)進行平面高程測量具有效率高、成本低及精度高的優(yōu)點[1]。但因受制于大地水準面誤差及布設(shè)平面控制網(wǎng)的高程異常值，GPS 測量技術(shù)常常會有一定的精度誤差，且在不同的地形測量區(qū)域，不同的GPS 高程轉(zhuǎn)換模型結(jié)果會有較大差異。常用的GPS 高程轉(zhuǎn)換模型主要有線性擬合、曲面擬合等[2]。結(jié)合實際地形，選用合理的GPS 水準高程轉(zhuǎn)換模型[3-4]，可以保證測量結(jié)果精度。因而，研究GPS 水 準 高程轉(zhuǎn)換模型在工程測量領(lǐng)域的應(yīng)用具有重要意義。

2 工程概況

工程位于某市工業(yè)園區(qū)建設(shè)開發(fā)用地，距城區(qū)4km，總面積14.25 平方公里，建成后可滿足8 萬人口入住。項目區(qū)域平面衛(wèi)星地形圖如圖1 所示。項目區(qū)域地勢較為平坦，前期已布置工程三、四等GPS 平面控制網(wǎng)，總計有28 個控制點，其中三等控制點11 個、四等控制點17 個。三等GPS 平面控制網(wǎng)利用GPS 接收信號機開展GPS 高程測量，采用1954 北京坐標系；四等平面控制網(wǎng)同步觀測，結(jié)合部分三等控制點，獲得項目區(qū)域測量點數(shù)據(jù)，將其作為工程數(shù)據(jù)校核。

圖1 項目區(qū)域衛(wèi)星地形圖

采用平差計算軟件計算出水準網(wǎng)內(nèi)各個控制點的平差值，保證平差值的中誤差為±4.7mm，各個數(shù)據(jù)點精度符合規(guī)范要求，可采用這一系列數(shù)據(jù)點作為高程擬合模型的基礎(chǔ)點。

3 GPS 水準高程擬合模型的應(yīng)用分析

根據(jù)28 個控制點的數(shù)據(jù)，計算出各個控制點的高程異常值，并繪出高程異常值趨勢圖（如圖2 所示）。從圖2 可以看出，測量區(qū)域28 個控制點的高程異常值范圍為11.33 ～11.79m。

圖2 高程異常值

獲得28 個控制點數(shù)據(jù)后，選定平面擬合、二次曲面、三次曲面、距離加權(quán)、移動二次曲面、抗差二次曲面、移動抗差二次曲面及多面函數(shù)模型共8 個模型進行GPS 水準高程擬合，并利用Matlab 編程技術(shù)計算出各個模型的精度參數(shù)值，綜合評估確定適合項目區(qū)域的擬合模型。

3.1 內(nèi)符合精度

內(nèi)符合精度包含擬合模型計算的中誤差、最大殘差、最小殘差。8 個模型的內(nèi)符合精度如表1 所示。

表1 各個模型精度參數(shù)表

根據(jù)表1 可知，多面函數(shù)模型的中誤差最大，為0.183m；二次曲面的中誤差最小，僅有0.0133m。多面函數(shù)模型的最大殘差值最大，為0.36m；二次曲面的最大殘差值最小，僅有0.0317m。距離加權(quán)模型的最小殘差值最大，為0.006m；二次曲面的最小殘差值最小，僅有0.0003m。綜合中誤差、最大殘差及最小殘差來看，二次曲面模型精度較高，適用于工程區(qū)域GPS 水準高程轉(zhuǎn)換。

3.2 殘差分布統(tǒng)計

各個模型殘差分布差異較大，根據(jù)不同模型下所有測量控制點的殘差值分布統(tǒng)計（如表2 所示），了解8個轉(zhuǎn)換模型的殘差分布情況，確定轉(zhuǎn)換模型精度較高者。

表2 各個模型殘差分布統(tǒng)計

從殘差值分布來看，殘差＞5cm 的控制點最多的是多面函數(shù)模型，達19 個，其中殘差超過4cm 的控制點有20 個，占比超過71.4%。殘差在2cm 以內(nèi)的控制點最多的是移動二次曲面，達26 個；二次曲面模型中殘差在2cm 以內(nèi)的控制點有24 個，占比超過85.7%。在GPS 水準高程擬合過程中，高程異常值在一定程度上有未知性，因而采用多面函數(shù)模型進行擬合轉(zhuǎn)換。當控制點分布不均勻時，可能會出現(xiàn)“插值振蕩”，造成多面函數(shù)模型殘差值分布出現(xiàn)急劇增大或減小的情況。

根據(jù)殘差分布統(tǒng)計，繪出各個模型的殘差分布直方圖（如圖3 所示）。從圖中可看出，在殘差低于1cm的控制點中，除距離加權(quán)與多面函數(shù)模型外，其他6 個模型擬合精度相符。其中二次曲面模型殘差分布更為穩(wěn)定，殘差在3cm 以內(nèi)的控制點占比超過90%。經(jīng)內(nèi)業(yè)復核證明，中誤差絕對值僅有1.36cm，滿足《城市測量規(guī)范》起算點要求。因而，在這個8 個轉(zhuǎn)換模型中，采用二次曲面模型進行GPS 水準高程轉(zhuǎn)換較為合適。

圖3 殘差分布直方圖

3.3 外業(yè)精度檢測

為確保選用的GPS 水準高程模型的合理性，采用外業(yè)檢測手段進行校核，獲得外業(yè)精度特征值，分析GPS 水準高程模型的適用性。為保證外業(yè)精度的獨立性，選用的控制點要與原測量區(qū)域的控制點不同，且重新采用一套設(shè)備及測量人員進行外業(yè)精度檢測，本項目采用了江西省某測繪院工作人員及GPS 設(shè)備。

選用6 個檢測點與4 個未知點獨立測量，構(gòu)成同步觀測GPS 平面網(wǎng)，檢測點與未知點的分布如圖4 所示。按照四等水準測量要求，結(jié)合GPS 觀測數(shù)據(jù)與檢測點水準測量數(shù)據(jù)，獲得水準測量觀測成果（如表3 所示）。

圖4 外業(yè)精度檢測點分布圖

表3 水準測量觀測成果

6 個外業(yè)精度檢測點的成果如表4 所示。從表4 可知，二次曲面模型轉(zhuǎn)換后的高程誤差較低，最大值為0.042m。根據(jù)公式（1）～（3）計算出外業(yè)精度檢測的誤差容許值為±5.44cm，與二次曲面模型擬合誤差最大值相比，最大差值小于容許值，表明外業(yè)檢測結(jié)果符合相關(guān)要求。為進一步得到二次曲面模型的表達式，采用Matlab 編程算法對參數(shù)進行反演，獲得參數(shù)最優(yōu)結(jié)果（如表5 所示），二次曲面函數(shù)為公式（3）。經(jīng)過對各個模型的內(nèi)符合精度驗算及外業(yè)精度檢測，可以得出結(jié)論：在項目區(qū)域，采用二次曲面模型進行GPS 水準高程轉(zhuǎn)換科學合理。

表4 外業(yè)精度檢測表

表5 二次曲面函數(shù)各參數(shù)值

4 結(jié)論

通過對項目區(qū)域8 個GPS 水準高程模型的內(nèi)符合精度參數(shù)對比分析，得出結(jié)論 ：二次曲面模型的中誤差、最大殘差及最小殘差值均為最小，二次曲面模型殘差分布范圍更為穩(wěn)定，殘差值在3cm 以內(nèi)的控制點占比超過90%?；谕鈽I(yè)精度檢測，獲得二次曲面模型的最大誤差值為0.042m，低于容許值，最終確定項目區(qū)域采用二次曲面模型進行高程轉(zhuǎn)換科學合理。