張興福 沈云中

(1)廣東工業(yè)大學(xué)測(cè)繪工程系，廣州510006) 2)同濟(jì)大學(xué)測(cè)量與國(guó)土信息工程系，上海200092

一種實(shí)用的GPS坐標(biāo)及高程同步轉(zhuǎn)換方法*

張興福1)沈云中2)

(
1)廣東工業(yè)大學(xué)測(cè)繪工程系，廣州510006) 2)同濟(jì)大學(xué)測(cè)量與國(guó)土信息工程系，上海200092

基于GPS控制網(wǎng)三維約束平差基本原理，提出一種GPS坐標(biāo)及高程同步轉(zhuǎn)換方法，即先由GPS點(diǎn)的空間直角坐標(biāo)反算點(diǎn)間基線向量，再將旋轉(zhuǎn)參數(shù)、尺度參數(shù)以及國(guó)家坐標(biāo)系中的坐標(biāo)作為待估參數(shù)，通過(guò)附加公共點(diǎn)約束進(jìn)行坐標(biāo)及高程轉(zhuǎn)換。利用工程實(shí)例數(shù)據(jù)對(duì)該方法進(jìn)行了試驗(yàn)，結(jié)果表明:該方法非常靈活、有效，其坐標(biāo)及高程轉(zhuǎn)換精度良好。

GPS基線反算;三維約束平差;坐標(biāo)轉(zhuǎn)換;高程轉(zhuǎn)換;精度分析

AbstractThe synchronous transformation method of GPS coordinate and height based on GPS control network three-dimension constraint adjustment theory is given out.The steps of the method are as follows.At first，the baseline vector between two points is calculated from the three-dimension rectangular coordinate of GPS points，and the rotation parameters，scale parameter and the coordinate parameters in state coordinate system are also regarded as the unknown parameters，the coordinate and height transformation of GPS points are completed through the common points constraint.Finally，a case of GPS control network is used to test the method.The results show that the method is flexible，effective，and the transformation accuracy of coordinate and height are also well.

Key words:GPS baseline calculation;three－dimension constraint adjustment;coordinate transformation;height transformation;accuracy analysis

1 引言

GPS測(cè)量技術(shù)已經(jīng)廣泛應(yīng)用到日常測(cè)繪工作中，極大地提高了測(cè)量工作效率。雖然我國(guó)已經(jīng)建立了高精度的三維地心坐標(biāo)系——2000國(guó)家大地坐標(biāo)系(CGCS2000)，但很多工程項(xiàng)目還涉及GPS測(cè)量成果與1980西安坐標(biāo)系、1954北京坐標(biāo)系以及獨(dú)立坐標(biāo)系的轉(zhuǎn)換。坐標(biāo)轉(zhuǎn)換一般可分為坐標(biāo)系轉(zhuǎn)換和坐標(biāo)基準(zhǔn)轉(zhuǎn)換［1－4］。另外由于GPS測(cè)量獲得的高程為大地高，而實(shí)際工程更需要點(diǎn)的正常高信息，常規(guī)的方法是GPS坐標(biāo)與高程分別進(jìn)行轉(zhuǎn)換，因此如何在保證一定轉(zhuǎn)換精度的前提下實(shí)現(xiàn)GPS坐標(biāo)及高程的同步轉(zhuǎn)換將具有重要的現(xiàn)實(shí)意義。

我國(guó)國(guó)家平面控制網(wǎng)和高程控制網(wǎng)是分開(kāi)布設(shè)的，有些控制點(diǎn)平面坐標(biāo)精度很高，而大地高精度較低，甚至根本就沒(méi)有大地高數(shù)據(jù);而有些控制點(diǎn)只有很高精度的正常高，而無(wú)大地高及平面坐標(biāo)。而要完成三維坐標(biāo)基準(zhǔn)轉(zhuǎn)換，就必須知道至少3個(gè)公共點(diǎn)的空間直角坐標(biāo)，這就意味著同一個(gè)點(diǎn)的平面坐標(biāo)和大地高都要已知，且具有足夠的精度，這在我國(guó)可能較難實(shí)現(xiàn)。在這種情況下，很難用常規(guī)的空間七參數(shù)法實(shí)現(xiàn)坐標(biāo)基準(zhǔn)轉(zhuǎn)換。另外基于點(diǎn)坐標(biāo)的常規(guī)空間七參數(shù)坐標(biāo)轉(zhuǎn)換方法也存在一些不足:1)在小區(qū)域范圍內(nèi)，旋轉(zhuǎn)參數(shù)以及尺度參數(shù)和平移參數(shù)是強(qiáng)相關(guān)的，但旋轉(zhuǎn)參數(shù)和尺度參數(shù)是不相關(guān)的［5］;2)未將點(diǎn)在某一坐標(biāo)系中的坐標(biāo)作為待估參數(shù)與坐標(biāo)轉(zhuǎn)換參數(shù)一起進(jìn)行嚴(yán)密平差;3)公共點(diǎn)的誤差無(wú)法有效地分配到其它待求點(diǎn);4)公共點(diǎn)的大地高精度低，甚至根本就沒(méi)有大地高。

