謝偉杰

（汕頭市潮南區(qū)規(guī)劃測(cè)繪大隊(duì)，廣東 汕頭515144）

1 背景介紹

坐標(biāo)轉(zhuǎn)換的途徑較多，最核心的是經(jīng)典七參數(shù)法，與之相對(duì)應(yīng)的數(shù)學(xué)模型類型豐富，如Bursa 模型等，各類模型均對(duì)坐標(biāo)精度、高程系統(tǒng)等提出較高的要求。但從實(shí)際工程狀況來看，資料缺失的問題尤為普遍，致使計(jì)算結(jié)果缺乏可靠性。相較之下，平面四參數(shù)法的限制性條件偏少，在坐標(biāo)轉(zhuǎn)換中有著更強(qiáng)的適用性。

1954 北京坐標(biāo)系對(duì)于我國(guó)測(cè)繪行業(yè)而言具有引導(dǎo)性意義，后續(xù)多數(shù)城市則創(chuàng)建了獨(dú)立坐標(biāo)系，但多數(shù)成果均是在傳統(tǒng)坐標(biāo)系技術(shù)的基礎(chǔ)上而衍生出的全新形式。社會(huì)持續(xù)發(fā)展之下，傳統(tǒng)坐標(biāo)系的適用性逐步下降，因此，我國(guó)對(duì)2000 國(guó)家大地坐標(biāo)系的應(yīng)用力度不斷加大，此時(shí)則存在新舊坐標(biāo)系的數(shù)據(jù)轉(zhuǎn)換問題【1】。為解決該問題，測(cè)繪行業(yè)中應(yīng)用最為典型的有坐標(biāo)轉(zhuǎn)換方式，本文引入平面四參數(shù)模型，分析其在坐標(biāo)轉(zhuǎn)換中的應(yīng)用效果。從轉(zhuǎn)換公式與算例計(jì)算2個(gè)角度切入，分析四參數(shù)模型下實(shí)現(xiàn)數(shù)據(jù)轉(zhuǎn)換的具體操作機(jī)制。

2 坐標(biāo)轉(zhuǎn)換的基本內(nèi)容

坐標(biāo)轉(zhuǎn)換可細(xì)分為以下2 類：

1）相同基準(zhǔn)下的坐標(biāo)轉(zhuǎn)換，即建立在同一坐標(biāo)參照系的基礎(chǔ)上。此時(shí)，坐標(biāo)表現(xiàn)形式可細(xì)分3 類，即空間直角坐標(biāo)（X，Y，Z），經(jīng)緯度+高程（B，L，H），平面坐標(biāo)+高程（x，y，h）。根據(jù)此特點(diǎn)，若要實(shí)現(xiàn)坐標(biāo)轉(zhuǎn)換，最為關(guān)鍵的是完成三相坐標(biāo)系之間的靈活轉(zhuǎn)換。

2）不同基準(zhǔn)下的坐標(biāo)轉(zhuǎn)換，其基本特點(diǎn)在于各自對(duì)應(yīng)的坐標(biāo)參考系存在差異。從該特點(diǎn)來看，坐標(biāo)轉(zhuǎn)換的流程更為復(fù)雜，統(tǒng)一坐標(biāo)表現(xiàn)形式是最為基礎(chǔ)的環(huán)節(jié)，隨后再實(shí)現(xiàn)坐標(biāo)間的轉(zhuǎn)換，其核心在于確定與各基準(zhǔn)相適應(yīng)的轉(zhuǎn)換參數(shù)。

由于基準(zhǔn)的不同，坐標(biāo)表示形式也存在差異，因此，要達(dá)到坐標(biāo)表示形式統(tǒng)一的效果，創(chuàng)建相同的基準(zhǔn)條件，最終實(shí)現(xiàn)坐標(biāo)轉(zhuǎn)換。應(yīng)根據(jù)待轉(zhuǎn)換區(qū)域的特點(diǎn)，挑選與之相適應(yīng)的轉(zhuǎn)換模型，而根據(jù)等級(jí)的不同，所對(duì)應(yīng)的模型選擇結(jié)果也存在差異：針對(duì)全國(guó)及省級(jí)的坐標(biāo)轉(zhuǎn)換，較為合適的是七參數(shù)轉(zhuǎn)換模型；針對(duì)省級(jí)以下的坐標(biāo)轉(zhuǎn)換，如平面四參數(shù)模型則具有可行性，其在應(yīng)用中限制條件相對(duì)較少。具體到本文的坐標(biāo)轉(zhuǎn)換分析工作中，選擇了平面四參數(shù)模型。

