劉洪曉 張國君

（招金礦業(yè)股份有限公司河?xùn)|金礦，山東 招遠265402）

0 引入

在利用手工數(shù)字化儀采集GIS數(shù)據(jù)中，往往由于數(shù)字化儀坐標(biāo)系與地面坐標(biāo)系不一致及圖紙變形而產(chǎn)生系統(tǒng)誤差，為了消除此誤差，通常是根據(jù)已知地面坐標(biāo)的控制點和格網(wǎng)點采用平面相似變換法進行處理。在本文中，筆者基于Excel實現(xiàn)了坐標(biāo)轉(zhuǎn)換的過程。

1 原理

為進行平面坐標(biāo)轉(zhuǎn)換，需要有一定數(shù)量的公共控制點，這些公共點應(yīng)具有兩個坐標(biāo)系中的雙重坐標(biāo)。

已知新舊坐標(biāo)系的坐標(biāo)變換方程為：

式中m為尺度比因子，α為旋轉(zhuǎn)因子。

圖1 坐標(biāo)轉(zhuǎn)換示意圖

如圖 1 所示，令 a=x0，b=y0，c＝mcosα，d＝msinα，則（1）式可寫成：

式中a、b、c、d為所求的未知量，即平差參數(shù)。

設(shè)兩坐標(biāo)系中有 n 個公共點（xi,yi）和（xi′,yi′），i＝1,2,…,n，令新坐標(biāo)系的坐標(biāo)為觀測值，舊坐標(biāo)系中坐標(biāo)設(shè)為無誤差，當(dāng)n＞3時，則可列出誤差方程為：

2 基于Excel的過程實現(xiàn)

2.1 已知坐標(biāo)數(shù)據(jù)

如圖2所示，在Excel中輸入已知點號為1、2、…、7的點在新、舊坐標(biāo)系中的數(shù)據(jù)：

圖2 輸入已知數(shù)據(jù)

2.2 列立誤差方程

依據(jù)式（3），取n＝7，則如圖3所示，誤差方程式如下

圖3 列立誤差方程

2.3 解算誤差方程

如圖4所示，誤差方程的解算如下：

圖4 誤差方程解算

于是可得平差參數(shù) a＝1907.0027，b＝-640.13873，c＝0.37107202，d＝0.9286021。

2.4 坐標(biāo)的轉(zhuǎn)換

為了驗證所計算的平差參數(shù)，如圖5所示，將點號為1、2、…、7的點坐標(biāo)代入，并同時計算了點號為8、9、10的點坐標(biāo)：

圖5 坐標(biāo)轉(zhuǎn)換

通過進行坐標(biāo)轉(zhuǎn)換，參照圖2與圖5的坐標(biāo)數(shù)據(jù)，可以比較點號為1、2、…、7的點在轉(zhuǎn)換前后是一樣的（忽略小數(shù)位的取舍影響）。依此，可以進一步對其他的點進行相應(yīng)的轉(zhuǎn)換。

3 結(jié)論

利用Excel，依據(jù)間接平差原理，本文很好地實現(xiàn)了坐標(biāo)的轉(zhuǎn)換。但在使用時需注意以下幾點：

（1）需要注意公共點的個數(shù)n必須大于3，理論上越大越好；

（2）由于小數(shù)位的取舍影響，使得參與計算的已知點在轉(zhuǎn)換前后的坐標(biāo)略有微小差異；

（3）本文重復(fù)利用了Excel的函數(shù)功能，在進行運算時要明白各種函數(shù)的意思。

［1］武漢大學(xué)測繪學(xué)院測量平差學(xué)科組.誤差理論與測量平差基礎(chǔ)[M].武漢大學(xué)出版社,2003.

［2］王永,等.利用Excel繪制誤差橢圓的方法[J].礦山測量,2008,12.

［3］泥立麗,等.基于Excel的繪制誤差曲線的方法[J].礦山測量,2010,6.

［4］王永,等.基于分形的線要素綜合數(shù)據(jù)處理方法[J].礦山測量,2009,10.