現(xiàn)有測繪成果向CGCS2000成果轉(zhuǎn)換的實踐*

黃發(fā)秀，鐘連琨

(廣西壯族自治區(qū)國土測繪院，廣西 南寧 530023)

現(xiàn)有測繪成果向CGCS2000成果轉(zhuǎn)換的實踐*

黃發(fā)秀，鐘連琨

(廣西壯族自治區(qū)國土測繪院，廣西 南寧 530023)

在比較各種坐標(biāo)系之間轉(zhuǎn)換模型適用性的基礎(chǔ)上，介紹現(xiàn)有各種坐標(biāo)系成果向CGCS2000成果轉(zhuǎn)換的方法。重點論述回歸分析模型參數(shù)的計算方法、模型的檢驗及不同軟件格式的數(shù)據(jù)成果在轉(zhuǎn)換中所遇到的技術(shù)問題及解決方法。

測繪成果；坐標(biāo)轉(zhuǎn)換；GNSS控制網(wǎng)；坐標(biāo)修正模型；CGCS2000

Abstract：This thesis is based on the applicability of the conversion model between various coordinate systems，having introduced the methods of converting the results of the existing various coordinate systems into that of CGCS2000.The emphases of this thesis are the calculation method of regression analysis model parameters，the test of model，and the technical problems and their solutions encountered in the conversion of data in different software formats.

Keywords：surveying and mapping results；coordinate conversion；GNSS contral network；coordinate correct model；CGCS2000

0 引言

根據(jù)國家測繪地理信息局《關(guān)于加快2000國家大地坐標(biāo)系推廣使用的通知》(國測國發(fā)[2013]11號)的要求，2016年需完成縣級現(xiàn)有數(shù)字化基礎(chǔ)測繪成果向2000國家大地坐標(biāo)系的轉(zhuǎn)換工作。為此，需要將20世紀(jì)50年代和80年代基于參考橢球的國家大地坐標(biāo)系統(tǒng)(參心坐標(biāo)系)——1954年北京坐標(biāo)系和1980西安坐標(biāo)系的測繪成果資料向CGCS2000進(jìn)行轉(zhuǎn)換[1-4]。廣西壯族自治區(qū)國土測繪院經(jīng)過競標(biāo)承擔(dān)了廣西賓陽縣原有基礎(chǔ)測繪數(shù)據(jù)轉(zhuǎn)換到CGCS2000的項目。由于這些資料時間久遠(yuǎn)、種類繁多、精度不一，給轉(zhuǎn)換工作帶來了較大困難。為此，本文闡述了轉(zhuǎn)換工作中遇到的技術(shù)難點及處理方法，供同行們參考。

1 項目概況

廣西賓陽縣地處北回歸線以南，東與貴港市交界，西、南連接南寧市，北與來賓市為鄰。東西最大橫距75.2 km，南北最大縱距60.6 km，面積2 314 km2。地勢南高北低，由西南向東北傾斜?？h境東、南、西三面邊緣土山環(huán)亙，北及東北面邊緣石山群立，中部為不閉合盆地，屬大片沖積平原，海拔高度在100 m～120 m之間。賓陽縣境內(nèi)，鐵路、高鐵、各等級公路相與交錯，四通八達(dá)，交通十分便利?？h城南距廣西首府南寧市70 km，東距華南第三大內(nèi)河港口——貴港港口僅70 km。

廣西賓陽縣國土資源局自20世紀(jì)末設(shè)立至今，為適應(yīng)各個時期國民經(jīng)濟(jì)和社會發(fā)展的需要，基于1954年北京坐標(biāo)系和1980西安坐標(biāo)系進(jìn)行了大量的測繪工作，積累了大量的成果資料。這些測繪成果從表現(xiàn)形式上分有數(shù)據(jù)、信息、圖件及相關(guān)技術(shù)資料；從種類上分有各等級的控制測量成果、各種比例尺的地形圖、專題圖；從成果的載體上分有紙、光、電磁3類；從數(shù)據(jù)格式上分有xlsx、doc、dwg、sdd、sdb、WP、WT、WL、tif、pmf、mxd等格式；從采用的軟件上分有Office、AutoCAD、SuperMap、ArcGIS、MapGIS等。這些成果是不同時期不同單位為完成當(dāng)時的工作，在不同坐標(biāo)系或同一坐標(biāo)系下利用不同時期的起算成果完成的局部成果。其特點是技術(shù)方案參差不齊，精度不一，資料不全。其中，有些控制點已被嚴(yán)重?fù)p毀，名存實亡無法恢復(fù)，這給轉(zhuǎn)換工作帶來了極大困難。

