蘭志武

(福州市勘測院　福建福州　350001)

?

基于CGCS2000橢球的城市坐標(biāo)系建立及應(yīng)用

蘭志武

(福州市勘測院福建福州350001)

簡要介紹了建立獨立坐標(biāo)系的常用方法，提出了采用“平移、旋轉(zhuǎn)、無縮放”方法，建立基于CGCS2000橢球的城市地方坐標(biāo)系。同時，描述了新坐標(biāo)系統(tǒng)的應(yīng)用以及新坐標(biāo)系統(tǒng)與國家(地方)坐標(biāo)系的相互轉(zhuǎn)換關(guān)系，實現(xiàn)現(xiàn)有測繪成果與新坐標(biāo)系統(tǒng)的準(zhǔn)確轉(zhuǎn)換。

CGCS2000；地方坐標(biāo)系；轉(zhuǎn)換精度

0　引言

為了統(tǒng)一城市規(guī)劃建設(shè)布局，需要統(tǒng)一的測繪基準(zhǔn)，特別是平面坐標(biāo)系統(tǒng)。根據(jù)文獻[1]規(guī)定：城市平面坐標(biāo)系統(tǒng)投影長度變形值不應(yīng)大于25mm/km，即相對誤差小于1/40 000，才能滿足大比例尺測圖、精密工程建設(shè)等需求。本文將以福州市永泰縣為例討論城市地方坐標(biāo)系的建立方法。

永泰縣位于福建省東部，介于北緯25°39′～26°05′，東經(jīng)118°23′～119°12′，東西長84km，南北寬46km，總面積2 243km2。城區(qū)中心距離120°中央子午線約117km，最東邊距120°中央子午線約80km,最西邊距120°中央子午線約162km；若直接采用120°中央經(jīng)線作為永泰縣區(qū)域3°帶投影的中央經(jīng)線，長度變形最大值約32cm/km，最小變形值約8cm/km，超出文獻[1]中規(guī)定的投影變形不大于2.5cm/km。因此，為了減小投影變形影響、滿足國土報批需要等要求，有必要建立新的地方獨立坐標(biāo)系。根據(jù)國測國發(fā)[2013]11號《關(guān)于加快2000國家大地坐標(biāo)系推廣使用的通知》。本文將討論基于CGCS2000參考橢球的城市地方獨立坐標(biāo)系建立方法及其應(yīng)用。

1　建立獨立坐標(biāo)系的方法

由文獻[2]中相關(guān)內(nèi)容可知，獨立坐標(biāo)系的建立方法大致可以歸納為以下3種：

(1)采用具有高程抵償面的國家統(tǒng)一3°帶高斯投影平面直角坐標(biāo)系統(tǒng)，其中央子午線與國家坐標(biāo)系采用的3°帶相同，投影面采用高程抵償面，高程抵償面一般選擇為當(dāng)?shù)氐钠骄叱堂?。此種方法保持中央子午線不變，選擇某一高程面作為歸化高程面，使高程歸化改正和高斯投影變形改正相互抵消，使得測區(qū)中央的兩項投影改正接近于零；

(2)采用自定義中央子午線高斯投影平面直角坐標(biāo)系統(tǒng)，其投影面采用國家坐標(biāo)系參考橢球面。此種方法保持高程歸化面不變，只移動中央子午線的方法，使得測區(qū)中央的兩項投影改正接近于零；

(3)采用具有高程抵償面的自定義中央子午線高斯投影平面直角坐標(biāo)系統(tǒng)，其投影面一般選擇為當(dāng)?shù)氐钠骄叱堂妫醒胱游缇€位于城市中心。此種方法以測區(qū)中央為中央子午線，歸化高程面提高到測區(qū)的平均高程面上，建立任意帶高斯正形投影平面直角坐標(biāo)系，這樣可以保證測區(qū)的兩項改正在測區(qū)中央幾乎為零。

