MapSource軟件在物探前期工作中的運用

劉彥華，楊青峰，王正亮

(江蘇省水文地質海洋地質勘查院，江蘇 淮安 223400)

0 引言

手持GPS是目前項目前期勘察過程中較為常見的定位工具，涉及的領域廣泛，例如工區范圍調查、地質調繪、物探測線清理，鉆孔孔位初定等。物探測線上的物理點(例如地震炮點、檢波點、電極點、記錄點等)的坐標可以通過手持GPS實現，野外數據采集過程中遇到的明顯的地質現象也可以通過手持GPS定位，為以后資料解譯提供相對精確的定位幫助(李鵬等，2008)。MapSource作為Garmin系列手持GPS的配套軟件，不僅能提供點位、導航等信息，也可以作為物探測線導航，還可以在Google Earth軟件中直觀地顯示航線和航跡，為物探解譯發揮更大的作用。

1 坐標轉換問題

在工程物探測量中，常見的是業主提供的平面直角坐標，例如北京54和西安80(國外常用UTM坐標)。而手持機系統默認的是WGS－84坐標，所以2種不同坐標系之間存在著平移和旋轉關系(WGS－84坐標系與我國應用的坐標系之間的誤差約為80 m)。通過手持GPS進行絕對定位必須進行坐標轉換，轉換后的絕對定位精度可由80 m提高到5～10 m，通過區域控制點校正定位精度可達到0～3 m，足以滿足工程物探的精度要求。

WGS－84坐標和北京54(或西安80)坐標系屬于兩個不同的模擬地球形狀的橢球模型。北京54坐標系采用的是克拉索夫斯基橢球模型，而WGS－84坐標系有1個WGS－84橢球，其常數采用 IUGG第17屆大會大地測量常數的推薦值。兩個橢球模型之間轉換首先要考慮兩者不同的橢球基準，即DA、DF。DA、DF分別是2個橢球基準之間半長軸的差和扁率倒數的差。WGS－84橢球2個最常用的幾何常數為：長半軸：(6 378 137±2)m;扁率：1∶298.257 223 563。在進行坐標轉換時，這2個參數是固定的，即：WGS－84→北京54：DA：－108;DF：0.000 000 5、WGS－84→西安80：DA：－3;DF：0。

橢球的基準轉換后，建立橢球的原點仍不一致，還需要在DX、DY、DZ這3個空間平移參量，以便將2個不同的橢球原點重合，這樣才能使2個坐標系的橢球完全轉換。由于各地的地理位置不同，因此在各個地方的這3個坐標軸平移參量也不同。

2 運用公式求取DX、DY、DZ三參數

傳統求三參數的方法即通過公式求取，具體計算方法如下。