2 技術(shù)方案

通過對所收集的成果資料進(jìn)行認(rèn)真分析后，設(shè)計了如圖1所示的技術(shù)方案。在全縣范圍內(nèi)新布設(shè)一個D級GNSS控制網(wǎng)，該網(wǎng)既要滿足本次坐標(biāo)轉(zhuǎn)換，又要顧及未來測繪需要。要求D級GNSS網(wǎng)覆蓋全縣區(qū)域，同時新、舊兩期要有分布均勻、數(shù)量適當(dāng)?shù)闹睾峡刂泣c，便于坐標(biāo)轉(zhuǎn)換模型的建立。利用自主開發(fā)的坐標(biāo)轉(zhuǎn)換軟件，將賓陽縣國土資源局提供的原有數(shù)據(jù)成果轉(zhuǎn)換為CGCS2000國家大地坐標(biāo)系成果。

圖1 坐標(biāo)轉(zhuǎn)換技術(shù)方案流程圖Fig.1 The technical flow of coordinate conversion

3 D級GNSS控制網(wǎng)的建立及成果分析

以5個C級以上GPS控制點為起算，按邊連接或邊網(wǎng)連接的方式新布設(shè)62點(其中重合原有等級GPS點19個)共67個點組成D級GNSS靜態(tài)平面控制網(wǎng)。該網(wǎng)共305條邊，最短邊為2 651.49 m，最長邊為21 895.21 m，平均邊長為9 424.81 m。采用3～8臺套中海達(dá)V30多功能全球定位系統(tǒng)雙頻接收機(jī)采集數(shù)據(jù)，時段觀測長度大于150 min，觀測時段數(shù)為1.94?；€解算采用Trimble Business Center-Version 3.50軟件處理，經(jīng)處理后該網(wǎng)同步環(huán)最大閉合差2.038×10-6·D，異步環(huán)最大閉合差5.996 0×10-6·D，重復(fù)邊最弱相對中誤差為0.669×10-6·D。網(wǎng)平差采用武漢測繪科技大學(xué)研發(fā)的PowerADJ4.0v軟件進(jìn)行，平差后該網(wǎng)最弱邊相對誤差1/695 095，最弱點位中誤差3.070 cm。高程采用廣西壯族自治區(qū)測繪地理信息局檔案館的似大地水準(zhǔn)面精化成果。平差后新、舊兩期19個重合點的成果對比情況，如表1所示。

由表1可見：兩期成果兩個坐標(biāo)系之間的誤差都較大，尤其是1954年北京坐標(biāo)系的成果，誤差高達(dá)2 m以上、符號基本上為“+”，具有明顯的系統(tǒng)特征；在1980西安坐標(biāo)系中，兩期成果誤差小一些，具有偶然誤差的一些特點，但是否具有同一分布仍需深入分析。初步分析后判斷，可能是局部舊網(wǎng)所用的起算數(shù)據(jù)為廣西賓陽縣國土資源局布網(wǎng)的成果及交叉使用相互成果數(shù)據(jù)。

表1 新、舊兩期重合點坐標(biāo)的比較

4 坐標(biāo)修正模型的參數(shù)計算

經(jīng)上述分析得出：賓陽縣國土資源局所有由上述控制成果所衍生的測繪資料并不算嚴(yán)格意義上的1954年北京坐標(biāo)系或1980西安坐標(biāo)系成果。因此，要對其進(jìn)行必要的修正，使之成為真正意義上的1954年北京坐標(biāo)系或1980西安坐標(biāo)系成果。然而由于這些成果施測年代久遠(yuǎn)，資料不齊全，測量報告及所用的起算點等基本資料都沒有，因此無法使用原有的測量數(shù)據(jù)對其重新進(jìn)行改算。

