CASS 8.0在繪制地質(zhì)工作部署圖坐標(biāo)轉(zhuǎn)換中的應(yīng)用

徐 陽

(山東省煤田地質(zhì)局第一勘探隊，山東青島 266555)

0 引言

在多金屬礦勘探過程中，往往采用磁法勘探和激電中梯的方法進行工作區(qū)域的勘探(崔東鄆和王偉建，2009；管志寧，2015；王猛和梁中朋，2015)，或運用綜合物探方法分區(qū)域探測多金屬礦在測區(qū)的分布狀況(李榮亮等，2017)。多金屬礦往往分布在地理環(huán)境較為復(fù)雜的地域，因此尋找創(chuàng)新型的找礦方法成為了重要的研究對象(任為望等，2012；郭世強，2013；王曉軍等，2013；林斌，2018；邰文星等，2019；李俊杰等，2019)。而地質(zhì)勘探往往處在人煙稀少的地段，環(huán)境比較惡劣，嚴(yán)重威脅到野外工作人員的生命安全，因此需要尋找合適的方法以大大降低野外工作環(huán)境的危險程度(苑立清和趙金凱，2009)。因而，施工前的工作部署圖的合理規(guī)劃以及坐標(biāo)點的布置成為了重中之重(曾群望，1985；劉勵慎和汪津華，1987；臧忠淑和蔣樹一，2003；艾斯卡爾和李華，2010；王寶江，2010；張建恩等，2010)，而將測繪遙感等技術(shù)合理應(yīng)用到地質(zhì)勘探工作部署圖的規(guī)劃中可以大大提高工作效率(包世泰等，2004；康世英和張宏偉，2006；胡志文等，2012；張維寬等，2012；何輝等，2014；徐力峰，2014；狄廣禮和劉劍英，2015；周朝憲等，2015；李飛等，2018；王浩峰等，2019)。在眾多測繪軟件中，CASS具有非常強大的地質(zhì)圖編輯功能，及坐標(biāo)點導(dǎo)入導(dǎo)出的功能，因此運用CASS軟件對多金屬礦的工作部署圖進行轉(zhuǎn)換具有十分強大的優(yōu)點(李得基等，2006；韋晶珠，2010；李丙法，2017)。

結(jié)合以往的研究發(fā)現(xiàn)，本文首次提出運用測繪程序CASS 8.0來直接繪制地質(zhì)勘查類的工作部署圖，并利用測繪學(xué)科軟件與地質(zhì)學(xué)科知識的交叉應(yīng)用優(yōu)勢，使數(shù)據(jù)導(dǎo)入導(dǎo)出操作簡便化，易于掌握，并且繪圖精確度提高，使點坐標(biāo)直接呈現(xiàn)在圖形上，便于后期修改調(diào)整。

針對大面積比例尺坐標(biāo)點導(dǎo)出導(dǎo)入的困難問題，本文嘗試直接運用測繪軟件CASS 8.0中對高程點的繪制和輸出，直接轉(zhuǎn)換為帶有測點線號的點坐標(biāo)。該點坐標(biāo)無須經(jīng)過再次變換，經(jīng)Excel調(diào)整變換即可按線號批量導(dǎo)入手持GPS定位儀中。

CASS 8.0在地貌、地物繪制，點坐標(biāo)的繪制、提取等方面具有明顯的優(yōu)勢。為了有效解決坐標(biāo)點批量生成及導(dǎo)入問題，本文首次提出運用測繪軟件CASS 8.0的點陣排列工具及有編碼高程點的批量輸出的功能，輸出帶有線號以及序號的點坐標(biāo)。此方法操作簡單快捷，無需通過編程或應(yīng)用Excel表格方程式，即可得出直觀并且方便測繪儀器使用的數(shù)據(jù)及工作部署圖，可極大提高坐標(biāo)點提取的可操作性及簡便性。

1 技術(shù)方法

1.1 測線坐標(biāo)點的生成

