GPS高程擬合方法

1、精選優(yōu)質(zhì)文檔-傾情為你奉上GPS高程擬合方法3.1等值線圖示法等值線圖示法是最直接的求算高程異常的方法。這種方法的核心思想就是內(nèi)插的思想，繪制高程異常的等值線圖，然后采用內(nèi)插法來確定未知點的高程異常值。具體操作十分的簡單，在測區(qū)內(nèi)制定分布均勻的GPS點，用水準測量的方法來測定這些點的水準高，根據(jù)公式=H-Hr求出這些點的高程異常，選擇適當?shù)谋壤甙凑找阎c的平面坐標展會在圖紙內(nèi)，對已知點標注出高程異常值，再確定等高距，繪制出高程異常值的等值線圖。之后就可以內(nèi)插出待測點的高程異常值，進而求出待測點的正常高。這種方法只適用地形相對平坦的地方，在此種測區(qū)內(nèi)采用這種方法擬合的高程精度可達到厘米級。測區(qū)

2、的地形相對復雜內(nèi)插出的高程異常值就不準確，而且這種內(nèi)插法的精度往往取決于兩個方面，分別是測區(qū)內(nèi)GPS點的分布密度和已知點大地高的精確度。首先GPS點的分布比較密集，那么內(nèi)插精度就相對較高，如果比較稀疏這時候就要借助于此測區(qū)的重力測量資料，提高內(nèi)插精度。且還要注意GPS點間高程異常的非線性變化。另外就是水準點的精度，聯(lián)測時盡量選取高精度的正常高，盡可能使得出的高程異常值準確，進而才能內(nèi)插出待測點高精度的高程異常值。這種方法雖然簡單易操作，但是有其弱點，就是精度不高，只有當對擬合精度要求不高的時候才使用此種方法(注：等值線法不需構造數(shù)學模型)。3.2狹長帶狀區(qū)域線性擬合解析內(nèi)插法作為擬合高程最常用

3、的方法，主要思想是把似大地水準面用數(shù)學曲面近似擬合，建立所在測區(qū)內(nèi)最為接近似大地水準面的數(shù)學模型，以此來計算測區(qū)內(nèi)任意點的高程異常值，從而計算出正常高。這種方法計算出的高程異常值的精度是由所采用的數(shù)學模型和似大地水準面的擬合程度所決定的。解析內(nèi)插法在選擇數(shù)學模型時，首先要考慮的就是GPS點的分布情況。GPS點的分布情況可分為帶狀分布和面狀分布。若GPS點是呈線狀布設，而且是以沿線似大地水準面為一條連續(xù)且光滑的曲線，這時就可以采用相對于狹長帶狀區(qū)域的解析內(nèi)插法來內(nèi)插出待定點的高程異常值，從而求出待定點的正常高。這種線狀分布的內(nèi)插原理是：測區(qū)內(nèi)已知水準點，用GPS測出其GPS高程，計算出已知水準點

4、高程異常值，根據(jù)已知點的平面坐標和計算得出的高程異常值，構造出一個插值函數(shù)，這個函數(shù)是用來擬合GPS分布線上的似大地水準面的。用這個函數(shù)內(nèi)插出位置點的高程異常值。下面是兩種用來擬合線狀分布的GPS高程的內(nèi)插法。3.2.1 多項式曲線擬合法多項式曲線擬合是線狀分布擬合的主要方法。多項式擬合顧名思義其插值函數(shù)是一個m次的代數(shù)多項式，若高程控制點的高程異常為，坐標為xl(或yi或di或擬合坐標或xi-x0或yi-y0)的函數(shù)關系為下式：l=a0+a1xl1+a2xl2+a3xl3+amxlm3-1各高程控制點的已知高程異常與其擬合值之差為下式所示：ri=l(x)-i(i=0，1，2n)3-2上式我們

