賀俊凱，徐東升，王明遠，李明峰

(1.南京工業(yè)大學 測繪科學與技術(shù)學院,南京 211800；2.江蘇省地質(zhì)礦產(chǎn)局第二地質(zhì)大隊,江蘇 常州213022)

0 引言

不同基準下的坐標轉(zhuǎn)換是基于兩種坐標系下的公共控制點(以下簡稱公共點)，通過坐標轉(zhuǎn)換模型計算轉(zhuǎn)換參數(shù)而實現(xiàn)的.影響坐標轉(zhuǎn)換精度的因素不僅包括轉(zhuǎn)換模型和控制點坐標誤差，還包括公共點的空間分布等因素.關(guān)于轉(zhuǎn)換模型與控制點精度影響的研究已相對比較成熟，形成了較為完善的理論和方法[1-4].

近年來，關(guān)于公共點選取的研究取得了一定進展.王玉成等[5]、焦立芬[6]等研究發(fā)現(xiàn)，適當增加公共點個數(shù)能有效地提高轉(zhuǎn)換精度；公共點分布越均勻，所得轉(zhuǎn)換精度越穩(wěn)定.周躍寅等[7]研究了公共點圖形外圍和內(nèi)部點位精度與公共點空間分布的關(guān)系.但控制點分布的均勻性尚缺乏數(shù)學描述，影響了坐標轉(zhuǎn)換模型的精度和效率.因此，分析公共點空間分布對坐標轉(zhuǎn)換模型精度的影響，構(gòu)建顧及控制點空間分布的坐標轉(zhuǎn)換模型，對當前廣泛推行的原有坐標系向CGCS2000 的轉(zhuǎn)換具有重要的理論意義和實際應用價值.

本文根據(jù)獨占圓[8]思想提出控制點均勻度(UCP)，利用平均最鄰近距離表征控制點密度，分析控制點空間分布對坐標轉(zhuǎn)換模型精度的影響，構(gòu)建公共點空間分布約束下的坐標轉(zhuǎn)換模型，并結(jié)合實例驗證模型的有效性.

1 控制點空間分布

1.1 控制點均勻度

表征點集空間分布均勻程度的指標有多種，主要分為基于密度的方法與基于最鄰近點間距離的方法.已有研究表明，基于密度的方法并不能有效地表征點的離散程度[9].本文將基于獨占圓思想，提出UCP 概念并建立表達模型.

在n個控制點組成的集合中，以某控制點為圓心，以其到最臨近控制點歐氏距離的一半為半徑所畫的圓為該控制點的獨占圓.平面內(nèi)各獨占圓皆不重疊，最近的空間關(guān)系是相切.

如圖1 所示，d12是點1、點2 間的歐氏距離；D1為點1 至所有點的距離中最小值，稱為點1 的最鄰近距離；r1為D1的一半，定義為點1 獨占圓的半徑r1.兩點間歐氏距離dij與最鄰近距離Di的計算公式如式(1)～(2)所示：

圖1 獨占圓示意圖

式中：xi、yi為i點的坐標；xj、yj為j點的坐標.

對各獨占圓位于坐標轉(zhuǎn)換區(qū)域內(nèi)的面積求和，該面積和可以在一定程度上表示控制點均勻程度；面積和越大，控制點集越均勻.設區(qū)域內(nèi)面積和為a，其計算公式為

式中：ai為第i個獨占圓在轉(zhuǎn)換區(qū)域內(nèi)的有效面積；si為其超出區(qū)域的面積；ri為其獨占圓半徑.

顯然，a相同并不代表其均勻程度一致.結(jié)合圖1可以推斷，在a一定的情況下，各獨占圓面積越相近，控制點分布越均勻.設為與面積ai對應的有效獨占圓半徑，則有

在控制點均勻度表達中應考慮獨占圓面積大小的離散程度，其程度可以用有效獨占圓半徑中誤差σr表示

為了使獨占圓面積和能更加準確地表示均勻程度，通過離散程度對面積和a賦以權(quán)重.設權(quán)為pr，經(jīng)實驗可知，pr取值范圍為2/3～1 時，均勻度與坐標轉(zhuǎn)換精度擬合效果最佳，因此有

式中，A為坐標轉(zhuǎn)換區(qū)域面積.pr越接近1，說明各獨占圓的面積越相近，均勻度越高.