在勘探野外工作中，采用手持GPS定位儀所需要的測點坐標(biāo)包括四個部分：橫坐標(biāo)、縱坐標(biāo)、線號、點號。文件格式為(李榮亮等，2017)：

橫坐標(biāo) 縱坐標(biāo) 線號 點號

4654445.415 48860.885 100 001

4654485.957 48860.609 100 002

4654524.640 48861.241 100 003

……

首先使用CASS 8.0繪制出需要的測點布陣圖，打開CASS軟件，點擊“繪圖處理”中的“展野外測點點號”，利用GPS定位儀測出的工區(qū)邊界點坐標(biāo)繪制出工區(qū)邊界。運用CASS中的對象捕捉功能，將邊界點連成線，繪制首條測線，并運用陣列復(fù)制功能繪制工區(qū)的所有測線，在每條測線端點處輸入測線號，輸入命令g，捕捉每條測線與邊界點的交點，并輸入相應(yīng)的測線號來作為高程點高程，選中所繪制高程點，點擊“陣列復(fù)制”輸入角度、間距、數(shù)目等信息，刪除工區(qū)邊界外高程點。繪制一條閉合的邊界線將所有高程點包含在內(nèi)，點擊“工程應(yīng)用”——“高程點生成數(shù)據(jù)文件”——“有編碼高程點”，此時操作平臺上會自動彈出“輸入數(shù)據(jù)文件名”的對話框，選擇相應(yīng)文件并點擊“打開”按鈕，輸入數(shù)據(jù)文件名后出現(xiàn)請選擇(1)選取高程點范圍(2)直接選取高程點或控制點，輸入1，點擊回車按鈕，此時命令欄出現(xiàn)“請選取建模區(qū)域邊界”的命令，點擊繪制的閉合邊界，即可生成所需坐標(biāo)文件(*.dat格式)(林斌，2018)。文件格式為：

測點序號 縱坐標(biāo) 橫坐標(biāo) 高程點

1， ，48860.885 4654445.415 1

2， ，48860.609 4654485.957 1

3， ，48861.241 4654524.640 1

……

對于CASS 8.0導(dǎo)出的坐標(biāo)文件，利用Microsoft Excel對數(shù)據(jù)進行處理，即可得到野外勘探工作的坐標(biāo)文件。具體步驟如下：①打開生成的測點坐標(biāo)；②對數(shù)據(jù)進行分裂；③按高程點的數(shù)值對整個表格進行重新排序，并改為測線序號；④將測點序號調(diào)整到最后一列，并重新按照測線的線號更改對應(yīng)的序號；將縱坐標(biāo)列與橫坐標(biāo)列互換。具體操作流程圖如圖1所示：

圖1 CASS 8.0生成坐標(biāo)流程Fig.1 Flowchart of CASS 8.0 coordinate generation

1.2 工作部署圖的繪制

將生成的坐標(biāo)文件按照以下格式重新修改并保存，文件格式為：

測點序號 縱坐標(biāo) 橫坐標(biāo) 高程點

1， ，48860.885 4654445.415 1

2， ，48860.609 4654485.957 1

3， ，48861.241 4654524.640 1

……

具體步驟如下：①利用Microsoft Excel打開生成的測點坐標(biāo)；②將測點序號調(diào)到第一列，并在后面插入新的一列；③將橫坐標(biāo)與縱坐標(biāo)對調(diào)；④將文件保存為*.csv格式；⑤將*.csv格式的文件修改為*.dat格式文件。

在CASS 8.0中打開相應(yīng)的地形地質(zhì)圖，點擊“繪圖處理”——“展野外測點點號”，此時操作平臺上會自動彈出“輸入坐標(biāo)數(shù)據(jù)文件名”的對話框，選擇以上生成的*.dat格式文件，并點擊“打開”，即可將野外測點點號展到圖上。繼續(xù)點擊“繪圖處理”——“展高程點”，此時操作平臺上會自動彈出“輸入坐標(biāo)數(shù)據(jù)文件名”的對話框，再次選擇以上生成的*.dat格式文件并點擊“打開”，即可將測線的線號展到地形地質(zhì)圖上。具體操作流程圖如圖2所示。