5、稱之為離差。(3-1)中xl是擬合點到參考點(x0，y0)的直線距離，x0，y0為設定的常數(shù)值。在一般情況下都認為，x0，y0就是測區(qū)內(nèi)已知點坐標的均值。多項式曲線擬合使用起來非常方便，但是它有自身的局限性，即是使用這種方法的時候，所測路線不能太長，要限制控制點到測點的距離不能太遠，通常把距離控制在300米以內(nèi)。這個要求是因為使用多項式曲線方法擬合似大地水準面，如果它擬合的范圍太大，點位的高程異常變化就越復雜，削高補低的方法不能滿足我們所要求的精度。隨著多項式階數(shù)的增大，也會使擬合出的曲線振蕩的更厲害，從而造成擬合的誤差增大。這些造成了多項式曲線擬合的缺陷，但是在路線較短的情況下，這種方法有足

6、夠的精度來擬合GPS點的正常高程。在式(3-1)中用m次多項式擬合似大地水準面，這個m的值如何取定，一般情況下如果測區(qū)不是很長，地形相對平坦，那么我們通常取m取為3。也就是說多項式為三次多項式。若測區(qū)比較長或者是測區(qū)地形比較復雜就要依情況而定，增加多項式的次數(shù)，提高擬合精度。依上述分析m的取值主要和測區(qū)長度以及測區(qū)的復雜程度有關。3.2.2三次樣條曲線擬合法三次樣條曲線擬合法針對測線長，已知點多的測區(qū)GPS高程擬合問題。由上述可以知道，當測線比較長已知點較多的時候，就需要構造高次的擬合多項式，當m值比較高的情況下，會出現(xiàn)不穩(wěn)定的現(xiàn)象，對求解高程異常值會有比較大的影響，并且最小二乘法在求多項式系

7、數(shù)中也會增大削高補低的誤差，因此為了避免測線長、已知點多這種情況下所出現(xiàn)的問題，通常采用分段擬合的方法，采用三次樣條函數(shù)擬合數(shù)學模型。這種方法很好的解決了因測線長而引起的問題。三次樣條曲線的實質(zhì)就是一個拼接而成的連續(xù)函數(shù)，在把測線分為多段的情況下，每段設為三次多項式函數(shù)，然后將這些多項式函數(shù)組成三次樣條函數(shù)。為了計算準確，應用中要求這種構造的曲線不僅在連接點處函數(shù)要連續(xù)。而且還要求這個函數(shù)的一級導數(shù)還有二級導數(shù)全部要是連續(xù)的，才能保證在分段之后構造的三次樣條函數(shù)后期運算中能夠計算出準確的高程異常值。設過n個已知點，i和xi(或yi或擬合坐標)在區(qū)間xi,xi+1(i=1,2,n-1)上有三次樣

8、條函數(shù)關系:(x)=(xi)+(x-xi)(xi,xi+1)+(x-xi)(x- xi+1)(x,xi,xi+1)3-3式中，x為待定點坐標,(xi,xi+1)為一階差商,(xi,xi+1)=(i+1-i)/(xi+1-xi)；(x，xi,xi+1)為二階差商,(x,xi,xi+1)=16,(xi)+,(x)+,(xi+1),而,(x)(i=1,2,n-1)，滿足系數(shù)矩陣為對稱三角陣的線性方程組(xi-xi+1),(xi-1)+2(xi+1,xi-1),(xi)+(xi+1,xi),(xi+1)=6,(x,xi+1)-(x,xi)3-4(x0)=,(xn)=0用追趕法解上面方程組，可求出,(x

9、i)和(xi,xi+1)，而,(x)=,(xi)+(x-xi),(xi,xi+1) 3-5這種做法有諸多好處，其中優(yōu)點有三點：其一計算簡便，其二保留了多項式的優(yōu)點，其三克服了多項式的缺點。多項式的缺點是單個多項式會有不靈活不穩(wěn)定的現(xiàn)象。由于三次樣條曲線的種種優(yōu)點，往往在實際中當遇到測線長已知點多的情況下采用此方法擬合高程。3.3曲面擬合法曲面擬合法是用于GPS點的分布在一定區(qū)域的時候，且可以選擇數(shù)學曲面擬合該區(qū)域的似大地水準面，構造適當?shù)臄?shù)學模型，計算該區(qū)域內(nèi)的高程異常值，然后求出正常高。這種擬合法的主體思想和曲線擬合法異曲同工的。具體思想是：已知測區(qū)的若干已知水準點，并且用GPS測定這些點的

