譚文專， 熊忠招，賈化明， 邵澤興

（1.湖北省國(guó)土測(cè)繪院，湖北 武漢 430010; 2.義烏市國(guó)土資源局，浙江 義烏 322000）

農(nóng)村地籍調(diào)查中長(zhǎng)度投影變形超限研究

譚文專1， 熊忠招1，賈化明2， 邵澤興1

（1.湖北省國(guó)土測(cè)繪院，湖北 武漢 430010; 2.義烏市國(guó)土資源局，浙江 義烏 322000）

利用坐標(biāo)投影邊長(zhǎng)歸化關(guān)系，結(jié)合Excel表格，采用數(shù)學(xué)轉(zhuǎn)換關(guān)系解決投影長(zhǎng)度變形值超限的問題，確保確定的坐標(biāo)系滿足地籍調(diào)查成果和實(shí)現(xiàn)調(diào)查成果在統(tǒng)一坐標(biāo)系下不同比例尺數(shù)據(jù)及圖件融合的精度要求，滿足生產(chǎn)需求，提高作業(yè)效率。

地籍調(diào)查；投影變形；坐標(biāo)系；GPS-RTK測(cè)量

《湖北省宅基地使用權(quán)和集體建設(shè)用地使用權(quán)確權(quán)登記發(fā)證工作實(shí)施細(xì)則》中對(duì)坐標(biāo)系統(tǒng)的要求是：平面坐標(biāo)系統(tǒng)采用“1980西安坐標(biāo)系統(tǒng)”，高程系統(tǒng)采用“1985國(guó)家高程基準(zhǔn)”；中央子午線按照地圖投影3°帶的標(biāo)準(zhǔn)方法選定；當(dāng)長(zhǎng)度變形值大于±2．5 cm/km時(shí)，界址點(diǎn)測(cè)量和地籍圖測(cè)繪的投影面為各城鎮(zhèn)村莊坐標(biāo)起算點(diǎn)高程值的高程面，但不進(jìn)行高程抵償計(jì)算。提出利用坐標(biāo)投影邊長(zhǎng)歸化關(guān)系，結(jié)合EXCEL表格編輯公式程序設(shè)計(jì)，采用數(shù)學(xué)轉(zhuǎn)換關(guān)系解決投影長(zhǎng)度變形值超限問題，確保確定的獨(dú)立坐標(biāo)系滿足地籍調(diào)查成果和實(shí)現(xiàn)調(diào)查成果在統(tǒng)一坐標(biāo)系下不同比例尺數(shù)據(jù)及圖件融合的精度要求。

1 基本思路

按《湖北省宅基地使用權(quán)和集體建設(shè)用地使用權(quán)確權(quán)登記發(fā)證工作實(shí)施細(xì)則》中對(duì)“當(dāng)長(zhǎng)度變形值大于±2．5 cm/km時(shí)，界址點(diǎn)測(cè)量和地籍圖測(cè)繪的投影面為各城鎮(zhèn)村莊坐標(biāo)起算點(diǎn)高程值的高程面，但不進(jìn)行高程抵償計(jì)算”的要求，提出以下思路。

1.1 歸化轉(zhuǎn)換獨(dú)立坐標(biāo)系的范圍與坐標(biāo)起算點(diǎn)及高程面的確定

利用現(xiàn)有的1∶1萬地形圖和已經(jīng)完成的第二次土地調(diào)查和農(nóng)村集體土地所有權(quán)確權(quán)登記發(fā)證工作成果資料，計(jì)算投影長(zhǎng)度變形值和調(diào)查區(qū)的變形情況，來確定進(jìn)行高斯投影長(zhǎng)度變形歸化和參考橢球體邊長(zhǎng)歸化到相對(duì)于坐標(biāo)起算點(diǎn)高程值的高程面平距獨(dú)立坐標(biāo)系的范圍。范圍依精度要求而定，通常以離起算點(diǎn)最遠(yuǎn)點(diǎn)（最長(zhǎng)邊）歸化后的坐標(biāo)與歸化前的坐標(biāo)差值△S小于工作底圖上的0．1 mm，來保證不同比例尺的圖件和數(shù)據(jù)融合的精度；起算點(diǎn)以范圍中心點(diǎn)和范圍內(nèi)平均高程來確定坐標(biāo)值，實(shí)地?zé)o需設(shè)點(diǎn)，只需相應(yīng)點(diǎn)坐標(biāo)統(tǒng)一向起算點(diǎn)進(jìn)行歸化計(jì)算即可。采用歸化后高程面的平距邊長(zhǎng)，方位角不變，計(jì)算相應(yīng)點(diǎn)坐標(biāo)。