目前，不同坐標(biāo)系(近似1980西安坐標(biāo)系也算一個坐標(biāo)系)間坐標(biāo)變換的常用方法有相似變換法、3參數(shù)法、7參數(shù)法等。這些方法都是通過對空間圖形進(jìn)行平移、旋轉(zhuǎn)以及縮放等手段來實現(xiàn)，它們都具有幾何關(guān)系明確、解求參數(shù)方便、變換公式規(guī)則、具有適當(dāng)外推的優(yōu)點，但對坐標(biāo)系的局部變形顧及較差[4]。由于測量誤差不可避免，如折光誤差、尺長誤差、信號傳播誤差等，因此系統(tǒng)中各點的坐標(biāo)值并不嚴(yán)格處于確定的數(shù)學(xué)模型下，從而使得各坐標(biāo)系統(tǒng)中產(chǎn)生局部性或整體性偏差，于是就存在變換模型的適度問題?；貧w分析是統(tǒng)計學(xué)中的一個重要分支，在自然科學(xué)領(lǐng)域和社會科學(xué)領(lǐng)域有著十分廣泛的應(yīng)用，是處理多個變量之間相關(guān)關(guān)系的一種數(shù)學(xué)方法。它是建立在對客觀事物進(jìn)行大量試驗和觀測的基礎(chǔ)上，用于尋找隱藏在哪些看上去是不確定的現(xiàn)象中的統(tǒng)計規(guī)律的數(shù)理統(tǒng)計方法。按照控制論的觀點，一個系統(tǒng)的內(nèi)部結(jié)構(gòu)并不清楚，但是只要掌握了它的外部聯(lián)系的狀態(tài)方式，闡明它的輸入和輸出的因果關(guān)系，就可以認(rèn)為充分掌握了這個系統(tǒng)。在坐標(biāo)轉(zhuǎn)換中，將不同坐標(biāo)系之間的變換過程視為黑箱，輸入的是變換前坐標(biāo)系的舊坐標(biāo)，輸出的是新坐標(biāo)系的新坐標(biāo)。應(yīng)用回歸分析方法建立輸入和輸出之間的數(shù)學(xué)關(guān)系，從而確定它們之間的變換關(guān)系。

鑒于上述原因，本文設(shè)計了賓陽縣1954年北京坐標(biāo)系、1980西安坐標(biāo)系的測繪成果資料向2000國家大地坐標(biāo)系之間的轉(zhuǎn)換的總體方案：

1)通過各種技術(shù)手段對現(xiàn)有不是嚴(yán)格意義上的1980西安坐標(biāo)系測繪成果進(jìn)行修正，使其成為嚴(yán)格意義的1980西安坐標(biāo)系成果。

2)所有測繪成果資料坐標(biāo)系之間的轉(zhuǎn)換以1980西安坐標(biāo)系的成果為基礎(chǔ)，分別向2000國家大地坐標(biāo)系、1954年北京坐標(biāo)系轉(zhuǎn)換。而原有的1954年北京坐標(biāo)系的成果不作任何處理，作為歷史資料保留。

以上述19個1980西安坐標(biāo)系下的重合點為樣本[3]，利用統(tǒng)計分析方法排除3個不具有同一誤差分布的點后，采用布爾莎(Bursa)7參數(shù)模型和回歸分析模型分別求取修正模型，其樣本的內(nèi)、外符合殘差的中誤差分別為9.4 cm與2.1 cm。結(jié)果表明，回歸分析模型在此樣本中精度更優(yōu)，因此在坐標(biāo)修正模型上采用回歸分析模型。

坐標(biāo)修正模型的方程式[2]為：

(1)

設(shè)Xi=X舊i-X0，Yi=Y舊i-Y0，ωi=X新i-X舊i，x1i=X舊i-X0，x2i=Y舊i-Y0，x3i=(X舊i-X0)2，x4i=(X舊i-X0)(Y54i-Y0)，x5i=(Y舊i-Y0)2，x6i=(X舊i-X0)3，x7i=(X舊i-X0)2(Y舊i-Y0)，x8i=(X舊i-X0)(Y舊i-Y0)2，x9i=(Y舊i-Y0)3。