此外，還有采用橢球變換的方法即運用橢球膨脹法、橢球平移法、橢球變形法等建立獨立坐標(biāo)系[5]。各種獨立坐標(biāo)系的建立方法各有其優(yōu)點、缺點，在實際的應(yīng)用過程中有其特殊要求。

本文將采用“平移、旋轉(zhuǎn)、無縮放”的方法建立城市地方坐標(biāo)系，并對其在實際的生產(chǎn)應(yīng)用中加以介紹。

2　城市地方坐標(biāo)系的建立

2.1坐標(biāo)系統(tǒng)的定義

由文獻[2]相關(guān)規(guī)范可知：建立城市平面坐標(biāo)系統(tǒng)涉及到參考橢球的選擇、中央子午線的選擇、投影面高程的選擇以及坐標(biāo)系原點和定向的選擇等內(nèi)容。

根據(jù)此規(guī)范要求，通過對永泰縣地理位置、城市建設(shè)中長期規(guī)劃發(fā)展?fàn)顩r、以及經(jīng)濟建設(shè)狀況等因素，綜合考慮采用自定義中央子午線高斯正形投影方式建立城市地方獨立坐標(biāo)系，以在城區(qū)中央的118°50′為坐標(biāo)系的中央子午線，投影后平移旋轉(zhuǎn)至中央子午線為120°的1980西安坐標(biāo)系(尺度比為1)，參考橢球采用CGCS2000坐標(biāo)橢球，投影面采用CGCS2000坐標(biāo)系參考橢球面，坐標(biāo)系原點和定向同CGCS2000坐標(biāo)系，高程系統(tǒng)采用1985國家高程基準(zhǔn)。此方法可以確保長度投影變形不大于2.5cm/km，滿足城市測量和城鎮(zhèn)地籍調(diào)查精度要求等；此外，新建立的城市地方坐標(biāo)系統(tǒng)即CGCS2000(Y)坐標(biāo)系采用CGCS2000國家坐標(biāo)系，坐標(biāo)形式為中央經(jīng)線120°的，真正的投影方式為118°50′。

2.2CGCS2000(Y)坐標(biāo)系控制點成果的獲取

在坐標(biāo)系定義完成后，需要獲取在CGCS2000(Y)坐標(biāo)系下相應(yīng)的控制點成果，為后續(xù)的相關(guān)工作提供基礎(chǔ)數(shù)據(jù)(如計算與1980西安坐標(biāo)系、WGS-84坐標(biāo)系、CGCS2000坐標(biāo)系等之間的轉(zhuǎn)換參數(shù))。

按照如下流程獲取控制點在新坐標(biāo)系下的成果：

(1)獲取控制點的CGCS2000-118.50和XA80-120坐標(biāo)，求取CGCS2000-118.50至XA80-120的首次四參數(shù)，并解算擬建坐標(biāo)系初次成果；

(2)將四參數(shù)中尺度比設(shè)置為1后反復(fù)迭代計算擇優(yōu)選取轉(zhuǎn)換參數(shù)，解算擬建坐標(biāo)系CGCS2000(Y)坐標(biāo)系的最終成果；

(3)計算其他控制點的CGCS2000(Y)坐標(biāo)系成果用于與其他坐標(biāo)系轉(zhuǎn)換參數(shù)的求解。

2.3CGCS2000(Y)坐標(biāo)系下投影變形分析

根據(jù)上述方法獲取了已有控制點成果的CGCS2000(Y)坐標(biāo)系坐標(biāo)成果，為了檢核投影長度變形是否滿足不大于2.5cm/km的要求，采用RTK測量方式實地采集部分控制點的坐標(biāo)，比較已知控制點間的距離與實測已知控制點間的距離，判斷新建坐標(biāo)系的測量成果是否滿足精度要求。投影長度變形實地檢核較差統(tǒng)計見表1。

表1　實測(RTK)邊長檢核較差統(tǒng)計表

由表1可知：邊長較差最大值為2.80cm，最小值為0.05cm，相對誤差小于1/14000,符合文獻[3]中RTK平面控制點檢核測量技術(shù)要求。