針對(duì)上述情況，本文基于GPS網(wǎng)的三維約束平差原理，利用已知GPS點(diǎn)的空間直角坐標(biāo)反算出各點(diǎn)間基線向量，并組成一定的閉合圖形，通過(guò)附加公共點(diǎn)約束實(shí)現(xiàn)坐標(biāo)及高程轉(zhuǎn)換。該方法屬嚴(yán)密平差方法，誤差方程式是基于坐標(biāo)差轉(zhuǎn)換模型，剛好能消除平移參數(shù)，避免參數(shù)強(qiáng)相關(guān)對(duì)坐標(biāo)轉(zhuǎn)換的影響。另外可綜合平面坐標(biāo)(或大地經(jīng)緯度)與高程(平面坐標(biāo)和高程可不是同一個(gè)點(diǎn)，大地高用正常高代替)，在實(shí)現(xiàn)三維空間七參數(shù)法坐標(biāo)轉(zhuǎn)換的同時(shí)，將GPS大地高轉(zhuǎn)換為正常高(高程約束點(diǎn)要等于或大于2個(gè))，避免平面坐標(biāo)與高程不匹配對(duì)坐標(biāo)轉(zhuǎn)換的影響，又可對(duì)點(diǎn)在國(guó)家坐標(biāo)或獨(dú)立坐標(biāo)系中的坐標(biāo)進(jìn)行估計(jì)，公共點(diǎn)的誤差會(huì)通過(guò)基線傳遞到待求點(diǎn)。最后利用某一GPS工程數(shù)據(jù)對(duì)本文方法進(jìn)行了實(shí)驗(yàn)研究，效果良好。

2 原理與方法

GPS控制網(wǎng)的三維約束平差的基本原理是將GPS基線向量作為直接觀測(cè)量，將控制點(diǎn)在國(guó)家或地方獨(dú)立坐標(biāo)系中的空間直角坐標(biāo)作為未知參數(shù)，同時(shí)考慮GPS測(cè)量獲得的WGS-84坐標(biāo)和國(guó)家坐標(biāo)或地方獨(dú)立坐標(biāo)的轉(zhuǎn)換關(guān)系，通過(guò)附加地面控制點(diǎn)約束的方法，將GPS測(cè)量成果直接轉(zhuǎn)換為國(guó)家坐標(biāo)系或地方獨(dú)立坐標(biāo)系［2］。本文利用GPS點(diǎn)已知的WGS-84坐標(biāo)反算相鄰點(diǎn)間基線向量，借助GPS控制網(wǎng)的三維約束平差原理進(jìn)行坐標(biāo)基準(zhǔn)轉(zhuǎn)換及高程轉(zhuǎn)換。圖1中A、B、C、1、2、3和4為7個(gè)GPS點(diǎn)，其中A、B、C為GPS和地面控制點(diǎn)的公共點(diǎn)，1、2、3、4為GPS點(diǎn)，基于三維約束平差理論的GPS坐標(biāo)轉(zhuǎn)換方法的具體實(shí)現(xiàn)過(guò)程為:先由7個(gè)點(diǎn)WGS-84坐標(biāo)反算1－2，1－A，A－B等基線組成如圖1所示的基線網(wǎng)，然后通過(guò)固定A，B，C 3點(diǎn)地面坐標(biāo)進(jìn)行三維約束平差即可將其他4個(gè)GPS點(diǎn)坐標(biāo)轉(zhuǎn)換到國(guó)家或獨(dú)立坐標(biāo)系。

圖1 GPS點(diǎn)反算基線Fig.1Baseline network inversed from GPS points

