某市地方獨立坐標系與獨立工程坐標系的轉換

摘 要:坐標系轉換是測量中經常需要解決的問題。在小區域范圍內，可將地方獨立坐標系和工程坐標系視為在平面上建立的坐標系，而不需考慮投影產生的距離和方向變形，從而采用平面上的坐標系相似變換方式實現兩個坐標系的想互轉換。本文詳細敘述了以工程坐標系控制點為重合點，采用坐標系相似變換方法實現某市地方獨立坐標系與某企業獨立工程坐標系的相互轉換。

關鍵詞:地方獨立坐標系 獨立工程坐標系 坐標系轉換 相似變換

中圖分類號:TB22文獻標識碼:A文章編號:1674-098X(2011)04(b)-0096-02

1 引言

在工程測量中，對于投影的長度變形都有一定的要求，變形不能過大。由于高斯投影邊長存在變形，會使高斯投影計算邊長與實測邊長產生差異。因此，在工程測量中，為了工程測量和施工的方便，經常采用獨立的工程坐標系。[1]這就必然涉及到國家坐標系和地方坐標系與工程坐標系的相互轉換。本文論述了某市地方獨立坐標系與獨立工程坐標系的轉換。

2 坐標系簡介

在測量中，平面坐標系有國家統一坐標系，地方獨立坐標系和工程坐標系。

國家統一坐標系，投影面采取國家大地基準所確定的國家參考橢球面，按照6度帶或3度帶分帶進行投影，投影中央子午線為投影分帶的中央子午線。按照6度帶投影則投影中央子午線為l=6n-3;按照3度帶投影則投影中央子午線為l=3n。在投影面上，投影中央子午線的投影為X軸，赤道投影為Y軸。為避免在同一帶內Y坐標出現負值，將X軸向西移動500Km，即所有Y坐標值加500km。[2]采用國家統一坐標系，有利于測繪成果的統一和成果一測多用，有利于地球空間數據的交流和共享。

但是，由于許多城市和工程或是遠離國家3度帶中央子午線，或是地勢較高，若唯一采用全國統一的3分帶并基于國家參考橢球的高斯平面坐標系，邊長將產生甚大的投影變形，實難滿足城市規劃、工程建設和管理等各方面的應用。為滿足施工放樣和測圖控制的要求，由控制點平面坐標反算的邊長應與實測平距相接近。只有當3度帶中央子午線在測區附近，且測區的平均海拔高度又不太大時，才可采用國家3度帶高斯平面直角坐標系。當不滿足上述基本原則時，須采用不同于國家統一3度帶坐標的地方獨立坐標系。[3]因此，基于限制變形、方便、實用和科學的目的，在許多城市和工程測量中，常常會建立適合本地區的地方獨立坐標系。

地方獨立坐標系有兩種:一種是以測區的平均高程面為投影面的任意高斯平面直角坐標系，即任意帶坐標系。其中央子午線經過測區東西向的中心附近，或是經過某一控制點，或是其經度取某一整度;另一種是以抵償高程面為投影面的3帶高斯平面直角坐標系，即抵償坐標系，其中央子午線仍取最靠近的統一3帶的中央子午線。[2]

在某些工程建設和管理中，出于施工方便和保密等原則，既不采用國家統一坐標系，也不采用地方獨立坐標系，而是自定一個工程坐標系。工程坐標系一般是獨立平面直角坐標系，即選一個自定義投影帶，采用與測區平均高程面相切且與參考橢球面相平行的橢球面，通過測區中部的子午線作為中央子午線。當測區面積大于25平方公里，須采用按高斯正形投影的平面直角坐標系。對于面積小于25平方公里的區域，只須使中央子午線大致通過測區中央部分，邊長和方向觀測值的高斯投影改化均甚小，高斯正形投影任意帶的平面直角坐標系幾乎等價于一般的平面直角坐標系。[3]因此，當測區的面積小于25km2，可不進行方向和距離改正，直接把局部地球表面作為平面建立獨立的平面直角坐標系。

3 坐標系轉換的理論分析

相對獨立的兩個坐標系存在聯系時必須確定它們轉換關系。通常將某個坐標系中的控制點作為重合點，利用重合點上兩個坐標系的兩套平面坐標，根據一定的轉換模型，求解轉換參數，實現兩套坐標系的相互轉換。

由于高斯坐標系并非是尺度處處相同的平面坐標系，不同的橢球、不同的橢球定位和不同的中央子午線位置都會造成兩種坐標系之間的尺度差異，采用平面上的相似變換會造成較大的模型誤差。[3]因此，當區域較大時，首先須進行橢球變換，通過高斯投影反算公式，計算平面坐標所對應的大地坐標經緯度，然后通過橢球模型轉換將該大地經緯度轉換為另一坐標系所采用的橢球模型的經緯度，通過高斯正算公式，計算高斯平面坐標。通過以上步驟，將不同橢球模型投影的高斯平面坐標轉換為同一橢球模型下的高斯投影平面坐標。最后在同一高斯平面上進行坐標系的轉換。