3　CGCS2000(Y)坐標(biāo)系與國家(地方)坐標(biāo)系之間的聯(lián)系

根據(jù)文獻[2]中相關(guān)規(guī)范要求，新建立的坐標(biāo)系應(yīng)實現(xiàn)與現(xiàn)行的國家坐標(biāo)系(如CGCS2000坐標(biāo)系、1980西安坐標(biāo)系、1954北京坐標(biāo)系)之間的相互轉(zhuǎn)換，以及城市平面坐標(biāo)系新、舊網(wǎng)之間的相互轉(zhuǎn)換。

坐標(biāo)轉(zhuǎn)換工作非常復(fù)雜又重要，應(yīng)本著“最匹配本區(qū)域原有平面控制點坐標(biāo)系統(tǒng)和轉(zhuǎn)換后成果精度損失最小”為原則，對各階段轉(zhuǎn)換結(jié)果進行嚴(yán)格檢查和反復(fù)驗算，每一階段成果都合格后，成果才能提交使用。

3.1坐標(biāo)轉(zhuǎn)換數(shù)學(xué)模型

常用的坐標(biāo)轉(zhuǎn)換一般包括各種空間直角坐標(biāo)系與大地坐標(biāo)系，地心空間直角坐標(biāo)系與參心空間直角坐標(biāo)系，以及不同參心空間直角坐標(biāo)系之間的相互轉(zhuǎn)換。當(dāng)不同坐標(biāo)系之間存在嚴(yán)密的數(shù)學(xué)轉(zhuǎn)換模型時，可以采用相應(yīng)的模型之間進行坐標(biāo)轉(zhuǎn)換，目前常采用轉(zhuǎn)換模型有布爾沙模型、三維七參數(shù)轉(zhuǎn)換模型、二維七參數(shù)轉(zhuǎn)換模型、二維四參數(shù)轉(zhuǎn)換模型、三維四參數(shù)轉(zhuǎn)換模型等眾多轉(zhuǎn)換模型。本文以采用二維四參數(shù)轉(zhuǎn)換模型求解CGCS2000(Y)坐標(biāo)系與CGCS2000坐標(biāo)系之間的轉(zhuǎn)換參數(shù)。

將CGCS2000坐標(biāo)系坐標(biāo)按118°50′3度帶高斯投影后與CGCS2000(Y)坐標(biāo)系通過二維四參數(shù)轉(zhuǎn)換模型進行轉(zhuǎn)換。

二維四參數(shù)轉(zhuǎn)換的模型如下：

式中：

x1，x2——源坐標(biāo)系坐標(biāo)，單位為m

x2,y2——目標(biāo)坐標(biāo)系坐標(biāo)，單位為m

Δx，Δy——為平移參數(shù)，單位為m

α——為旋轉(zhuǎn)參數(shù)，單位為rad

m——尺度參數(shù)，無量綱

1980西安坐標(biāo)系、WGS-84坐標(biāo)系、1954北京坐標(biāo)系與CGCS2000(Y)坐標(biāo)系轉(zhuǎn)換四參數(shù)均可按照圖2計算流程求解。

3.2精度評定

為了科學(xué)、客觀地評價計算轉(zhuǎn)換參數(shù)的可靠性和準(zhǔn)確性，將用于參數(shù)計算的重合點代入轉(zhuǎn)換模型，用得到的轉(zhuǎn)換參數(shù)計算重合點坐標(biāo)殘差，剔除殘差大于3倍點位中誤差的重合點，再以剩下的點重新計算轉(zhuǎn)換參數(shù)，直到所有參與計算轉(zhuǎn)換參數(shù)的點都滿足精度要求為止[6]。根據(jù)已有點位分布情況，最終選擇了等級較高、精度可靠的28個控制點作為求解參數(shù)重合點(其中GPSA級點2個、GPSC級點8個、城市三等點18個)，如圖4.2所示；以及參數(shù)驗算重合點分布略圖，如圖4.3所示。根據(jù)圖2計算流程，計算CGCS2000(Y)坐標(biāo)系至CGCS2000坐標(biāo)系二維轉(zhuǎn)換四參數(shù)。鑒于轉(zhuǎn)換參數(shù)的保密性，本文沒有列舉4個轉(zhuǎn)換參數(shù)的具體數(shù)值。