10、高程，利用公式求得這些點的高程異常，有了已知點的高程異常，已知點的平面坐標是已知的，所以利用其平面坐標(x，y)和高程異常值構造出來的數(shù)學模型擬合最為接近于該測區(qū)的似大地水準面，然后內(nèi)插出未知點的高程異常值，進而求出正常高。3.3.1多項式曲面擬合法測站點的大地高H與正常高h之間有如下關系：h=H-3-6多項式函數(shù)擬合法的基本思想是在小區(qū)域GPS網(wǎng)內(nèi)，將似大地水準面看成曲面(或平面)，將高程異常表示為平面坐標(x，y)的函數(shù)，通過網(wǎng)中起算點(既進行了GPS測量又進行了幾何水準聯(lián)測的點)已知的高程異常確定測區(qū)的似大地水準面形狀，求出其余各點的高程異常，然后根據(jù)式(3-6)求出其他點的正常高，其數(shù)

11、學模型為：=f(x，y)+3-7式中f(x，y)是擬合的似大地水準面；是擬合誤差，f(x，y)=a0+a1x +a2y+a3x2+a4xy+a5y2+3-8x=B-B0 y=L-L0B0=1nBL0=1nL其中：n為GPS網(wǎng)中點的數(shù)量，(B,L)為已知點的大地坐標，a0，a1，a2，a3，a4，a5為擬合待定參數(shù)；x，y為各GPS點的平面坐標的近似值，一般取起算點的平面坐標減去網(wǎng)中全部點平面坐標的均值。1二次曲面擬合取(3-8)式中的一、二次項后將大地水準面擬合為：f(x，y)=a0+a1x+a2y+a3x2+a4xy+a5y23-9即得二次曲面擬合模型：=a0a1a2a3a4a5 1 x y

12、 x2 x y y2 T+3-10每一個起算點可以組成一個上式，若共存在m個這樣的起算點，則可列出m個方程：1=a0+a1x1+a2y1+a3x12+a4x1y1+a5y11+12=a0+a1x2+a2y2+a3x22+a4x2y2+a5y21+2m=a0+a1xm+a2ym+a3xm2+a4xmym+a5ym1+m3-11從而組成誤差方程：V=-Bx+L3-12上式中，B=1 x1 y1x12y12 x1y11 x2y2x22 y22 x2y2 1 xm ymxm2 ym2xmym，x=a0 a1 a2 a3 a4 a5T，L=1 2 mT，解得x=(BTPB)-1BTPL3-13解算出ai

13、即可求出網(wǎng)中其余點的高程異常，并利用式(3-6)求出各未知點的正常高h。2多項式平面擬合在小范圍或平原地區(qū)，可以認為大地水準面趨近于平面。此時，可選用公式（3-8）的前三項，將大地水準面擬合為：f(x，y)=a0+a1x+a2y 3-14擬合模型為：=a0a1a2 1 x y T+3-15其中，ai(i=0,1,2)為未知參數(shù)，此時要求公共點至少3個。3多項式相關平面擬合：也叫做四參數(shù)曲面擬合，若選用公式(3-8)的前三項和第五項進行擬合，則擬合曲面的表達式變?yōu)椋篺(x，y)=a0+a1x+a2y+a3xy 3-16擬合模型為：=a0a1a2a5 1 x y xy T+3-17其中，ai(i=

14、0,1,2,3)為未知參數(shù)，此時需要公共點至少4個。3.3.2移動曲面擬合法移動曲面擬合法是一種局部逼近法，其基本思想是以每一個內(nèi)差點為中心，利用內(nèi)差點周圍數(shù)據(jù)點的值，建立一個擬合曲面，使其到各個數(shù)據(jù)點的距離之加權平方和為極小，而這個曲面在內(nèi)插點上的值就是所求的內(nèi)插值。設P為內(nèi)插的點，下面對P構造相應的曲面。本文取如下的二次多項式曲面為例：f(x，y)=a0+a1x +a2y+a3x2+a4xy+a5y23-18設選取數(shù)據(jù)點的坐標為(xi,yi)，i=1,2,n；n6且設內(nèi)插點P的坐標為(xp,yp)。將(xi,yi)改化到以P為原點的局部坐標系中，即： xi,=xi-xp yi,=yi-yp

