趙會會

摘要：在進(jìn)行GPS水準(zhǔn)高程測量時，水準(zhǔn)聯(lián)測點(diǎn)對高程擬合有較大的影響，本文首先對水準(zhǔn)聯(lián)合測點(diǎn)的優(yōu)化選擇問題進(jìn)行了分析，并對擬合誤差最小目標(biāo)函數(shù)水準(zhǔn)聯(lián)測點(diǎn)的優(yōu)化措施進(jìn)行了研究，提出了GPS水準(zhǔn)聯(lián)測點(diǎn)的兩種方法，最后經(jīng)過計(jì)算證明，使用這兩種方法可得到相同的GPS水準(zhǔn)點(diǎn)，并得到了最小擬合誤差。

關(guān)鍵詞：GPS;水準(zhǔn)擬合;優(yōu)化

GPS水準(zhǔn)擬合的應(yīng)用很廣，大到似大地水準(zhǔn)面的精化，小到一個工程網(wǎng)中用來代替三、四等水準(zhǔn)測量。而GPS水準(zhǔn)擬合的模型也因區(qū)域大小而異，如全國或省市級似大地水準(zhǔn)面精化時可用多面函數(shù)、樣條函數(shù)或調(diào)和函數(shù)等;在小范圍工程測量中可用多面函數(shù)和多項(xiàng)式等;當(dāng)高程控制網(wǎng)布設(shè)成線狀或帶狀時（如道路等），也可采用曲線形式來擬合;在測區(qū)面積較大，高程異常變化復(fù)雜的區(qū)域內(nèi)，可采用分區(qū)擬合的方法。

1 GPS水準(zhǔn)擬合的原理

GPS測量得到的大地高是相對于WGS-84橢球的，而正常高是以似大地水準(zhǔn)面為高程起算面的，二者的其轉(zhuǎn)換關(guān)系為：Hr=H+ξ（1）式中Hr為正常高;H為大地高;ξ為高程異常。由上式可知，由大地高轉(zhuǎn)化得到高精度的正常高的關(guān)鍵是求定高精度的高程異常，GPS水準(zhǔn)擬合是常用的方法之一。GPS水準(zhǔn)擬合就是通過GPS與水準(zhǔn)點(diǎn)的重合求出重合點(diǎn)的高程異常，運(yùn)用重合點(diǎn)的高程異常擬合出某一區(qū)域的似大地水準(zhǔn)面，經(jīng)過坐標(biāo)內(nèi)插得到該區(qū)域任意點(diǎn)的高程異常，從而得到其他GPS點(diǎn)的正常高。對于似大地水準(zhǔn)面的擬合，為了不失一般性，也顧及下面的實(shí)例，這里僅采用二次曲面函數(shù)作為擬合函數(shù)。

2 GPS水準(zhǔn)擬合模型

2.1曲線擬合法

當(dāng)GPS點(diǎn)呈線狀布設(shè)（未構(gòu)成網(wǎng)狀，如鐵路、公路的初測導(dǎo)線等），在認(rèn)定沿線似大地水準(zhǔn)面為一條連續(xù)而光滑的曲線的前提下，可應(yīng)用線性擬合法，求待定點(diǎn)的正常高。其原理是：根據(jù)高程控制點(diǎn)的平面坐標(biāo)及其高程異常值，通過構(gòu)造一個插值函數(shù)來擬合測線方向上的似大地水準(zhǔn)面曲線，然后據(jù)此內(nèi)插其它點(diǎn)的高程異常。以下是本文研究的幾種典型的曲線擬合法。

2.2逐點(diǎn)剔除法

這里介紹沈云中教授提出的逐點(diǎn)剔除法，逐點(diǎn)剔除法是與全組合方法相對而提出的，因?yàn)樵贕PS水準(zhǔn)優(yōu)化中關(guān)于最佳重合點(diǎn)數(shù)及其分布的確定是要在某一準(zhǔn)則的基礎(chǔ)上進(jìn)行逐一比較才能確定的。假設(shè)某一GPS網(wǎng)中聯(lián)測了n個GPS水準(zhǔn)點(diǎn)，現(xiàn)選取m個點(diǎn)參與擬合計(jì)算，客觀上存在一組m個點(diǎn)的組合使最終擬合計(jì)算的結(jié)果最優(yōu)，在此我們稱這m個GPS水準(zhǔn)點(diǎn)為最佳擬合方案。全組合法就是將每組擬合方案的擬合結(jié)果都算出來，精度最高的一組即為最優(yōu)方案，實(shí)際上這種算法往往是不可能實(shí)現(xiàn)的，而逐點(diǎn)剔除法的基本思想是從n個重合點(diǎn)中逐點(diǎn)剔除對擬合精度貢獻(xiàn)最小的點(diǎn)，最后剩下的m個重合點(diǎn)就是GPS水準(zhǔn)擬合計(jì)算的最優(yōu)擬合方案。從n個點(diǎn)中剔除一點(diǎn)，共有n種剔除方案，若其中某種方案求得的擬合結(jié)果精度最高，則采用該方案，剔除該方案所剔除的點(diǎn)，使得總點(diǎn)數(shù)降為n-1。依此類推，直到剩下m個點(diǎn)，即為最優(yōu)的GPS水準(zhǔn)點(diǎn)。該算法總共需要計(jì)算n+（n-1）+Λ+（m+1）種方案，這是普通計(jì)算機(jī)可以勝任的。接下來的問題就是如何確定重合點(diǎn)數(shù)m，只需將逐點(diǎn)剔除法略加拓展，因?yàn)楫?dāng)給定某一擬合函數(shù)后，為能采用參數(shù)平差法進(jìn)行求解，最少的重合點(diǎn)數(shù)是已經(jīng)確定的，如本文采用二次曲面函數(shù)，最少重合點(diǎn)數(shù)為6，因此只需采用逐點(diǎn)剔除法將重合總點(diǎn)數(shù)n逐點(diǎn)剔除直至為6，然后比較重合點(diǎn)為6～n中精度最高的即為最佳重合點(diǎn)數(shù)，相應(yīng)點(diǎn)的分布即為最佳分布。

