基于重心基準的平面坐標轉換

范方標

（深圳市巖土綜合勘察設計有限公司，廣東深圳518100）

基于重心基準的平面坐標轉換

范方標

（深圳市巖土綜合勘察設計有限公司，廣東深圳518100）

在工程測量領域，實現新舊坐標系之間的順利轉換是充分利用以前測繪成果的關鍵，是繪制各種比例尺地圖的基礎。本文以平面坐標轉換為例，介紹平面直角坐標轉換模型，并根據實際生產實踐分析該模型中存在的病態現象，據此提出一種重心基準的平面坐標轉換模型，同時將這兩種模型進行對比分析，提出了修正平面直角坐標轉換模型的方法，并通過工程實例分析轉換結果和精度。

坐標轉換；重心；參數；矩陣

0 引 言

我國大地坐標系主要有北京54坐標和西安80坐標系，后者只對一等點和二等點進行布測，因此對三、四等點的布測只能利用北京54坐標。然而在實際測繪工程中，當需要對三、四等點進行布測時，要將西安80坐標轉換為北京54坐標，因此選擇正確的轉換模型是保證轉換可靠性的基礎，關系到測量結果的精確性。

1 坐標轉換概述

坐標轉換包括參心坐標轉換、不同地心坐標轉換，參心坐標與地心坐標之間的轉換、相同坐標系的直角坐標與大地坐標轉換、大地坐標與高斯平面坐標轉換等[1］。對二維坐標系轉換來說：如果空間中的兩個直角坐標系原點相同，在經過三次旋轉后，可實現兩個坐標系的重合；若兩個坐標系的原點處于不同位置，那么通過坐標軸的平移和旋轉可實現兩個坐標系的重合；若兩個坐標系之間還存在尺度的不同，那么則需要進行尺度縮放。因此二維坐標系的轉換參數包括平移、旋轉、尺度等。

由于轉換模型非常嚴密，故而在實際生產實踐中，利用數學轉換模型轉換北京54、西安80、WGS－84等坐標系而得到的結果與實際存在較大差異。出現這一現象的原因是：由于自由網平差缺少起算數據，因此往往通過附加條件、約束條件解決平差基準問題。通過對比分析坐標轉換模型和重心基準條件方程發現，二者在形式上非常類似。因此，如何將重心基準條件應用到坐標轉換中，是生產實踐中需要重點解決的課題。

2 轉換矩陣狀態

現假設某點在原坐標系下的坐標為（x1，y1），對應到新坐標系下的坐標為（X，Y），而原坐標系的原點對應到新坐標系的坐標為（x0，y0），同時原坐標系與新坐標系的軸夾角為α，尺度參數為m，因此轉換模型為

X＝x0＋x1mcos a－y1msin a，

Y＝y0＋y1mcos a＋x1msin a.（1）

為簡化矩陣的形式，現假設x0＝a，y0＝b，mcosα＝c，msinα＝d，則

X＝a＋x1c－y1d，

Y＝b＋y1c＋x1d.（2）

可將轉換模型寫成矩陣形式

式中：a、b表示平移參數；c、d 表示縮放參數和旋轉參數綜合作用下的參數，而要想解出這4個未知參數，至少需要這兩個坐標系的2個公共點。假設這兩個坐標系中有k（k＞2）的公共點，那么從前面的矩陣中可以得出2k個誤差方程，利用最小二乘原理求解4個轉換參數?，F假設：

寫成

AX＝L＋V，（4）

式中：A表示系數矩陣；X表示需要求解的未知參數；V表示新坐標系坐標的改正數。可以看出；A為列滿軼，根據VTPV＝min原則可求解出未知參數為

X＝（ATPA）－1ATPL＝N－1W，（5）

式中：

W＝ATPL；N為四階滿軼正定矩陣；P為一般單位矩陣，則計算出來的未知參數是唯一的；但是假如N是病態的，則解出的未知參數不可靠。

在實際工作中，平面直角坐標值x、y比較大，當以m為單位時，一般有6～7位數；而在系列矩陣A中，其他元素為0或1，因此當（x1，y1）的絕對值較大時，其列向量的長度差別非常大，使得A的矩陣條件數過大，進而使得N條件數更大，這時候的矩陣呈現出病態現象，解出的未知參數不準確，平面坐標轉換不可靠。而在病態的矩陣中，即使是非常微小的擾動，都可能產生較大的誤差，最終導致得出的參數不準確，使得兩個坐標系難以通過旋轉、位移、縮放等重合到一起[2］。

3 重心基準轉換模型分析

基于重心基準的平面坐標轉換矩陣為ST：

矩陣中，（x0i，y0i）表示近似坐標，其具體的計算條件為

在這個計算條件中：前面兩個說的是兩個坐標系的重心坐標不變；第三個說的是重心到其他各點的帶權向徑方位角平均數不變，這里的權是重心到各點距離的平方；第四個說的是重心到其他各點的帶權向徑長度的平均值不變。從這四個條件中可看出，該平面坐標的轉換模型是以重心為基準的，將重心作為不變量，因此中心是固定的。

將這個重心基準的轉換矩陣與前文得到的坐標轉換矩陣對比分析發現，二者的區別僅在于行向量順序的變化。如果不考慮轉換模型參數的物理性質，那么這兩個矩陣在形式上是一致的。

在前文提到的矩陣之所以造成轉換不準確是因為：在實際工作中，x和y的取值整數位數在6～7位之間，而坐標之間的平移參數整數位數只有2～3位，因此平移量非常微小，容易被人忽視。同時，平移參數對應的公共點數量非常少，與旋轉和尺度參數對應的公共點相差甚遠。研究表明：重心基準坐標轉換模型在平差前后，其新坐標系和原坐標系的重心都沒有任何變化，只不過平面直角坐標轉換模型公共點轉換前后的平均值之間相差一個常數，即平移參數[3］。因此，平面直角坐標轉換模型解決的主要問題還是基準問題。

從上述平面直角坐標轉換模型出現病態的原因來分析，解決病態的根本方法是消除掉平移參數。由于在坐標系的某一個特定區域中，轉換前后的公共點重心與相對位置保持不變，因此可以建立約束條件：

通過方程的聯合求解，消除其中的平移參數（a、b），這樣平面直角坐標轉換模型中就只有公共點的坐標向量。此時，根據約束條件判斷該模型是否處于良態，進而求解出縮放和旋轉綜合作用下的c、d兩個參數，從而實現兩個坐標系的可靠轉換。

在當前的測繪領域，坐標轉換多采用平面直角坐標轉換模型，又被稱為平面四參數模型，通過求解出a、b、c、d這四個未知參數，實現兩個坐標系的順利轉換重疊。而本文提出的重心基準平面坐標轉換模型將重心作為固定的不變量，通過約束條件對平面四參數模型進行相應的修正，消除a、b兩個參數，使得轉換模型處于良態，保證轉換結果的精確度，確保工程測量工作的順利開展[4］。

4 實例計算

某地測繪的數據原先采用的是BJ－54坐標系，將數據成果轉換為80坐標系，規定比例尺為1∶10 000，精度達到1m以上。具體的轉換計算可以參考下表，表1示出了以六個公共點為基準，進行參數成果轉換的結果。

表1 原始數據

根據式（3）和式（5）采用最小二乘法計算；根據式（3）、式（8）、式（5），利用附加重心基準對平移參數予以消去，進而計算。結果如表2所示。

表2 矩陣條件數和轉換參數結果

由表2可知，根據最小二乘方法進行計算，設計矩陣條件參數很大，該矩陣嚴重病態，根據式（2）計算出尺度系數，得到：m＝8 134.9＞1，說明BJ－54的坐標和80坐標之間的尺度變化率較小。而通過第二種方案進行計算，設計矩陣與方程條件數均為1，屬于良態矩陣，能夠予以求逆計算，得到旋轉參數之后，就可以對公共點轉換殘差予以計算，得到最大點位誤差是0.388m，滿足應用需求。

除此之外，從數據精化需求角度而言，其整體殘差較大，特別是第3、第6點位，誤差超過了0.3 m，此兩點間可能存在粗差，需要予以剔除，重新計算殘差，利用剔除較大殘差后計算的轉換參數進行轉換精度。

5 結束語

隨著測繪行業的發展，社會生活中的平面坐標轉換將會越來越多，研究平面坐標轉換模型將成為該領域重要發展方向。利用經典最小二乘理論計算坐標轉換參數在理論非常嚴密，但是設計矩陣的病態情況對解算結果有著一定的影響。通過對設計矩陣病態的分析，根據經典秩虧自由網平差，對比其前后重心的相對性，合理解釋設計矩陣病態消除的附加條件，保證結果的可靠性。對此，本文通過對平面四參量轉換模型病態的分析，利用重心基準平面坐標轉換模型對病態進行修正，通過增加約束條件以消除a、b參數，從而確保坐標轉換轉換的可靠性。

[1]韋 鋮，馬文雙，李明君，等.工程測量中平面坐標轉換軟件設計及應用[J］.測繪工程，2013，22（4）：76－79.

[2]陳紹杰，秦永洋，孫興平，等.顧橋礦區GPS基準網高程系統轉換研究[J］.測繪，2010，33（1）：34－38.

[3]彭愛文，曹佩瑤.平面坐標轉換方法探討及轉換軟件的設計思路[J］.測繪與空間地理信息，2007，30（3）：189－191.

[4]劉陶勝，黃聲享，羅 力，等.基于重心基準的平面坐標轉換研究[J］.大地測量與地球動力學，2011，31（2）：102－106.

On Plane Coordinate Transformation Based on Center of Gravity Datum

FAN Fangbiao
（Shenzhen Integrated Geological Exoloration ＆Design Co.Led，Shenzhen518100，China）

in the field of engineering measurement，to achieve a smooth transition between the old and the new coordinate system is to make full use of the key before the surveying and mapping results，is the basis for mapping various scale map.This paper takes the plane coordinate transformation as an example，introduces the plane right angle coordinate transformation model，and according to the analysis of morbid phenomena exist in the model of the actual production practice，plane coordinate transformation model proposes a benchmark，at the same time，by a comparative analysis of these two kinds of models，and put forward the method of correcting plane right angle coordinate transformation model，and through the the analysis of engineering example the conversion results and precision.

Coordinate transformation；center of gravity；parameters；matrix

P226.3

A

1008－9268（2015）01－0079－04

10.13442/j.gnss.1008－9268.2015.01.018

范方標（1978－），男，工程師，研究方向為GPS技術。

2014－06－10

聯系人：范方標E－mail：1330843245＠qq.com