天津市陸海坐標(biāo)系統(tǒng)一轉(zhuǎn)換研究

張奇，王文旭，張志全，王文江

(天津市測繪院，天津 西青 300381)

?

天津市陸海坐標(biāo)系統(tǒng)一轉(zhuǎn)換研究

張奇*，王文旭，張志全，王文江

(天津市測繪院，天津 西青 300381)

依托國家首批測繪地理信息公益性行業(yè)科研專項《天津市陸海一體化地理信息服務(wù)平臺建設(shè)研究》，收集了天津市陸域、海域各坐標(biāo)系及陸域、海域控制網(wǎng)資料，選取了坐標(biāo)轉(zhuǎn)換模型和坐標(biāo)重合點，進行了坐標(biāo)轉(zhuǎn)換參數(shù)計算，確定了天津市陸、海坐標(biāo)系統(tǒng)一轉(zhuǎn)換關(guān)系并評定了轉(zhuǎn)換精度，選取了檢測點對坐標(biāo)轉(zhuǎn)換精度進行了檢測，實現(xiàn)了天津市陸海坐標(biāo)系的統(tǒng)一轉(zhuǎn)換。

陸海一體化；平面基準(zhǔn)；坐標(biāo)轉(zhuǎn)換；轉(zhuǎn)換參數(shù)

1 研究背景

本研究依托國家首批測繪地理信息公益性行業(yè)科研專項《天津市陸海一體化地理信息服務(wù)平臺建設(shè)研究》項目開展，目的是提供陸海地理空間信息平面基準(zhǔn)轉(zhuǎn)換與統(tǒng)一的理論和關(guān)鍵技術(shù)支持，通過對陸、?，F(xiàn)用以及歷史成果資料所用坐標(biāo)系資料的搜集與調(diào)查、坐標(biāo)系間轉(zhuǎn)換關(guān)系的確定來實現(xiàn)陸、海多種平面坐標(biāo)系之間的相互轉(zhuǎn)換，為陸海地理信息數(shù)據(jù)融合提供支持[1]，目前國內(nèi)尚無天津地區(qū)陸海坐標(biāo)系統(tǒng)一轉(zhuǎn)換的研究。

2 研究內(nèi)容

天津市陸域現(xiàn)用坐標(biāo)系為2000國家大地坐標(biāo)系[2]以及1990年天津市任意直角坐標(biāo)系，陸域歷史成果用到1980西安坐標(biāo)系；天津市海域現(xiàn)用的坐標(biāo)系為2000國家大地坐標(biāo)系，歷史成果用到1980西安坐標(biāo)系。由于陸域、海域現(xiàn)用的通用坐標(biāo)系都是2000國家大地坐標(biāo)系，因此本研究將重點放在涉及2000國家大地坐標(biāo)系的坐標(biāo)轉(zhuǎn)換，包括2000國家大地坐標(biāo)系與1980西安大地坐標(biāo)系的轉(zhuǎn)換，2000國家大地坐標(biāo)系與1990年天津市任意直角坐標(biāo)系的轉(zhuǎn)換。本研究的主要內(nèi)容包括收集天津市陸域、海域各坐標(biāo)系資料及控制網(wǎng)資料，選取坐標(biāo)轉(zhuǎn)換模型和坐標(biāo)重合點進行坐標(biāo)轉(zhuǎn)換參數(shù)計算，確定各坐標(biāo)系之間的坐標(biāo)轉(zhuǎn)換關(guān)系，評定坐標(biāo)轉(zhuǎn)換精度，選取檢測點對坐標(biāo)轉(zhuǎn)換精度進行檢測[3]。

3 天津市陸、海坐標(biāo)系統(tǒng)一轉(zhuǎn)換

3.1 陸、?？刂凭W(wǎng)

天津陸、?？刂凭W(wǎng)統(tǒng)計如表1所示，分布如圖1所示，陸域共有各等級控制點248個，海域共有各等級控制點34個。

陸、海平面控制點統(tǒng)計 表1

3.2 轉(zhuǎn)換模型

常用的坐標(biāo)轉(zhuǎn)換模型包括Bursa七參數(shù)轉(zhuǎn)換模型[4]以及二維高斯平面坐標(biāo)轉(zhuǎn)換模型[5]。本研究選用理論上更加嚴密的Bursa七參數(shù)轉(zhuǎn)換模型進行陸海坐標(biāo)系統(tǒng)一轉(zhuǎn)換。

(1)

其中，三個平移參數(shù)[dxdydz]T，三個旋轉(zhuǎn)參數(shù)[ωXωYωZ]T和一個尺度參數(shù)m。

圖1 Bursa七參數(shù)轉(zhuǎn)換示意圖

3.3 重合點及檢測點選取