3 坐標(biāo)轉(zhuǎn)換過程

3.1 技術(shù)路線

坐標(biāo)轉(zhuǎn)換應(yīng)嚴(yán)格遵循特定的流程依次執(zhí)行，最為基礎(chǔ)的是挑選轉(zhuǎn)換模型，再通過對(duì)坐標(biāo)重合點(diǎn)的分析，明確可行的模型轉(zhuǎn)換參數(shù)，通過對(duì)模型殘差的分析，可掌握其精度水平。根據(jù)實(shí)際結(jié)果，若精度評(píng)估未達(dá)到設(shè)計(jì)標(biāo)準(zhǔn)，表明所選擇的模型缺乏可行性，需再次挑選模型，重復(fù)執(zhí)行此流程直至滿足精度要求，最后再完成成果的轉(zhuǎn)換。

3.2 模型選擇

本文選擇的是平面四參數(shù)模型，選取某一P 點(diǎn)，為實(shí)現(xiàn)該點(diǎn)在O1、O2坐標(biāo)系中的轉(zhuǎn)換，提出了轉(zhuǎn)換模型公式，坐標(biāo)轉(zhuǎn)換示意圖如圖1 所示。具體分析如下。

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

將P 點(diǎn)對(duì)應(yīng)至O1坐標(biāo)系中，將其表示為P1（x1，y2）；類似的，對(duì)應(yīng)至O2坐標(biāo)系中，可得到P2（x2，y2）。根據(jù)2 坐標(biāo)系的特點(diǎn)，提出尺度變化因子，即1+m。分析P1向P2轉(zhuǎn)換的具體流程：

經(jīng)過轉(zhuǎn)換后產(chǎn)生的坐標(biāo)可按式（1）、式（2）計(jì)算：

式中，x2、y2為經(jīng)過轉(zhuǎn)換后產(chǎn)生的坐標(biāo)；x0、y0為平移量；m 為尺度參數(shù)；x1、y1為轉(zhuǎn)換前的坐標(biāo)；α 為旋轉(zhuǎn)參數(shù)（因旋轉(zhuǎn)方向的不同而發(fā)生變化，即順時(shí)針為正，逆時(shí)針為負(fù)）。將式（1）、式（2）轉(zhuǎn)變?yōu)榫仃囆问?，則有：

選取O1坐標(biāo)系，將其進(jìn)行旋轉(zhuǎn)操作，采取順時(shí)針轉(zhuǎn)α 的方式，使其與O2坐標(biāo)系保持相平行的關(guān)系，得到平面四參數(shù)計(jì)算公式，即為式（3）。

可以得知，該轉(zhuǎn)換適應(yīng)的是二維平面條件；若要實(shí)現(xiàn)三維坐標(biāo)的轉(zhuǎn)換，則要?jiǎng)?chuàng)建同一基準(zhǔn)，在此條件下獲得平面坐標(biāo)，再確定合適的轉(zhuǎn)換參數(shù)。換言之，如坐標(biāo)轉(zhuǎn)換發(fā)生在不同基準(zhǔn)的條件下，此時(shí)必須做到轉(zhuǎn)換條件的一致性，而關(guān)鍵方法在于創(chuàng)建相同的基準(zhǔn)。

3.3 重合點(diǎn)的選取

坐標(biāo)轉(zhuǎn)換過程中，應(yīng)對(duì)比分析2個(gè)坐標(biāo)系的特點(diǎn)，確定的坐標(biāo)成果點(diǎn)應(yīng)同時(shí)分布在各坐標(biāo)系中，將滿足此條件的點(diǎn)視為重合點(diǎn)。但該點(diǎn)并非絕對(duì)，需選擇轉(zhuǎn)換參數(shù)，根據(jù)該值求得重合點(diǎn)坐標(biāo)殘差，通過所得結(jié)果分析殘差值與中誤差的關(guān)系，若出現(xiàn)了殘差值超過3 倍中誤差的情況，表明該轉(zhuǎn)換參數(shù)缺乏可行性，需分析后重新確定【2】。同時(shí)，轉(zhuǎn)換區(qū)域大小將直接決定重合點(diǎn)數(shù)量，為提升轉(zhuǎn)換結(jié)果精確性，要求該點(diǎn)數(shù)量至少為5個(gè)。

3.4 計(jì)算轉(zhuǎn)換參數(shù)

經(jīng)上述流程后可得到重合點(diǎn)坐標(biāo)，在此基礎(chǔ)上引入平面四參數(shù)模型，經(jīng)分析后求得模型參數(shù)，此處選擇的是最小二乘法。

