蘭溪市國土數據資料坐標轉換工作的思考

章永斌

（蘭溪市測繪事務所，浙江蘭溪321100）

1 引言

依據國家測繪地理信息局發布的公告，2000國家坐標系從2008年7月1日開始啟用，2000國家坐標系轉換與現行國家坐標系、銜接的需要8-10年過渡期；在過渡期內現有地理信息系統應逐步轉換到2000國家坐標系，2008年7月1日后新建的地理信息系統，必須采用2000國家坐標系。由于1980西安坐標系已經采用20多年，蘭溪市國土資源局從1998年開始啟用1980西安坐標系測繪數據，工作中形成了大量的國土數據資料和相關的信息系統采用了1980西安坐標系測繪數據成果，于是面臨著大量的國土數據資料的坐標轉換問題。

2 我國不同時期所采用的坐標系

在上世紀50年代和80年代分別創建了1954年北京坐標系與1980西安坐標系，測制了各種比例尺地形圖，由于受當時技術條件限制，傳統技術是實現中國大地坐標系的主要依賴手段。1954北京坐標系采用的是克拉索夫斯基橢球體，該橢球長軸為6378245M，短軸為6356863.0188M，扁率為1/298.3，坐標起算原點在前蘇聯的普爾科沃，該系統在定位和計算過程中，沒有采用中國的數據，它在中國范圍內不適合，不能滿足高精度定位以及空間科學、地球科學和戰略武器發展的需要。

在上世紀70年代，經過我國大地測量工作者們二十多年的艱巨努力，終于完成了全國一、二等天文大地網的布測。經過整體平差，采用1975年IUGG第十六屆大會推薦的參考橢球，建立了1980西安坐標系，該橢球長軸為6378140M，短軸為6356755.2882M，扁率為1/298.257，坐標起算原點在我國中部陜西省涇陽縣永樂鎮。它在定位和計算過程中，采用了中國的數據，該系統在中國范圍內非常符合，在中國建設發展和科學研究領域發揮了巨大作用。

20088年4月國務院正式批準2000國家坐標系是我國當前最新的國家坐標系，它的原點為包括大氣、海洋的整個地球的質量中心，Z軸由原點指向歷元2000.0的地球參考極的方向，X軸由原點指向格林尼治參考子午面與歷元2000.0地球參考橢球的赤道面的交點，Y軸與Z軸、X軸構成右手正交坐標系。該橢球長軸為6378137M，短軸為6356752.3141M，扁率為1/298.257。

當前GPS發展過程中，國際上以及手機GPS用戶采用WGS-84坐標系，該坐標系由于框架、歷元混亂，因此其精度較底。它與2000國家坐標系相比實質是一樣的，屬于地心坐標系，采用的地球參考橢球非常接近，扁率差異所引起的橢球面上的緯度和高度變化最大達0.1MM，當前測量精度范圍內，其差異可以忽略，也可以說該坐標系與2000國家坐標系相容至CM級水平。

3 坐標投影

地球上的點用經緯度表示的是地理坐標系，也稱大地坐標系。由于地理坐標系不夠方便，所以人們比較習慣于使用平面坐標系，平面坐標系分別用xy表示。把球體表面的坐標轉成平面坐標需要一定的手段，這個手段稱為投影。由于投影方法也不是唯一的，但是投影就是為了一個目的，務求使當地坐標最準確。目前就有好多投影方法，比如高斯投影、墨卡托投影等。

我國的投影坐標系主要采用高斯-克呂格投影，分為6度和3度分帶投影，1：2.5萬-1：50萬比例尺地形圖采用經差6度分帶，1：1萬比例尺的地形圖采用經差3度分帶。具體分帶法是：6度分帶從本初子午線（格林尼治天文臺所在的子午線）開始，按經差6度為一個投影帶自西向東劃分，全球共分60個投影帶，中國跨13-23帶;3度投影帶是從東經1度30分經線(1.5°)開始，按經差3度為一個投影帶自西向東劃分，全球共分120個投影帶，中國跨25-45帶。該投影的特點：在投影面上，中央子午線和赤道的投影都是直線，并以中央子午線和赤道的交點O作為坐標原點，以中央子午線的投影為縱坐標軸，以赤道的投影為橫坐標軸，形成高斯平面直角坐標系。因此這種投影，將中央經線投影為直線，其長度沒有變形，其余經線為凹向中央子午線的曲線，距離中央經線愈遠，變形愈大；赤道線投影后是直線，但其長度有變形，除赤道外的其余緯線，投影后為凸向赤道的曲線，并以赤道為對稱軸；經線與緯線投影后保持正交。所以高斯-克呂格投影采用分帶投影，目的就是使投影帶邊緣的變形不致過大。