重合點選取的原則是盡量選取數(shù)量足夠的高等級、高精度且分布均勻的點作為坐標(biāo)轉(zhuǎn)換參數(shù)計算的重合點。在2000國家大地坐標(biāo)系與1990年天津市任意直角坐標(biāo)系的轉(zhuǎn)換中，經(jīng)過試算，均勻選取了48個重合點進行坐標(biāo)轉(zhuǎn)換參數(shù)的計算。在2000國家大地坐標(biāo)系與1980西安坐標(biāo)系轉(zhuǎn)換參數(shù)計算中，除選取前述的48個陸域重合點外，還選取了19個海域坐標(biāo)重合點，包括北方海區(qū)GPS控制網(wǎng)B級點6個以及C級點13個。

坐標(biāo)轉(zhuǎn)換參數(shù)計算完成、坐標(biāo)轉(zhuǎn)換關(guān)系確定后，須對坐標(biāo)轉(zhuǎn)換精度進行檢測，本研究在天津陸、海域一共選取均勻分布的65個重合點進行坐標(biāo)轉(zhuǎn)換精度驗證。

圖2 坐標(biāo)轉(zhuǎn)換重合點及檢測點分布

檢測點分布如圖2右中綠色圓點所示，65個檢驗點包括陸域50個(37個GPS C級點和13個連續(xù)運行基準(zhǔn)站點)重合點，海域15個重合點(天津港附近GPS D級網(wǎng)控制點)。

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

雖然陸域、海域控制點都有2000國家大地坐標(biāo)系及1980西安坐標(biāo)系成果，但因陸域、海域控制網(wǎng)施測時間不同，且未經(jīng)過聯(lián)合平差，其成果之間由于地表位移以及平差基準(zhǔn)不完全相同而可能存在偏差。因此，必須聯(lián)合陸域、海域控制點進行坐標(biāo)轉(zhuǎn)換參數(shù)計算。

利用前述重合點和坐標(biāo)轉(zhuǎn)換模型，采用最小二乘法[6]計算出2000國家大地坐標(biāo)系與1980西安坐標(biāo)系、1990年天津市任意直角坐標(biāo)系之間的轉(zhuǎn)換參數(shù)。具體計算過程如下，首先，將重合點坐標(biāo)換算為空間直角坐標(biāo)；其次，使用全部的重合點求取轉(zhuǎn)換參數(shù)，再利用該轉(zhuǎn)換參數(shù)回代，求得已知點的坐標(biāo)，計算已知點的殘差；再次，分析殘差值大于 0.1 m點的可靠性，檢查點位和成果的正確性，重新確定重合點；最后，重復(fù)上述步驟，直到回代殘差到一定范圍內(nèi)為止。計算得到各重合點轉(zhuǎn)換殘差(平面)如圖3所示。從圖3可以看出，2000國家大地坐標(biāo)系到1980西安坐標(biāo)系轉(zhuǎn)換中，在天津市東部沿海以及西部市界處，陸域重合點殘差較大，個別點如唐家河以及西北斜殘差(平面)達到了 0.3 m。2000國家大地坐標(biāo)系到1990年天津市任意直角坐標(biāo)系轉(zhuǎn)換中，重合點轉(zhuǎn)換殘差(平面)小于 0.11 m，在武清西部、漢沽、大港重合點坐標(biāo)轉(zhuǎn)換殘差較大。

表2統(tǒng)計了重合點轉(zhuǎn)換殘差絕對值在x方向、y方向、平面上的最小值、最大值、平均值以及中誤差。從表中可以看出，2000國家大地坐標(biāo)系到1980年西安坐標(biāo)系轉(zhuǎn)換中，重合點坐標(biāo)轉(zhuǎn)換殘差在x方向上為 0.104 m，y方向為 0.102 m，平面上為 0.146 m。2000國家大地坐標(biāo)系到1990年天津市任意直角坐標(biāo)系轉(zhuǎn)換中，重合點坐標(biāo)轉(zhuǎn)換殘差在x方向上為 0.042 m，y方向為 0.043 m，平面上為 0.060 m。

圖3 坐標(biāo)系轉(zhuǎn)換殘差 重合點殘差絕對值統(tǒng)計表 表2

圖4 重合點2000國家大地坐標(biāo)系與1980西安坐標(biāo)系轉(zhuǎn)換殘差分布

各重合點殘差分布如圖4所示。從圖4可以看出，2000國家大地坐標(biāo)系到1980年西安坐標(biāo)系轉(zhuǎn)換中重合點殘差x方向和y方向絕大多數(shù)在 ±0.1 m之間。2000國家大地坐標(biāo)系到1990年天津市任意直角坐標(biāo)系轉(zhuǎn)換中重合點殘差x方向上在 ±5 cm以內(nèi)，在y方向殘差在在 ±0.11 m之內(nèi)，在平面上小于 0.11 m。