3.5 精度評(píng)估與檢核

數(shù)據(jù)準(zhǔn)確性是坐標(biāo)轉(zhuǎn)換最為基礎(chǔ)的目標(biāo)，因此要確保數(shù)據(jù)精度，坐標(biāo)轉(zhuǎn)換時(shí)必須加大力度做好精度評(píng)估工作。確定模型參數(shù)后需進(jìn)一步驗(yàn)證，評(píng)定其準(zhǔn)確性，此處需引入重合點(diǎn)殘差的誤差這一指標(biāo)，并選擇適量的外部檢驗(yàn)點(diǎn)，根據(jù)轉(zhuǎn)換參數(shù)求得各點(diǎn)的實(shí)際坐標(biāo)情況，并將其與已知坐標(biāo)對(duì)比，實(shí)現(xiàn)外部檢核。此時(shí)，重合點(diǎn)的選擇尤為關(guān)鍵，要求點(diǎn)數(shù)至少達(dá)6個(gè)，且彼此間具備分布均勻的特點(diǎn)，否則將對(duì)外部檢核結(jié)果的可靠性造成影響【3】。若精度滿足要求，進(jìn)一步完成坐標(biāo)成果的改算處理；若精度偏低，則需要在原坐標(biāo)轉(zhuǎn)換方案的基礎(chǔ)上做出調(diào)整，選擇更為合適的模型，并改進(jìn)重合點(diǎn)的選擇方式，此后再計(jì)算轉(zhuǎn)換參數(shù)，分析其與精度要求的差異，滿足要求后才可用于坐標(biāo)轉(zhuǎn)換工作中。

3.6 坐標(biāo)成果改算

根據(jù)外部重合點(diǎn)計(jì)算精度情況，若滿足要求，此時(shí)計(jì)算所得的4 項(xiàng)參數(shù)具有可行性，可用于待定點(diǎn)的坐標(biāo)轉(zhuǎn)換。此處需全面檢查數(shù)據(jù)，確保無誤后方可提交。根據(jù)數(shù)據(jù)轉(zhuǎn)換需求，利用VB 編寫程序，將其作為數(shù)據(jù)轉(zhuǎn)換的基本工具。

4 算例

確定公共點(diǎn)坐標(biāo)A1、A2、A3的信息，詳細(xì)內(nèi)容見表1。

表1 公共點(diǎn)坐標(biāo)

根據(jù)式（1）、式（2），經(jīng)計(jì)算后可求得4 項(xiàng)參數(shù)，具體有：x0=40.986，y0=-95.817，α=0°0′0.2″，m=1.000。

在待轉(zhuǎn)換的O1與O2坐標(biāo)系中，若得知其中的任意坐標(biāo)，在本文所提出的四參數(shù)模型的支持下，即可實(shí)現(xiàn)在2 坐標(biāo)系間的轉(zhuǎn)換。舉例而言，若O1坐標(biāo)系中存在某點(diǎn)，其坐標(biāo)為x1=4056960.835，y1=536443.590，根據(jù)上述提到的程序，運(yùn)行后即可獲得其在O2下的坐標(biāo)，具體結(jié)果為x2=4056920.168，y2=536499.498。

5 結(jié)語

大數(shù)據(jù)背景下，基于科學(xué)的方式實(shí)現(xiàn)數(shù)據(jù)的改算，可為各行業(yè)的發(fā)展提供充足的支持，于測(cè)繪行業(yè)而言亦是如此。當(dāng)前，基于坐標(biāo)的轉(zhuǎn)換模型類型較為豐富，但各自在實(shí)際應(yīng)用中均具有局限性。坐標(biāo)轉(zhuǎn)換可分為同一基準(zhǔn)和不同基準(zhǔn)2 種條件，彼此間的坐標(biāo)轉(zhuǎn)換并非獨(dú)立存在，而是關(guān)聯(lián)緊密、相輔相成。因此，在后續(xù)的研究工作中，相關(guān)人員可從算法模型的角度入手，實(shí)現(xiàn)二者的統(tǒng)一。此外，影響轉(zhuǎn)化成果精度的因素較多，如橢球的扁率等，而現(xiàn)階段的模型尚未在此方面實(shí)現(xiàn)統(tǒng)一的整合，缺乏固定算法模型的支持，因此，在此方面的研究較為空白。鑒于此，希望學(xué)者在后續(xù)研究中可深化算法模型，使其融合各項(xiàng)影響因素，提升算法模型的適用性，為坐標(biāo)轉(zhuǎn)換提供可靠的指導(dǎo)。