李圣明，史俊波，李 林，董新瑩

（1.中交第二公路勘察設計研究院有限公司，湖北 武漢 430050；2.武漢大學 測繪學院，湖北 武漢 430079）

經典GNSS高程擬合基于“移去—擬合—恢復”的思想[1]，引入EGM2008地球重力場模型和DEM格網數(shù)據，分別確定長波項和短波項，最后選擇合適的擬合模型計算高程異常殘余項。多項式擬合法是最常用的高程擬合模型之一，構建多項式擬合模型后，利用已知點真值，基于最小二乘準則確定未知參數(shù)，外推得到未知點的高程異常大小。唐詩華[2]等系統(tǒng)地研究了不同估計方法的原理和特點，得出最小二乘估計在高程擬合中精度最高的結論；陳安平[3]等比較分析了二次曲線擬合法和三次曲線擬合法等模型的適用性和可靠性；王雪林[4]和楊坤[5]等通過實際測量數(shù)據證明了多項式曲線擬合法可滿足四等水準測量要求，但涉及的測區(qū)均處于地勢平坦的區(qū)域，相對高差較小。多項式擬合法模型簡單、應用廣泛，但不適用于起伏較大的測區(qū)。因此，針對大高差區(qū)域，本文利用反距離加權法和支持向量機（SVM）模型進行高程異常殘余項的擬合；再通過實際算例，比較分析了不同模型的擬合效果，給出大高差區(qū)域最優(yōu)GNSS高程異常擬合模型的建議。

1 研究原理與方法

1.1 多項式擬合模型

多項式擬合模型是將呈帶狀分布的聯(lián)測點視為一個曲面，認為曲面上各控制點的高程異常與其坐標xi（或yi）之間存在函數(shù)關系。本文以xi為自變量，高程異常模型可表示為：

式中，a0,a1,…,a m為待定系數(shù)。

對于每個(xi,yi,ζi)，聯(lián)立上述方程進行擬合后，再利用最小二乘法求解式（1）中的各項系數(shù)并代回，即可求出曲面上任意一點的高程異常值。

1.2 反距離加權模型

反距離加權法[6]常見于散點數(shù)據的插值應用中，適用于散點分布較均勻的情況。待求點P處的高程異常插值結果ζP，即為測區(qū)內所有已知點高程異常值ζi(i=1,2,…,n)的加權和。

式中，di為第i個已知點到待求點P距離；u決定距離增加時權重減小程度。

1.3 SVM模型

SVM[7]模型，即某區(qū)域的N個訓練樣本點可表示為{(xk,y k)|k=1,2,…,N}，其中xk∈Rn為n維兩列輸入向量，兩列分別為坐標x和y；y k∈R為n維1列輸出向量，代表已知點的高程異常殘余項。利用非線性映射φ(x k):Rn→Rnk將輸入空間的非線性回歸問題轉換為高維特征空間中的線性回歸問題，可表示為求解以下約束優(yōu)化問題：

式中，w∈Rnk為權向量；γ為調整參數(shù)；b為偏差量；e k∈R為誤差變量。

2 研究結果與分析

本文選取兩個算例驗證上述高程異常擬合模型。算例1：數(shù)據點分布在東西5.35 km、南北15.19 km的范圍內，最高最低海拔差距約為445.9 m，將其中9個聯(lián)測點作為已知點，28個聯(lián)測點作為檢核點。算例2：數(shù)據點分布在東西20.22 km、南北12.95 km的范圍內，最高最低海拔差距約為705.4 m，將其中9個聯(lián)測點作為已知點，38個聯(lián)測點作為檢核點。算例點位分布情況如圖1所示。

圖1 算例1和算例2的點位分布

根據參與建模的已知點的高程異常真值ζi和計算得到的高程異常值ζiY，求得殘差Vi Y=ζi-ζiY，作為內符合精度。根據參與檢核的各點高程異常真值ζi和計算得到的高程異常值ζiC，求得殘差Vi C=ζi-ζiC，作為外符合精度。利用多項式擬合法對算例1和算例2的9個已知點進行高程異常殘余項建模，并統(tǒng)計內符合殘差，如圖2、3所示。

圖2 算例1多項式擬合法建模內符合殘差

圖3 算例2多項式擬合法建模內符合殘差

利用多項式擬合法、反距離加權法和SVM模型對算例1的28個檢核點和算例2的38個檢核點進行檢驗計算，外符合殘差對比如圖4、5所示。

圖5 3種模型外符合殘差對比（算例2）

兩個算例的區(qū)域特征以及檢核點采用3種模型計算得到的外符合殘差平均值如表1所示，繪出的對比圖如圖6所示，可以看出，多項式擬合法總體上能達到高程擬合的精度要求，但精度略差于另外兩種模型；反距離加權法和SVM模型體現(xiàn)出良好的適用性，外符合殘差的差距分別為0.52 cm和0.27 cm；綜合兩個算例結果判斷，在地勢起伏比較大的崎嶇地區(qū)，SVM模型精度優(yōu)于反距離加權法。

圖6 3種模型精度對比圖

表1 高程擬合算例綜合分析

3 結語

針對高程轉換時常用的多項式擬合法在大高差區(qū)域適用性較差的問題，本文將反距離加權法和SVM模型替代了多項式擬合法；并通過兩個最大高差分別為445.9 m和705.4 m的實際算例比較分析了3種模型在大高差區(qū)域對于高程異常殘差項的擬合效果。結果表明，利用SVM模型的外符合精度分別為1.73 cm和1.13 cm，優(yōu)于反距離加權法的2.25 cm和1.40 cm和多項式擬合法的3.17 cm和2.18 cm。因此，SVM模型可作為大高差區(qū)域GNSS高程異常擬合模型的最優(yōu)選擇。