根據(jù)轉(zhuǎn)換前后的坐標(biāo)數(shù)據(jù)，進行轉(zhuǎn)換參數(shù)的精度評定。具體方法如下：

V(轉(zhuǎn)換殘差)=重合點轉(zhuǎn)換坐標(biāo)值-重合點已知坐標(biāo)值

式(1)

利用式(4.2)計算各平面點轉(zhuǎn)換誤差M點：

式(2)

同時，利用式(4.3)評定轉(zhuǎn)換參數(shù)轉(zhuǎn)換精度：

式(3)

式中，MX為X方向轉(zhuǎn)換中誤差，MY為Y方向轉(zhuǎn)換中誤差，

則可以利用式(4.4)計算轉(zhuǎn)換中誤差M：

式(4)

根據(jù)上式計算、統(tǒng)計采用二維轉(zhuǎn)換四參數(shù)的X方向轉(zhuǎn)換中誤差、Y方向轉(zhuǎn)換中誤差以及轉(zhuǎn)換中誤差等精度指標(biāo)。限于篇幅，文中只列舉了部分參數(shù)驗算重合點的轉(zhuǎn)換誤差，具體見表2和表3。

表2　二維四參數(shù)轉(zhuǎn)換誤差(部分)

表3　二維四參數(shù)轉(zhuǎn)換精度統(tǒng)計表

由表2、表3統(tǒng)計結(jié)果可知，二維轉(zhuǎn)換四參數(shù)轉(zhuǎn)換結(jié)果中，求解參數(shù)重合點點位誤差最小值為0mm，點位誤差最大值為1.1mm，點位中誤差為0.5mm，滿足規(guī)范要求的參數(shù)計算重合點的點位誤差小于3倍點位中誤差；對于參數(shù)驗算重合點點位誤差最小值為0.3mm，點位誤差最大值為1.1mm，點位中誤差為0.4mm，轉(zhuǎn)換精度較高。

此外，結(jié)合表2和圖5、圖6曲線走勢可以看出，求解參數(shù)重合點在x方向轉(zhuǎn)換誤差均在1mm以內(nèi)，96%的點在y方向的轉(zhuǎn)換誤差及轉(zhuǎn)換中誤差均在1mm以內(nèi)；參數(shù)驗算重合點在x方向轉(zhuǎn)換誤差均在1mm以內(nèi)，99%的點在y方向的轉(zhuǎn)換誤差均在1mm內(nèi)，98%的點轉(zhuǎn)換中誤差均在1mm以內(nèi)，轉(zhuǎn)換精度比較高。

3.3CGCS2000(Y)坐標(biāo)系的應(yīng)用

在新坐標(biāo)系統(tǒng)建立完成后，不但要建立新坐標(biāo)系統(tǒng)與原有坐標(biāo)系統(tǒng)的轉(zhuǎn)換關(guān)系，實現(xiàn)現(xiàn)有測繪成果與新坐標(biāo)系統(tǒng)之間相互精確轉(zhuǎn)換，而且還要規(guī)范新坐標(biāo)系統(tǒng)的使用。因此，在新坐標(biāo)系使用過程中要注意以下四點：

(1)除國土用地報批特殊要求外，城市規(guī)劃、工程建設(shè)、勘測定界等工作均在新的坐標(biāo)系統(tǒng)下實施完成；