圖2 CASS 8.0繪制工作部署圖流程Fig.2 Flowchart of drawing work deployment map by CASS 8.0

2 應(yīng)用實例

2.1 1∶1萬高精度磁測工作中坐標(biāo)點生成及工作部署

在內(nèi)蒙古自治區(qū)額濟納旗黑鷹山南區(qū)鉛多金屬礦普查項目中，勘查區(qū)位于內(nèi)蒙古自治區(qū)西部額濟納旗的達來呼布鎮(zhèn)西偏北約200 km處，地理極值坐標(biāo)為：東經(jīng)98°51′56″—98°59′36″，北緯42°01′31″—42°04′01″，勘查區(qū)面積20.14 km2。

勘查過程中，首先在該區(qū)進行了1∶1萬地質(zhì)草測，1∶1萬高精度磁法測量，然后選擇位于F1、F2斷層附近的成礦有利部位，進行了1∶1萬激電測量、槽探、激電測深工作。

(1)GPS衛(wèi)星定位儀參數(shù)的設(shè)置

磁測掃面測網(wǎng)以本地區(qū)國家三角點為起算點，施測GPS(E級)控制網(wǎng)作為該工區(qū)的首級控制。

測網(wǎng)基線在首級控制網(wǎng)的基礎(chǔ)上發(fā)展，采用動態(tài)GPS(RTK)進行，其相應(yīng)的參數(shù)設(shè)置如下：

坐標(biāo)系統(tǒng)：選擇西安80坐標(biāo)系，高程采用1985高程系。

(2)測網(wǎng)布置

測區(qū)內(nèi)測線按照0°方位角布置。

點、線號的編排方法原則是：測線號與測點號從南到、北從西向東由小至大。

(3)測量網(wǎng)度(100×40)

每條測線的測點插衛(wèi)生筷子標(biāo)記，并用記號筆清楚寫明點線號，以便后續(xù)工作時能在野外現(xiàn)場識別。測地工作依據(jù)執(zhí)行GB/T18314-2009《全球定位系統(tǒng)(GPS)測量規(guī)范》、GB/T18314-2001《地質(zhì)礦產(chǎn)勘查測量規(guī)范》、DZ/T0153-95《物化探工程測量規(guī)范》、《物化探測量規(guī)范》，嚴(yán)格遵守有關(guān)規(guī)范和設(shè)計要求。對外業(yè)測量采用不定時的自檢、互檢、確保測量質(zhì)量。為保證點位相對誤差小于5 m，工作中衛(wèi)星數(shù)量不少于5顆。

普查區(qū)高精度磁測測線方位角0°點、線號的編排方法原則是：測線號與測點號從南到北、從西向東由小至大。測網(wǎng)網(wǎng)度100 m×40 m，每條測線的測點插衛(wèi)生筷子標(biāo)記，并用記號筆清楚寫明點線號，以便后續(xù)工作時能在野外現(xiàn)場識別。該工區(qū)邊界4個點的坐標(biāo)分別為：

4657687.170，488859.530

4657685.430，499360.740

4654445.420，488860.890

4654445.280，499361.040

按照1.1的方法得到的前10個點的坐標(biāo)如下：

橫坐標(biāo) 縱坐標(biāo) 線號 點號

4654445.42, 488860.89, 1, 1

4654485.96, 488860.61, 1, 2

4654524.64, 488861.24, 1, 3

4654565.46, 488860.98, 1, 4

4654605.28, 488860.93, 1, 5

4654645.43, 488861.03, 1, 6

4654685.39, 488861.1, 1, 7

4654725.23, 488861.01, 1, 8

4654765.97, 488861.09, 1, 9