若設(shè)i、j為國(guó)家或獨(dú)立坐標(biāo)系中的兩點(diǎn)，其空間直角坐標(biāo)分別為(XiYiZi)和(XjYjZj)，相應(yīng)的坐標(biāo)差為(ΔXijΔYijΔZij)，兩點(diǎn)對(duì)應(yīng)WGS-84坐標(biāo)系坐標(biāo)差為(ΔˉXijΔˉYijΔˉZij)，現(xiàn)通過(guò)坐標(biāo)差形式計(jì)算兩坐標(biāo)系轉(zhuǎn)換參數(shù)(3個(gè)旋轉(zhuǎn)參數(shù)，1個(gè)尺度參數(shù))，同時(shí)將國(guó)家或獨(dú)立坐標(biāo)系的空間直角坐標(biāo)作為未知參數(shù)一起求解，則誤差方程式為

其中，δ^Xi、δ^Yi、δ^Zi、δ^Xj、δ^Yj、δ^Zj為國(guó)家坐標(biāo)系下i、j點(diǎn)的空間坐標(biāo)改正數(shù)，^εx、^εy、^εz為兩坐標(biāo)系旋轉(zhuǎn)參數(shù)，且為微小量，δ^μ為尺度參數(shù)，上標(biāo)“0”表示變量的近似值。

若將點(diǎn)的站心地平坐標(biāo)系作為過(guò)渡坐標(biāo)系［6］，取與局部橢球的定位定向有聯(lián)系的繞地平北、東方向和天頂方向的旋轉(zhuǎn)角^εη、^εξ和^εA為旋轉(zhuǎn)參數(shù)，這3個(gè)旋轉(zhuǎn)參數(shù)中前兩個(gè)旋轉(zhuǎn)參數(shù)用于實(shí)現(xiàn)橢球變換，第三個(gè)旋轉(zhuǎn)參數(shù)用于坐標(biāo)方位旋轉(zhuǎn)，則公式(1)可改寫為，

（5）菌苗基礎(chǔ)研究較為落后。從目前的具體實(shí)際分析，菌苗生產(chǎn)和使用的菌種都是20世紀(jì)60～70年代分離的菌株，由于時(shí)間過(guò)長(zhǎng)，在免疫特性方面與現(xiàn)代流行性的菌株存在一定的差距，會(huì)在一定程度上造成免疫失敗。

空間直角坐標(biāo)系和站心地平坐標(biāo)系的旋轉(zhuǎn)矩陣，B0和L0為大地經(jīng)緯度近似值，其他參數(shù)含義同公式(1)。

若只考慮將WGS-84坐標(biāo)系的三維空間直角坐標(biāo)轉(zhuǎn)換為國(guó)家坐標(biāo)系或獨(dú)立坐標(biāo)系中的兩維坐標(biāo)，則式(2)中的基準(zhǔn)轉(zhuǎn)換參數(shù)可只取尺度參數(shù)δμ和繞天頂方向的旋轉(zhuǎn)角εA。

利用式(2)列立坐標(biāo)基準(zhǔn)轉(zhuǎn)換誤差方程式后，其法方程秩虧，無(wú)法解算轉(zhuǎn)換參數(shù)及待求的國(guó)家或獨(dú)立坐標(biāo)系中的未知坐標(biāo)參數(shù)，因此需附加基準(zhǔn)條件，若K點(diǎn)在國(guó)家或獨(dú)立坐標(biāo)系中的大地緯度、大地經(jīng)度和大地高已知，則可列立約束條件，

其中，BK、LK、HK為已知值，B0K、L0K、H0K為概略值，而δB、δL、δH為改正數(shù)，若取B0K=BK，L0K=LK，H0K=HK，則δB=δL=δH=0。

由于將空間直角坐標(biāo)作為待求參數(shù)，為了建立大地經(jīng)緯度改正數(shù)以及大地高改正數(shù)與空間直角坐標(biāo)改正數(shù)的關(guān)系式，可對(duì)大地坐標(biāo)與空間直角坐標(biāo)轉(zhuǎn)換公式進(jìn)行全微分，整理可得，其中，P為WGS-84坐標(biāo)系中的基線向量權(quán)矩陣，若GPS點(diǎn)坐標(biāo)方差已知?jiǎng)t可計(jì)算相應(yīng)基線的權(quán)陣，若GPS點(diǎn)方差未知?jiǎng)t可用單位陣代替，而K為聯(lián)系系數(shù)，利用式(5)即可以將GPS測(cè)量獲得的WGS－84坐標(biāo)轉(zhuǎn)換為國(guó)家或獨(dú)立坐標(biāo)系。實(shí)際計(jì)算中，若公共點(diǎn)的平面坐標(biāo)已知，則可以用公式(3)、(4)中的前兩項(xiàng)列立條件方程，若公共點(diǎn)的高程已知，則可用公式(3)、(4)中的最后一項(xiàng)列立條件方程。該方法可通過(guò)綜合公共點(diǎn)的平面坐標(biāo)和高程信息實(shí)現(xiàn)GPS坐標(biāo)轉(zhuǎn)換及高程轉(zhuǎn)換。

