CORS測量成果轉換的一步法及其精度分析

杜向鋒，張興福，張永毅

(1.廣東工貿職業技術學院，廣東 廣州 510510； 2.廣東工業大學，廣東 廣州 510006)

?

CORS測量成果轉換的一步法及其精度分析

杜向鋒1，張興福2，張永毅2

(1.廣東工貿職業技術學院，廣東 廣州 510510； 2.廣東工業大學，廣東 廣州 510006)

針對目前 CORS測量成果轉換中平面和高程獨立轉換這一情況，提出了適用于CORS測量成果轉換的一步法，其本質是將區域橢球平移至區域似大地水準面，即將點的正常高作為大地高使用，并通過空間七參數轉換模型實現坐標和高程的同步轉換。兩個工程實例計算結果表明，該方法簡單有效，其平面和高程轉換精度均較為理想，可滿足一般工程測量的要求，具有較好的應用價值。

CORS；成果轉換；一步法；精度分析；區域似大地水準面模型

一、引 言

隨著CORS技術的不斷發展和完善，普通工程測量的手段發生了根本性的變化，利用CORS網絡，用戶很容易獲得測量點的WGS-84坐標，然后分別利用坐標轉換程序和似大地水準面內插程序即可將該點的WGS-84坐標轉換為國家(城市)坐標及國家(城市)高程，該方法具有測量速度快(一般僅需數十秒至數分鐘)、可實時定位、精度較高等優點，目前已經被廣泛應用于工程測量[1-3]。CORS測量成果轉換常采用坐標轉換和高程轉換分離的方法，即坐標轉換通過平面四參數或空間七參數實現[1-3]，而高程轉換可通過區域似大地水準面模型或采用高程擬合的方法實現。針對這種情況[4-10]，本文提出了一種適用于CORS成果轉換的一步法。試驗證明，本文方法只需要少量控制點(≥3點)算出7個坐標轉換參數，即可實現CORS成果的一步轉換，轉換精度較為理想，可用于普通工程測量，無須復雜的區域似大地水準面支持或單獨進行高程擬合計算等，具有很好的應用價值。

二、原理與方法

CORS 成果轉換一步法的基本原理是將區域橢球平移至區域似大地水準面，此時控制點的正常高可作為大地高使用；然后再選取一定量的公共點，計算坐標轉換的空間七參數，進而實現CORS測量成果的一步轉換，具體過程為：

1) 選取控制網中適量且分布合理的國家(城市)控制點數據，利用控制點的平面坐標，通過高斯反算獲得這些點在對應橢球上的大地緯度B和大地經度L，然后將該點的大地高用正常高代替，并將其轉換為空間直角坐標，公式為[1-3]

(1)

式中，(X,Y,Z)為測站在WGS-84橢球下空間直角坐標；B和L分別為點的大地緯度和大地經度；h為該點的正常高；N為卯酉圈的半徑；e為參考橢球第一偏心率；a和b分別為對應橢球的長半軸和短半軸。且

2) 將步驟1)中控制點在WGS-84坐標系下的大地緯度B、大地經度L和大地高H也轉換為對應參考基準下的空間直角坐標，計算公式見式(1)。

3) 將步驟1)和步驟2)獲得的控制點組成公共點對，并利用布爾莎-沃爾夫模型求取兩個基準間的七參數，公式為[1-3]

(2)

式中，X0、Y0、Z0為兩個空間直角坐標系之間的3個平移參數；R1、R2、R3為3個旋轉參數；D為尺度改正量。利用3個以上公共點按照最小二乘原理即可計算出以上7個轉換參數。

4) 利用步驟(3)獲得的七參數將CORS測量成果轉換為國家(城市)大地坐標，此時的大地高即為正常高，利用高斯正算公式將大地經緯度轉換為平面坐標，具體流程如圖1所示。該方法的本質是在區域范圍內將國家(城市)參考橢球平移到區域平均似大地水準面，此時正常高可作為大地高使用。

圖1 CORS測量成果轉換的一步法流程

三、實例分析

為了分析本文提出的一步法的效果，分別對兩個控制網數據進行了分析，并同時考慮公共點數量及位置分布等因素對計算結果的影響。為了檢驗轉換方法的可靠性及精度，采用內符合和外符合進行精度評定。內符合即為利用參與轉換參數計算的公共點的坐標和高程殘差進行精度統計；而外符合是指利用非參與(待轉點)轉換參數計算點的坐標差值和高程差值進行精度統計。

1.A市城市控制網

A市城市D級控制網由44個網點組成，平均點間距離4.6 km，控制范圍東西方向約25 km，南北方向約30 km，每個控制點均進行了三等水準聯測，平面坐標為1980西安坐標系，高程系統為1985國家高程基準，其控制點分布如圖2和圖3所示。其中，“▲”代表公共點，“+”代表待定點，下同。

A市控制網計算采用兩種方案，具體如圖2和圖3所示。方案Ⅰ只選取了最外圍3個點作為公共點，其余點均作為待定點進行計算；方案Ⅱ將公共點均勻地增加到了13個點。將各計算方案獲得的結果與已知結果進行對比統計，其結果見表1。

圖2 A市計算方案Ⅰ

圖3 A市計算方案Ⅱ

表1 A市控制網不同方案計算結果精度統計 mm