控制點等間距行列排布時控制點分布最均勻.為了更好地表述控制點均勻度，顧及獨占圓間的間隙，控制點均勻度L計算公式為

L的大小直觀反映了控制點的均勻分布程度.控制點分布越均勻，各獨占圓的面積越相近，則獨占圓面積總和a越大，L越趨近于1.即有

1.2 控制點密度

控制點密度描述轉(zhuǎn)換區(qū)域內(nèi)控制點的疏密程度，即坐標轉(zhuǎn)換時研究區(qū)范圍不變，其密度指標是指控制點的個數(shù).坐標轉(zhuǎn)換中公共點個數(shù)應大于轉(zhuǎn)換模型要求的最低個數(shù)，同時與坐標轉(zhuǎn)換計算量有關(guān).已有研究表明，坐標轉(zhuǎn)換區(qū)域內(nèi)公共點個數(shù)存在最佳值[6].

在控制點均勻度相近的情況下，控制點的個數(shù)直接影響控制點間平均最臨近距離，距離越大則控制點集分布越稀疏，反之控制點集分布越密集.因此，控制點間平均最鄰近距離dmin可以同等反映控制點分布的疏密程度，同樣存在最優(yōu)值.坐標轉(zhuǎn)換區(qū)域大小不同，公共點最優(yōu)個數(shù)也不同.設該種地形下某等級控制點的最優(yōu)平均最鄰近距離為，則最佳公共點個數(shù)t可根據(jù)下式計算：

2 控制點均勻度與密度約束下的坐標轉(zhuǎn)換模型

基于上述分析，構(gòu)建UCP 與密度約束下的坐標轉(zhuǎn)換模型，探討地方坐標系到CGCS2000 的坐標轉(zhuǎn)換流程.常用的坐標轉(zhuǎn)換模型有布爾莎模型、莫洛金斯基模型、二維四參數(shù)模型等模型.在三維坐標轉(zhuǎn)換中，布爾莎模型應用最廣泛，其通過對坐標系的平移、縮放與旋轉(zhuǎn)完成不同基準下的空間直角坐標轉(zhuǎn)換.因此，本文基于布爾莎模型展開討論，該模型共有7 個轉(zhuǎn)換參數(shù)，至少需要3 個公共已知點，其坐標轉(zhuǎn)換公式為

式中：X1、Y1、Z1為原坐標系下坐標；X2、Y2、Z2為待轉(zhuǎn)換坐標系下坐標；TX、TY、TZ為平移參數(shù)；εX、εY、εZ為旋轉(zhuǎn)參數(shù)；k為尺度參數(shù).

顧及高程異常值精度較低對大地高的影響，結(jié)合王解先等[10]研究得出大地高誤差對轉(zhuǎn)換得到的平面坐標影響很小的結(jié)論，本模型在轉(zhuǎn)換后將坐標投影到平面，并評定轉(zhuǎn)換精度.

控制點密度和均勻度從兩個角度反映了控制點的空間分布特征，兩者間存在一定聯(lián)系.t較小時，隨機選取的不同組公共點間L差異較大；隨著t值增加，L離散程度有所下降.實際轉(zhuǎn)換中，應首先根據(jù)式(9)確定公共點個數(shù)t，再尋找L能符合條件的公共點組合，因此兩者并不會相互影響.

模型從平面內(nèi)外符合精度來評估坐標轉(zhuǎn)換精度.內(nèi)符合精度 σin由公共點平差的驗后單位權(quán)中誤差評定.外符合精度σout由檢核點轉(zhuǎn)換前后坐標殘差的均方根值表示，其計算方法如式(11)所示：

式中：Δxi、Δyi為檢驗點轉(zhuǎn)換前后x、y方向上的坐標差；n為檢驗點個數(shù).

基于上述討論，建立顧及控制點空間分布約束的坐標轉(zhuǎn)換模型，轉(zhuǎn)換流程如圖2 所示.

圖2 均勻度與密度約束下的坐標轉(zhuǎn)換流程

坐標轉(zhuǎn)換具體步驟如下：

2)隨機選取t個控制點為公共點，計算均勻度L.設閾值為L0，若L＞L0，則使用其計算轉(zhuǎn)換參數(shù)；否則繼續(xù)隨機選取公共點.若循環(huán)次數(shù)m高于限定值，則會選取其中L最大的一組公共點.

