劉思遠(yuǎn) 李小勇 馬壯 王永生

（1.河南省遙感測繪院,河南 鄭州 450003；2.河南省科源測繪中心,河南 鄭州 450003；3. 河南省測繪地理信息局信息中心,河南 鄭州 450003）

2000國家大地坐標(biāo)系是我國自主建立、適應(yīng)現(xiàn)代空間技術(shù)發(fā)展趨勢的地心坐標(biāo)系，自2008年7月1日起啟用，到2018年全面完成國土資源系統(tǒng)各類空間數(shù)據(jù)向2000國家大地坐標(biāo)系轉(zhuǎn)換工作[1]，這關(guān)系到國土資源規(guī)劃、用地、地籍（不動產(chǎn)登記）、礦業(yè)權(quán)審批等工作的開展，有利于省、市、縣三級國土業(yè)務(wù)數(shù)據(jù)的共建共享，推動國土資源管理、空間規(guī)劃編制等工作有序開展。面對種類繁雜的國土資源空間數(shù)據(jù)和巨大的數(shù)據(jù)量，對地市級國土資源空間數(shù)據(jù)2000國家大地坐標(biāo)系轉(zhuǎn)換技術(shù)流程進(jìn)行研究，可為國土管理部門提供技術(shù)支撐。區(qū)工礦倉儲用地級別與基準(zhǔn)地價圖、土地利用現(xiàn)狀等數(shù)據(jù)，主要通過不同的改正量進(jìn)行坐標(biāo)轉(zhuǎn)換。

（3）1∶2000～1∶500大比例尺數(shù)據(jù)主要為農(nóng)村集體土地所有權(quán)（界址點(diǎn)精度不超過5cm）、農(nóng)村集體土地使用權(quán)、城鎮(zhèn)地籍、不動產(chǎn)登記等種類數(shù)據(jù)。該類數(shù)據(jù)主要為1980西安坐標(biāo)系，轉(zhuǎn)換前應(yīng)選取數(shù)據(jù)對應(yīng)的起算控制點(diǎn)進(jìn)行分析，開展必要的普查和聯(lián)測，統(tǒng)一轉(zhuǎn)換或歸算至省衛(wèi)星導(dǎo)航定位基準(zhǔn)站網(wǎng)2000國家大地坐標(biāo)系框架下[2]。利用起算控制點(diǎn)分析轉(zhuǎn)換或聯(lián)測平差結(jié)果，按照點(diǎn)位坐標(biāo)的轉(zhuǎn)換方法計算轉(zhuǎn)換參數(shù)，開展數(shù)據(jù)轉(zhuǎn)換。數(shù)據(jù)轉(zhuǎn)換技術(shù)流程如圖1所示。

圖1 數(shù)據(jù)轉(zhuǎn)換技術(shù)流程圖

1 總體思路

國土資源空間數(shù)據(jù)類型豐富，不同基準(zhǔn)、不同比例尺數(shù)據(jù)進(jìn)行坐標(biāo)轉(zhuǎn)換時，涉及的數(shù)據(jù)分析、轉(zhuǎn)換流程等不同。國土資源空間數(shù)據(jù)原始坐標(biāo)系主要包括1980西安坐標(biāo)系、1954北京坐標(biāo)系和城市獨(dú)立坐標(biāo)系。原始坐標(biāo)為1954北京坐標(biāo)系的數(shù)據(jù)，需先行轉(zhuǎn)換為1980西安坐標(biāo)系，再轉(zhuǎn)換為2000國家大地坐標(biāo)系。城市獨(dú)立坐標(biāo)系（主要為大比例尺數(shù)據(jù)），包括為建設(shè)管理各類工程的勘察設(shè)計、施工和管理、地形圖測繪、城鎮(zhèn)地籍測量等需要而建立的高精度城市坐標(biāo)系，這類坐標(biāo)系的數(shù)據(jù)可采用參數(shù)轉(zhuǎn)換方式進(jìn)行轉(zhuǎn)換。

根據(jù)國土資源空間數(shù)據(jù)的不同比例尺，可采用參數(shù)轉(zhuǎn)換或格網(wǎng)改正量進(jìn)行坐標(biāo)系轉(zhuǎn)換，并檢核轉(zhuǎn)換精度，確保數(shù)據(jù)轉(zhuǎn)換成果滿足要求。具體轉(zhuǎn)換方式如下：

（1）比例尺小于1∶250000的數(shù)據(jù)無需轉(zhuǎn)換。