1）當(dāng)調(diào)查區(qū)范圍內(nèi)高差變化?。ㄩL(zhǎng)度變形值不大于±2．5 cm/km）時(shí)，一般以行政村（地籍子區(qū)）為單位確定歸化轉(zhuǎn)換獨(dú)立坐標(biāo)系的范圍，便于數(shù)據(jù)匯總統(tǒng)計(jì)等工作。

2）山區(qū)高差變化大（長(zhǎng)度變形值大于±2．5 cm/ km）時(shí)，以村莊為單位確定歸化轉(zhuǎn)換獨(dú)立坐標(biāo)系范圍、再劃小歸化轉(zhuǎn)換獨(dú)立坐標(biāo)系范圍。

3）歸化轉(zhuǎn)換獨(dú)立坐標(biāo)系范圍的確定，以符合不同比例尺的圖件和數(shù)據(jù)融合精度為準(zhǔn)，但要注意一個(gè)宗地必須在一個(gè)歸化轉(zhuǎn)換獨(dú)立坐標(biāo)系區(qū)內(nèi)，確保宗地內(nèi)界址點(diǎn)間的精度。

1.2 歸化計(jì)算的作業(yè)方法

歸化轉(zhuǎn)換獨(dú)立坐標(biāo)系計(jì)算利用EXCEL編制程序進(jìn)行，起算點(diǎn)是設(shè)計(jì)方案中預(yù)定的，所以程序編制后，只要把需歸化數(shù)據(jù)導(dǎo)入，運(yùn)用EXCEL進(jìn)行計(jì)算和檢核。數(shù)據(jù)歸化形式根據(jù)不同的作業(yè)方法和資源配置而定。

1）如果以行政村（地籍子區(qū)）為單位確定歸化轉(zhuǎn)換獨(dú)立坐標(biāo)系范圍時(shí)，可首先在全調(diào)查區(qū)布設(shè)首級(jí)控制（一級(jí)點(diǎn)），一般利用HBCORS-RTK測(cè)量方法按規(guī)范要求施測(cè)統(tǒng)一坐標(biāo)系成果， 再按確定歸化轉(zhuǎn)換獨(dú)立坐標(biāo)系范圍的控制點(diǎn)（確保控制點(diǎn)分布及控制范圍）計(jì)算獨(dú)立坐標(biāo)，然后在此基礎(chǔ)上采用RTK（RTK測(cè)量時(shí)應(yīng)直接利用首級(jí)控制原測(cè)量成果與獨(dú)立坐標(biāo)數(shù)據(jù)求其轉(zhuǎn)換參數(shù)）或?qū)Ь€法加密圖根點(diǎn)，進(jìn)行界址點(diǎn)測(cè)量和地籍圖測(cè)繪。

2）圖根控制（包括RTK測(cè)量方法測(cè)定的界址點(diǎn)）均施測(cè)統(tǒng)一坐標(biāo)系成果，再按所需范圍對(duì)相關(guān)點(diǎn)計(jì)算獨(dú)立坐標(biāo)，進(jìn)行界址點(diǎn)測(cè)量和地籍圖測(cè)繪。

2 坐標(biāo)轉(zhuǎn)換數(shù)學(xué)模型建立

1）地面兩點(diǎn)之間的水平長(zhǎng)度歸算到參考橢球面所產(chǎn)生變形的近似值，用式（1）計(jì)算：

式中，Hm為范圍內(nèi)的平均高程；R為平均曲率半徑；S0為2點(diǎn)投影到參考橢球面上的弦長(zhǎng)。

2）參考橢球面上的長(zhǎng)度投影到高斯平面上所產(chǎn)生變形的近似值，用式（2）計(jì)算：