其中，X0、Y0為項目范圍坐標(biāo)的平均值。則多項式(1)可以變換成多元線性方程式：

y=β0+β1x1i+β2x2i+β3x3i+β4x4i+β5x5i+β6x6i+β7x7i+β8x8i+β9x9i+…+ωi

(2)

設(shè)因變量為Y(預(yù)測對象)，自變量為X，(xi，yi)為X、Y的觀測值，i=1，2，...，n，則

Y=BX+L

(3)

式中：

YT=

XT=

根據(jù)最小二乘法(OLS)原理有：

N=BTB

(4)

U=BTW

(5)

回歸模型的系數(shù)為：

(6)

以回歸分析模型設(shè)定的模型是正確的為前提條件來開展模型參數(shù)估計和假設(shè)檢驗。這種假設(shè)一旦模型設(shè)定有問題，那么據(jù)此進(jìn)行參數(shù)估計和假設(shè)檢驗也會有問題。由于多項式(1)中β3的自變量是由xi、yi這兩個基本自變量派生出來的，因此不能同時完全滿足線性回歸模型的線性、正交、獨立、正態(tài)分布4個基本假定。若用經(jīng)典線性回歸模型中F、t檢驗判斷其顯著性，往往會得出錯誤的結(jié)論。

從式(1)中可以看出，增加一個系數(shù)與不增加之前的兩個模型相比，其實質(zhì)就是對新增系數(shù)是否被限制為“0”之間的差別。利用模型之間的這種嵌套關(guān)系，可以對交互項的存在進(jìn)行統(tǒng)計檢驗。

H0∶β3=0?H1∶β3≠0

構(gòu)造F檢驗統(tǒng)計：

(7)

如果F>Fα((df2-df1)，df1)，拒絕零假設(shè)H0，認(rèn)為受限制模型與不受限制模型在數(shù)據(jù)擬合上是有顯著差異的，即交互項XY的偏回歸系數(shù)β3顯著區(qū)別于0。

其中：

(8)

(9)

(10)

(11)

而點位中誤差的SSR、SSE、SST計算公式為：

(12)

(13)

(14)

同理，通過式(7)可以檢驗判斷模型對點位殘差是否顯著、有效，從而達(dá)到對模型的檢驗。

在本項目中，按此方法進(jìn)行檢驗，修正模型采用最高至二次項，即β6起后面的系數(shù)項全部省略。這樣既能滿足坐標(biāo)修正的精度要求，又能簡化工作。

5 坐標(biāo)系之間轉(zhuǎn)換模型及參數(shù)計算

以新布設(shè)的D級GNSS網(wǎng)點及項目區(qū)附近的B、C級GPS控制點共77個點作為起算數(shù)據(jù)。這些點覆蓋全縣縣域，共同具有1954年北京坐標(biāo)系、1980西安坐標(biāo)系、2000國家大地坐標(biāo)3套成果。坐標(biāo)轉(zhuǎn)換模型參數(shù)的求解方法是：分兩次分別從77個起算數(shù)據(jù)中隨機(jī)選取6個不重復(fù)點作為外檢核點，其余點作為參數(shù)求取的起算點。采用布爾莎(Bursa)模型各自求取一套轉(zhuǎn)換參數(shù)，并用所選的外部檢核點檢驗其參數(shù)可靠性；在二次轉(zhuǎn)換參數(shù)都達(dá)到規(guī)范規(guī)定的精度后，再將檢核點也作為參數(shù)計算的起算點，重新計算一套轉(zhuǎn)換參數(shù)；計算的結(jié)果經(jīng)內(nèi)部檢核合格后，作為本項目坐標(biāo)系間轉(zhuǎn)換參數(shù)的最終成果。

為了檢驗修正模型和轉(zhuǎn)換模型的正確性，本項目專門成立了一個檢核小組對模型的正確性進(jìn)行檢驗。檢核的方法是采用廣西COSR站提供的數(shù)據(jù)鏈，利用中海達(dá)V30儀器按一級GPS控制測量點的精度要求實施，在全縣范圍內(nèi)，分布均勻兼顧各個區(qū)域，一共檢查了26個等級控制點或圖根埋石點，檢測結(jié)果與轉(zhuǎn)換結(jié)果的誤差都小于70.7 mm。說明本項目確定的坐標(biāo)修正模型和坐標(biāo)系之間的轉(zhuǎn)換模型理論依據(jù)是正確的，方法是合理的，結(jié)果也是可靠的。