從表1的結果中可以看出，兩個計算方案獲得的結果與已知結果比較互差均在厘米級，其中方案Ⅰ的結果相對較差。方案Ⅰ中內符合精度為：平面X、Y方向標準差分別為5.6 mm和18.2 mm，高程標準差為0 mm(原因點數少)，外符合精度為：平面X、Y方向標準差分別為5.3 mm和11.2 mm，高程標準差為31.6 mm。方案Ⅱ中內符合精度為：平面X、Y方向標準差分別為9.3 mm和11.3 mm，高程標準差為26.7 mm，外符合精度為：平面X、Y方向標準差分別為4.0 mm和7.6 mm，高程標準差為23.5 mm。說明隨著公共點數量的增加和點位分布更加均勻，計算結果的精度得到一定提高，由此可見本文提出的一步成果轉換法是有效的，可一次性實現平面與高程的同步轉換，且轉換結果的精度可滿足普通工程測量的要求。

2.B市某工程控制網

B市因市政工程需要，利用該市CORS網絡及似大地水準面建立了一個施工控制網，網點數99個，平均點間距約4 km，控制范圍南北方向跨度約35 km，東西方向跨度約50 km，平面坐標基準及高程基準分別為城市獨立坐標系和城市獨立高程基準，選擇了4種計算方案，各計算方案選點情況如圖4—圖7所示。

圖4 B市計算方案Ⅰ

圖5 B市計算方案Ⅱ

圖6 B市計算方案Ⅲ

圖7 B市計算方案Ⅳ

由于B市坐標成果為城市獨立坐標系，其數值大小和國家坐標系相差較大，故在坐標轉換時首先選擇控制網某一點計算該點兩坐標系平面坐標的近似平移參數，然后再進行成果轉換，計算結果見表2。

表2 B市某工程控制網不同方案計算結果精度統計 mm

從表2可以看到，方案Ⅰ中內符合精度：平面X、Y方向標準差分別為10.9 mm和4.7 mm，高程標準差為6.4 mm；外符合精度：平面X、Y方向標準差分別為4.1 mm和3.5 mm，高程標準差為18.8 mm。方案Ⅱ中內符合精度：平面X、Y方向標準差分別為9.5 mm和4.4 mm，高程標準差為5.6 mm；外符合精度：平面X、Y方向標準差分別為4.1 mm和3.5 mm，高程標準差為18.9 mm。方案Ⅲ中內符合精度：平面X、Y方向標準差分別為7.1 mm和4.6 mm，高程標準差為18.2 mm；外符合精度：平面X、Y方向標準差分別為4.1 mm和3.3 mm，高程標準差為17.3 mm。方案Ⅳ中內符合精度：平面X、Y方向標準差分別為5.1 mm和3.6 mm，高程標準差為17.8 mm；外符合精度：平面X、Y方向標準差分別為4.2 mm和3.3 mm，高程標準差為16.2 mm。所有計算方案獲得內外符合精度在X、Y方向的標準差均不超過11 mm，高程標準差不超過20 mm，隨著公共點數的增加，高程轉換精度有一定提高，平面精度變化不大，高程轉換的外符合精度由18.8 mm提高到16.2 mm。

四、結 論

通過對以上兩個控制網數據的計算分析，可以得出以下結論：

1) 本文提出的用于CORS成果轉換的一步法實現方法簡單有效，其轉換結果精度較為理想，平面精度可達到1 cm，高程精度優于3 cm，能夠滿足普通工程測量的精度需求。

2) 隨著公共點數量的增加和分布均勻化，平面坐標轉換的精度提高有限，而高程轉換的精度有一定的提高，但精度提高速度逐漸放緩，若精度要求不高，也可認為在區域不大情況下，本文方法對公共點數要求不高。

3) 對于數百平方千米甚至1000多平方千米的測量區域(如中小城市等)，一套七參數即可滿足生產的需要，在實際生產中，將該參數固化到測量設備中即可實現實時、快速高精度的三維測量，方法具有很好的應用價值。

[1] 李征航，黃勁松．GPS測量與數據處理[M]．武漢：武漢大學出版社，2005．

[2] 陳俊勇.中國現代大地基準——中國大地坐標系統2000(CGCS 2000)及其框架[J].測繪學報,2008,37(3):269-271.

[3] 劉基余.GPS衛星導航定位原理與方法[M].北京：科學出版社, 2003.

[4] 張興福,魏德宏.基于似大地水準面格網的插值方法及精度分析[J].大地測量與地球動力學,2011,31(1): 113-116，122.

[5] 楊光,林鴻,歐海平,等.廣州市亞厘米級高精度似大地水準面的確定[J].測繪通報，2007(1)：24-25，32.

[6] 張瓊,胡亞明,蘇貝,等.基于超高階地球重力場模型的哈爾濱市區域似大地水準面確定[J].測繪科學,2011,36(4):56-57，82.

[7] 寧津生,羅志才,楊沾吉,等.深圳市1 km高分辨率厘米級高精度大地水準面的確定[J].測繪學報,2003,32(2):102-107.

[8] 陳俊勇,李建成,寧津生,等.全國及部分省市地區高精度高分辨率似大地水準面的研究和實施[J].測繪通報,2005(5):1-5.

[9] 張興福,李博峰,魏德宏,等.多類重力場模型的精度分析及聯合確定GPS點正常高的方法[J].測繪學報,2013,42(1):6-12.

[10] 張興福，劉成.綜合EGM2008模型和SRTM/DTM2006.0剩余地形模型的GPS高程轉換方法[J].測繪學報,2012,41(1): 25-32.

One-step Method and Accuracy Analysis of CORS Coordinate and Height Transformation

DU Xiangfeng，ZHANG Xingfu，ZHANG Yongyi

杜向鋒，張興福，張永毅.CORS測量成果轉換的一步法及其精度分析[J].測繪通報，2015(7)：23-26.

10.13474/j.cnki.11-2246.2015.0203

2015-03-04

國家自然科學基金(41104002)

杜向鋒(1980—)，男，碩士，講師，主要從事GPS和工程測量研究。E-mail：dunac@163.com

P228.4

：B

：0494-0911(2015)07-0023-04