4654805.61, 488861.25, 1, 10

將該方法產(chǎn)生的坐標(biāo)文件生成*.csv格式后導(dǎo)入GPS高精度手持儀器(合眾思壯MG838)即可進行野外定點。該方法的優(yōu)點是能利用CASS 8.0的專業(yè)性繪圖工具直觀、精準(zhǔn)、簡潔地繪制工區(qū)坐標(biāo)點的位置。在勘探工區(qū)中若存在不可逾越的障礙物時，如河流、湖泊、陡崖等，亦可根據(jù)地形地質(zhì)圖及時調(diào)整測點位置，合理規(guī)劃路線，提高工作效率，降低野外工作風(fēng)險，迅速準(zhǔn)確輸入輸出測點坐標(biāo)。

將以上文件轉(zhuǎn)變?yōu)镃ASS 8.0使用的*.dat文件如圖3所示。

圖3 1∶1萬磁測測點轉(zhuǎn)換后的dat文件格式Fig.3 Dat file of 1∶10000 magnetic measurement points after transformation

將以上文件展到CASS打開的地形地質(zhì)圖中，調(diào)整坐標(biāo)點的格式、大小等特性，便可得到相應(yīng)的工作部署圖，如圖4所示。本次應(yīng)用共生成測線106條，測點坐標(biāo)8693個，使用天寶R4 GPS差分定位系統(tǒng)復(fù)測，精確度達到mm級別。

2.2 1∶1萬高精度激電中梯工作中坐標(biāo)點生成及工作部署

普查區(qū)激電區(qū)測線方位角45°，測網(wǎng)網(wǎng)度100 m×20 m，線號、點號均以1 m為單位編號，線號由南向北、點號由西向東增大。每條測線的測點插衛(wèi)生筷子標(biāo)記，并用記號筆清楚寫明點線號，以便后續(xù)工作時能在野外現(xiàn)場識別。該工區(qū)邊界6個點的坐標(biāo)分別為：

圖4 黑鷹山南區(qū)鉛多金屬礦普查1∶1萬磁測測點坐標(biāo)投影到地圖上效果Fig.4 Map projection of coordinates of 1∶10000 magnetic survey points for lead polymetallic ore survey in the south of Heiyingshan district

4656292.976，490727.507

4656294.945，494406.776

4657708.181，494406.058

4657708.181，494406.058

4654455.492，495071.401

4654455.191，490729.890

按照1.1的方法得到的前10個點的坐標(biāo)如下：

橫坐標(biāo)， 縱坐標(biāo)， 線號， 點號

4655799.580， 490730.022, 1300.000， 1

4655813.722， 490744.164, 1300.000， 2

4655827.864， 490758.306， 1300.000， 3

4655842.006， 490772.448， 1300.000， 4

4655856.148， 490786.590， 1300.000， 5

4655870.290， 490800.732， 1300.000， 6

4655884.433， 490814.875， 1300.000， 7

4655898.575， 490829.017， 1300.000， 8

4655912.717， 490843.159， 1300.000， 9

4655926.859， 490857.301， 1300.000， 10

將該方法產(chǎn)生的坐標(biāo)文件生成*.csv格式后導(dǎo)入GPS高精度手持儀器(合眾思壯MG838)即可進行野外定點，將以上文件轉(zhuǎn)變?yōu)镃ASS 8.0使用的*.dat文件,如圖5所示。

圖5 1∶1萬電測測點轉(zhuǎn)換后的dat文件格式Fig.5 dat file format of 1∶10000 electric measuring points after transformation

將以上文件展到CASS打開的地形地質(zhì)圖中，

調(diào)整坐標(biāo)點的格式、大小等特性，便可得到相應(yīng)的工作部署圖，如圖6所示。本次應(yīng)用共生成測線44條，測點坐標(biāo)4511個，使用天寶R4 GPS差分定位系統(tǒng)復(fù)測，精確度達到mm級別。