6 圖形數(shù)據(jù)文件的轉(zhuǎn)換

本項目需要進(jìn)行轉(zhuǎn)換的圖形數(shù)據(jù)格式有4類，分別為：CAD、MapGIS、ArcGIS和SuperMap格式。轉(zhuǎn)換的原則是按不同格式采用原有相應(yīng)平臺，并選擇相應(yīng)功能逐一進(jìn)行轉(zhuǎn)換。

6.1 CAD格式數(shù)據(jù)文件的轉(zhuǎn)換

CAD格式數(shù)據(jù)的轉(zhuǎn)換是本項目中工作量最大、數(shù)據(jù)情形最復(fù)雜的一項工作。它既有1954年北京坐標(biāo)系和1980年西安坐標(biāo)系的成果，又有近似1980西安坐標(biāo)系成果。這些成果資料分布區(qū)域零散，在工作中得一一區(qū)別對待。在CAD格式數(shù)據(jù)轉(zhuǎn)換工作中遇到的技術(shù)問題和處理方法如下：

1)根據(jù)該項目特點，以Microsoft Visual Studio 2012.NetFramework 版本為平臺，利用C#語言開發(fā)了獨立于AutoCAD系統(tǒng)外運(yùn)行的GTY-Coord坐標(biāo)轉(zhuǎn)換軟件。該軟件可支持三維3參數(shù)、布爾莎7參數(shù)、平面4參數(shù)、平面多項式擬合模型的坐標(biāo)轉(zhuǎn)換。該軟件是基于直接讀取CAD文件數(shù)據(jù)庫中相關(guān)坐標(biāo)數(shù)據(jù)的屬性，然后用內(nèi)嵌的修正模型、坐標(biāo)系轉(zhuǎn)換模型逐一修正并置換其坐標(biāo)值的原理開發(fā)。對其它圖元屬性等均不更改，保證轉(zhuǎn)換后的數(shù)據(jù)文件與原有數(shù)據(jù)除坐標(biāo)值外均保持一致，不存在圖元內(nèi)容和屬性的更改、丟失等情況發(fā)生。

2)選擇GTY-Coord軟件對原有近似1980西安坐標(biāo)系的圖形文件進(jìn)行坐標(biāo)修正，使之轉(zhuǎn)為正確的1980西安坐標(biāo)系成果。修正后的成果及不需修正的成果，選擇GTY-Coord軟件中坐標(biāo)系之間的轉(zhuǎn)換功能，經(jīng)二道程序，分別得到2000國家大地坐標(biāo)系、1954年北京坐標(biāo)系的兩套成果資料。在處理1∶10 000地形圖數(shù)據(jù)文件時，由于坐標(biāo)系之間的差異，實地上同一點的經(jīng)、緯度各不相同，而3個坐標(biāo)系下的地形圖分幅經(jīng)、緯度是相同的。為使轉(zhuǎn)換后的圖幅嚴(yán)格與國家分幅一致，在數(shù)據(jù)轉(zhuǎn)換之前應(yīng)先將這些圖幅拼接成一整塊。坐標(biāo)系之間轉(zhuǎn)換完成后，再用CASS 9.0 軟件中的小比例尺分幅功能按國家標(biāo)準(zhǔn)切塊分幅，分幅后的圖廓仍然按原同名圖幅整飾。

6.2 MapGIS數(shù)據(jù)文件的轉(zhuǎn)換

本項目需要對新、舊版本MapGIS平臺生成的縣、鄉(xiāng)兩級多個土地利用總體規(guī)劃圖數(shù)據(jù)庫進(jìn)行轉(zhuǎn)換。這項工作主要采用MapGIS 6.7和MapGIS K9兩個平臺來完成。在轉(zhuǎn)換過程中，所遇到的技術(shù)問題及解決方法如下：