15、3-19形成新的坐標( xi,，yi,)，為移動坐標。任一點數(shù)據(jù)(xi,yi)假設距離d的遞減函數(shù)為：（d）=xi-xp2+yi-yp2 3-20將(d)作為權函數(shù)，對每個數(shù)據(jù)點賦予權i，這里i不是代表數(shù)據(jù)點的觀測精度，而是反映該點與內(nèi)插點的相關程度的大小。因此，權i確定的原則應該與該數(shù)據(jù)點和內(nèi)插點的距離di有關，di越小，它對內(nèi)插點的影響越大，則權應越大。相反，di越小，它對內(nèi)插點的影響越小，則權應越小。最后，由最小二乘法解如下帶權的極小值問題：Min丨i=1nivi2丨=R3-21為了給出二次多項式曲面，式(3-18)的系數(shù)，那么這時需要選取P點周圍的數(shù)據(jù)點。當點數(shù)不夠多時，則應擴大R的取

16、值?，F(xiàn)在這里由n個數(shù)據(jù)點的值，可得到如下的方程式：vi=a0+a1xi,+a2yi,+a3x,i2+a4xi,yi,+a5y,i2-fi(i=1，2，n) 3-22由此得系數(shù)ai(i=1,2,n)，從而得到所對應的二次曲面方程，進而得到所求內(nèi)插點的高程異常值。3.3.3多面函數(shù)曲面擬合法多面函數(shù)擬合法的本質(zhì)是數(shù)學曲面逼近的方法。其基本思想是用數(shù)學表面逼近所測區(qū)域的大地水準面，通常認為任何表面，無論這個表面是否是有規(guī)則的，都能通過一定的方法構造出來一個有規(guī)則的數(shù)學表面逼近其表面。通過構造數(shù)學表面，用數(shù)學表達式高精度的逼近并且代替其真實表面。也就是說每個插值點都可以和已知點建立起來相應的函數(shù)關系式

17、，然后將這些函數(shù)關系式迭加在一起，組成一個全新的函數(shù)關系式，那么稱這個迭加函數(shù)為多面函數(shù)，由于這是每個插值點與已知數(shù)據(jù)建立的函數(shù)關系，因此多面函數(shù)具有計算最佳擬合值的特點，正因如此，多面函數(shù)曲面擬合法就能夠更準確的擬合出未知點的高程值。多面函數(shù)的數(shù)學表達式為：=f(x，y)=i=1naiQ(x，y，xi，yi)3-23多面函數(shù)式中包含了待定系數(shù)ai，核函數(shù)Q(x，y，xi，yi)，其中核函數(shù)是關于x，y的函數(shù)，核函數(shù)的中心在(xi，yi)處。理論上講核函數(shù)是可以任意構造的，在實際應用中，通常用以下幾種函數(shù)來充當核函數(shù)。1錐面Q(x，y,xi,yi)=C+x-xi2+y-yi2123-242雙曲