蘭溪建成區位于東經119°18，按照國家高斯-克呂格投影3度分帶，其中央子午線為東經120°，由于蘭溪建成區東經119°18，落在東經120°中央子午線投影帶邊緣，其國家控制點坐標應用到工程測量上，投影變形過大，歷史上蘭溪采用將東經120°中央子午線投影帶的國家控制點坐標平移轉換到東經119°18，作為中央子午線的投影帶坐標，建立蘭溪獨立坐標系，以消除因投影帶邊緣變形過大對工程測量的影響。因此蘭溪市國土數據資料所涉及的坐標系一般都是在東經119°18，作為中央子午線的投影帶下，將測量數據的坐標進行整體平差后，再把坐標數據平移轉換到東經120°中央子午線投影帶下的數據坐標，提交上述數據資料。

4 Bursa七參數坐標轉換方法和流程

坐標轉換就是利用擁有兩種坐標系坐標的重合點，采用坐標轉換模型，利用最小二乘法原理，得到相應的目標坐標系的坐標。由于Bursa七參數轉換模型理論比較嚴密，一般測量界坐標轉換大多采用Bursa七參數為三維模型，其轉換模型有3個平移參數，3個旋轉參數和1個尺度參數。本次蘭溪市國土數據資料坐標轉換工作也采用Bursa七參數轉換模型。

重合點的選取一般遵循高等級、高精度、分布均勻、覆蓋整個轉換區域等原則，并盡量避免選取變形或沉降較大區域的點，因此本次蘭溪市國土數據資料坐標轉換工作主要利用城鎮地籍數字化測量時上交的地籍一級控制點，并對所有的地籍一級控制點重新委托測量2000國家大地坐標系的坐標。按照坐標轉換選取坐標的重合點的要求，每批6個重合點，6個重合點盡量均勻分布在蘭溪市市域范圍，采用Bursa七參數為三維模型，采用最小二乘法原理求解轉換七參數；運用求解的轉換七參數對其他的坐標的重合點進行轉換，得到求算的2000國家大地坐標系的坐標，并將求算的2000國家大地坐標系的坐標與實測的2000國家大地坐標系的坐標進行比較，計算重合點坐標殘差，剔除殘差大于3倍單位權中誤差的重合點，直到重合點坐標殘差小于3倍單位權中誤差的重合點，確定重合點計算轉換七參數；使用其它地籍一級控制點未參與計算轉換參數的重合點作為外部檢核點，利用已確定的轉換七參數計算這些外部檢核點2000國家大地坐標系的坐標，與該外部檢核點的實測2000國家大地坐標系的坐標進行比較，并求算坐標轉換殘差。殘差檢核公式是否符合精度要求，當絕大多數外部檢核點坐標轉換殘差符合精度要求，經過大量的試算與分析剔除一些變化大的重合點，最后確定轉換七參數，并計算轉換中誤差，并上報市地理信息管理處存檔，便于今后其他兄弟單位坐標轉換使用；利用最后確定的轉換七參數，將所有的地籍二級控制點和地籍圖根點坐標轉換成2000國家大地坐標系的坐標。

4.1 轉換參數的計算步驟

①按照轉換區域，選取合適的轉換模型；

②通過選取重合點，計算轉換參數；

③計算重合點坐標殘差，剔除殘差大于3倍單位權重誤差的重合點；

④重復上述②和③，使重合點坐標殘差均小于3倍單位權重誤差，并且計算轉換參數的重合點數量根據轉換模型確定。本次采用七參數轉換模型時，重合點不少于6個；

⑤由確定的重合點計算轉換參數。

4.2 坐標轉換

通過利用轉換參數進行坐標轉換，而求得各控制點在2000國家大地坐標系坐標。