參與2000國家大地坐標(biāo)系與1980西安坐標(biāo)系間轉(zhuǎn)換參數(shù)計算的國家一、二等三角點大多數(shù)于上世紀50年代～70年代布設(shè)施測，由于我國地殼運動變化復(fù)雜，經(jīng)過多年，這些三角點點位發(fā)生了不同程度的變化。為了使坐標(biāo)轉(zhuǎn)換結(jié)果和1980西安坐標(biāo)系下的歷史成果進一步符合以達到項目設(shè)計書中坐標(biāo)轉(zhuǎn)換精度優(yōu)于 10 cm的要求，同時兼顧不對坐標(biāo)系的均勻性產(chǎn)生過大的影響，本研究依據(jù)最小二乘原理，采用二次多項式對2000國家大地坐標(biāo)系到1980西安坐標(biāo)系轉(zhuǎn)換參數(shù)計算中的重合點的殘差進行了擬合，采用擬合得到的二次曲面進行殘差分配，圖5為x方向和y方向的殘差擬合曲面。

圖5 2000國家大地坐標(biāo)系到1980年西安坐標(biāo)系轉(zhuǎn)換殘差擬合

采用二次曲面進行殘差分配后，得到2000國家大地坐標(biāo)系到1980年西安坐標(biāo)系轉(zhuǎn)換殘差中誤差在x方向和y方向上分別為 0.076 cm和 0.073 cm，大大優(yōu)于不進行殘差擬合的結(jié)果。

3.5 精度推算

坐標(biāo)轉(zhuǎn)換關(guān)系確定后，需要對坐標(biāo)轉(zhuǎn)換精度進行推算[7]。

(1)V(殘差)=重合點已知坐標(biāo)-重合點轉(zhuǎn)換坐標(biāo)

(2)殘差中誤差

坐標(biāo)轉(zhuǎn)換精度是通過求取轉(zhuǎn)換參數(shù)的重合點的殘差中誤差體現(xiàn)的。2000國家大地坐標(biāo)系坐標(biāo)與1980西安坐標(biāo)、1990年天津市任意直角坐標(biāo)系坐標(biāo)轉(zhuǎn)換精度統(tǒng)計如表3所示。

坐標(biāo)轉(zhuǎn)換殘差中誤差統(tǒng)計表 表3

從表3可以看出，2000國家大地坐標(biāo)系與1980西安坐標(biāo)系轉(zhuǎn)換精度較差，x方向和y方向上分別為 0.071 4m和 0.068 8m，2000國家大地坐標(biāo)系到1990年天津市任意直角坐標(biāo)系坐標(biāo)轉(zhuǎn)換的精度在x方向和y方向分別為 0.041 6m和 0.042 6m，優(yōu)于2000國家大地坐標(biāo)系與1980西安坐標(biāo)系轉(zhuǎn)換精度，原因是參與其轉(zhuǎn)換參數(shù)計算的重合點于2005年進行過復(fù)測，成果現(xiàn)勢性較好。

3.6 外符合精度評估

確定坐標(biāo)轉(zhuǎn)換關(guān)系后，進行了坐標(biāo)轉(zhuǎn)換精度檢測。以檢測點通過坐標(biāo)轉(zhuǎn)換得到的坐標(biāo)與其原始坐標(biāo)的較差來檢測坐標(biāo)轉(zhuǎn)換精度。精度檢測點轉(zhuǎn)換后坐標(biāo)與原始坐標(biāo)較差分布如圖6所示，統(tǒng)計表如表4所示。

圖6 精度檢測點轉(zhuǎn)換經(jīng)轉(zhuǎn)換得到坐標(biāo)與原始坐標(biāo)較差分布

檢測點偏差統(tǒng)計表 表4

由表4及圖6可以看出，在2000國家大地坐標(biāo)系到1980年西安坐標(biāo)系轉(zhuǎn)換中，檢測點偏差在x方向、y方向、平面中誤差分別為 0.061m、 0.022m、和 0.065m，在2000國家大地坐標(biāo)系到1990年天津市任意直角坐標(biāo)系轉(zhuǎn)換中，檢測點偏差在x方向、y方向、平面中誤差分別為 0.005m、 0.003m、和 0.005m。

聯(lián)合坐標(biāo)轉(zhuǎn)換重合點殘差與精度檢測點偏差[8]，可以得到2000國家大地坐標(biāo)系與1980西安坐標(biāo)系、1990年天津市任意直角坐標(biāo)系之間轉(zhuǎn)換精度空間分布如圖7所示。