3 算例與精度分析

選擇某一城市GPS控制網(wǎng)數(shù)據(jù)作為算例，該控制網(wǎng)控制范圍約1 000 km2，平均點(diǎn)間距約4 km，從控制網(wǎng)中選擇42個(gè)點(diǎn)，其中含已知平面點(diǎn)5個(gè)，坐標(biāo)基準(zhǔn)為1980西安坐標(biāo)系，34個(gè)GPS點(diǎn)聯(lián)測(cè)了正常高，高程基準(zhǔn)為1985國(guó)家高程基準(zhǔn)。坐標(biāo)基準(zhǔn)轉(zhuǎn)換過(guò)程為:

1)利用42個(gè)GPS點(diǎn)在WGS-84坐標(biāo)系中的空間直角坐標(biāo)反算兩點(diǎn)間的基線向量，并取各基線向量權(quán)陣為單位陣，并組成一定的閉合圖形，GPS點(diǎn)的WGS-84坐標(biāo)為無(wú)約束平差結(jié)果，組成的GPS基線網(wǎng)的示意圖見(jiàn)圖2。圖中實(shí)方框表示平面已知點(diǎn)，實(shí)三角形表示三維坐標(biāo)已知點(diǎn)，實(shí)圓點(diǎn)表示三維坐標(biāo)都未知的點(diǎn)，十字絲表示正常高已知的點(diǎn)。無(wú)論是否采用無(wú)約束平差坐標(biāo)，由于基線是由坐標(biāo)推算出來(lái)的，因此組成的閉合圖形的閉合差均為零。

圖2 GPS點(diǎn)基線組網(wǎng)Fig.2GPS baseline network

3)采用自編程序進(jìn)行坐標(biāo)基準(zhǔn)轉(zhuǎn)換，求取兩坐標(biāo)系3個(gè)旋轉(zhuǎn)參數(shù)(僅考慮平面坐標(biāo)轉(zhuǎn)換時(shí)，旋轉(zhuǎn)參數(shù)為一個(gè))，1個(gè)尺度參數(shù)以及點(diǎn)在國(guó)家或獨(dú)立坐標(biāo)系中的坐標(biāo)。此過(guò)程需要迭代計(jì)算，直到待估坐標(biāo)參數(shù)改正數(shù)小于設(shè)置的限差為止，國(guó)家坐標(biāo)系或獨(dú)立坐標(biāo)系中點(diǎn)的待估坐標(biāo)的初值可用WGS-84坐標(biāo)代替，此過(guò)程一般只需要迭代3次即可收斂。

3.1 GPS三維坐標(biāo)轉(zhuǎn)換為平面坐標(biāo)及精度分析

當(dāng)將三維的GPS網(wǎng)轉(zhuǎn)換到由常規(guī)技術(shù)構(gòu)建的國(guó)家或獨(dú)立坐標(biāo)系時(shí)，實(shí)際上這種控制網(wǎng)只是兩維網(wǎng)，因此很難獲得可靠且相互獨(dú)立的3個(gè)旋轉(zhuǎn)參數(shù)，對(duì)于地面點(diǎn)僅能固定一個(gè)點(diǎn)的三維坐標(biāo)，其余地面點(diǎn)充其量?jī)H能固定兩維坐標(biāo)(大地經(jīng)緯度)?，F(xiàn)選擇CB點(diǎn)的三維坐標(biāo)(大地經(jīng)緯度，大地高用正常高代替或取概略值)，BSZ點(diǎn)，WCL點(diǎn)以及EHG點(diǎn)的大地經(jīng)緯度作為約束條件，利用推算出來(lái)的86條基線進(jìn)行坐標(biāo)基準(zhǔn)轉(zhuǎn)換，由于只固定了一個(gè)點(diǎn)的高程，該高程起高程基準(zhǔn)作用，因此基準(zhǔn)轉(zhuǎn)換參數(shù)只取旋轉(zhuǎn)參數(shù)εA以及尺度參數(shù)δμ，此過(guò)程獲得的點(diǎn)的高程成果不可用，將其基準(zhǔn)轉(zhuǎn)換獲得的平面坐標(biāo)成果稱為CASE 1，相應(yīng)的基線改正數(shù)統(tǒng)計(jì)結(jié)果見(jiàn)表1，旋轉(zhuǎn)參數(shù)εA=－0.19 s，δμ=2.65 ppm，而將用原始GPS數(shù)據(jù)解算得到的基線及方差信息和同樣控制點(diǎn)進(jìn)行約束平差獲得平面坐標(biāo)成果稱為CASE 0，現(xiàn)對(duì)兩者進(jìn)行比較，結(jié)果見(jiàn)圖3。