2.3 工作部署圖生成后測網(wǎng)敷設(shè)及質(zhì)量評述

2.3.1 測網(wǎng)敷設(shè)及質(zhì)量評述

物探網(wǎng)測線放樣全部采用科力達風(fēng)云K-9雙頻接收機進行實時定位放樣，投入一個基準(zhǔn)站、四個移動站。這樣既保證了工作效率的提高，又保證了測量精度。

動態(tài)RTK測量精度指標(biāo)：

平面精度：±2 cm+2 ppm

高程精度：±3 cm+1 ppm

儀器通過山東省計量科學(xué)研究院質(zhì)量檢測，其精度指標(biāo)符合作業(yè)要求，可以投入生產(chǎn)。

(1)作業(yè)方法

本次物探網(wǎng)的布設(shè)采用網(wǎng)絡(luò)RTK模式，即在工區(qū)內(nèi)只架設(shè)一個基準(zhǔn)站，通過GPRS網(wǎng)絡(luò)與通訊控制中心連接，流動站也通過GPRS將自己的初始化位置傳輸給控制中心，并接受中心發(fā)送的差分信號，生成厘米級的位置信息。網(wǎng)絡(luò)RTK的作業(yè)半徑可達20 km，能覆蓋整個工區(qū)。

圖6 黑鷹山南區(qū)鉛多金屬礦普查1∶1萬電測測點坐標(biāo)投影到地圖上效果Fig.6 Map projection of coordinates of 1∶10000 electric survey points for lead polymetallic ore survey in the south of Heiyingshan district

選取通視開闊，設(shè)置相對制高點作為基準(zhǔn)站(基站架設(shè)在未知點上)，在其上安置一臺接收機，連續(xù)跟蹤所有的可見衛(wèi)星，連接通訊控制中心，發(fā)送實時改正信息；另四臺有接收GPRS網(wǎng)絡(luò)信號的GPS接收機作為流動站，流動站的移動手簿通過在已知控制點上進行點校正，校正后便可在事先輸入好的測線上進行測點的布設(shè)。在測點位置，壓有寫有點線號的紅紙條以作標(biāo)記，并且每到達一個點都在電子手薄中記錄下當(dāng)前的坐標(biāo)位置。

(2)RTK作業(yè)要求

①至少同步觀測4顆以上分布良好的衛(wèi)星，并在運動過程中保持連續(xù)的跟蹤。

②在實時動態(tài)狀態(tài)下顯示移動站的三維坐標(biāo)，以及其相應(yīng)的坐標(biāo)分量的中誤差，通過設(shè)置最大點位中誤差實現(xiàn)對測量精度的限制，符合要求的記錄，否則拒絕記錄。

2.3.2 高精度磁測工作及其質(zhì)量評述

普查區(qū)高精度磁測測線方位角0°點、線號的編排方法原則是：測線號與測點號從南到北、從西向東由小至大。測網(wǎng)網(wǎng)度100 m×40 m，每條測線的測點插衛(wèi)生筷子標(biāo)記，并用記號筆清楚寫明點線號，以便后續(xù)工作時能在野外現(xiàn)場識別。

磁法測量使用5臺GSM19T型質(zhì)子磁力儀。在開工前和工作結(jié)束后對儀器性儀器一致性進行校驗。

測點使用GSM19T型質(zhì)子磁力儀作總場觀測，每次觀測時探頭的高度保持一致，日變和野外磁測儀器時間調(diào)成一致。觀測人員嚴(yán)格“去磁”，未攜帶任何鐵磁性物品。在生產(chǎn)期間，避開了干擾。通過每天對日變數(shù)據(jù)的分析，未遇到磁暴。變觀測點位于磁測總基點上。

日變觀測使用了性能較穩(wěn)定的GSM19T質(zhì)子磁力儀，進行日變觀測，每20秒觀測一次。每天的日變觀測早于各儀器的早校正點觀測之前開始，遲于晚校正點觀測之后結(jié)束。

