PAGEPAGE17化探數(shù)據(jù)處理成圖的過程畢武1、2段新力1、2黃顯義1、2袁小龍1、2彭仲秋1、2李永華1、21.烏魯木齊金維圖文信息科技有限公司，新疆，烏魯木齊，8300912.新疆地礦局物化探大隊(duì)計(jì)算中心，新疆，昌吉，8311000前言原測(cè)數(shù)據(jù)水文地質(zhì)數(shù)據(jù)處理數(shù)據(jù)預(yù)處理重磁數(shù)據(jù)處理化探數(shù)據(jù)處理數(shù)學(xué)地質(zhì)數(shù)據(jù)處理點(diǎn)位數(shù)據(jù)圖打印輸出原測(cè)數(shù)據(jù)水文地質(zhì)數(shù)據(jù)處理數(shù)據(jù)預(yù)處理重磁數(shù)據(jù)處理化探數(shù)據(jù)處理數(shù)學(xué)地質(zhì)數(shù)據(jù)處理點(diǎn)位數(shù)據(jù)圖打印輸出數(shù)據(jù)檢查數(shù)據(jù)整理數(shù)據(jù)網(wǎng)格化投影變換格式轉(zhuǎn)換電法數(shù)據(jù)處理數(shù)據(jù)處理圖件繪制輸出數(shù)據(jù)輸入圖1金維地學(xué)信息處理研究應(yīng)用系統(tǒng)(GeoIPAS)圖件制作流程圖圖像輸出文本文件圖形文件剖面圖等值線圖場(chǎng)強(qiáng)圖三元圖三維立體圖點(diǎn)位符號(hào)圖……1處理步驟化探處理的成果包括：(1)參數(shù)統(tǒng)計(jì)表；(2)R型聚類分析-譜系圖；(3)重復(fù)樣三層套合方差分析或者重復(fù)樣合格率計(jì)算結(jié)果；(4)點(diǎn)位數(shù)據(jù)圖；(5)地球化學(xué)圖；(6)直方圖；(7)組合異常圖；(8)綜合異常圖；(9)遠(yuǎn)景區(qū)劃圖；(10)單元素異常參數(shù)統(tǒng)計(jì)(附表冊(cè))；(11)異常剖析(附圖冊(cè))；(12)綜合異常登記卡(附表冊(cè))。在GeoIPAS系統(tǒng)中，化探數(shù)據(jù)處理分為以下幾個(gè)主要步驟：1.1數(shù)據(jù)檢查數(shù)值檢查，坐標(biāo)檢查，重復(fù)樣坐標(biāo)檢查。1.2分析處理重復(fù)樣三層套合方差分析、重復(fù)樣合格率計(jì)算、化探特征參數(shù)統(tǒng)計(jì)、化探背景值分析、R型聚類分析、因子分析。1.3數(shù)據(jù)分析數(shù)據(jù)變換；襯值、累加襯值；數(shù)據(jù)累加、累乘、比值；異常歸一化。1.4網(wǎng)格化離散數(shù)據(jù)網(wǎng)格化、XYZ數(shù)據(jù)轉(zhuǎn)網(wǎng)格數(shù)據(jù)1.5成圖點(diǎn)位數(shù)據(jù)圖、彩色等量線圖、直方圖、組合異常圖、單元素異常圖、綜合異常圖、剖析圖。1.6單元素異常參數(shù)統(tǒng)計(jì)1.7綜合異常登記卡開始開始數(shù)據(jù)輸入數(shù)據(jù)預(yù)處理網(wǎng)格化成圖異常統(tǒng)計(jì)數(shù)據(jù)分析處理結(jié)束分析樣品與坐標(biāo)對(duì)應(yīng)數(shù)值檢查坐標(biāo)檢查重復(fù)樣提取采樣編碼統(tǒng)計(jì)三層套合方差分析或重復(fù)樣合格率計(jì)算化探特征參數(shù)統(tǒng)計(jì)R型聚類分析離散數(shù)據(jù)網(wǎng)格化XYZ數(shù)據(jù)轉(zhuǎn)網(wǎng)格數(shù)據(jù)點(diǎn)位數(shù)據(jù)圖直方圖地球化學(xué)圖組合異常圖單元素異常圖綜合異常圖剖析圖TXTGRD單元素異常參數(shù)統(tǒng)計(jì)圖2化探數(shù)據(jù)處理成圖流程TXTTXTTXTTXTTXT矢量GRD矢量TXT矢量CSV2具體處理過程2.1數(shù)據(jù)檢查我們的數(shù)據(jù)處理工作從化驗(yàn)室提供的樣品分析報(bào)告開始，項(xiàng)目要提供坐標(biāo)和樣品對(duì)應(yīng)的分析數(shù)據(jù)，坐標(biāo)我們一般取實(shí)際工作中的米單位，系統(tǒng)中默認(rèn)東西向橫坐標(biāo)為X坐標(biāo)，不加帶號(hào)，南北向縱坐標(biāo)為Y坐標(biāo)，需要時(shí)還要提供樣品對(duì)應(yīng)的地質(zhì)編碼，我們拿到這個(gè)數(shù)據(jù)后首先進(jìn)行數(shù)據(jù)檢查，以確保數(shù)據(jù)中不出現(xiàn)寫錯(cuò)、漏填、負(fù)數(shù)、0、>等字符，如果有這樣的情況要找實(shí)驗(yàn)室給予糾正。數(shù)據(jù)準(zhǔn)備好后，我們要把數(shù)據(jù)轉(zhuǎn)換成TXT后綴的文本文件，這就做好了處理前的準(zhǔn)備數(shù)據(jù)工作。2.2分析處理2.2.1重復(fù)樣三層套合方差分析一般是從分析樣品的結(jié)果中挑出來重復(fù)樣的分析值，每組四個(gè)樣品，按如下順序排列：11第一次采樣第一次分析12第一次采樣第二次分析21第二次采樣第一次分析22第二次采樣第二次分析結(jié)果：三層套合方差分析成果--元素:Cu方差來源離差平方和自由度平均離差平方和F-檢驗(yàn)值地點(diǎn)間13.224180.73511.152樣品間1.252190.0661.482分析間1.689380.044三層套合方差分析成果--元素:Au方差來源離差平方和自由度平均離差平方和F-檢驗(yàn)值地點(diǎn)間0.573180.03210.781樣品間0.056190.0030.483分析間0.232380.006查表：Fa0.05(18,19)=1.35Fb0.05(19,38)=1.22Cu：11.152>1.351.482>1.22即：F1>FaF2>FbAu：10.781>1.350.483<1.22即：F1>FaF2<Fb說明：一般是F1＞臨界值，F(xiàn)2＜臨界值。Fa>f0.05()理論值,說明采樣與分析誤差不顯著，此次分析結(jié)果可靠，F(xiàn)b>f0.05()理論值,說明采樣與分析誤差顯著，且以采樣誤差為主。例如報(bào)告中的描述：先列一個(gè)表格說明所有元素的F1和F2；重復(fù)樣三層套合方差分析，39種元素F1值遠(yuǎn)超過臨界值，說明采樣與分析誤差變化均低于區(qū)內(nèi)地球化學(xué)自然變化，數(shù)據(jù)是可用的。F2值表現(xiàn)為兩種情況，Ag、Pb、P、U、MgO、La等6個(gè)元素F2值小于臨界值，說明采樣誤差不大于分析誤差，以分析誤差為主；其余33種元素F2值均大于臨界值，說明誤差應(yīng)以采樣誤差為主。2.2.2地層編碼統(tǒng)計(jì)從樣品相對(duì)應(yīng)的地質(zhì)圖中自動(dòng)讀取地層的顏色號(hào)作為編碼，實(shí)現(xiàn)每個(gè)樣品的地層自動(dòng)編碼。在GeoIPAS系統(tǒng)中的化探->化探特征參數(shù)統(tǒng)計(jì)->地層編碼統(tǒng)計(jì)其處理結(jié)果是在輸入的數(shù)據(jù)文件中增加一列，過程是將坐標(biāo)投影到地質(zhì)圖相應(yīng)的投影參數(shù)上，以相同地層顏色來套取數(shù)據(jù)，取地層的顏色號(hào)為落入此地層的采樣點(diǎn)的地層編碼，所以要想一次做好，要先對(duì)地質(zhì)圖做個(gè)預(yù)處理，最好是按編碼順序，將顏色號(hào)改成相應(yīng)的順序號(hào)，比如Q的顏色改為1，N的顏色改為2，E的顏色改為3，等等。相對(duì)應(yīng)的出來的結(jié)果就是按順序的地層編碼了，可以直接用了，繪制的直方圖也是這個(gè)順序了。2.2.3化探特征參數(shù)統(tǒng)計(jì)化探參數(shù)統(tǒng)計(jì)成果--元素:Cu地質(zhì)單元樣品數(shù)(N)平均值(X)中位數(shù)(M)標(biāo)準(zhǔn)離差(S)變化系數(shù)(Cv)最大值最小值濃集克拉克值(C)致礦系數(shù)(Z)1.算術(shù)全區(qū)210541.91934.7328.6280.683250.492.491.1038.6162.對(duì)數(shù)全區(qū)21051.5241.5410.3070.2012.3990.3960.042.2553.算術(shù)剔除高值全區(qū)203939.01234.0222.9730.589108.722.491.02710.7344.對(duì)數(shù)剔除高值全區(qū)21051.5241.5410.3070.2012.3990.3960.042.255注：濃集克拉克值可以選用系統(tǒng)提供默認(rèn)的，也可以直接輸入用戶所在不同地區(qū)的濃集克拉克值。這個(gè)值可以保存為文件，以后導(dǎo)入使用。2.2.4R型聚類分析從譜系圖中可以看到，Cr、Ni、Co相關(guān)性最好，我們?cè)谧鼋M合異常圖時(shí)可以參考譜系圖來做元素的組合，但這個(gè)也不能完全作為組合參考，我們通常還是以常用組合方式來做組合異常圖，如：Cu-Pb-Zn；W-Sn-Mo；Cr-Ni-Co；Au-Ag-As-Sb。2.3網(wǎng)格化2.3.1工作比例尺與測(cè)網(wǎng)密度DZ／T0145-94土壤地球化學(xué)測(cè)量規(guī)范表1工作比例尺與測(cè)網(wǎng)密度工作階段比例尺矩形網(wǎng)格正方形網(wǎng)格點(diǎn)／km2線距(m)×點(diǎn)距(m)點(diǎn)線距(m)區(qū)域調(diào)查1：2000002000×500100011．1000001000×250～500500～10001～4普查1：50000500×100～250250～5004～201：25000250×50～100125～25016～80詳查1：10000100×20～5050～100100～5001：500050×10～2525～501：200020×5～10從上表可以看到，化探在實(shí)際工作中一般和物探工作一起工作的常采用矩形網(wǎng)格，但在實(shí)際成圖中卻是用正方形網(wǎng)格的好，不用經(jīng)過網(wǎng)格化處理，沒有數(shù)據(jù)損失，也符合繪圖的正方形網(wǎng)。1:20萬化探在全國范圍已基本完成，只在新疆等高山區(qū)有少量的沒有覆蓋全，最常采用的是組合樣，即：1點(diǎn)／4km2，少量的有單點(diǎn)樣的，這主要用轉(zhuǎn)網(wǎng)格數(shù)據(jù)的方式來將數(shù)據(jù)做網(wǎng)格化處理。單點(diǎn)樣的網(wǎng)格化一般采用500×500到1000×1000的網(wǎng)格距，搜索半徑3000×3000到5000×5000。組合樣做轉(zhuǎn)網(wǎng)格數(shù)據(jù)，注意起始坐標(biāo)為奇數(shù)公里，間距為2000，然后要進(jìn)行一次補(bǔ)空擴(kuò)邊處理，這主要是處理數(shù)據(jù)成圖到圖框邊。1:5萬化探是現(xiàn)在化探普查的主要工作，一般在山區(qū)采用自由網(wǎng)的水系樣，土壤測(cè)量一般用正方形網(wǎng)格采樣，1:1萬化探作為1:5萬化探異常檢查的主要工作手段，常采用100×40的矩形網(wǎng)格，少量有100×20、200×40、150×40等加密或者抽稀采樣。也有采用物探布設(shè)的測(cè)網(wǎng)非南北方向的測(cè)線。2.3.2下面用幾個(gè)問題來具體說明。問題1：1:5萬土壤測(cè)量，網(wǎng)距330×300的怎么設(shè)置搜索半徑和間距？5萬化探，網(wǎng)距500*500的，用GeoIPAS做地球化學(xué)圖的時(shí)候，怎么設(shè)置搜索半徑?答：a、關(guān)于網(wǎng)格距，最好采用正方形網(wǎng)格，網(wǎng)格距為野外采樣的點(diǎn)線距的一半，比如野外采樣為500×500，網(wǎng)格化的網(wǎng)格距選250×250就可以了，100×40的選為50×40或者50×50。b、網(wǎng)格化的搜索半徑選擇需要一定的經(jīng)驗(yàn)，建議你首先看看點(diǎn)位數(shù)據(jù)圖，選擇一個(gè)既能將所有數(shù)據(jù)聯(lián)系起來，又不至于網(wǎng)格后把未采樣區(qū)填充的參數(shù)，一般是網(wǎng)格距的2-3倍。問題2：在進(jìn)行1:5萬化探采樣時(shí)，乙類異常區(qū)采樣密度設(shè)計(jì)為12個(gè)點(diǎn)／km2，丙、丁類異常區(qū)及其它基巖區(qū)采樣密度設(shè)計(jì)為個(gè)8個(gè)點(diǎn)／km2，網(wǎng)格化時(shí)XY坐標(biāo)的間距、行列數(shù)如何選？答：一般1:5萬化探點(diǎn)線距250×250米或者500×500米；搜索半徑以滿足丙、丁類異常區(qū)網(wǎng)格后只有個(gè)別空白區(qū)或者無空區(qū)，搜索半徑一般1000×1000到2000×2000米；如果對(duì)網(wǎng)格化參數(shù)設(shè)置不太清楚，可以選擇“計(jì)算網(wǎng)格化參數(shù)”；網(wǎng)格化后如果異常比較細(xì)碎零亂，可用數(shù)據(jù)預(yù)處理→二維數(shù)據(jù)濾波處理→周圍點(diǎn)平均圓滑處理，X和Y方向點(diǎn)數(shù)選擇3，原點(diǎn)倍數(shù)選擇2-4倍。問題3：請(qǐng)教如下問題．1:2.5萬土壤測(cè)量，線距250m，點(diǎn)距40怎么設(shè)置搜索半徑和間距？答：搜索半徑500*500，間距40*40或者50*40，注意最好與原始點(diǎn)位重疊.問題4：1:1萬化探，南北向測(cè)線，100×40網(wǎng)格距。答：這樣的數(shù)據(jù)本身符合網(wǎng)格數(shù)據(jù)的格式，只是以TXT文本格式保存，所以我們可以借助GeoIPAS軟件的“數(shù)據(jù)預(yù)處理”→“數(shù)據(jù)格式轉(zhuǎn)換”→“XYZ數(shù)據(jù)轉(zhuǎn)為網(wǎng)格數(shù)據(jù)”功能，這里注意數(shù)據(jù)坐標(biāo)的起始值、終止值和網(wǎng)格間距，網(wǎng)格間距就選擇100×40，起始值和終止值要保證數(shù)據(jù)的原始點(diǎn)位都與網(wǎng)格點(diǎn)重合。這里有一個(gè)問題，就是有些項(xiàng)目用的是GPS坐標(biāo)，不一定是與網(wǎng)格點(diǎn)位完全重合，一般我們建議用規(guī)則網(wǎng)的坐標(biāo)。圓滑方法：周圍點(diǎn)平均圓滑；圓滑參數(shù)：X方向點(diǎn)數(shù)：3，Y方向點(diǎn)數(shù)：3，原點(diǎn)倍數(shù)：3問題5：我的測(cè)網(wǎng)不是正南北向是45度方向，用咱們的軟件進(jìn)行網(wǎng)格化后變成正南北的，而且圖幅形狀不是原來測(cè)網(wǎng)的形狀。請(qǐng)問這樣的問題怎樣解決?1:1萬化探，北西向測(cè)線，100×40網(wǎng)格距。答：這樣的數(shù)據(jù)雖然也是網(wǎng)格數(shù)據(jù)，但是經(jīng)過旋轉(zhuǎn)的測(cè)線方向與我們GRD所定義的網(wǎng)格間距是有夾角的，最常采用的是離散數(shù)據(jù)網(wǎng)格化方法。離散數(shù)據(jù)網(wǎng)格化選擇冪指數(shù)加權(quán)，指數(shù)因子：10；搜索范圍：圓域，300米，網(wǎng)格化間距：40×40。要用XYZ數(shù)據(jù)轉(zhuǎn)網(wǎng)格數(shù)據(jù)的方式，就要知道測(cè)線的方位角，采樣的點(diǎn)線號(hào)，起始坐標(biāo)；利用點(diǎn)線號(hào)轉(zhuǎn)換為實(shí)際點(diǎn)線距，然后再把點(diǎn)線距轉(zhuǎn)為網(wǎng)格數(shù)據(jù)，成圖后，統(tǒng)一處理，包括旋轉(zhuǎn)方位角，然后平移到圖上的起始坐標(biāo)位置套圖框。2.3.3網(wǎng)格化成圖的比較：三角剖分繪制的地化圖，有被測(cè)線方向拉長的趨勢(shì)，雖然其能保留原占坐標(biāo)位的原始值，但成圖結(jié)果不是很美觀實(shí)用。距離平方導(dǎo)數(shù)加權(quán)和冪指數(shù)加權(quán)網(wǎng)格化，對(duì)數(shù)據(jù)的壓制較大，異常線的形態(tài)較刻板。Kriging網(wǎng)格化選用了兩種搜索范圍，一種是圓域：300×300；一種是橢圓域：300×100，并且旋轉(zhuǎn)45度。這兩種方法整體上差異不大，只在局部的空區(qū)上表現(xiàn)出來不同，對(duì)數(shù)據(jù)的壓制較小，異常的形態(tài)也保持較好，也是常選用的網(wǎng)格化方法，但其也有一個(gè)問題，就是對(duì)特高值邊上的空區(qū)會(huì)處理出負(fù)值，這個(gè)不是我們希望得到的結(jié)果。2.3.4經(jīng)驗(yàn)：在處理1:1萬化探數(shù)據(jù)時(shí)，可以根據(jù)測(cè)線方向和數(shù)據(jù)坐標(biāo)點(diǎn)位的不同選擇不同的數(shù)據(jù)處理方法，要遵循的一個(gè)原則是盡量使網(wǎng)格化后的坐標(biāo)點(diǎn)位與原始數(shù)據(jù)坐標(biāo)點(diǎn)位接近，網(wǎng)格化搜索的范圍越小數(shù)據(jù)損失越小。對(duì)于網(wǎng)格間距以小于等于最小的原始點(diǎn)距為準(zhǔn)，如果要保留更多的細(xì)節(jié)，也可以將網(wǎng)格距減小一倍，如果圖面太過零亂，可以用GeoIPAS的二維濾波處理的周圍點(diǎn)平均圓滑處理，這里可以加大原點(diǎn)倍數(shù)，默認(rèn)為2，可以選擇3-5倍。所有這些處理方法可以先選擇一個(gè)元素用不同的方法做出結(jié)果來比較一下確定效果最好的方法，記住處理參數(shù)然后再統(tǒng)一批處理。2.4成圖2.4.1圖件編制的要求：a、土壤地球化學(xué)測(cè)量圖件分兩部分：基礎(chǔ)圖與推斷解釋圖。b、區(qū)域調(diào)查和普查工作的圖件編制按化探區(qū)域調(diào)查有關(guān)規(guī)范執(zhí)行。c、圖件編制必須符合地球化學(xué)勘查圖式圖例及用色標(biāo)準(zhǔn)規(guī)定。d、成果報(bào)告需做交通測(cè)區(qū)位置圖、實(shí)際材料圖、等值線圖、綜合異常圖及其它推斷解釋圖。2.4.2彩色等量線圖問題1：有些異常為什么不標(biāo)不上高值點(diǎn)?答：這是因?yàn)榫W(wǎng)格化后成圖時(shí)，最高的那一個(gè)等值線沒有套住高值點(diǎn)的坐標(biāo)位置。對(duì)這種情況一般我們只能用手動(dòng)加注極值點(diǎn)，目前還沒有太好的辦法解決。程序只搜索那個(gè)最高的或者最低的等值線范圍內(nèi)的坐標(biāo)點(diǎn)，然后從搜索的點(diǎn)中選擇最大或者最小的進(jìn)行標(biāo)注。加注后是一般刪線，保證極值點(diǎn)在最高的或者最低的等值線范圍內(nèi)。有時(shí)點(diǎn)剛好在線邊上也就不做處理了。加標(biāo)注時(shí)，在MAPGIS工程中，加入點(diǎn)位數(shù)據(jù)圖文件，將其設(shè)為關(guān)閉狀態(tài)，先在等值線圖上找到?jīng)]標(biāo)的地方，然后打開點(diǎn)位數(shù)據(jù)圖，把相應(yīng)的極值拷貝粘貼過來，然后修改子圖+，-。問題2：對(duì)于等值線高值的圈圈很密怎么處理?答：對(duì)于等值線高值的圈圈很密可以減少或者帛稀特高值的畫線間隔，不改變常量范圍的等值線。另外，有時(shí)為了圖面好看也可以處理一下特高值數(shù)據(jù)，我們?cè)贕eoIPAS系統(tǒng)中有這樣的處理模塊，就是壓制一下特高值數(shù)據(jù)，然后網(wǎng)格化成圖，可以視情況選用。使用注意事項(xiàng)：（1）數(shù)據(jù)橫坐標(biāo)(東向距)不帶投影帶號(hào)。（2）在繪制投影經(jīng)緯網(wǎng)圖框時(shí)，在原數(shù)據(jù)投影參數(shù)設(shè)置中輸入正確的數(shù)據(jù)坐標(biāo)單位和投影參數(shù)，通常高斯坐標(biāo)數(shù)據(jù)的單位為米，比例尺為1。（3）在選擇標(biāo)注原始極值時(shí)，必須選擇正確的橫坐標(biāo)、縱坐標(biāo)和元素名稱或變量名稱，否則極值標(biāo)注將會(huì)出錯(cuò)。2.4.3直方圖生成直方圖的流程為：首先選擇數(shù)據(jù)文件，程序使用TXT格式數(shù)據(jù)，其次選定元素和組間隔，當(dāng)選擇繪制全區(qū)及各地層單元直方圖時(shí)需要數(shù)據(jù)中有編碼列，指定編碼列，有地層編碼文件則裝入，沒有就手動(dòng)設(shè)置，最后設(shè)置剖面繪制起始值等設(shè)置，完成直方圖繪制。a、數(shù)據(jù)文件格式為TXT文本格式，數(shù)據(jù)要包含列標(biāo)識(shí)，數(shù)據(jù)中必須包含地質(zhì)編碼。b、直方圖的最大分組數(shù)和起始值可以用戶自己選擇，并且每個(gè)地層編碼可以分別進(jìn)行設(shè)置。c、地質(zhì)編碼相當(dāng)于數(shù)字的地質(zhì)圖，將地質(zhì)體用數(shù)字代表，每個(gè)樣品都對(duì)應(yīng)一個(gè)地質(zhì)時(shí)代，當(dāng)某個(gè)時(shí)代的樣品個(gè)數(shù)少于30個(gè)時(shí)可以將其合并到其它地質(zhì)單元中。d、當(dāng)用戶沒有編制地層編碼文件時(shí)，可臨時(shí)在編碼窗口中挑選編碼所在的地層。方法是雙擊“地層代號(hào)”下的空白行。e、使用注意事項(xiàng)：（1）編碼列必須選擇正確。（2）當(dāng)修改直方圖起始值時(shí)，小數(shù)點(diǎn)后第二次的正值的尾數(shù)為7，負(fù)值的尾數(shù)為3。f、直方圖繪圖參數(shù)設(shè)置（1）選擇繪制對(duì)應(yīng)形態(tài)曲線：即程序自動(dòng)根據(jù)直方圖勾繪是否屬于正態(tài)分布的圓滑曲線；（2）柱體內(nèi)填充顏色：即每個(gè)柱體內(nèi)填充相應(yīng)的顏色；（3）對(duì)于X刻度線，用戶可以根據(jù)實(shí)際需要進(jìn)行設(shè)置。2.4.4組合異常圖a、我們推薦在做組合異常圖時(shí)，用單個(gè)元素來做，在MAPGIS工程中加入不同單個(gè)元素異常線文件來任意組合。b、另外如果要做元素異常參數(shù)統(tǒng)計(jì)，推薦在MAPGIS中先對(duì)單元素異常進(jìn)行編輯整理壓縮保存線文件后用單元素異常參數(shù)統(tǒng)計(jì)來做。c、生成組合異常圖的流程為：首先選擇數(shù)據(jù)文件所在的路徑，其次選擇元素組合，然后進(jìn)行異常統(tǒng)計(jì)設(shè)置，最后繪制圖框完成組合異常圖。d、組合異常圖使用的數(shù)據(jù)文件為Grd格式的網(wǎng)格數(shù)據(jù)，并且要求各元素網(wǎng)格化時(shí)的參數(shù)一致。e、在數(shù)據(jù)列選擇框中選擇需組合的元素，每選擇一個(gè)元素，在異常參數(shù)設(shè)置表中就會(huì)增加該元素，并提供系統(tǒng)缺省的異常下限值和線寬度、顏色等繪制參數(shù)，缺省的異常下限計(jì)算方法有累頻和均值標(biāo)準(zhǔn)差兩種方式，兩種方式分別提供了弱、中、強(qiáng)3種選擇，對(duì)累計(jì)頻率而言弱、中、強(qiáng)分別對(duì)應(yīng)了累計(jì)頻率的85%、90%和95%，對(duì)均值標(biāo)準(zhǔn)差而言弱、中、強(qiáng)分別對(duì)應(yīng)了平均值+2倍的離差、+2.3倍的離差和+2.7倍的離差，如果對(duì)系統(tǒng)提供的缺省下限值不滿意，可以輸入確定的異常下限值；每種元素異常對(duì)應(yīng)的線顏色、線寬、線型均可在異常參數(shù)設(shè)置表中進(jìn)行設(shè)置。2.4.5單元素異常圖單元素異常圖使用的數(shù)據(jù)文件為Grd格式的網(wǎng)格數(shù)據(jù)，繪制方式有3種：直接用異常下限勾繪異常、采用異常襯度值勾繪異常和采用剩余異常值勾繪異常：a、直接用異常下限值勾繪異常，并以一定的含量等級(jí)劃分出異常濃度外帶、中帶和內(nèi)帶。當(dāng)數(shù)據(jù)為單一地質(zhì)單元或數(shù)據(jù)基本符合正態(tài)分布時(shí)，按照全區(qū)異常下限值的1、2～4、3～8倍劃分三個(gè)濃度帶，勾繪異常外、中、內(nèi)帶。b、采用異常襯度值勾繪異常，當(dāng)數(shù)據(jù)為兩個(gè)以上的多個(gè)母體分布時(shí)，可依多個(gè)母體分別統(tǒng)計(jì)出的襯值，統(tǒng)一確定襯度異常下限，勾繪異常。并以襯度異常下限的1、2～4、3～8倍劃分三個(gè)異常濃度分帶。c、采用剩余異常值勾繪異常，以求得普查區(qū)區(qū)域背景趨勢(shì)面，然后用原始數(shù)據(jù)減去該點(diǎn)的背景趨勢(shì)值所獲得的剩余異常值作圖。當(dāng)求得剩余異常值后，以正剩余異常值確定剩余異常下限，并以一定的倍數(shù)勾繪正剩余異常及濃度分區(qū)。在應(yīng)用剩余異常值制作異常圖時(shí)，應(yīng)對(duì)求取地球化學(xué)背景趨勢(shì)面的方法進(jìn)行試驗(yàn)、研究，以確保圈出的異常具有明確的找礦指示意義。2.4.6剖析圖GeoIPAS做異常剖析圖有兩種方法：地化圖切剖析圖和GRD數(shù)據(jù)生剖析圖。首先都要準(zhǔn)備好切剖析圖的矩形框，這個(gè)可以用MAPGIS中線編輯繪制矩形線，將異常包括進(jìn)去。a、地化圖切剖析圖在GeoIPAS中做，利用地化圖切剖析圖模塊，添加地質(zhì)圖和地化圖的MAPGIS文件，然后設(shè)置地質(zhì)圖的投影參數(shù)，導(dǎo)入裁剪框，設(shè)置裁剪框的投影參數(shù)，選擇文件時(shí)，將地質(zhì)圖放在第一個(gè)位置，選擇結(jié)果文件目錄，進(jìn)行排版設(shè)置，確定后進(jìn)行裁剪，裁剪結(jié)果就組成一個(gè)完整的剖析圖MAPGIS工程。推薦在MAPGIS中做，可以用剖析圖的裁剪框，做成一個(gè)區(qū)，然后在MAPGIS地質(zhì)圖編輯中用工程裁剪，把地質(zhì)圖和地化圖的裁剪結(jié)果文件加入到剖析圖工程中，進(jìn)行排版，組成一個(gè)完整的剖析圖MAPGIS工程，第一個(gè)位置放地質(zhì)圖。b、GRD數(shù)據(jù)生剖析圖推薦在GeoIPAS中用GRD數(shù)據(jù)生剖析圖，地質(zhì)圖可在MAPGIS中用工程裁剪得到，裁剪框也用MAPGIS中線編輯繪制矩形線得到。要注意的是有三處設(shè)置投影參數(shù)的地方，源數(shù)據(jù)投影參數(shù)就是GRD文件的，生成的結(jié)果剖析圖的投影參數(shù)，以及裁剪框的投影參數(shù)，選擇排列好剖析的圖順序，可以生成規(guī)范要求的三級(jí)濃度分帶的剖析圖。2.5單元素異常參數(shù)統(tǒng)計(jì)GeoIPAS中有兩種方式：a、在組合異常圖時(shí)，選擇上“異常參數(shù)統(tǒng)計(jì)”輸入原始數(shù)據(jù)文件名以及坐標(biāo)，選擇一個(gè)或者多個(gè)網(wǎng)格化后的GRD文件，選擇異常下限，以及最后一列的統(tǒng)計(jì)元素名，然后計(jì)算出明碼文件的同時(shí)還會(huì)產(chǎn)生一個(gè)結(jié)果文件+pSta.csv的屬性文件，如結(jié)果文件名是Ag(.wal.wat),那么屬性文件名AgpSta.csv，這個(gè)文件就是和異常對(duì)應(yīng)的參數(shù)統(tǒng)計(jì)，GeoIPAS同時(shí)也在結(jié)果線文件上掛接好屬性。b、由異常線來做異常參數(shù)統(tǒng)計(jì)做元素異常參數(shù)統(tǒng)計(jì)時(shí)，我們推薦用組合異常圖把每一個(gè)元素的異常生成出來，在生圖時(shí)選擇不套圖框，不繪比例尺，然后在MAPGIS中編輯對(duì)異常線做處理。具體包括：刪除小的異常線，閉合被圖框或工區(qū)邊角截?cái)嗟漠惓＞€，連接并閉合被空區(qū)截?cái)嗟木€，使元素異常個(gè)數(shù)與實(shí)際統(tǒng)計(jì)的物理數(shù)一致，然后按編排的異常順序號(hào)修改線的ID號(hào)，并壓縮保存線文件(MAPGIS編輯子系統(tǒng)中，在左邊工程文件欄內(nèi)，點(diǎn)右鍵，選擇壓縮保存工程，就把所有編輯狀態(tài)的文件都?jí)嚎s保存了)，然后將這個(gè)異常線文件作為輸入，在化探->單元素異常參數(shù)統(tǒng)計(jì)中統(tǒng)計(jì)異常參數(shù)，注意設(shè)置線文件投影參數(shù)，這個(gè)就是我們生成異常線時(shí)的參數(shù)，還需要設(shè)置數(shù)據(jù)文件的投影參數(shù)。運(yùn)行后就會(huì)有一個(gè)CSV的結(jié)果，其中包括“元素名-編號(hào)”列，如Au-1，就和我們的異常編號(hào)對(duì)應(yīng)一致，可以直接當(dāng)作單元素異常參數(shù)統(tǒng)計(jì)附表了，同時(shí)元素異常線上也掛接好屬性了。單元素異常的參數(shù)，包括：異常點(diǎn)數(shù)、異常面積、異常的平均值、異常的對(duì)數(shù)平均值（幾何平均值）、異常內(nèi)樣品最大值、異常標(biāo)準(zhǔn)離差、異常襯度、異常規(guī)模、異常NAP，同時(shí)進(jìn)行了排序。異常編號(hào)點(diǎn)數(shù)序數(shù)異常面積序數(shù)異常均值序數(shù)幾何平均值序數(shù)異常極大值序數(shù)異常均方差序數(shù)異常襯度序數(shù)異常規(guī)模序數(shù)異常NAP序數(shù)帶數(shù)序數(shù)排序結(jié)果Cu-011453.81983460.5557191.779431969.62196.88729171.513891978.519475.7828261314Cu-022936.40139363.18181.791821892.911213.224151.579518148.384410.1113229Cu-032140.7054951273.41121.8197114105.51045.382121.835251223.5706131.29476132211Cu-04971.93215890.516751.918896195.7244.343832.26292597.605964.37237312Cu-052140.56871652.405251.709522663.442315.6058121.31012257.05472210.745068221322Cu-085101.13022989.12661.921125130936.564862.22815655.523492.518319227Cu-102140.6018811454.09231.72112466.642117.7484111.35225238.4805200.