(2)工程放樣：利用相應(yīng)的轉(zhuǎn)換參數(shù)(如1980西安坐標(biāo)系至CGCS2000(Y)坐標(biāo)系轉(zhuǎn)換四參數(shù))將設(shè)計坐標(biāo)轉(zhuǎn)換至CGCS2000(Y)坐標(biāo)系，采用CGCS2000(Y)坐標(biāo)進行放樣；

(3)日?？刂茰y量平差計算(用于GPS靜態(tài)解算)：平差計算應(yīng)放在投影經(jīng)線為118度50分下進行，投影高為0；無約束平差后，輸入已知平面控制點時，已知點東西方向加上*米后參與計算，投影設(shè)置選擇118度50分，待平差后東西方向坐標(biāo)成果減去*米；

(4)針對國土用地報批的地塊，建立新坐標(biāo)系統(tǒng)與西安80坐標(biāo)系的轉(zhuǎn)換關(guān)系，實現(xiàn)用地紅線由新坐標(biāo)系統(tǒng)向西安80坐標(biāo)系的轉(zhuǎn)換(面積不變)，滿足國土用地報批需求。

4　結(jié)語

本文以永泰縣為例，討論了基于CGCS2000橢球建立城市地方坐標(biāo)系統(tǒng)的過程，根據(jù)永泰縣的具體情況確定了相應(yīng)的建立方案，并建立了新坐標(biāo)系統(tǒng)與國家(地方)坐標(biāo)系統(tǒng)聯(lián)系。根據(jù)轉(zhuǎn)換精度可知，轉(zhuǎn)換結(jié)果準(zhǔn)確、可靠，可以用來實現(xiàn)已有測繪成果與新坐標(biāo)系統(tǒng)之間相互轉(zhuǎn)換，而且不會造成精度的損失。

此外，新坐標(biāo)系統(tǒng)的建立解決了長度變形問題帶來的困擾，有利于提高測繪成果的精度，更加簡單方便地獲取可靠的高精度成果。同時，將已有各種坐標(biāo)系統(tǒng)的測繪成果統(tǒng)一到新坐標(biāo)系統(tǒng)下，解決了多種坐標(biāo)成果資料混亂問題，為日常工作中測繪成果的使用提供了方便。

[1]CJJ/T8-2011 城市測量規(guī)范 [S]

[2]GB/T28584-2012 城市坐標(biāo)系統(tǒng)建設(shè)規(guī)范 [S]

[3]CJJ/T73-2010衛(wèi)星定位城市測量技術(shù)規(guī)范[S]

[4]李東,毛之琳.基于CGCS2000的地方坐標(biāo)系統(tǒng)建立方法的研究[J].測繪技術(shù)裝備，2009(4):3～5

[5]海清.通過橢球變換建立區(qū)域獨立坐標(biāo)系的方法[J].海洋測繪，2007(9)：31～34

[6]大地測量控制點坐標(biāo)轉(zhuǎn)換技術(shù)規(guī)程[M].北京：國家測繪地理信息局，2013.

蘭志武(1967.1-)，高級工程師，主要從事測繪與地理信息工程。

OntheEstablishmentandApplicationofCityCoordinateSystemBasedonCGCS2000Ellipsoid

LAN Zhiwu

(FuZhouInvestigationandSurveyingInstitute,Fuzhou350001)

Inviewoftheexistingproblemsofurbancoordinatesystem,thispaperbrieflyintroducesthecommonmethodsofestablishingindependentcoordinatesystemandputforwardthe"translation,rotation,noscaling"methods,andestablishesurbanlocalcoordinatesystembasedonCGCS2000ellipsoid.Atthesametime,thispaperdescribestheapplicationofthenewcoordinatesystem,andthemutualchangerelationofnewcoordinatesystemwiththenational(local)coordinatesystem,achievesaccurateconversionofexistingsurveyingandmappingresultsandthenewcoordinatesystem.

CGCS2000；localcoordinatesystem；TransformationAccuracy

蘭志武(1967.1-)，高級工程師。

E-mail:610428934@qq.com

2016-02-10

TU198

A

1004-6135(2016)03-0109-05