關(guān)于ArcGIS坐標(biāo)轉(zhuǎn)換方法及精度評(píng)估研究

楊騫

關(guān)于ArcGIS坐標(biāo)轉(zhuǎn)換方法及精度評(píng)估研究

楊騫

（甘肅煤炭地質(zhì)勘查院，甘肅蘭州 730000）

在構(gòu)建和運(yùn)用地理信息體系時(shí)，常常會(huì)需要借助坐標(biāo)轉(zhuǎn)換來(lái)完成。在當(dāng)前國(guó)內(nèi)的測(cè)繪數(shù)據(jù)成果資料中，大部分1∶50000地形圖是采取西安80坐標(biāo)系建立的，還有部分地形圖采用北京54坐標(biāo)系完成的，我國(guó)已實(shí)現(xiàn)將國(guó)家現(xiàn)行的大地坐標(biāo)系向國(guó)家2000大地坐標(biāo)系的轉(zhuǎn)變過(guò)程。ArcGIS是建立地理信息體系的主流平臺(tái)，在國(guó)內(nèi)地理數(shù)據(jù)信息的產(chǎn)出、數(shù)據(jù)庫(kù)建立及系統(tǒng)應(yīng)用開(kāi)發(fā)之中都獲得了廣泛采用，從而形成大量根據(jù)ArcGIS建立的矢量信息數(shù)據(jù)。文中深入分析了ArcGIS的坐標(biāo)轉(zhuǎn)換方法，并且基于ArcGIS的坐標(biāo)轉(zhuǎn)換形式，針對(duì)國(guó)內(nèi)常用坐標(biāo)系轉(zhuǎn)換系統(tǒng)，對(duì)點(diǎn)位坐標(biāo)轉(zhuǎn)換精度實(shí)施了評(píng)估。

ArcGIS；坐標(biāo)轉(zhuǎn)換；精度評(píng)估

近年來(lái)，伴隨信息技術(shù)、電子技術(shù)、數(shù)據(jù)庫(kù)技術(shù)、空間信息技術(shù)以及地理信息體系等技術(shù)的不斷發(fā)展與逐漸完善，國(guó)內(nèi)外針對(duì)地理信息體系的開(kāi)發(fā)和應(yīng)用均已邁進(jìn)入了新的發(fā)展階段。當(dāng)前，隨著地理信息數(shù)據(jù)的作用越發(fā)重要，各個(gè)行業(yè)對(duì)于地理信息數(shù)據(jù)的響應(yīng)時(shí)間與精確度也都提出了越來(lái)越高的需求。而圖像圖形數(shù)據(jù)處理是地理信息數(shù)據(jù)之中最重要的一部分，把地理坐標(biāo)向平面體系中進(jìn)行映射，則必須要利用坐標(biāo)轉(zhuǎn)換來(lái)完成。以往常用的坐標(biāo)轉(zhuǎn)換方式精確度較低，且轉(zhuǎn)化的誤差也較大。ArcGIS是GIS的主流系統(tǒng)，被廣泛運(yùn)用于系統(tǒng)開(kāi)發(fā)、數(shù)據(jù)產(chǎn)出及數(shù)據(jù)庫(kù)建立之中，產(chǎn)生了基于ArcGIS系統(tǒng)的大量矢量數(shù)據(jù)信息。在當(dāng)前共同存在多類坐標(biāo)系的條件下，利用ArcGIS進(jìn)行矢量坐標(biāo)數(shù)據(jù)轉(zhuǎn)換的需要越來(lái)越多。文中分析了ArcGIS系統(tǒng)矢量數(shù)據(jù)的坐標(biāo)轉(zhuǎn)換方法，可以達(dá)到ArcGIS系統(tǒng)矢量數(shù)據(jù)無(wú)損、高效、高精確度的轉(zhuǎn)換效果。

1　坐標(biāo)系分析

若要實(shí)施坐標(biāo)系的轉(zhuǎn)換，則首先必須要對(duì)圖形轉(zhuǎn)換坐標(biāo)系及目標(biāo)坐標(biāo)系進(jìn)行了解，且坐標(biāo)系間最主要的差別就在地球橢球參數(shù)與坐標(biāo)投影上，該參數(shù)通常在進(jìn)行坐標(biāo)系定義時(shí)便已給定。

1.1地球橢球參數(shù)分析