18、面Q(x，y,xi,yi）=C+x-xi2+y-yi2+2123-253倒曲面Q(x,y,xi,yi)=C+x-xi2+y-yi2+2-123-264三次曲面Q(x，y,xi,yi）=C+x-xi2+y-yi2+2323-27在上述各式中，(x，y)表示內(nèi)插點坐標，(xi,yi)表示的是已知點的坐標，那么核函數(shù)中的x-xi2+y-yi2表示的是內(nèi)插點到已知點的水平距離，式中的參數(shù)為光滑系數(shù)。其具體求解過程為：以核函數(shù)為雙曲面為例，說明多項式曲面擬合法的具體求解過程，設測區(qū)內(nèi)的已知點個數(shù)為n，求解(3-23)中的系數(shù)(a1,a2,a3,an)，其矩陣形為下式所示：=XB3-28其中，=1 2 n

19、T，B=a0 a1 a2 a3 anT由此方程組可解得系數(shù)(a1,a2,a3,an)的唯一解：B=X-13-29求解未知點的高程異常值，根據(jù)公式(3-28)和(3-29)即可得到求解p=XPB=XPX-13-30根據(jù)以上求解過程可知，(3-30)式中的已知點的高程異常值直接關系到未知點的高程異常值的計算結(jié)果，因此，如果想要更好的解算出未知點的高程異常值，必須認真選取已知點，并且使所選的已知點的高程異常值相差比較大，因為這些點能最好的描述地形變化特征，即高程異常值的分布特征。這些特征點的選擇一般在地勢高和地勢較低的地方。在選擇多面函數(shù)求解測區(qū)內(nèi)的點的高程異常值的時候，需要注意的是以及核函數(shù)的選取

20、的問題，由于其取值是自主取值，為了能達到擬合最佳效果，就要逐步的試驗進而改進，然后選定一個最佳取值。3.4地球重力場模型擬合法所謂的地球重力場模型擬合法的關鍵是要收集相關的重力場信息，這些數(shù)據(jù)包括衛(wèi)星跟蹤數(shù)據(jù)、衛(wèi)星測高數(shù)據(jù)以及地球重力數(shù)據(jù)等。收集到足夠的數(shù)據(jù)后利用地球撓動位的球諧函數(shù)級數(shù)展開式求算測區(qū)內(nèi)點的高程異常值進而求得點位的正常高。求取測區(qū)內(nèi)一點P高程異常之前要先計算出該點的撓動位，由物理學知識可以知道地面一點P的撓動位的計算公式是由該點引力位V減去該點的正常引力位U求得的。其公式如下：T=V-U3-31那么地面點P的高程異常值就可以利用撓動位求得：=Tr3-32上式中，r為是P點的正常

21、重力值。由上述兩式中可以看出，如果要求出P點高程異常就要知道該點的引力位、正常引力位以及和正常重力值。其中正常重力值和正常引力位都是可以精確的計算出來的。如果P點的引力位知道，那么就可以計算出該點的高程異常值。提高引位的精度就是提高高程異常值的精度。這種方法的缺點就是要收集比較多的數(shù)據(jù)，有時候在測區(qū)內(nèi)會缺少某些數(shù)據(jù)，采用這種方法就會受到限制。而且這種方法的精度受到收集到的數(shù)據(jù)的精度的限制往往比不上前面所述方法的精度高。3.5地球重力場結(jié)合GPS水準擬合法從前面可以看出無論是GPS水準擬合法或是利用地球重力場計算高程異常值，其分別都會有優(yōu)點或缺點，在實際應用中，往往希望突出優(yōu)點避免缺點，所以如果

22、把兩種方法結(jié)合起來，這是一個提高高程異常值精度的新思路?？梢栽谝院蟮膶嵺`中應用。該方法的基本思路是：首先在已知水準點用GPS測出大地高，利用大高和正常高的差值求出高程異常，然后再利用地球重力場模型法求出已知點的高程異常m，兩種方法的高程異常值求出后，由于所用方法的不同，所以會有差值，計算出兩者的差值。如下所示：=-m3-33已知點的兩種高程異常的插值計算出后，用其平面坐標和差值在構造出的曲面數(shù)學模型中推算出未知點的，由于m是由地球重力場模型法求出的，所以就可以計算未知點的正常高。其計算公式如下：Hr=H-m-3-343.6 本章小結(jié)當前在實際應用的方法可以大致分為四類，線性擬合模型、曲面擬合模

23、型、地球重力場模型擬合法、地球重力模型與GPS水準相結(jié)合法。這些方法都要根據(jù)地形的情況、數(shù)據(jù)收集的情況等選擇一種擬合效果最佳的擬合方法。只有選擇對的擬合方法，才能在后面的過程中得到理想的擬合結(jié)果。當GPS點成線狀分布的時候，本文介紹了二種模型可供選擇，這二種模型是多項式曲線擬合法和三次樣條曲線擬合法。其擬合的精度與自身的條件和擬合條件相關。一般情況下，如果選點比較合適且已知聯(lián)測點分布比較均勻，都能達到厘米級的精度。多項式曲線擬合法簡單易理解且要求的已知點數(shù)少，缺點是其只考慮一個方向，求出的高程異常值不是非常準確。三次樣條曲線擬合法適用測線長控制點多的情況。此種方法把測區(qū)分為幾段，每段用多項式函

24、數(shù)表示，然后拼接為三次樣條函數(shù)，這樣做的目的是為了避免最小二乘原理的弊端。由于最小二乘原理在擬合比較多的已知點的時候會出現(xiàn)誤差較大的現(xiàn)象，對擬合結(jié)果產(chǎn)生不良的影響。當GPS點成面狀分布的時候，通常用曲面擬合法進行擬合，曲面擬合的方法有很多，其中主要有多項式曲面擬合法、平面及平面相關擬合法、移動曲面擬合法及多面函數(shù)曲面擬合法。和前面帶狀曲線擬合一樣，這四種方法對于曲面擬合來說各自有優(yōu)缺點，各自適用于不同情況。其中多項式曲面擬合最為普遍與基礎，這種方法就是用多項式函數(shù)擬合測區(qū)的似大地水準面，與多項式曲線擬合不同的是它考慮兩個方向的高程異常值，所以它的項數(shù)是比較多的，因此通常用二次或三次多項式。這種

25、方法一般用于地勢比較平坦的地區(qū)，當?shù)貏荼容^復雜的時候，就要求加入地形改正值消除地形的影響。平面及平面相關擬合法主要用于地勢近似于平面，或所測測區(qū)面積較小的時候，用這種模型擬合高程計算簡單，高程精度基本上可以滿足一般工程的需要。移動曲面擬合法是一種局部逼近法，這種方法適合地形較復雜且對擬合結(jié)果精度要求較高的情況下。當?shù)匦伪容^復雜的地區(qū)時，用面函數(shù)曲面擬合法無限逼近測區(qū)的似大地水準面。其中要分清楚多項式曲面擬合法和多面函數(shù)曲面擬合法的區(qū)別，多項式擬合適合地形平坦且高程異常變化值單一，表現(xiàn)在地形上就是其地勢向一個方向突出，且擬合區(qū)域不能太大。但是多面函數(shù)曲面擬合法就和其恰恰相反，它適合按網(wǎng)狀分布的GPS聯(lián)測點，地形比較復雜，控制點比較多的測區(qū)。而且當已知聯(lián)測點及顯著點越多的時候，其擬合出的高程異常值就會越高，其缺點就是沒有辦法先計算出合適光滑因子。只有經(jīng)過多次的驗證改進才能確定出合適的光滑因子，然后計算出精度比較高的高程異常值，進而求出高精度的正常高。相對于其他擬合法，重力場模型擬和法就相對用的較少，主要是對數(shù)據(jù)的要求比較苛刻，要得到高程異常值，首先要收集相關的重力場資料，這對于一般的工程是比較困難的，有時測區(qū)內(nèi)因缺少某些資料就不能采用此種方法，而且這種擬合方法與收集的數(shù)據(jù)的精度有關，若受到數(shù)據(jù)誤差影響，那么擬合的結(jié)果往往