1)由于MapGIS6.7軟件不支持2000國家大地坐標(biāo)的橢球參數(shù)且沒有批量轉(zhuǎn)換功能。因此，MapGIS6.7版本的數(shù)據(jù)只能利用MapGIS K9軟件進(jìn)行批量轉(zhuǎn)換。

2)由于MapGIS6.7與MapGISK9軟件在開發(fā)上存在著根本性的區(qū)別，即兩個版本軟件的底層數(shù)據(jù)庫的開發(fā)平臺發(fā)生了根本性的改變，造成區(qū)、線、點、注記這4類文件的結(jié)構(gòu)完全不同，軟件之間的數(shù)據(jù)互用需要相互轉(zhuǎn)換。這樣數(shù)據(jù)轉(zhuǎn)入、轉(zhuǎn)出后，部分?jǐn)?shù)據(jù)會發(fā)生改變。區(qū)數(shù)據(jù)在拓?fù)涔δ苤?，兩者的?guī)則不同，如MapGISK9軟件對部分在MapGIS6.7軟件中拓?fù)錂z查報錯的數(shù)據(jù)不報錯。對這類問題的處理方法是：先將數(shù)據(jù)在MapGISK9中進(jìn)行轉(zhuǎn)換并導(dǎo)出后，再在MapGIS6.7軟件中拓?fù)渲匦略靺^(qū)，然后將Lable點屬性重新賦值于區(qū)中。點和注記文件導(dǎo)出后，在MapGISK9軟件中進(jìn)行追加合并(MapGIS 6.7軟件的此功能會使注記的字體、大小和子圖大小發(fā)生改變，而線文件則不受影響，不需更改)。

3)原數(shù)據(jù)庫存在拓?fù)溴e誤、區(qū)重疊、弧段相交重疊、小于限差范圍的小區(qū)、多部件組合區(qū)等問題。如果不對這些問題進(jìn)行處理會造成數(shù)據(jù)轉(zhuǎn)換后區(qū)部分?jǐn)?shù)據(jù)缺失、變形、拆分等現(xiàn)象。為此，在數(shù)據(jù)轉(zhuǎn)換前應(yīng)重新拓?fù)湓靺^(qū)。

4)原數(shù)據(jù)文件日期BGRQ字段為空時，MapGISK9軟件轉(zhuǎn)換后對日期等屬性默認(rèn)賦值為1900-1-1，解決的方法是使用MapGIS6.7軟件屬性庫管理功能清空BGRQ字段。

5)MapGISK9軟件不兼容MapGIS6.7軟件中版面類注記，此類數(shù)據(jù)在MapGISK9軟件轉(zhuǎn)換后會丟失。解決方法是：先利用MapGISK9軟件進(jìn)行轉(zhuǎn)換，再用MapGIS6.7軟件對缺失內(nèi)容逐一補(bǔ)回。

6.3 ArcGIS格式、SuperMap格式數(shù)據(jù)文件的轉(zhuǎn)換

需要轉(zhuǎn)換的數(shù)據(jù)庫有：耕地質(zhì)量等別成果庫、耕地質(zhì)量等別圖數(shù)據(jù)庫、集體土地所有權(quán)地籍?dāng)?shù)據(jù)庫、農(nóng)村宅基地及集體建設(shè)用地登記發(fā)證地籍?dāng)?shù)據(jù)庫。這些數(shù)據(jù)的坐標(biāo)系間轉(zhuǎn)換主要采用超圖6.0.5軟件平臺中7參數(shù)Coordinate Frame算法投影轉(zhuǎn)換功能完成。在本項目數(shù)據(jù)轉(zhuǎn)換過程中，所遇到的技術(shù)問題及處理方法如下：

1)耕地質(zhì)量等別數(shù)據(jù)庫坐標(biāo)系成果為1980西安坐標(biāo)系成果，可直接采用超圖6.0.5軟件平臺中的相應(yīng)功能完成。

2)集體土地所有權(quán)地籍?dāng)?shù)據(jù)庫包含兩方面內(nèi)容：一是耕地權(quán)屬數(shù)據(jù)，二是農(nóng)村宅基地權(quán)屬數(shù)據(jù)。這兩類數(shù)據(jù)的情形不盡相同，具體處理方式如下：