圖3和表1結(jié)果顯示:利用本文方法將GPS網(wǎng)三維坐標(biāo)轉(zhuǎn)換為國(guó)家坐標(biāo)系的成果與嚴(yán)密平差方法轉(zhuǎn)換的結(jié)果相差不大，除3個(gè)點(diǎn)坐標(biāo)的y方向差值超過(guò)10mm外，其他點(diǎn)x、y方向差值均不超過(guò)10 mm，x、y方向差值的中誤差分別為3.3 mm和5.2 mm，說(shuō)明本文基準(zhǔn)轉(zhuǎn)換方法是可靠的。在基準(zhǔn)轉(zhuǎn)換過(guò)程中基線改正數(shù)絕大部分都小于1 cm，最大為2.35 cm，由于推算獲得的GPS基線網(wǎng)的閉合差為零，基線改正數(shù)應(yīng)當(dāng)主要是由點(diǎn)在國(guó)家坐標(biāo)系中的誤差分配造成的。

表1 坐標(biāo)轉(zhuǎn)換基線改正數(shù)統(tǒng)計(jì)結(jié)果(CASE 1)Tab.1Statistic results of corrected baseline value(CASE 1)

圖3 CASE 1與CASE 0坐標(biāo)比較Fig.3Comparison of the coordinates between CASE 1 and CASE 0

3.2 GPS三維坐標(biāo)同步轉(zhuǎn)換為平面坐標(biāo)和高程及精度分析

由于我國(guó)平面控制和高程控制通常是分開(kāi)布設(shè)的，因此很難獲得地面點(diǎn)的真正三維坐標(biāo)，有些地面點(diǎn)可能只有大地經(jīng)緯度(平面坐標(biāo)可反算出大地經(jīng)緯度)，而無(wú)準(zhǔn)確的大地高信息;有些點(diǎn)可能只有高程信息，特別是正常高信息，而無(wú)大地經(jīng)緯度。如何綜合利用地面點(diǎn)分離的大地經(jīng)緯度和正常高信息，通過(guò)三維坐標(biāo)轉(zhuǎn)換的方法，實(shí)現(xiàn)將GPS測(cè)量成果轉(zhuǎn)換為國(guó)家坐標(biāo)系中的平面坐標(biāo)，同時(shí)將GPS大地高轉(zhuǎn)換為正常高，這正是本文研究的內(nèi)容。本文在將CB點(diǎn)的三維坐標(biāo)(大地經(jīng)緯度及正常高)，BSZ點(diǎn)、WCL點(diǎn)以及EHG點(diǎn)大地經(jīng)緯度作為約束條件外，另將點(diǎn)G02、G58、G07、G37和G41的正常高作約束條件，其平面坐標(biāo)轉(zhuǎn)換結(jié)果稱為CASE 2，相應(yīng)的基線改正數(shù)統(tǒng)計(jì)結(jié)果見(jiàn)表2，與CSAE 1相應(yīng)的旋轉(zhuǎn)參數(shù)，εA=－0.20 s，δμ=2.55 ppm，由于采用三維坐標(biāo)轉(zhuǎn)換方法，轉(zhuǎn)換參數(shù)除εA和δμ外還有和兩個(gè)轉(zhuǎn)換參數(shù)，其與CASE 0坐標(biāo)轉(zhuǎn)換的比較結(jié)果見(jiàn)圖4，CASE 2和CASE 1轉(zhuǎn)換坐標(biāo)比較結(jié)果見(jiàn)圖5。從兩種方法計(jì)算獲得的對(duì)應(yīng)的旋轉(zhuǎn)參數(shù)以及尺度參數(shù)可以看出:兩者差別非常小，方位旋轉(zhuǎn)參數(shù)差了0.01 s，而尺度參數(shù)差了0.1 ppm，因此增加高程約束點(diǎn)對(duì)方位旋轉(zhuǎn)參數(shù)以及尺度參數(shù)影響可忽略。