當(dāng)天原始記錄傳入計算機之后進行100%自檢驗收，并填寫磁法資料日驗收記錄表,確保記錄數(shù)據(jù)完整齊全，合格后保存原始觀測文件。內(nèi)業(yè)計算主要包括儀器一致性精度計算、測點ΔT值計算、質(zhì)量檢查精度計算等。

質(zhì)量檢查采用“一同三不同”(同點位、不同儀器、不同日期、不同操作員)的質(zhì)量檢查原則進行。本區(qū)1∶10000高精度磁測工作共完成34.55 km2，全區(qū)掃面共完成物理點8693個，質(zhì)檢點490個，檢查率5.6%，野外觀測精度：±1.3 nT。野外觀測質(zhì)量滿足《地面高精度磁測勘查技術(shù)規(guī)程》和設(shè)計的要求，質(zhì)量合格。

ΔDi—第i點觀測值與放樣坐標(biāo)的點位差；n—觀測總點數(shù)，i=1,2,3,……n。

平面點位中誤差：m位(剖)=±0.57 m;

m允(剖)=±6.25 m

高程點位中誤差：m高(剖)=±0.41 m;

m允(剖)=±5 m

從檢查結(jié)果來看，本區(qū)的測量工作精度可靠，精度指標(biāo)滿足設(shè)計及規(guī)范要求，為野外物探數(shù)據(jù)的采集提供了可靠的保障。

2.3.3 激電中梯工作及其質(zhì)量評述

普查區(qū)激電區(qū)測線方位角45°，測網(wǎng)網(wǎng)度100 m×20 m，線號、點號均以1 m為單位編號，線號由南向北、點號由西向東增大。

使用儀器為重慶地質(zhì)儀器廠生產(chǎn)的DFZ10-2型大功率智能發(fā)射系統(tǒng)和與之配套的DJS-8型數(shù)字直流激電接收機，供電電源選擇為本田發(fā)電機ET12000交流發(fā)電機，配以直流穩(wěn)壓裝置。

在野外施工過程中嚴(yán)格按照《時間域激發(fā)極化法技術(shù)規(guī)程(DZ/T0070-93)》、(DZ/T0153-2014)《物化探工程測量規(guī)范》執(zhí)行。

激電異常是明顯高于背景值的ηa變化。為減少背景值本身的波動和觀測誤差給異常圈定帶來的影響，一般通過異常下限來劃定異常，異常下限取決于背景值的大小、穩(wěn)定性和ηa的觀測精度。在激電中梯所測量的全部極化率數(shù)據(jù)中，剔除畸變、突變數(shù)據(jù)后去掉大于3.0%的數(shù)據(jù)，將剩余數(shù)據(jù)取算術(shù)平均值作為該測區(qū)極化率背景值。

經(jīng)計算本區(qū)極化率背景值為1.45%。根據(jù)公式ηa下限=(1+3Mηa)×計算異常下限，其中Mηa=3.34%，ηa背景=1.45%，計算得ηa下限=1.60%。

在電腦上對測區(qū)所有測線起始點坐標(biāo)和重復(fù)觀測點坐標(biāo)等資料數(shù)據(jù)進行統(tǒng)一整理記錄。由技術(shù)人員進行100%自檢及互檢驗收，確保記錄數(shù)據(jù)完整齊全。對已檢查野外數(shù)據(jù)進行室內(nèi)綜合處理，保證每個數(shù)據(jù)采集對畸變點進行分析并剔除，保證處理結(jié)果的高保真度和電性分布特征的客觀真實性。結(jié)合數(shù)據(jù)整理激電中梯視電阻率和視極化率等值線平面圖和激電測深成果圖件，根據(jù)所得資料對工作區(qū)電性特征進行分析解釋，并編寫物探成果報告點與實際坐標(biāo)相對應(yīng)，計算視電阻率數(shù)值，根據(jù)實際數(shù)據(jù)采集點位整理物探工作實際材料圖。