3)將公共點平面坐標經(jīng)高斯反算得到大地坐標，再將大地坐標(B，L，H)換算為空間直角坐標.

4)使用布爾莎模型計算轉(zhuǎn)換參數(shù)，將檢驗點轉(zhuǎn)換至目標坐標系下.

5)將檢驗點轉(zhuǎn)換后空間直角坐標投影至平面，并依據(jù)式(11)評定精度.

3 算例分析

某區(qū)域面積5 079.2 km2，整體地形較為平坦，共有GPS 控制網(wǎng)成果點235 個，包括D 級點81 個與E級點154 個.以該區(qū)域地方坐標轉(zhuǎn)換至CGCS2000坐標為例，分析空間分布對轉(zhuǎn)換精度的影響，求解該區(qū)域的坐標轉(zhuǎn)換模型并驗證其有效性.

3.1 控制點密度對坐標轉(zhuǎn)換精度影響

為便于求解坐標轉(zhuǎn)換參數(shù)與分析討論，首先選取3 個控制點作為公共點，并逐漸增加至15 個.對于不同的t值，隨機選取若干組 0.4＜L＜0.5 的D 級控制點集計算轉(zhuǎn)換參數(shù)，使用均勻分布的E 級點作為檢核點，記錄各組內(nèi)外符合精度的精度平均值.實驗結(jié)果如表1 所示.

表1 控制點密度與轉(zhuǎn)換精度關(guān)系

由表1 可以得出以下結(jié)論：

1)隨著t值增加，基本上呈不斷增大趨勢，內(nèi)符合精度降低；

2)t=3 時，控制點組合dmin值過大，最大，＞3mm 的公共點集所占比例最高，此時坐標轉(zhuǎn)換的精度最低；

4)t＞7時，坐標轉(zhuǎn)換精度提升減緩；當t=11 時，坐標轉(zhuǎn)換精度最高，即該地形下D 級公共點值為24 km；當t＞13 時，坐標轉(zhuǎn)換精度反而降低.

3.2 公共點均勻性對坐標轉(zhuǎn)換精度影響

基于3.1 節(jié)中的結(jié)論，公共點數(shù)t取11.從D 級點中隨機選取11 個作為公共點，統(tǒng)計控制點均勻度L和外符合精度 σout.為充分反映規(guī)律，進行多組實驗，并對結(jié)果進行冪函數(shù)擬合.以上述方法計算得到的L和σout的關(guān)系，如圖3 所示，圖中每個點表示每組轉(zhuǎn)換的 σout.

圖3 L 對σout 的影響擬合曲線

由圖3 可知：

1)L較低時的公共點計算得到的 σout離散且基本偏大；

2)隨著L的提高，σout逐漸向2 mm 以下集中，且不再出現(xiàn)較大誤差；

3)隨L的提高，σout降幅越來越??；在L超過0.40 后，擬合曲線趨于平緩，L增加對精度的影響作用有限.

因此，在實際應用中選取L大于0.40 的組合，可有效保證坐標轉(zhuǎn)換精度.

4 結(jié)語

本文通過研究控制點均勻度與密度對坐標轉(zhuǎn)換精度的影響，構(gòu)建了顧及控制點空間分布的坐標轉(zhuǎn)換模型，解決了坐標轉(zhuǎn)換中公共點選取缺乏標準的問題.具體結(jié)論如下：

1)公共點選取應兼顧轉(zhuǎn)換區(qū)域控制點的均勻度與密度，先確定公共點密度，再決定均勻度；顧及控制點空間分布的坐標轉(zhuǎn)換模型可有效提高轉(zhuǎn)換精度；

2)坐標轉(zhuǎn)換精度隨公共點均勻度增加呈提高趨勢，L達到0.4 后可滿足坐標轉(zhuǎn)換的精度要求；

3)公共點密度對應坐標轉(zhuǎn)換精度具有最佳值；以平坦地區(qū)D 級控制點作為公共點時，dmin取24 km所得轉(zhuǎn)換精度最高.

本文只考慮了均勻度與密度指標，地形起伏較大區(qū)域內(nèi)坐標轉(zhuǎn)換精度可能與所選公共點的豎直方向分布有關(guān)，后續(xù)研究中應考慮控制點高程的影響.