圖4、圖5和表2結(jié)果顯示:雖然增加了5個(gè)高程約束點(diǎn)，采用3個(gè)旋轉(zhuǎn)參數(shù)和一個(gè)尺度參數(shù)實(shí)現(xiàn)坐標(biāo)和高程的同步轉(zhuǎn)換，坐標(biāo)轉(zhuǎn)換結(jié)果和CASE 0結(jié)果比較，同樣除3個(gè)點(diǎn)坐標(biāo)的y方向差值超過(guò)10 mm外，其他點(diǎn)x、y方向差值均不超過(guò)10 mm，x、y方向差值的中誤差分別為3.6 mm和6.3 mm。將CASE 2和CASE 1結(jié)果比較，x、y方向差值的中誤差分別為1.6 mm和3.6 mm，在該區(qū)域，增加高程點(diǎn)約束對(duì)y方向的影響要稍微大于x方向影響，整體看，增加高程約束點(diǎn)對(duì)平面坐標(biāo)轉(zhuǎn)換結(jié)果的影響比較小，但對(duì)基線改正還是有一定的影響，特別是對(duì)y和z方向影響較大。

圖4 CASE 2與CASE 0坐標(biāo)比較Fig.4Comparison of the coordinates between CASE 2 and CASE 0

圖5 CASE 1與CASE 2坐標(biāo)比較Fig.5Comparison of the coordinates between CASE 1 and CASE2

表2 坐標(biāo)轉(zhuǎn)換基線改正數(shù)統(tǒng)計(jì)結(jié)果(CASE 2)Tab.2Statistic results of corrected baseline values(CASE 2)

當(dāng)作為高程約束點(diǎn)超過(guò)1個(gè)時(shí)，若將選擇的正常高為高程約束量，則通過(guò)三維坐標(biāo)轉(zhuǎn)換方法，在將GPS點(diǎn)三維坐標(biāo)轉(zhuǎn)換到國(guó)家坐標(biāo)系大地經(jīng)緯度的同時(shí)，也可以將GPS大地高轉(zhuǎn)換為正常高，圖6為CSAE 2轉(zhuǎn)換成果中的高程結(jié)果與實(shí)測(cè)水準(zhǔn)高程的比較，圖6結(jié)果顯示:其高程轉(zhuǎn)換結(jié)果和實(shí)測(cè)水準(zhǔn)高的差值絕對(duì)值不超過(guò)50 mm，差值中誤差為23.8 mm，高程轉(zhuǎn)換的精度非常高。

圖6 CASE 2高程與真實(shí)高程比較Fig.6Comparison of the heights between CASE 2 and true one

4 結(jié)論

1)基于GPS三維約束平差理論，推導(dǎo)了顧及旋轉(zhuǎn)參數(shù)(3個(gè)或一個(gè))、尺度參數(shù)以及國(guó)家坐標(biāo)系或獨(dú)立坐標(biāo)系坐標(biāo)參數(shù)為待求參數(shù)的坐標(biāo)轉(zhuǎn)換模型，該方法利用GPS點(diǎn)空間直角坐標(biāo)反算點(diǎn)間基線向量，并通過(guò)附加大地經(jīng)緯度和正常高的方法實(shí)現(xiàn)GPS坐標(biāo)及高程的同步轉(zhuǎn)換，地面約束點(diǎn)無(wú)需具有三維坐標(biāo)，可根據(jù)情況分別選擇大地經(jīng)緯度和高程約束點(diǎn)，非常靈活。

2)實(shí)際應(yīng)用中，GPS點(diǎn)的空間直角坐標(biāo)可以采用無(wú)約束平差成果，也可以采用其他方法獲得的成果，如GPS-RTK等，對(duì)于有點(diǎn)的方差信息，可通過(guò)點(diǎn)的方差信息計(jì)算對(duì)應(yīng)基線權(quán)矩陣;若沒(méi)有精度信息，可將反算基線權(quán)矩陣定為單位權(quán)。無(wú)論采用那種方法獲得的成果，若形成閉合圖形，則閉合差均為零。