本區(qū)1∶10000激電中梯工作共完成9.00 km2，共計4644個有效物理點。質(zhì)量檢查按照“同點位、不同接收機、不同操作員、不同電流”的“一同三不同”原則進行，共計檢查150個點，質(zhì)檢比例為3.1%(設(shè)計及規(guī)范要求不低于3.0%)。經(jīng)統(tǒng)計計算，質(zhì)量檢查視極化率實達均方相對誤差為±3.34%(設(shè)計要求低于4.0%)；視電阻率實達均方相對誤差為±4.40%(設(shè)計要求低于7.0%)，表明質(zhì)量滿足設(shè)計及規(guī)范要求，整體工作完成情況良好。

3 結(jié)論

基于對工作中測線點坐標(biāo)批量提取及輸入的需求，本文提出了一種運用CASS 8.0更加簡便得出工作部署圖以及測線點坐標(biāo)的方法，并批量輸入到GPS手持定位儀中。對比結(jié)合VB6.0語言或Visual C++編程等得出坐標(biāo)點數(shù)據(jù)的方法，此操作方式更加簡單、直觀。通過直接繪制點和提取點坐標(biāo)，無需坐標(biāo)轉(zhuǎn)換即可獲得能夠批量導(dǎo)入GPS手持儀中的數(shù)據(jù)。同時，在繪制點陣的過程中，運用CASS 8.0帶來的優(yōu)勢，可以在施工前結(jié)合地形地質(zhì)圖合理刪除一些無效點，甚至可以繪制曲線型的測線來規(guī)避危險區(qū)域或者障礙物(懸崖、峭壁等)，以提高工作效率，保障工作人員的生命安全。將這些技術(shù)運用到內(nèi)蒙古自治區(qū)額濟納旗黑鷹山南區(qū)鉛多金屬礦普查工作中，得到了較好的效果。

[附中文參考文獻]

艾斯卡爾，李華.2010.關(guān)于磁法勘探測點坐標(biāo)的生成與工作部署圖繪制[J].物探與化探，34(2):250-252.

包世泰，夏斌，蔣鵬，付文生.2004.基于GIS的地質(zhì)勘察信息系統(tǒng)設(shè)計與實現(xiàn)[J].地理與地理信息科學(xué)，20(4):31-35.

崔東鄆，王偉建.2009.激電中梯測量工作在多金屬找礦中的應(yīng)用[J].貴州地質(zhì)，26(3):218-220.

狄廣禮，劉劍英.2015.基于GPS-RTK技術(shù)在地質(zhì)勘探工程測量中的應(yīng)用研究[J].科技資訊，13(5):51-52.

管志寧.2005.地磁場與磁力勘探[M].北京:地質(zhì)出版社:1-309.

郭世強.2013.淺談在地質(zhì)找礦中地質(zhì)勘探的應(yīng)用[J].科技創(chuàng)業(yè)家，(6):198,200.

何輝，薛典軍，王林飛.2014.地球物理數(shù)據(jù)處理解釋系統(tǒng)中坐標(biāo)轉(zhuǎn)換的設(shè)計及實現(xiàn)[J].物探與化探，38(5):1060-1063.

胡志文，歐陽燕，羅湘.2012.水工環(huán)地質(zhì)勘察及遙感技術(shù)在地質(zhì)工作中的應(yīng)用[J].江西建材，(5):187-188.

康世英，張宏偉.2006.GPS測量在石油物探中的應(yīng)用探討[J].礦產(chǎn)與地質(zhì)，20(4-5):552-555.

李榮亮，田建榮，劉洋，任杰，白順寶，劉曉峰.2017.綜合物探方法在甘肅梧桐井鐵銅多金屬礦勘查中的應(yīng)用[J].地質(zhì)與勘探，53(4):755-764.