第一步：搜集應用區域內GPS網3個以上網點WGS－84坐標系B、L、H值及我國坐標系(北京54或西安80)B、L、h、x值(B、L、H分別為大地坐標系中的大地緯度、大地經度及大地高，h、x分別為大地坐標系中的高程及高程異常。

第二步：計算不同坐標系三維直角坐標值，對應橢球的有關常數見表1，計算公式如下：

表1 不同坐標系對應橢球的有關常數Table 1 Relevant constants of corresponding ellipsoidsfor different coordinate systems

其中，X、Y、Z為大地坐標系中的三維直角坐標;A為大地坐標系對應橢球之長半軸;e2為大地坐標系對應橢球第一偏心率;N為該點的卯酉圈曲率半徑，N=A/(1－e2sin2B)/2;H=h+x，H為北京54或西安80坐標系中的大地高(劉大杰等，1997;杜大彬等，2007)

第三步：求出 DX，DY，DZ。即利用 WGS－84 坐標系的X、Y、Z值，減去我國坐標系的對應值，得出實現坐標系統轉換的3個參數(應算出WGS－84與北京和西安坐標系2套參數)。

第四步：參數驗證。參數計算之后必須對其進行驗證。驗證的方法是在應用區域內選擇5個以上控制點進行實測，實測值與測繪部門提供的理論值對比，若平均誤差≤10 m，則計算出的參數可以使用，否則要重新計算或查找出現問題的原因。

3 用MapSource求DX、DY、DZ三參數

傳統求三參數的方法仍然需要現場驗證，通過不斷修改參數來達到定位效果。通過MapSource軟件反復推算三參數可以達到現場驗證的目的，節省大量的時間。下面通過1組實例說明MapSource軟件的優越性。

選取能控制工程區域的3個以上三角點的WGS－84坐標系B、L及我國坐標系(北京54或西安80)X、Y值。以3點407P、408P和424P為例進行計算。

表2 福建某工程三角點數據Table 2 Triangulation point data of a project in Fujian

3.1 MapSource軟件預設

MapSource V6.15.6版本為中文版本，雙擊圖標打開軟件，單擊菜單欄的“編輯”，在彈出的下拉菜單中點擊“首選項”命令，彈出“首選項”對話框。

單擊“位置”標簽，對下列選項進行設置：網格(用戶定義的網格)—從下拉列表框中可選擇“緯度/經度 hddd°mm'ss.s″”(度、分、秒格式);基準(用戶定義的基準)—從下拉列表框中可選擇“WGS－84”(Garmin Corp，2004)(圖1)。

3.2 輸入3點的WGS－84坐標

單擊菜單欄的“編輯”，在彈出的下拉菜單中點擊“新航點”，在彈出的“新航點”對話框中輸入對應點的經緯度坐標，并保存設置。重復該步驟可以得到3個點的經緯度坐標(Garmin Corp，2006)。

圖1 三角點的WGS－84坐標Fig.1 WGS-84 coordinates of the triangulation points

3.3 輸入DA、DF數值

在“位置”對話框中將坐標系改成西安80坐標系，輸入 DA、DF 數值(Garmin Corp，2004)。

單擊“位置”標簽，對應的選項進行下列設置：網格—用戶定義的坐標，這時后面出現一個“屬性”按鈕;基準—用戶定義的基準，后面也出現“屬性”按鈕。下面對這些屬性按鈕進行設置。

首先對“用戶定義的網格”進行設置，點擊“屬性”按鈕，彈出的對話框中，分別對中央經線、比例系數、東西偏差假定值、南北偏差假定值進行賦值。

中央經線：通常坐標X是7位，Y是8位，其中Y坐標前2位是帶號，表明是中央子午線(通常指經度L)。一般來說，20左右的應為6度分帶，帶號為n則6n－3;40左右的為3度分帶則3n。

比例系數為1，東西偏差假定值為500 000 m，南北偏差假定值為0 m。

再次對“用戶定義的基準屬性”進行設置，點擊“屬性”按鈕，彈出的對話框中分別有DeltaX、DeltaY、DeltaZ、Delta半長軸、Delta平面度等參數。Delta半長軸、Delta平面度2個常數即DA、DF。DA、DF分別填寫－3、0(圖2)。

圖2 MapSource下網格參數設置Fig.2 Grid parameter setting in MapSource

3.4 運用MapSource快速求取三參數

假定DX、DY、DZ為0，可以得出在西安80坐標系統下三角點的假定數值，將假定數值與實際數值相減，并不斷地修正DX、DY、DZ，可以得到最接近實際數值的假定數值，即擬合程度最高的DX、DY、DZ(圖 3)(馬翠玲，2010)。

通過編寫的程序反復計算，得出DX=－111、DY=－46、DZ=－5。該參數在3個點的假定數值與實際數值對比見表3。

通過MapSource軟件可以在室內計算出DX、DY、DZ三參數，并與實際數值對比分析，誤差極小。經大量實際驗證，該方法簡單適用。

4 MapSource與Google Earth軟件結合布置測線

4.1 物探測線編輯

物探野外工作主要以物探測線為主，將布置好的測線按照里程結合線位弧度縮小誤差，以物探點的形式輸入到MapSource軟件中，并編輯好航線(圖4)。野外操作過程中可以導航此航線，隨時提醒偏航距，調整清理線位，達到精確勘探的目的(圖5)(劉凡珍等，2007;王會峰等，2008)。

圖3 MapSource下三參數快速設置Fig.3 Rapid setting of three parameters in MapSource

表3 三角點平面直角坐標實際值和假定值關系Table 3 Relationship between actual values and assumed values of plane rectangular coordinates of triangulation points

4.2 Google Earth顯示的物探線位

在MapSource軟件的菜單欄“視圖”最下一個命令欄，有與Google Earth軟件接口，可以在 Google Earth軟件中查看測線(圖6)。這項功能可以大大提高野外工作的效率，提前做好分段清理的計劃，節省大量時間。

線路途徑的小路清晰可見，里程樁號位置可定，可以預先做好分段數據采集處理，計算好排列長度，盡最大可能采集有效數據(圖7)。

4.3 運用Google Earth軟件進行航跡記錄

物探野外施工過程中，難免遇到陡坎、公路、河流等無法施工地段，通過手持GPS的航跡記錄，還可以給這些位置定位，對數據處理有極大幫助(圖8)。

航跡記錄不僅能記錄施工經過的位置，也保留了對應位置上的時間，有助于回憶施工的細節。運用此項功能還可以為地質調繪提供很好的幫助。

5 結論

(1)通過MapSource軟件可以快速計算GPS定位的空間平移參量DX、DY、DZ 3個參數，為前期工作節省大量時間。除了測量部門提供的三角點的坐標之外，還可以通過轉換線位圖，設置相對坐標系統來布置測線。

(2)通過Google Earth軟件自帶坐標轉換功能，可以將相對坐標系統下布置的測線轉到Google Earth軟件中查看。

圖4 MapSource下物探航點編輯Fig.4 Edition of geophysical navigation points in MapSource

圖5 MapSource下物探航線編輯Fig.5 Edition of geophysical navigation lines in MapSource

圖6 Google Earth軟件下物探航線Fig.6 Geophysical exploration routes in the software Google Earth

圖7 Google Earth軟件下物探航線的地形Fig.7 Terrain of geophysical exploration routes in the software Google Earth

圖8 Google Earth軟件下野外物探航跡Fig.8 Field geophysical exploration paths in the software Google Earth

(3)結合Google Earth軟件，MapSource軟件不僅可以對物探前期測線清理、后期數據處理提供幫助，還可以運用到鉆孔位置初定、地質調繪點初定等初步勘察中。

杜大彬，張寬房，張開盾，等.2007.手持GPS坐標系轉換方法［J］.陜西地質，25(1)：96－102.

GARMIN Corp.2004.MapSource使用說明書［R］.北京：Garmin Corp.

GARMIN Corp.2006.eTrex Venture－奇遇 GPS使用手冊［R］.北京：Garmin Corp.

李鵬，楊婷婷，王秋菊，等.2008.手持GPS的功能及其在野外調查中的應用［J］.內蒙古農業科技，(5)：98－99.

劉大杰，施一民，過靜珺.1997.全球定位系統(GPS)的原理與數據處理［M］.上海：同濟大學出版社.

劉凡珍，郭玉軍，孫萍，等.2007.GPS與 Excel、MapGIS相結合在化探工作中的應用［J］.吉林地質，26(1)：61－65.

馬翠玲.2010.GPS手持機坐標系統轉換參數的一種求法［J］.中國電子商務，(6)：72－73.

王會峰，彭立華，安興，等.2008.GPS技術在化探工作中的應用［J］.物探與化探，32(5)：477－479.