v(轉換殘差)=已知點2000國家大地坐標系的平面坐標-轉換求算的該點2000國家大地坐標系平面坐標

X坐標轉換中誤差：MX=∑△vx2/n

Y標轉換中誤差：MY=∑△vy2/n

轉換中誤差：M=(MX2+MY2)/2

4.3 外部符合精度檢核的方法

①把未參與計算轉換參數的重合點作為外部檢核點，其中點數一般不少于6個并且均勻分布；

②利用轉換參數計算外部檢核點的坐標與該外部檢核點的已知坐標進行比較，檢核公式同坐標轉換殘差計算公式相同。

5 蘭溪市國土歷史數據資料坐標的轉換

實現空間矢量數據、柵格數據、文本文件數據、Oracle數據庫空間矢量數及切片影像數據坐標轉換。依據實際工作需要將全系統的數據分出層次，分步驟向2000國家大地坐標系進行坐標轉換。首先優先轉換國土資源管理基礎平臺數據，滿足國家規定時間使用2000國家大地坐標后的國土資源管理要求。其次對使用頻率較低的歷史檔案數據在后續的工作中逐步轉換。同時對已有地籍圖根點以上的控制點，用Bursa七參數坐標轉換方法和流程求解的參數，一次性轉換成國家2000坐標系，便于今后工程測量的需要，同時對Bursa七參數坐標轉換方法和流程所求解的參數報備測繪管理部門，便于其他測量單位在蘭溪區域范圍內測量坐標轉換使用。最后，一直以來，蘭溪市各部門北京54坐標系轉換西安80坐標系已有固定的轉換參數，為避免坐標系多次轉換造成數據銜接誤差，在蘭溪市各部門北京54坐標系轉換2000國家大地坐標系的實際工作中，應該按照固定的轉換參數將蘭溪市各部門北京54坐標系轉換成西安80坐標系，再將數據按80系與2000國家大地坐標系的新參數進行轉換。

矢量數據管理單元的轉換：

①首先取得80系各要素的坐標，逐一計算其2000系下各要素的坐標。

②并且將2000系下各要素的坐標寫回原要素。

③添加2000系下新的方里格網層及標注，刪除原方里格網數據層及標注。

該數據為CGCS200，其范圍為原80系范圍。

④數據編輯、鄰帶之間接邊、拓撲重建。

⑤先進行數據相互拼接，再按2000系標準圖幅范圍進行裁切。

⑥更新相關字段屬性值。

⑦在元數據中，增加2000系下4個圖廓角點坐標，橢球長半徑、橢球扁率、所采用的大地基準以及產品更新日期。

矢量數據空間數據庫的轉換：連接后臺數據庫，先把每個要素類加載空間數據庫中，再讀取各要素80系坐標，然后逐一計算2000系下各要素坐標，將2000系下的要素存儲到空間數據庫中新建的要素類，具體方法如下：

①新建一個新的要素類要與原要素類結構相同；

②通過先獲取各要素80系坐標，來逐點計算2000系下各要素的坐標；

③將2000系下要素寫入新建要素類；

④先添加2000系下新的方里格網層，再刪除原方里格網數據層；

⑤更新相關字段屬性值。

柵格數據的轉換：

①對1：1萬標準分幅DOM讀取80系頭文件中定位坐標，計算DOM中心點坐標“改正量”，按照像素關系計算移動量（像素數），避免圖幅之間接邊數據重新采樣；

②先進行數據相互拼接，再按2000系標準圖幅范圍及重疊像素進行裁切；

③DOM數據勻色處理；

④更改數據頭文件中定位坐標；

⑤在元數據中，增加2000系圖廓坐標范圍，橢球長半徑、橢球扁率、所采用的大地基準以及產品更新日期。對標準分幅DOM讀取80系頭文件中定位坐標，計算DOM中心點坐標“改正量”，按照像素關系計算移動量（像素數），避免圖幅之間接邊數據重新采樣；更改數據頭文件中定位坐標。

通過對歷年的國土調查矢量數據庫和地籍矢量數據庫以及歷年航片、衛片進行坐標轉換，以達到國家對國土數據上報必須使用國家2000坐標系的要求。