馬春生，趙華偉，臧興運(yùn)，朱小鋒

吉林省地質(zhì)調(diào)查院，吉林 長(zhǎng)春 130102

0 引言

數(shù)字地質(zhì)調(diào)查系統(tǒng)DGSS(Digital Geological Survey System)是中國(guó)地質(zhì)調(diào)查局在原數(shù)字填圖野外數(shù)據(jù)采集系統(tǒng)、數(shù)字剖面系統(tǒng)、固體礦產(chǎn)野外數(shù)據(jù)采集系統(tǒng)、礦產(chǎn)資源調(diào)查數(shù)據(jù)處理與綜合分析子系統(tǒng)、資源量估算系統(tǒng)和礦體三維顯示系統(tǒng)等6大系統(tǒng)集成為一體化的數(shù)字地質(zhì)調(diào)查系統(tǒng)軟件（2010年）。該系統(tǒng)自面世以來(lái)，在區(qū)域地質(zhì)調(diào)查、固體礦產(chǎn)勘查、礦體模擬、品位估計(jì)、資源儲(chǔ)量估算、礦山開(kāi)采系統(tǒng)優(yōu)化等方面得到了廣泛的應(yīng)用，深受廣大地質(zhì)工作者的喜愛(ài)。

利用DGSS系統(tǒng)，可以方便地對(duì)地球化學(xué)數(shù)據(jù)進(jìn)行數(shù)理統(tǒng)計(jì)、相關(guān)性計(jì)算等數(shù)據(jù)分析以及各類圖件的制圖，與專業(yè)的地球化學(xué)處理軟件相比，利用DGSS系統(tǒng)做圖最大的優(yōu)點(diǎn)是成圖既符合數(shù)據(jù)庫(kù)的要求，無(wú)需各種轉(zhuǎn)換即可以直接參與數(shù)據(jù)庫(kù)的整理工作[1]。

但是在利用DGSS系統(tǒng)進(jìn)行地球化學(xué)有關(guān)圖件制作時(shí)，由于采樣布點(diǎn)等原因，也會(huì)遇到成圖范圍與實(shí)際圖件范圍不一致的情況[2]。本文以XXXX幅1：5萬(wàn)水系沉積物測(cè)量地球化學(xué)圖為例，闡述在利用DGSS系統(tǒng)做圖時(shí)，配合擴(kuò)邊數(shù)據(jù)及載剪框的使用，有效地解決了這一問(wèn)題。

1 準(zhǔn)備工作

是正確成圖的基礎(chǔ)，包括數(shù)據(jù)準(zhǔn)備和圖件準(zhǔn)備。

1.1 數(shù)據(jù)準(zhǔn)備

打開(kāi)DGSS系統(tǒng)XXXX幅，點(diǎn)擊野外總圖庫(kù)→勾選左側(cè)“GEOWATER.WT”文件→編輯欄“工具”→圖層屬性導(dǎo)出EXCEL，導(dǎo)出圖幅內(nèi)采樣點(diǎn)屬性。

在采樣點(diǎn)屬性的EXCEL表格里，保留“ID”號(hào)、采樣點(diǎn)坐標(biāo)“X”、“Y”，樣品編號(hào)“PNT_CODE”，其余“采樣地區(qū)”等屬性刪掉，代之以各元素分析結(jié)果。需要注意的是坐標(biāo)“X”不能有高斯投影帶號(hào)[3，4]，分析結(jié)果直接標(biāo)注各元素既可，含量單位Au元素默認(rèn)為“10-9”，其余元素為“10-6”。最后整理結(jié)果如表1所示，當(dāng)然，實(shí)際樣品數(shù)及分析元素并不限于表1的數(shù)量。

將整理好的表1中數(shù)據(jù)由EXCEL表轉(zhuǎn)出TXT文本文件，保存?zhèn)溆谩?/p>

1.2 圖件準(zhǔn)備

包括樣品分布圖的制做及裁剪框的制作。

首先利用MAPGIS中的實(shí)用服務(wù)→投影變換功能，將上述TXT文本文件轉(zhuǎn)換為點(diǎn)文件DATA_GEOCHEMICAL_STREAM .WT，加入到該圖幅的PRB庫(kù)工程中備用。

利用圖幅的內(nèi)圖框制作載剪框。為防止圖形參數(shù)丟失及線屬性改變，建議將A_FRAME.WL文件復(fù)制，在復(fù)制文件中將除內(nèi)圖框的線全部刪除，執(zhí)行自動(dòng)剪斷線功能，刪除四角處的線頭，四條內(nèi)框直線聯(lián)接成一個(gè)載剪框，將文件重新命名為載剪框.WL并加入到該圖幅的圖幅PRB庫(kù)中備用（該文件為臨時(shí)文件，用后可刪除）。

表1 XXXX幅1：5萬(wàn)水系沉積物測(cè)量數(shù)據(jù)表Table 1 Measurement data of 1:50 000 stream sediment