在地球測(cè)量時(shí)，用以表示地球的橢球被稱作地球橢球，一般也可稱為橢球，是數(shù)字形式的地球表示方式。對(duì)于幾何參數(shù)特定，且定向、定位的，用來(lái)表示某區(qū)域大地水準(zhǔn)面的橢球被稱為參考橢球。地面上所有觀測(cè)參數(shù)均應(yīng)被歸算于參考橢球上，同時(shí)在該面上開(kāi)展計(jì)算分析過(guò)程。參考橢球面為計(jì)算測(cè)量大地的基準(zhǔn)面，并且也是地球形狀與地圖投影研究的參考面，通常坐標(biāo)體系不同，其所具備的橢球參數(shù)也會(huì)不同（表1中所示）。

表1　不同的坐標(biāo)系統(tǒng)的橢球參數(shù)表

1.2坐標(biāo)投影

一般來(lái)說(shuō)，大地橢球體表面為曲面，所以，為了可以把不能展開(kāi)的曲面繪制到平面上，通常需要采取某種數(shù)學(xué)方法以對(duì)平面坐標(biāo)和地理坐標(biāo)間關(guān)系進(jìn)行確定。在平面與曲面間創(chuàng)建點(diǎn)點(diǎn)之間的函數(shù)表達(dá)式的方法即被稱作是地圖投影。投影實(shí)際上就是把大地橢球面內(nèi)各經(jīng)緯度上的點(diǎn)坐標(biāo)，根據(jù)特定數(shù)學(xué)關(guān)系向平面對(duì)應(yīng)的坐標(biāo)上進(jìn)行轉(zhuǎn)換。當(dāng)前國(guó)內(nèi)采用的西安80、北京54坐標(biāo)體系都是利用Gauss-Kruger投影方式。我國(guó)地處地球北半部，為防止Y坐標(biāo)線產(chǎn)生負(fù)值，因此規(guī)定X軸西移500km，如此所有的坐標(biāo)均呈現(xiàn)為正值。

2　ArcGIS坐標(biāo)轉(zhuǎn)換方法

ArcGIS地理對(duì)照解決措施是把空間數(shù)據(jù)與坐標(biāo)系進(jìn)行分開(kāi)儲(chǔ)存，因此在轉(zhuǎn)換坐標(biāo)時(shí)，不但需要對(duì)地理參照系進(jìn)行定義，還應(yīng)了解空間坐標(biāo)數(shù)據(jù)的處理形式。通常而言，ArcGIS包括兩套坐標(biāo)體系，分別為投影坐標(biāo)系和地理坐標(biāo)系。后者是以經(jīng)緯度對(duì)橢球或球表面上的點(diǎn)位參照系進(jìn)行定義，前者是以三維、二維點(diǎn)線面等要素位置進(jìn)行（x，y，z）參照系的定位。同時(shí)，在ArcGIS中預(yù)裝了全球數(shù)百類地理參照系，這里國(guó)內(nèi)常會(huì)用到的包括北京54坐標(biāo)系、西安80坐標(biāo)系及2000大地坐標(biāo)系等。此外，在這基礎(chǔ)之上，ArcGIS也包括有動(dòng)態(tài)轉(zhuǎn)換、靜態(tài)轉(zhuǎn)換與即時(shí)轉(zhuǎn)換這三類轉(zhuǎn)換形式（如表2中所示）。

表2　ArcGIS3種坐標(biāo)轉(zhuǎn)換方式的比較

由上可得，ArcGIS的坐標(biāo)轉(zhuǎn)換方法，不論地理參照系的定義，或是地理空間坐標(biāo)的轉(zhuǎn)換形式，均存有較為顯著的區(qū)別。對(duì)于ArcGIS空間坐標(biāo)完成轉(zhuǎn)換的原理若無(wú)法做到深刻了解，則無(wú)法避免該坐標(biāo)轉(zhuǎn)換形式被錯(cuò)誤使用的可能性。

3　ArcGIS中坐標(biāo)轉(zhuǎn)換

在ArcGIS系統(tǒng)中，全部數(shù)據(jù)均是通過(guò)圖層方式進(jìn)行儲(chǔ)存的，所以在實(shí)施坐標(biāo)轉(zhuǎn)換過(guò)程中，也必須依據(jù)圖層形式，即按照每個(gè)圖層實(shí)施轉(zhuǎn)換。在轉(zhuǎn)換坐標(biāo)前，還應(yīng)該定義總體工程和單個(gè)圖層的坐標(biāo)體系，并且對(duì)坐標(biāo)投影也應(yīng)進(jìn)行定義。

3.1坐標(biāo)系的定義