對于小區域坐標系的轉換而言，當區域較小時，邊長和方向觀測值的高斯投影改化均甚小，高斯正形投影任意帶的平面直角坐標系幾乎等價于一般的平面直角坐標系，因此可以將坐標系視為假定的平面直角坐標系，在平面上直接進行轉換。

平面坐標系的轉換一般采用四參數的相似變換模型，即采用平移、旋轉、縮放實現坐標系的相互轉換。在數學直角坐標系下，相似變換公式如公式(1)所示:

(1)

相似變換模型的四個參數△X、△Y、m、α，其中△X和△Y為兩個坐標系的原點偏移量;m為兩坐標系的縮放因子，當采用同一長度單位時，m值為1，無需進行縮放;α為角度旋轉因子。

測量坐標系是左手坐標系，因此其變換公式與數學直角坐標系的轉換公式稍微有些不同，考慮到m值為1，在此省略，則測量坐標系下相似變換如公式(2)所示:

(2)

該方式在平面上進行坐標系的轉換，參數簡單，不需跌代。上述公式一共三個參數，而根據一個點的兩套坐標可以列出兩個式子，因此只需已知兩個點的兩套坐標即可求出該轉換模型的三個參數。

如圖1 所示，已知兩點P1，P2，假設在地方獨立坐標系中的坐標分別是:P1(X1，Y1)，P2(X2，Y2)，在獨立工程坐標系下的坐標分別是:P1(x1，y1)，P2(x2，y2)，如公式(3)所示，則有:

(3)

根據圖2所示，則有公式(4):

(4)

根據公式(3)、(4)，已知兩點坐標即可計算出所有轉換參數，從而實現兩個坐標系的相互轉換。

為了減弱點坐標誤差對轉換參數的影響，同時也考慮到轉換模型存在一定的誤差，因此一般利用多個點的兩套坐標，列出多余方程，采用最小二乘平差法計算轉換參數，以保證轉換參數精度。

4 應用分析

該市已建立完善的CORS系統，可以通過VRS快速、精確的獲得該市任一點的地方獨立坐標系的坐標。

該工程所處地勢平緩，面積較小(遠小于25平方公里)，因此其獨立工程坐標系未進行方向和距離改正，直接將局部地球表面作為平面建立了獨立的平面直角坐標系，即假定平面直角坐標系。而在該工程范圍內的地方獨立坐標，因高斯投影后的距離和方向改化較小而可以忽略不計，也可以看做是在一個平面上的直角坐標系。因此，地方獨立坐標系與獨立工程坐標系的轉換采用了在平面上的直角坐標系的相似變換方法。

(1)首先根據“長邊定向”原則，在該獨立工程坐標系內選取了已知獨立工程系坐標的四個控制點。

(2)利用VRS獲得控制點的WGS84坐標，然后轉換為地方獨立坐標系坐標。表1為選取的控制點的兩套坐標。

(3)根據上述公式(3)、(4)，每兩個點即可求出一套參數，然后對這些參數加權取平均。由表1四個控制點求得參數:

△X=403944.245;△Y=206484.924;

sinα=-0.901562119;

cosα=0.432649679

(4)根據所解參數和相似變換公式(2)可知兩個坐標系的轉換關系，實現兩個坐標系的相互轉換。

地方坐標系向工程坐標系的轉換關系如公式(5)所示:

(5)

工程坐標系向城市坐標系的轉換關系如公式(6)所示:

(6)

表2列出了部分檢驗點的工程坐標和通過由地方坐標轉換計算的工程坐標及其誤差。

從上表2轉換誤差可以看出，坐標轉換的精度都在3mm以內，轉換精度比較高，能夠滿足工程需要。

5 總結

以上分析表明，對于較小區域的坐標轉換，以控制點作為重合點，利用重合點上兩套平面坐標，根據相似變換原理，求得平移，旋轉和尺度四個參數，在平面上進行坐標轉換，精度可以滿足使用需要，且該轉換方法簡單易懂，易于實現。但是要保證轉換數據成果的質量，就必須保證已知坐標和所解轉換參數有一定的有效數字位數，且必須以“長邊定向”原則選擇控制點。

參考文獻

[1]孔祥元，郭際明，劉宗泉.大地測量學基礎[M].武漢大學出版社，2001.

[2]潘正風，楊正堯，等.數字測圖原理與方法[M].武漢大學出版社，2004.

[3]施一民.現代大地控制測量[M].測繪出版社，2008.

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

[5]暢開獅.建立城市獨立坐標系相關問題的探討[J].城市勘測，2008，Vol.01:86-90.

[6]韓劍虹，等.地方獨立坐標系與1980西安坐標系之間的轉換[J]，Vol.0.:78-80.

[7]張正祿，等.工程測量學[M].武漢大學出版社，2005.

[8]劉文強.GPS RTK技術在城市工程測量中的應用探討.科技創新導報，2010(1).

[9]楊景彬.深入探討GPS城市控制網現代化改造方法[J].科技創新導報，2010(2).

[10]齊江.深度探討城市工程測量中坐標系的選用思路[J].科技創新導報，2010(2).

“本文中所涉及到的圖表、公式、注解等請以PDF格式閱讀”