式中，Ym為2點(diǎn)橫坐標(biāo)(自然值)的平均值；R為平均曲率半徑；S為2點(diǎn)(長(zhǎng)度)歸算到參考橢球面上的長(zhǎng)度。

3）地面測(cè)量邊長(zhǎng)水平距離的高程改化和投影改化的近似改正數(shù)之和的計(jì)算式為：

地面上有2點(diǎn)A、B，它們?cè)诟咚雇队捌矫嫔系闹苯亲鴺?biāo)分別為A（XA，YA）、B（XB、YB），則可由式（4）計(jì)算出AB間的距離S：

式中，S表示在高斯投影平面上2點(diǎn)間的距離。4）由式（1）～（4）計(jì)算可得：

式中，SG表示在高斯投影平面上兩點(diǎn)間的距離,即GPS距離；ST表示參考橢球面上的長(zhǎng)度。

式（5）是參考橢球面距離與高斯投影面平面距離的關(guān)系。利用方位角不改變?cè)恚梢郧蟪鱿鄬?duì)于各城鎮(zhèn)村莊坐標(biāo)起算點(diǎn)高程值（范圍內(nèi)的平均高程）的高程面其他點(diǎn)的坐標(biāo)值。

3 坐標(biāo)數(shù)據(jù)批量處理

根據(jù)上面推導(dǎo)的公式，利用ECCEL的編程功能進(jìn)行設(shè)計(jì)坐標(biāo)的批量處理（見圖1）。

1）表中第5行，C5、D5、E5單元格內(nèi)為起算點(diǎn)坐標(biāo)和高程，即以范圍中心點(diǎn)和范圍內(nèi)平均高程來確定坐標(biāo)值；第2行，B2為平均曲率半徑。應(yīng)注意與程序代碼相對(duì)應(yīng)，如果在復(fù)制、粘貼中單元格不對(duì)應(yīng)，應(yīng)注意修正。

2）表中第6行，C6、D6、E6單元格內(nèi)為需歸化轉(zhuǎn)換的點(diǎn)坐標(biāo)和高程；F6單元格為經(jīng)程序計(jì)算得出的歸化轉(zhuǎn)換點(diǎn)與起算點(diǎn)的邊長(zhǎng)；G6單元格為經(jīng)程序計(jì)算得出的歸化轉(zhuǎn)換后的點(diǎn)X坐標(biāo)；H6單元格為經(jīng)程序計(jì)算得出的歸化轉(zhuǎn)換后的點(diǎn)Y坐標(biāo)；I6單元格內(nèi)原高程直接轉(zhuǎn)入；J6單元格內(nèi)經(jīng)程序計(jì)算得出的歸化轉(zhuǎn)換后的坐標(biāo)與歸化轉(zhuǎn)換前坐標(biāo)差值△S，也是歸化轉(zhuǎn)換點(diǎn)相對(duì)于起算點(diǎn)長(zhǎng)度變形值。以下各項(xiàng)相同。

3）相對(duì)應(yīng)的由程序公式計(jì)算的單元格內(nèi)程序代碼如下：

F6單元格為：

=SQRT(SUMSQ(C6-$C$5)+SUMSQ(D6-$D$5))

G6單元格為：

=IF(AND(C6-$C$5>0,D6-$D$5>0),($C$5+F6/ (1+((D6+$D$5)/2-500000)^2/2/$B$2^2-$E$5/$B$2)* COS(ATAN((D6-$D$5)/(C6-$C$5)))), IF(AND(C6-$C$5＜0,D6-$D$5>0),$C$5+F6/(1+((D6+$D$5)/2-500000)^2/2/$B$2^2-$E$5/$B$2)*COS(ATAN((D6-$D$5)/(C6-$C$5))+PI()),IF(AND(C6-$C$5＜0,D6-$D$5＜0),$C$5+F 6/(1+((D 6+$D$5)/2-500000)^2/2/$B$2^2-$E$5/$B$2)*COS(ATAN((D6-$D$5)/(C6-$C$5))+PI()),$C$5+F6/(1+((D6+$D$5)/2-500000)^2/2/$B$2^2-$E$5/$B$2)*COS(ATAN((D6-$D$5)/(C6-$C$5))+2*PI()))))

H6單元格為：