打開(kāi)DGSS系統(tǒng)，切換到XXXX幅實(shí)際材料圖界面，在左側(cè)文件欄中插入兩個(gè)上述文件“載剪 框 .WL”、“DATA_GEOCHEMICAL_STREAM.WT”并使之處于可編輯狀態(tài)，其余文件關(guān)閉，則得到XXXX幅水系樣品分布圖，如圖1所示（比例尺為1：5萬(wàn)，圖中略）。

2 地球化學(xué)圖（等值線圖）的制作

在實(shí)際材料圖界面勾選點(diǎn)文件“DATA_GEOCHEMICAL_STREAM .WT”，點(diǎn)擊上面編輯欄綜合數(shù)據(jù)處理→等值線圖，在彈出的對(duì)話框中，選擇點(diǎn)圖層一欄選取“DATA_GEOCHEMICAL_STREAM .WT”文件，字段名稱選取“Cu”元素，其余按提示要求即可做出地球化學(xué)圖（等值線圖），這里以Cu元素為例。本例網(wǎng)格化方法選用GRID方法，橫向網(wǎng)格數(shù)、縱向網(wǎng)格數(shù)均用50，可根據(jù)做為參考標(biāo)尺，就近將復(fù)制點(diǎn)移到載剪框端點(diǎn)之外，保存。

圖1 水系樣品分布圖Fig.1 Distribution of water system samples

擴(kuò)邊數(shù)據(jù)的選取，一般來(lái)說(shuō)選取最靠近載剪框（內(nèi)圖框）的點(diǎn)，參考邊部異常線的形態(tài)及趨勢(shì)。當(dāng)采樣點(diǎn)較多時(shí)（例如一個(gè)整5萬(wàn)圖幅水系采樣），每邊選取一個(gè)即可，不會(huì)對(duì)異常形態(tài)造成明顯影響；當(dāng)采樣點(diǎn)較少時(shí)，需根據(jù)圖面情況選取有代表性的點(diǎn)，選取不當(dāng)由于邊界效應(yīng)容易使異常線變形。本例選取5個(gè)，如圖3所示。

點(diǎn)擊上面編輯欄綜合數(shù)據(jù)處理→等值線圖，重做等值線過(guò)程，得到新的Cu元素地球化學(xué)圖，如圖3所示。與圖2相比，可見(jiàn)成圖范圍擴(kuò)大了，囊括整個(gè)實(shí)際圖件范圍。

由于本次成圖過(guò)程使用了邊框處部的擴(kuò)邊數(shù)實(shí)際情況加密或減稀，制圖幅面選取原始數(shù)據(jù)范圍。本例略去直方圖、色階圖例、極值點(diǎn)標(biāo)注等。

如圖2所示XXXX幅Cu元素地球化學(xué)圖，實(shí)際圖件范圍為一略向右傾斜的標(biāo)準(zhǔn)5萬(wàn)圖幅（傾斜方向與傾斜程度與該圖幅所處中央子午線位置有關(guān)），而成圖范圍則與采樣點(diǎn)分布范圍有關(guān)，為一規(guī)則的矩形。實(shí)際圖件范圍與成圖范圍在圖邊處不相吻合，有的地方超出，有的地方欠缺，圖面的完整性較差。

3 擴(kuò)邊數(shù)據(jù)與載剪框的應(yīng)用

回到實(shí)際材料圖界面，勾選點(diǎn)文件“DATA_GEOCHEMICAL_STREAM .WT”，將鼠標(biāo)設(shè)置為大十字，選取擴(kuò)邊數(shù)據(jù)進(jìn)行復(fù)制，利用鼠標(biāo)大十字據(jù)，因此成圖后需要以標(biāo)準(zhǔn)圖框?yàn)檩d切范圍，對(duì)生成的圖件進(jìn)行載切，以消除因擴(kuò)邊數(shù)據(jù)導(dǎo)致圖面范圍過(guò)大的影響[2]。關(guān)閉DGSS系統(tǒng)，打開(kāi)MAPGIS或SECTION系統(tǒng)，利用載剪框?qū)D3載剪（內(nèi)載），得到完整的Cu元素地球化學(xué)圖，如圖4所示，可以發(fā)現(xiàn)，擴(kuò)邊數(shù)據(jù)同時(shí)也已被載掉。

圖2 Cu元素地球化學(xué)圖（1）Fig.2 Cu element geochemical map (1)

4 圖件的檢驗(yàn)

直接對(duì)比圖2、圖4，或在MAPGIS中打開(kāi)圖2，在左邊欄中插入圖4點(diǎn)線文件，可以發(fā)現(xiàn)在載剪框（內(nèi)圖框）內(nèi)，二圖的異常曲線形態(tài)、大小基本相同；采樣點(diǎn)完全重合。由此可見(jiàn)該方法形成的地球化學(xué)圖件，既能使圖面完整，又確保了其原始性、真實(shí)性。

圖3 Cu元素地球化學(xué)圖（2）Fig.3 Cu element geochemical map (2)

圖4 Cu元素地球化學(xué)圖（3）Fig.4 Cu element geochemical map (3)