李丙法.2017.CASS 8.0應(yīng)用于大比例尺地形圖繪制及使用技巧[J].工程技術(shù)，(3):307.

李得基，梅貴福，田淑梅.2006.“南方測繪軟件CASS6.1”在工程測量中的應(yīng)用[J].青海國土經(jīng)略，(2):41-43.

李飛，丁志強，崔志強，胥值禮，李軍峰.2018.CH-3無人機航磁測量系統(tǒng)在我國新疆不同地形區(qū)的應(yīng)用示范[J].地質(zhì)與勘探,54(4):735-746.

李俊杰，張紅綱，王偉，何建設(shè)，李劍強，郭佳豪.2019.綜合物探技術(shù)在灰?guī)r地區(qū)隧洞超前預(yù)報中的應(yīng)用[J].地質(zhì)與勘探,55(6):1452-1462.

林斌.2018.關(guān)于找礦地質(zhì)勘探技術(shù)創(chuàng)新的思考[J].世界有色金屬，(3):84-85.

劉勵慎,汪津華.1987.PC-1500微型計算機處理鉆孔在空間坐標(biāo)位置的實用方法[J].地質(zhì)與勘探23(12):41-45.

任為望，宋子萱，耿博.2012.關(guān)于找礦地質(zhì)勘探技術(shù)創(chuàng)新的探討[J].科學(xué)技術(shù)創(chuàng)新，(12):17.

邰文星，金尚剛，何彥南，楊梅珍，劉萬龍.2019.綜合物探方法在湖北大冶雞冠咀銅金礦深部及外圍找礦中的應(yīng)用[J].地質(zhì)與勘探,55(4):1010-1025.

王寶江.2010.VB6.0物探線平面坐標(biāo)推算程序設(shè)計[J].山西建筑，36(9):358-359.

王浩鋒，劉波，陳霜，陳鳳興，薛文浩.2019.綜合物探測量在二連盆地沉積巖型鈾礦勘查中應(yīng)用[J].地質(zhì)與勘探,55(1):127-134.

王猛，梁中朋.2015.大功率激電在內(nèi)蒙古多金屬礦勘查的應(yīng)用[J].工程地球物理學(xué)報，12(2):176-182.

王曉軍，陳孝強，傅揚.2013.草原覆蓋區(qū)綜合物探找礦工作方法—以內(nèi)蒙古阿吉勒銀鉛鋅多金屬礦區(qū)為例[J].高校地質(zhì)學(xué)報，19(Z):376.

韋晶珠.2010.淺談如何利用CASS軟件繪制地形圖和斷面圖[J].現(xiàn)代企業(yè)文化，(8):119-120.

徐力峰.2014.物化探測點批量導(dǎo)入手持GPS的方法[J].地球，(1):141-142.

苑立清，趙金凱.2009.地質(zhì)勘探安全生產(chǎn)工作的思考與對策[J].安全與環(huán)境工程，26(3):87-91.

臧忠淑,蔣樹一.2003.ArcView GIS中圖件比例系數(shù)及坐標(biāo)位置的設(shè)置技巧[J].地質(zhì)與勘探,39(1):58-61.

曾群望.1985.繪制幾種地質(zhì)坐標(biāo)圖的BASIC程序―適用于SHARP PC-1500袖珍機[J].地質(zhì)與勘探，21(9):32-35.

張建恩，劉亞寧，李康虎，賈傳杰.2010.RTK物探測量放樣坐標(biāo)轉(zhuǎn)換參數(shù)求取方法的探討[J].物探裝備，20(4):264-268.

張維寬，張永勝，孟紅.2012.全站儀極坐標(biāo)法測設(shè)物探網(wǎng)的一種新方法[J].價值工程，(11):36-37.

周朝憲，董少波，和志軍，姜高珍，燕丹晨.2015.自定義坐標(biāo)系的建立及其坐標(biāo)變換實現(xiàn)[J].地質(zhì)與勘探,51(4):699-704.