（2）1∶100000～ 1∶25000小 比 例 尺 數(shù) 據(jù) 和1∶10000～1∶5000中比例尺數(shù)據(jù)主要為礦產(chǎn)資源總體規(guī)劃、永久基本農(nóng)田數(shù)據(jù)庫、土地利用總體規(guī)劃、城坐標(biāo)轉(zhuǎn)換需要注意：①轉(zhuǎn)換后成果的坐標(biāo)基準(zhǔn)需統(tǒng)一到衛(wèi)星導(dǎo)航定位基準(zhǔn)站網(wǎng)的2000國家大地坐標(biāo)系框架下；②轉(zhuǎn)換前后各類成果的精度、拓?fù)潢P(guān)系和屬性結(jié)構(gòu)不變，要素不丟漏；③分幅數(shù)據(jù)轉(zhuǎn)換后，重新生成方里格網(wǎng)，修改圖廓點(diǎn)及方里格網(wǎng)標(biāo)注坐標(biāo)，不進(jìn)行數(shù)據(jù)內(nèi)容的重新分幅裁切和整飾處理；④數(shù)據(jù)轉(zhuǎn)換后應(yīng)同步開展質(zhì)量檢查，并按照不同的比例尺類別統(tǒng)計轉(zhuǎn)換精度。

2 坐標(biāo)轉(zhuǎn)換技術(shù)流程

2.1 控制點(diǎn)選取與聯(lián)測

（1）國家標(biāo)準(zhǔn)坐標(biāo)系數(shù)據(jù)轉(zhuǎn)換重合點(diǎn)選取。從市級D級GPS控制網(wǎng)(該數(shù)據(jù)主要為1980西安坐標(biāo)系)均勻選擇覆蓋市域的若干控制點(diǎn)，通過GNSS連續(xù)運(yùn)行基準(zhǔn)站觀測，得到被選取控制點(diǎn)的2000國家大地坐標(biāo)系成果；1980西安坐標(biāo)系和2000國家大地坐標(biāo)系的控制點(diǎn)成果一起組成重合點(diǎn)，用于轉(zhuǎn)換參數(shù)解算和數(shù)據(jù)轉(zhuǎn)換。多余的控制點(diǎn)用于模型校驗(yàn)和精度評定。

（2）相對獨(dú)立的平面坐標(biāo)系數(shù)據(jù)轉(zhuǎn)換重合點(diǎn)選取。均勻選擇覆蓋地市的8～15個控制點(diǎn)，通過GNSS連續(xù)運(yùn)行基準(zhǔn)站觀測，得到被選取控制點(diǎn)的2000國家大地坐標(biāo)系成果；與原相對獨(dú)立的平面坐標(biāo)系成果組成重合點(diǎn)，用于轉(zhuǎn)換參數(shù)解算和數(shù)據(jù)的轉(zhuǎn)換。多余的控制點(diǎn)用于模型校驗(yàn)和精度評定。采用基于GNSS連續(xù)運(yùn)行基準(zhǔn)站的觀測模式?？刂泣c(diǎn)聯(lián)測流程如圖2所示。

圖2 控制點(diǎn)聯(lián)測流程圖

2.2 轉(zhuǎn)換參數(shù)模型

市級大地坐標(biāo)系數(shù)據(jù)成果的坐標(biāo)系統(tǒng)轉(zhuǎn)換選擇二維七參數(shù)轉(zhuǎn)換模型，市級平面直角坐標(biāo)系數(shù)據(jù)成果的坐標(biāo)系統(tǒng)轉(zhuǎn)換選擇“二維四參數(shù)＋換帶”模型。對于相對獨(dú)立的平面坐標(biāo)系統(tǒng)與2000國家大地坐標(biāo)系的聯(lián)系，可采用“二維四參數(shù)＋換帶”“逆過程＋二維四參數(shù)＋換帶”或二維多項(xiàng)式回歸模型。相同中央子午線、不同橢球和不同投影基準(zhǔn)面采用二維四參數(shù)轉(zhuǎn)換法；相同橢球、不同中央子午線采用換帶方法 ；不同橢球、不同投影基準(zhǔn)面、不同中央子午線采用“二維四參數(shù)＋換帶”方法。根據(jù)確定的重合點(diǎn)坐標(biāo)、坐標(biāo)系統(tǒng)轉(zhuǎn)換模型利用最小二乘法計算坐標(biāo)系統(tǒng)轉(zhuǎn)換參數(shù)。

2.3 數(shù)據(jù)坐標(biāo)轉(zhuǎn)換