2.3三次樣條曲線擬合法

當(dāng)測線長，己知點(diǎn)ζi多變化較大時，為避免高次插值的振蕩現(xiàn)象，同時又保證分段低次插值連接點(diǎn)上的光滑性，為此，通常采取分段計(jì)算，以三次樣條函數(shù)作為擬合模型。

2.4多項(xiàng)式曲面擬合

多項(xiàng)式曲面擬合法是GPS水準(zhǔn)中最常用的幾何方法之一。多項(xiàng)式曲面擬合法是認(rèn)為高程異常在一定范圍內(nèi)變化平緩的前提下，將高程異常（或高程異常差）近似地看作是一定范圍內(nèi)各點(diǎn)坐標(biāo)的曲面函數(shù)，用這一擬合函數(shù)來計(jì)算其它GPS點(diǎn)的高程異常和正常高。

2.5解析法

在對測區(qū)中大地水準(zhǔn)面進(jìn)行擬合時，解析法主要是使用數(shù)學(xué)模型來進(jìn)行擬合的，當(dāng)待測點(diǎn)的大地坐標(biāo)是已知的，通過使用數(shù)學(xué)模型就可將高程異常精度計(jì)算出來。其中精度主要是由測區(qū)大地水準(zhǔn)面和數(shù)學(xué)模型的擬合度來決定的。

3案例分析

3.1控制網(wǎng)和觀測

某地形平緩區(qū)域GPS控制網(wǎng)的平均變長為1km，見圖1。選用ASHTECHZ12GPS定位儀對各個控制點(diǎn)進(jìn)行觀測，觀測時間持續(xù)90min，每間隔20s進(jìn)行一次采樣，圖形強(qiáng)度因子GDOP<4。在進(jìn)行研究時，采用Ⅱ級標(biāo)準(zhǔn)對控制點(diǎn)進(jìn)行測量，使用精密星歷對GPS網(wǎng)單獨(dú)水平差進(jìn)行測量，平面中最弱點(diǎn)的平面點(diǎn)位誤差要控制在3mm以內(nèi)，經(jīng)過計(jì)算證明，差異的最大值為2.5mm以內(nèi)，GPS定位精確度比較高。

3.2分析擬合模型

分別建模對差值模型和二次模型進(jìn)行計(jì)算后證明，受粗差的影響，擬合誤差在23.71mm內(nèi)，模型精度較小。進(jìn)行粗差剔除，并檢測其顯著性后，只需要7個點(diǎn)就可建立差值模型，模型擬合誤差在0.55mm內(nèi)，測量準(zhǔn)確度顯著提升。利用此模型可將各點(diǎn)高程推求出來，相較于Ⅱ等水準(zhǔn)高程，中誤差為1.26mm，達(dá)到了Ⅱ等水準(zhǔn)的基本要求。

4結(jié)束語

從理論分析及實(shí)例計(jì)算可以看出，逐點(diǎn)剔除法在GPS水準(zhǔn)擬合中重合點(diǎn)數(shù)及其分布的確定方面是一種行之有效的辦法，可以改變傳統(tǒng)上憑經(jīng)驗(yàn)選取重合點(diǎn)的盲目性。當(dāng)然，GPS水準(zhǔn)擬合的優(yōu)化還需兼顧模型選取以及評定標(biāo)準(zhǔn)的合理判定等因素，另外應(yīng)對采取逐點(diǎn)剔除法所得到的最優(yōu)重合點(diǎn)進(jìn)行進(jìn)一步分析，看是否與地形的特征點(diǎn)相一致，這樣使得GPS水準(zhǔn)優(yōu)化不僅具有數(shù)學(xué)意義，而且兼顧物理特征，更具實(shí)際指導(dǎo)價(jià)值。

參考文獻(xiàn)：

[1]李明峰，馮寶紅，劉三枝.GPS定位技術(shù)及其應(yīng)用[M].北京：國防工業(yè)出版社，2006.

[2]簡程航.GPS高程擬合方法研究及其工程應(yīng)用[D].北京：中國地質(zhì)大學(xué)，2014.

（作者單位：河北省水利水電勘測設(shè)計(jì)研究院）