①耕地權(quán)屬數(shù)據(jù)庫的數(shù)據(jù)屬于1980西安坐標(biāo)系成果，可以按上述方法直接轉(zhuǎn)換。該數(shù)據(jù)庫中包含影像文件，此類文件較大，計算機(jī)無法支撐如此大的文件在同一數(shù)據(jù)庫中進(jìn)行坐標(biāo)轉(zhuǎn)換，需要單獨處理后再進(jìn)行合并。

②農(nóng)村宅基地權(quán)屬數(shù)據(jù)庫比較復(fù)雜，既有大比例尺圖層數(shù)據(jù)，又有小比例尺圖層數(shù)據(jù)。其中，大比例尺數(shù)據(jù)有測量控制點、宗地界線、界址點等地籍要素內(nèi)容，其成果因所在區(qū)域不同，既有1980西安坐標(biāo)系成果，又有近似1980西安坐標(biāo)系成果。而小比例尺數(shù)據(jù)有行政區(qū)劃界線、地類圖斑、注記等數(shù)據(jù)，其成果屬于1980西安坐標(biāo)系成果。此類數(shù)據(jù)的處理方法是：對大比例尺地形圖數(shù)據(jù)，根據(jù)數(shù)據(jù)所在區(qū)域和坐標(biāo)修正的方案對其相應(yīng)的坐標(biāo)數(shù)據(jù)進(jìn)行修正處理，經(jīng)坐標(biāo)修正后，全部大比例尺數(shù)據(jù)均統(tǒng)一為1980西安坐標(biāo)系成果；然后將此數(shù)據(jù)重新合并到小比例尺數(shù)據(jù)中，采用超圖6.0.5軟件平臺中的相應(yīng)功能完成。

7 結(jié)論

通過前文分析，筆者得出以下結(jié)論：

1)只有運(yùn)用數(shù)據(jù)統(tǒng)計的方法認(rèn)真分析現(xiàn)有測繪成果質(zhì)量情況，通過充分的論證，正確鑒別出哪些成果具有同一模型從而選擇不同的方法加以處理，才能使轉(zhuǎn)換后的成果達(dá)到最優(yōu)。

2)線性回歸模型是此類坐標(biāo)轉(zhuǎn)換的最佳選擇。在模型的選擇上，應(yīng)以簡約為原則，否則既浪費了自由度，有時又會導(dǎo)致估計精度下降。由于高階項具有極其不穩(wěn)定性的特點，具有指數(shù)式放大樣本噪音，模型極其不穩(wěn)定。因此，在低次方項滿足要求的前提下，建議使用低次模型。

3)在數(shù)據(jù)轉(zhuǎn)換中，應(yīng)盡量使用原有版本的軟件進(jìn)行轉(zhuǎn)換。最佳的方法是通過讀取各數(shù)據(jù)庫中坐標(biāo)屬性數(shù)據(jù)實行點對點置換，避免數(shù)據(jù)轉(zhuǎn)入、轉(zhuǎn)出，減少數(shù)據(jù)丟失，提高工作效率。

[1] 王福保.概率論及數(shù)理統(tǒng)計[M].上海：同濟(jì)大學(xué)出版社，1994.

[2] 朱華統(tǒng).大地坐標(biāo)系統(tǒng)的建立[M].北京：測繪出版社，1986.

[3] 鐘連琨，黃發(fā)秀.GPS水準(zhǔn)擬合方式的統(tǒng)計分析及擬合方式的選擇[J].工程勘察，2000(5)：61-63.

[4] 鐘連琨.回歸分析法在GPS測量坐標(biāo)變換的應(yīng)用[J].測繪通報，2012(S1)：55-57.

PracticeofConvertingExistingSurveyingandMappingCoordinatesintoCGCS2000Results

HUANG Fa-xiu，ZHONG Lian-kun

(GuangxiInstituteofLandandSurveyingandMapping，NanningGuangxi530023，China)

2017-06-02——以廣西賓陽縣為例

P 226+.3

A

1007-9394(2017)03-0001-04

黃發(fā)秀(1963～)，男，廣西扶綏縣人，高級工程師，現(xiàn)主要從事測繪生產(chǎn)和管理方面的工作。