3)利用常規(guī)的空間七參數(shù)方法進(jìn)行坐標(biāo)轉(zhuǎn)換，旋轉(zhuǎn)參數(shù)及尺度參數(shù)與平移參數(shù)強(qiáng)相關(guān)，而本方法能消除這種相關(guān)性;另外公共點(diǎn)會(huì)有殘差，其殘差一般很難進(jìn)行分配，而本方法將公共點(diǎn)作為強(qiáng)制約束點(diǎn)，轉(zhuǎn)換之后不存在殘差，那么公共點(diǎn)的誤差會(huì)通過(guò)反算的基線傳遞到其他待求點(diǎn)，并進(jìn)行整體分配。

4)GPS坐標(biāo)及高程的同步轉(zhuǎn)換方法既可將GPS坐標(biāo)轉(zhuǎn)換為國(guó)家坐標(biāo)系中的平面坐標(biāo)，又可實(shí)現(xiàn)GPS高程轉(zhuǎn)換，其平面坐標(biāo)轉(zhuǎn)換精度與GPS控制網(wǎng)嚴(yán)密平差方法相差很小，高程轉(zhuǎn)換精度也很高，更重要的是增加高程約束點(diǎn)對(duì)平面坐標(biāo)轉(zhuǎn)換的影響較小。因此本文方法比較適合區(qū)域不大且地形起伏相對(duì)較小測(cè)區(qū)的GPS坐標(biāo)及高程同步轉(zhuǎn)換。

1施一民.現(xiàn)代大地控制測(cè)量［M］.上海:同濟(jì)大學(xué)出版社，2003.(Shi Yimin.Modern geodesy control surveying［M］.Shanghai:Tongji University Press，2003)

2劉大杰，等.全球定位系統(tǒng)(GPS)的原理與數(shù)據(jù)處理［M］.上海:同濟(jì)大學(xué)出版社，1996.(Liu Dajie，et al.The theory and data processing of Global Positioning system［M］.Shanghai:Tongji University Press，1996)

3朱華統(tǒng)，等.GPS坐標(biāo)系統(tǒng)的變換［M］.北京:測(cè)繪出版社，1994.(Zhu Huatong，et al.Transformation of GPS coordinate system［M］.Beijing:Surveying and Mapping Press，1994)

4楊元喜，徐天河.不同坐標(biāo)系綜合變換法［J］.武漢大學(xué)學(xué)報(bào)(信息科學(xué)版)，2001，26(6):509－512.(Yang Yuanxi and Xu Tianhe.The combined method of datum transformation between different coordinate systems［J］.Geomatics and Information Science of Wuhan University，2001，26(6):509－512)

5王解先.七參數(shù)轉(zhuǎn)換中參數(shù)之間的相關(guān)性［J］.大地測(cè)量與地球動(dòng)力學(xué)，2007，(2):43－46.(Wang Jiexian.Correlations among parameters in seven parameter transformation model［J］.Journal of Geodesy and Geodynamics，2007，(2):43–46)

6沈云中，衛(wèi)剛.利用過(guò)渡坐標(biāo)系改進(jìn)3維坐標(biāo)變換模型［J］.測(cè)繪學(xué)報(bào)，1998，27(2):161－165.(Shen Yunzhong and Wei Gang.Improvement of three dimensional coordinate transformation model by use of interim coordinate system［J］.Acta Geodaetica et Cartographica Sinica，1998，27(2): 161－165 )

A PRACTICAL METHOD FOR SYNCHRONOUS TRANSFORMATION OF GPS COORDINATE AND HEIGHT

Zhang Xingfu1)and Shen Yunzhong2)(
1)Department of Surveying and Mapping，Guangdong University of Technology，Guangzhou510006) 2)Department of Surveying and Geomatics，Tongji University，Shanghai 200092

P226+.3

A

1671-5942(2011)03-0063-06

2010-11-13

中國(guó)科學(xué)院動(dòng)力大地測(cè)量學(xué)重點(diǎn)實(shí)驗(yàn)室開(kāi)放基金(L09－07)

張興福，男，1977年生，副教授，主要從事衛(wèi)星重力、GPS數(shù)據(jù)處理及應(yīng)用軟件開(kāi)發(fā).E－mail:xfzhang77@163.com