在通過(guò)ArcGIS進(jìn)行地圖的編制時(shí)，通常會(huì)產(chǎn)生多層數(shù)據(jù)，包括點(diǎn)、面、線等相關(guān)數(shù)據(jù)。在定義時(shí)，先應(yīng)對(duì)總體工程的坐標(biāo)系統(tǒng)開(kāi)展，點(diǎn)開(kāi)Data Frame Property選項(xiàng)卡之中的CoordinateSystem項(xiàng)，接著定義該工程所在處的坐標(biāo)體系。ArcGIS中帶著很多已完成定義過(guò)程的坐標(biāo)體系，包括：西安80、北京54、WGS-84等，并且對(duì)于西安80、北京54坐標(biāo)系定義根據(jù)3°帶與6°帶展開(kāi)劃分，僅需按照工程所在處范圍實(shí)施選取。若是獨(dú)立的地方性坐標(biāo)體系，則需定義此坐標(biāo)系，選新建選項(xiàng)展開(kāi)設(shè)置。在實(shí)現(xiàn)以上坐標(biāo)系的定義后，可以觀察到圖形之中坐標(biāo)顯出單位：degree或meter，然而此定義只是定義了工程的坐標(biāo)系統(tǒng)，對(duì)于各圖層并未產(chǎn)生影響，這個(gè)結(jié)果可利用查看各圖層Source選項(xiàng)卡得知（如圖1中所示）。

圖1　圖層屬性

進(jìn)行任何圖層坐標(biāo)系的定義都可利用工具箱ArcToolbox來(lái)展開(kāi)。在DataManagementTools工具欄中，有個(gè)工具ProjectionsandTransformations，可選用此工具中DefineProjection選項(xiàng)以定義圖層坐標(biāo)系。

3.2轉(zhuǎn)換參數(shù)的建立

坐標(biāo)體系間轉(zhuǎn)換主要包含在不同空間直角坐標(biāo)體系間轉(zhuǎn)換，在地球坐標(biāo)體系和空間直角坐標(biāo)間轉(zhuǎn)換，以及在高斯平面坐標(biāo)體系與大地坐標(biāo)體系間轉(zhuǎn)換等。ArcGIS的坐標(biāo)轉(zhuǎn)換方法包括七參數(shù)法和三參數(shù)法等。在定義轉(zhuǎn)換參數(shù)前，應(yīng)對(duì)需要進(jìn)行區(qū)域轉(zhuǎn)換的參數(shù)展開(kāi)計(jì)算。

3.3幾類坐標(biāo)體系的轉(zhuǎn)換辦法

對(duì)于不同種類、大小和格式的數(shù)據(jù)，在ArcGIS之中所對(duì)應(yīng)的坐標(biāo)轉(zhuǎn)換方式也各不相同。在ArcGIS的自帶Toolbox內(nèi)對(duì)于不同的數(shù)據(jù)坐標(biāo)體系，提供出兩個(gè)不同的轉(zhuǎn)換工具，即ProjectRaster和Project。

1）矢量圖的坐標(biāo)轉(zhuǎn)換。ArcGIS對(duì)于shp格式數(shù)據(jù)和矢量數(shù)據(jù)類型所設(shè)置的坐標(biāo)轉(zhuǎn)換方法為通過(guò)Toolbox工具箱中Project選項(xiàng)進(jìn)行（轉(zhuǎn)換方式如圖2所示）。

圖2　坐標(biāo)轉(zhuǎn)換的流程圖

2）影像圖的坐標(biāo)轉(zhuǎn)換。在測(cè)繪時(shí)所用影像圖大多是正射影像圖（DOM）。此處既有整幅拼接數(shù)據(jù)，也存在分幅數(shù)據(jù)。在對(duì)DOM實(shí)施坐標(biāo)轉(zhuǎn)換時(shí)，通常在DOM建立中便已完成，如需在應(yīng)用DOM中開(kāi)展坐標(biāo)轉(zhuǎn)換，則對(duì)于不同的數(shù)據(jù)量大小能用兩類方式：DOM分幅，對(duì)于小數(shù)據(jù)量，可以選擇Toolbox中Project Raster來(lái)實(shí)施轉(zhuǎn)換，且在轉(zhuǎn)換時(shí)要和矢量坐標(biāo)轉(zhuǎn)換相一致。完成拼接后的DOM數(shù)據(jù)，一般數(shù)據(jù)量很大。此時(shí)，選擇ProjectRaster工具實(shí)施轉(zhuǎn)換會(huì)造成系統(tǒng)崩潰，因此，可以采取數(shù)據(jù)輸出形式以轉(zhuǎn)換坐標(biāo)系。如將西安80轉(zhuǎn)換成北京54，做法包括：先定義DOM的坐標(biāo)系為西安80，之后再對(duì)Transformation的轉(zhuǎn)換參數(shù)進(jìn)行定義，然后利用結(jié)果輸出便能得出所需的北京54坐標(biāo)系DOM數(shù)據(jù)信息。該法最大可處理的影像坐標(biāo)轉(zhuǎn)換數(shù)據(jù)能達(dá)到30G。