1∶5萬、1∶1萬、1∶5000等中小比例尺數(shù)據(jù)，采用全國系列比例尺圖幅改正量或全國高精度高分辨率格網(wǎng)改正量（3″× 3″）進(jìn)行轉(zhuǎn)換。其中，1∶5萬比例尺數(shù)據(jù)的轉(zhuǎn)換采用全國1∶5萬比例尺圖幅改正量，1∶1萬比例尺數(shù)據(jù)的轉(zhuǎn)換采用全國1∶1萬比例尺圖幅改正量，1∶5000比例尺數(shù)據(jù)的轉(zhuǎn)換采用全國高精度高分辨率格網(wǎng)改正量（3″× 3″）?？蓪⒃擃悢?shù)據(jù)的改正量參數(shù)寫入定制的2000國家大地坐標(biāo)系轉(zhuǎn)換平臺中，通過轉(zhuǎn)換平臺進(jìn)行人機(jī)交互，自動進(jìn)行相關(guān)數(shù)據(jù)的坐標(biāo)轉(zhuǎn)換。

大比例尺矢量和柵格空間數(shù)據(jù)，利用兩個坐標(biāo)系下的重合點(diǎn)計算轉(zhuǎn)換參數(shù)，通過點(diǎn)位坐標(biāo)轉(zhuǎn)換法，運(yùn)用定制的2000國家大地坐標(biāo)系轉(zhuǎn)換平臺，實(shí)現(xiàn)向2000國家大地坐標(biāo)系或其相對獨(dú)立的平面坐標(biāo)系轉(zhuǎn)換。

3 精度評估與檢核

3.1 中小比例尺的精度評估與檢核

為有效進(jìn)行坐標(biāo)轉(zhuǎn)換精度評定，在轉(zhuǎn)換數(shù)據(jù)區(qū)域（以市域?yàn)閱挝唬?，將衛(wèi)星導(dǎo)航定位基準(zhǔn)站點(diǎn)、分析轉(zhuǎn)換或聯(lián)測平差后的數(shù)據(jù)起算控制點(diǎn)作為檢核點(diǎn)，采用數(shù)據(jù)對應(yīng)比例尺圖幅或格網(wǎng)改正量對其原坐標(biāo)系成果進(jìn)行2000國家大地坐標(biāo)系轉(zhuǎn)換，與其平差解算的2000國家大地坐標(biāo)系成果進(jìn)行對比統(tǒng)計，分別計算評定1∶5萬、1∶1萬、1∶5000比例尺數(shù)據(jù)的轉(zhuǎn)換精度。轉(zhuǎn)換點(diǎn)位的平均精度應(yīng)優(yōu)于相應(yīng)比例尺數(shù)據(jù)圖上的0.1mm。

（1）檢核點(diǎn)殘差

其中，X1、Y1為用圖幅或格網(wǎng)改正量轉(zhuǎn)換獲取的檢核點(diǎn)在2000國家大地坐標(biāo)系下的平面坐標(biāo)；X2、Y2為檢核點(diǎn)分析轉(zhuǎn)換或聯(lián)測平差后的2000國家大地坐標(biāo)系下的平面坐標(biāo)[3]。

（2）平面點(diǎn)位中誤差

Mp=，式中，平面坐標(biāo)x殘差中誤差，平面坐標(biāo)y殘差中誤差，n為檢核點(diǎn)個數(shù)。

3.2 大比例尺精度評估與檢核

選擇部分重合點(diǎn)作為外部檢核點(diǎn)，不參與坐標(biāo)系統(tǒng)轉(zhuǎn)換參數(shù)計算，用轉(zhuǎn)換參數(shù)計算其轉(zhuǎn)換后的坐標(biāo)，與已知坐標(biāo)進(jìn)行比對，統(tǒng)計轉(zhuǎn)換的外符合精度。區(qū)域內(nèi)應(yīng)選定至少6個均勻分布的重合點(diǎn)對坐標(biāo)系統(tǒng)轉(zhuǎn)換精度進(jìn)行檢核。轉(zhuǎn)換點(diǎn)位的平均精度應(yīng)優(yōu)于相應(yīng)比例尺數(shù)據(jù)圖上的0.1mm。

（1）重合點(diǎn)殘差

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

（2）平面點(diǎn)位中誤差

4 結(jié)論

開展國土資源空間數(shù)據(jù)的2000國家大地坐標(biāo)系轉(zhuǎn)換工作，統(tǒng)一城市空間信息成果坐標(biāo)系統(tǒng)和大地基準(zhǔn)，是有效保障空間規(guī)劃體系建設(shè)與監(jiān)督實(shí)施、自然資源開發(fā)利用與保護(hù)監(jiān)管的基礎(chǔ)，通過對國土資源空間數(shù)據(jù)2000國家大地坐標(biāo)系轉(zhuǎn)換的總體思路、技術(shù)流程、精度檢核等環(huán)節(jié)的研究，可為其他行業(yè)的2000國家大地坐標(biāo)系轉(zhuǎn)換工作提供解決思路。