=IF(AND(C6-$C$5>0,D6-$D$5>0),($D$5+F6/ (1+((D 6+$D$5)/2-5 0 0 0 0 0)^2/2/$B$2^2-$E$5/$B$2)*S I N(ATA N((D 6-$D$5)/(C 6-$C$5)))),IF(AND(C6-$C$5＜0,D6-$D$5>0),$D$5+F6/ (1+((D 6+$D$5)/2-5 0 0 0 0 0)^2/2/$B$2^2-$E$5/$B$2)*S I N(ATA N((D 6-$D$5)/(C 6-$C$5))+PI()),IF(AND(C6-$C$5＜0,D6-$D$5＜0),$D$5+F6/

圖1 村莊控制點(diǎn)歸算至起算點(diǎn)高程面計(jì)算表

(1+((D 6+$D$5)/2-5 0 0 0 0 0)^2/2/$B$2^2-$E$5/$B$2)*S I N(ATA N((D 6-$D$5)/(C 6-$C$5))+PI()),$D$5+F6/(1+((D6+$D$5)/2-500000)^2/2/$B$2^2-$E$5/$B$2)*SIN(ATAN((D6-$D$5)/ (C6-$C$5))+2*PI()))))

J6單元格為：

=SQRT((C6-G6)*(C6-G6)+(D6-H6)*(D6-H6))

程序代碼編制方式很多，以上程序代碼供參考。

4 結(jié) 語

農(nóng)村地籍調(diào)查中長(zhǎng)度投影變形超限時(shí)坐標(biāo)系確定方法很多，本文通過介紹農(nóng)村宅基地使用權(quán)和集體建設(shè)用地使用權(quán)確權(quán)登記發(fā)證工作試點(diǎn)中的應(yīng)用，總結(jié)出此方法的優(yōu)勢(shì)為：

1）便于作業(yè)員掌握，特別適用于CORS系統(tǒng)網(wǎng)絡(luò)RTK測(cè)量技術(shù)應(yīng)用和農(nóng)村散列式居民點(diǎn)的作業(yè)，可有效地提高作業(yè)效率。

2）可確保確定的獨(dú)立坐標(biāo)系滿足地籍調(diào)查成果需要和實(shí)現(xiàn)調(diào)查成果在統(tǒng)一坐標(biāo)系下不同比例尺數(shù)據(jù)及圖件融合的精度要求，形成“一張圖”。

3）便于成果使用，2套成果之間轉(zhuǎn)換方便，起算點(diǎn)實(shí)地不必設(shè)點(diǎn)，應(yīng)用時(shí)在統(tǒng)一坐標(biāo)系下作業(yè)后進(jìn)行簡(jiǎn)便計(jì)算即可得出需用成果。在實(shí)際工作中，村莊邊控制點(diǎn)難以保存，成果檢查驗(yàn)收時(shí)直接利用CORS系統(tǒng)網(wǎng)絡(luò)RTK測(cè)量技術(shù)獲得成果，直接進(jìn)行比對(duì)檢查，提高成果的正確性。

[1] 孫祥元，梅是義． 控制測(cè)量學(xué)[M]． 武漢：武漢大學(xué)出版社，2003

[2] TD/T 1001—2012．地籍調(diào)查規(guī)程[S]．

[3] CH/T 2009—2010．全球定位系統(tǒng)實(shí)時(shí)動(dòng)態(tài)（RTK）測(cè)量技術(shù)規(guī)范[S]．

[4] CJJ/T 8—2011．城市測(cè)量規(guī)范[S]．

[5] 許婭婭，黃文元．山區(qū)公路測(cè)量坐標(biāo)系的選擇方法研究[J]．測(cè)繪通報(bào)，2008(7)：26-28

[6] 王剛，童凌飛，關(guān)穎．CORS系統(tǒng)的主要技術(shù)及其應(yīng)用[J]．地理空間信息，2012，10(5)：96-98

[7] 潘國(guó)勝．工程控制網(wǎng)的投影變形分析與計(jì)算[J]．大眾科技，2009(6)：101-105

P282

B

1672-4623（2014）05-0068-02

10.3969/j.issn.1672-4623.2014.05.024

譚文專，高級(jí)工程師，主要從事測(cè)繪生產(chǎn)管理與技術(shù)工作。

2014-03-14。