4　坐標(biāo)轉(zhuǎn)換精度評(píng)估

為了評(píng)估坐標(biāo)轉(zhuǎn)換的精度質(zhì)量，在某區(qū)域選擇15個(gè)控制點(diǎn)坐標(biāo)數(shù)據(jù)，已知它的精準(zhǔn)西安80坐標(biāo)信息，根據(jù)以上方法以對(duì)此點(diǎn)的圖層實(shí)施坐標(biāo)轉(zhuǎn)換。先采取七參數(shù)法來(lái)計(jì)算出該區(qū)域的轉(zhuǎn)換坐標(biāo)參數(shù)，獲得坐標(biāo)的3個(gè)轉(zhuǎn)動(dòng)分量、3個(gè)偏移量以及1個(gè)尺度因子這七參數(shù)，針對(duì)ArcGIS中的圖層坐標(biāo)體系設(shè)置其轉(zhuǎn)換形式。由西安80坐標(biāo)轉(zhuǎn)換成北京54坐標(biāo)，可得最終的數(shù)據(jù)精度（如表3中坐標(biāo)轉(zhuǎn)換數(shù)據(jù)），表內(nèi)數(shù)據(jù)單位是米。

表3　利用ArcGIS轉(zhuǎn)換坐標(biāo)數(shù)據(jù)的精度表

由上表可看到，通過(guò)ArcGIS軟件對(duì)坐標(biāo)系按照七參數(shù)法進(jìn)行轉(zhuǎn)換，其數(shù)據(jù)精確度可以達(dá)到地形圖的精確度需要。現(xiàn)分析引起ArcGIS轉(zhuǎn)換坐標(biāo)誤差的因素為：北京54坐標(biāo)系和西安80坐標(biāo)系都屬于參心坐標(biāo)系，但是在定義地理坐標(biāo)體系時(shí)也存在差異；ArcGIS地理坐標(biāo)轉(zhuǎn)換的原理存在差別，然而ArcGIS并未向外部公布該轉(zhuǎn)換的原理；因轉(zhuǎn)換空間坐標(biāo)參數(shù)有著時(shí)空上的制約，因此計(jì)算轉(zhuǎn)換參數(shù)時(shí)的誤差也會(huì)來(lái)源于轉(zhuǎn)換坐標(biāo)的數(shù)據(jù)誤差。

5　結(jié)語(yǔ)

ArcGIS具有不同的坐標(biāo)轉(zhuǎn)換形式，因此，對(duì)于控制點(diǎn)坐標(biāo)轉(zhuǎn)換質(zhì)量的影響會(huì)不相同，并且在相同地理坐標(biāo)轉(zhuǎn)換法下，由于坐標(biāo)系選取不同，其轉(zhuǎn)換的精度也將不同。文中基于具體的坐標(biāo)轉(zhuǎn)換項(xiàng)目，針對(duì)實(shí)現(xiàn)ArcGIS系統(tǒng)地理坐標(biāo)轉(zhuǎn)換方法展開(kāi)了分析與研究，并對(duì)轉(zhuǎn)換項(xiàng)目中所遇問(wèn)題開(kāi)展總結(jié)與探討，從而希望可以為相關(guān)項(xiàng)目技術(shù)人員提供一些經(jīng)驗(yàn)上的參考。此外，因空間數(shù)據(jù)坐標(biāo)在進(jìn)行轉(zhuǎn)換時(shí)比較復(fù)雜，所以在轉(zhuǎn)換坐標(biāo)中，方案完成還必須針對(duì)實(shí)際數(shù)據(jù)狀況做好完善與改進(jìn)措施。

［1］ 趙慧慧，葛瑩.ArcGIS坐標(biāo)轉(zhuǎn)換方法及其精度評(píng)估［J］.地理空間信息.2016（03）.

［2］ 李平，盧立.ArcGIS中幾種坐標(biāo)系轉(zhuǎn)換方法的應(yīng)用研究［J］.城市勘測(cè).2012（01）.

［3］ 周衛(wèi)，張彥彥，龍毅.圖形坐標(biāo)轉(zhuǎn)換方法與實(shí)現(xiàn)［J］.地球信息科學(xué).2007（02）.

P226.3