從轉(zhuǎn)換精度分布圖中可以看出，2000國家大地坐標(biāo)系到1980西安坐標(biāo)系轉(zhuǎn)換精度在天津西部市界(武清西部)以及天津東南沿海(包括漢沽東部、塘沽、大港)低于市內(nèi)其他區(qū)域，在天津南部近海區(qū)域轉(zhuǎn)換精度低于東部近海。2000國家大地坐標(biāo)系到1990年天津市任意直角坐標(biāo)系轉(zhuǎn)換精度在市界處(薊縣東部、寧河?xùn)|部、漢沽東部、大港西南部、武清西南部與靜海西北部等區(qū)域)低于市內(nèi)區(qū)域。

4 結(jié) 論

本文收集了天津市陸、海各坐標(biāo)系資料及陸、?？刂凭W(wǎng)資料，選取了坐標(biāo)轉(zhuǎn)換模型和坐標(biāo)重合點，進行了坐標(biāo)轉(zhuǎn)換參數(shù)計算，確定了各坐標(biāo)系之間的坐標(biāo)轉(zhuǎn)換關(guān)系，評定了坐標(biāo)轉(zhuǎn)換精度，并選取檢測點對坐標(biāo)轉(zhuǎn)換進行了精度檢測。研究結(jié)果表明，采用本文所確定的坐標(biāo)轉(zhuǎn)換關(guān)系進行2000國家大地坐標(biāo)系到1980年西安坐標(biāo)坐標(biāo)系的轉(zhuǎn)換，x方向和y方向中誤差分別為 0.071 4m和 0.068 8m，2000國家大地坐標(biāo)系到坐標(biāo)1990年天津市任意直角坐標(biāo)系坐標(biāo)轉(zhuǎn)換，x方向和y方向中誤差分別為 0.041 6m和 0.042 6m。

坐標(biāo)轉(zhuǎn)換除本文所研究的正向轉(zhuǎn)換外，還包括逆向轉(zhuǎn)換，即1980西安坐標(biāo)系、1990年天津市任意直角坐標(biāo)系到2000國家大地坐標(biāo)系的轉(zhuǎn)換。反向轉(zhuǎn)換與正向轉(zhuǎn)換采用的轉(zhuǎn)換模型相同，重合點相同，轉(zhuǎn)換精度及其分布特性理論上也相同，本文不予贅述。

[1] 殷曉冬，胡家升，張立華等. 海岸帶多源數(shù)據(jù)三維無縫拼接技術(shù)[J].大連海事大學(xué)學(xué)報，2008(2)：19～23.

[2] 楊華忠，李軍，汪舟平等. 我國常用大地坐標(biāo)系與2000 中國大地坐標(biāo)系間的轉(zhuǎn)換[J]. 中國測繪學(xué)會大地測量專業(yè)委員會2007年綜合性學(xué)術(shù)年會，2007(11)：1～10.

[3] 焦立芬. 基于坐標(biāo)轉(zhuǎn)換重合點的分布、密度、精度與轉(zhuǎn)換精度分析[J]. 測繪技術(shù)裝備，2013(4)：25～28.

[4] 陳宇，白征東，羅滕. 基于改進的布爾沙模型的坐標(biāo)轉(zhuǎn)換方法[J]. 大地測量與地球動力學(xué)，2010(03)：71～73.

[5] 陳義，沈云中，劉大杰. 適用于大旋轉(zhuǎn)角的三維基準(zhǔn)轉(zhuǎn)換的一種簡便模型[J]. 武漢大學(xué)學(xué)報信息科學(xué)版，2004(12)：1101～1105.

[6] 邱衛(wèi)寧，陶本藻，姚宜斌. 測量數(shù)據(jù)處理理論與方法[M]. 武漢：武漢大學(xué)出版社，2008.

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

[8] 鐘業(yè)勛，童新華，王龍波. 從1980西安坐標(biāo)系到2000國家大地坐標(biāo)系的坐標(biāo)變換[J]. 海洋測繪，2010(30)：1～3.

Study on the Transformation of Coordinate System Integrating Land and Sea in Tianjin

Zhang Qi，Wang Wenxu，Zhang Zhiquan，Wang Wenjiang

(Tianjin Institute of Surveying and Mapping，Tianjin 300381，China)

This study is supported by the first batch of the national public sector special research of surveying and mapping geographic information "The study on the construction of the platform of geographic information service integrating sea and land in Tianjin”. In this paper，firstly the information of the coordinate system and the control network of sea and land of Tianjin are collected，then the coordinate transformation model and coordinate coincidence points are selected，thirdly the coordinate transformation parameters are calculated and the transformation relationship of coordinate system of Tianjin land and sea is determined，then the accuracy of the transformation is evaluated and detected，finally the unity of Tianjin land and sea coordinate system transformation has been realized.

integration of land and sea；datum of plane；coordinate transformation；transformation parameters

1672-8262(2016)05-91-05

P226.3

A

2016—04—18

張奇(1988—)，男，碩士，工程師，主要從事天津市地理信息空間基準(zhǔn)維護方面的工作。

國家測繪地理信息局2014年測繪地理信息公益性科研專項(201412012)