區(qū)化數(shù)據(jù)預(yù)處理

化探資料數(shù)據(jù)準(zhǔn)備與預(yù)處理向運(yùn)川區(qū)域地球化學(xué)研究課題組

2007年6月全國(guó)礦產(chǎn)資源潛力評(píng)價(jià)主要內(nèi)容資料與數(shù)據(jù)準(zhǔn)備數(shù)據(jù)預(yù)處理應(yīng)用目標(biāo)

應(yīng)用區(qū)域地球化學(xué)進(jìn)行礦產(chǎn)資源預(yù)測(cè)研究是以1:20萬(wàn)和1:50萬(wàn)水系沉積物資料為基礎(chǔ)，收集局部研究區(qū)域的1:5萬(wàn)水系沉積物和1:1萬(wàn)土壤測(cè)量數(shù)據(jù)，以及巖石地球化學(xué)資料等，以成礦地質(zhì)理論和地球化學(xué)理論為指導(dǎo)，對(duì)數(shù)據(jù)進(jìn)行處理、分析、成圖、推斷解釋、預(yù)測(cè)靶區(qū)圈定、及定量預(yù)測(cè)分析等。針對(duì)《全國(guó)礦產(chǎn)資源潛力評(píng)價(jià)》項(xiàng)目目的，結(jié)合區(qū)域地球化學(xué)的特點(diǎn)，重點(diǎn)對(duì)銅、鉛、鋅、鎳、鎢、錫、金、鉻、鉬、銻、稀土、銀、硼等礦種進(jìn)行預(yù)測(cè)分析、并初步估算單礦種資源量。數(shù)據(jù)源1:20萬(wàn)或1:50萬(wàn)水系沉積物資料與數(shù)據(jù)1:5萬(wàn)水系沉積物資料與數(shù)據(jù)1:5萬(wàn)－1:50萬(wàn)巖石測(cè)量資料與數(shù)據(jù)1:5000－1:2.5萬(wàn)土壤測(cè)量資料與數(shù)據(jù)化探異常檢查與查證報(bào)告與數(shù)據(jù)其他相關(guān)地質(zhì)、礦產(chǎn)、地理景觀等資料數(shù)據(jù)分類采樣信息數(shù)據(jù)元素分析數(shù)據(jù)質(zhì)量監(jiān)控?cái)?shù)據(jù)圖形數(shù)據(jù)文檔報(bào)告

數(shù)據(jù)準(zhǔn)備采樣信息數(shù)據(jù)▲采樣記錄卡->數(shù)據(jù)表（Access或Excel）數(shù)據(jù)結(jié)構(gòu)與內(nèi)容參考《區(qū)域地球化學(xué)勘查規(guī)范》及相關(guān)標(biāo)準(zhǔn)

▲組合樣信息->數(shù)據(jù)表▲采樣點(diǎn)位布置圖->數(shù)字化（MAPGIS）

水系沉積物采樣點(diǎn)分布信息組合樣信息數(shù)據(jù)準(zhǔn)備元素分析數(shù)據(jù)▲數(shù)據(jù)表存儲(chǔ)，采用Access或Excel

▲空間坐標(biāo)信息，經(jīng)度、緯度（十進(jìn)制度），直角坐標(biāo)（公里）

▲地質(zhì)等編碼信息（地層、巖石、景觀區(qū)、構(gòu)造單元、成礦單元等）

▲分析元素值，單位：Aung/g,氧化物％，其他ug/g

數(shù)據(jù)準(zhǔn)備質(zhì)量監(jiān)控?cái)?shù)據(jù)

▲分析元素信息表－記錄元素分析方法

▲標(biāo)樣推薦值

▲標(biāo)樣分析結(jié)果

▲標(biāo)樣控制樣質(zhì)量控制參數(shù)

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

▲基礎(chǔ)地質(zhì)、地理

▲采樣點(diǎn)位分布圖

▲異常圖

▲異常查證分布圖

采用MAPGIS或其他數(shù)字化軟件分層數(shù)字化

數(shù)據(jù)準(zhǔn)備文檔資料

▲地球化學(xué)勘查報(bào)告

▲異常查證成果報(bào)告

▲省內(nèi)有關(guān)地球化學(xué)與礦產(chǎn)研究報(bào)告及文章

▲填寫各類報(bào)告主要內(nèi)容信息表，包括：報(bào)告名稱，提交單位，提交時(shí)間，報(bào)告編寫主要人員，涉及范圍，內(nèi)容摘要，主要成果

條件許可的情況下將各類報(bào)告轉(zhuǎn)換成Word文檔保存

數(shù)據(jù)預(yù)處理空間坐標(biāo)轉(zhuǎn)換異點(diǎn)檢查分析系統(tǒng)誤差校正景觀區(qū)校正

數(shù)據(jù)表數(shù)據(jù)預(yù)處理空間坐標(biāo)轉(zhuǎn)換▲根據(jù)各省確定的地理投影模型，對(duì)準(zhǔn)備各類空間數(shù)據(jù)（數(shù)據(jù)表和圖形數(shù)據(jù)）轉(zhuǎn)換為統(tǒng)一的坐標(biāo)系統(tǒng)，橢球參數(shù)建議采用西安80/IUUG1975坐標(biāo)系(即長(zhǎng)半軸：6378140米，短半軸：6356755.3041米)。

▲采用專業(yè)軟件系統(tǒng)或GIS系統(tǒng)

1.表數(shù)據(jù)空間坐標(biāo)轉(zhuǎn)換

(1)確定源坐標(biāo)系統(tǒng)、參數(shù)，、單位，以及坐標(biāo)項(xiàng)

(2)確定目標(biāo)坐標(biāo)系統(tǒng)、參數(shù)、單位，給定新坐標(biāo)項(xiàng)

源坐標(biāo)系統(tǒng)目標(biāo)坐標(biāo)系統(tǒng)源坐標(biāo)項(xiàng)目標(biāo)坐標(biāo)項(xiàng)數(shù)據(jù)預(yù)處理空間坐標(biāo)轉(zhuǎn)換

2.圖形數(shù)據(jù)空間坐標(biāo)轉(zhuǎn)換建議在數(shù)字化完成后，進(jìn)行此項(xiàng)轉(zhuǎn)換

(1)確定源坐標(biāo)系統(tǒng)、參數(shù)、坐標(biāo)單位

(2)確定目標(biāo)坐標(biāo)系統(tǒng)、參數(shù)，坐標(biāo)單位

源坐標(biāo)系統(tǒng)目標(biāo)坐標(biāo)系統(tǒng)數(shù)據(jù)預(yù)處理異點(diǎn)檢查▲檢查空間坐標(biāo)和分析元素的奇異值

▲推薦采用AccessSQL查詢語(yǔ)言

如：SELECT索引編碼,經(jīng)度,緯度,AuFROMHHLLJJWHERE(Au<=0OrAu>2000）AndAu>-1E+30▲異點(diǎn)值替換，直接在表中替換或用相關(guān)軟件

低于檢出限的數(shù)據(jù)取檢出限的1/2,如：<3.0=1.5

高于檢出最高極限的數(shù)據(jù)取125%,如：>1000=1250

表名數(shù)據(jù)項(xiàng)條件數(shù)據(jù)預(yù)處理元素分析系統(tǒng)誤差校正

由于區(qū)域地球化學(xué)數(shù)據(jù)受地理景觀、采樣介質(zhì)、分析手段的影響，不可避免的產(chǎn)生明顯的系統(tǒng)誤差，尤其是涉及較大區(qū)域的化探數(shù)據(jù)，這種誤差更為突出。因此，在數(shù)據(jù)處理之前，有必要分別對(duì)各元素進(jìn)行系統(tǒng)誤差的處理?！Ｕ煞治龌虿蓸拥纫鸬南到y(tǒng)誤差

▲推薦采用GeoExpl

可采用系統(tǒng)的數(shù)據(jù)校正模塊或二維空間分析子系統(tǒng)結(jié)合SQL數(shù)據(jù)分析功能可分別進(jìn)行離散數(shù)據(jù)和網(wǎng)格數(shù)據(jù)的校正

▲建議直接對(duì)原始離散數(shù)據(jù)校正

數(shù)據(jù)預(yù)處理元素分析系統(tǒng)誤差校正系統(tǒng)誤差校正的原則：（1）分析誤差源，所展示的數(shù)據(jù)誤差與周邊數(shù)據(jù)值具有明顯的臺(tái)階狀；（2）數(shù)據(jù)誤差在空間上具有區(qū)域性特點(diǎn)，區(qū)域、圖幅或分析批次；（3）在數(shù)據(jù)值的分布上，掩蓋了地球化學(xué)特征和地質(zhì)特征展布的延續(xù)性和規(guī)律性。（4）在數(shù)據(jù)校正方法上，盡可能地選擇線性校正，通過(guò)簡(jiǎn)單的計(jì)算可以復(fù)原數(shù)據(jù)?！`差識(shí)別方法

數(shù)據(jù)表瀏覽圖形方式瀏覽數(shù)據(jù)點(diǎn)位圖數(shù)據(jù)預(yù)處理元素分析系統(tǒng)誤差校正▲校正步驟

(1)按原始點(diǎn)位采用符號(hào)分級(jí)的方式生成元素的符號(hào)圖，分級(jí)方法采用累計(jì)頻率方式。(2)通過(guò)校正圖示窗瀏覽原始數(shù)據(jù)全圖，確定具有明顯的數(shù)據(jù)臺(tái)階區(qū)域，區(qū)域的確定原則是由區(qū)域->圖幅->批次；采用圖形編輯工具，在圖上直接圈定要處理的區(qū)域（用面的方式表示）。(3)確定局部圖幅和分析批次范圍產(chǎn)生的系統(tǒng)誤差，校正單元由系統(tǒng)提供的工具直接在顯示窗中勾繪，確定完所有需要校正單元，各校正單元的ID需設(shè)定為唯一。

數(shù)據(jù)預(yù)處理元素分析系統(tǒng)誤差校正▲校正步驟

(4)建立校正單元與處理數(shù)據(jù)表空間位置索引關(guān)系。

(5)確定各單元的校正值或校正系數(shù)，主要方法是與單元周邊數(shù)據(jù)進(jìn)行對(duì)比分析，部分規(guī)律性較復(fù)雜的單元可以通過(guò)統(tǒng)計(jì)規(guī)律確定，同時(shí)還需考慮地球化學(xué)分布的整體空間分布趨勢(shì)和地質(zhì)背景；計(jì)算方法推薦采用：Vai=AVi+B

其中：Vai為校正點(diǎn)校正后數(shù)據(jù)，A為校正系數(shù)，Vi為校正點(diǎn)原始數(shù)據(jù)，B校正常數(shù)。A與B值的確定參照校正單元周邊數(shù)據(jù)單元（正常的數(shù)據(jù)單元）確定，可以由經(jīng)驗(yàn)判斷，也可通過(guò)統(tǒng)計(jì)規(guī)律確定。

(6)數(shù)據(jù)校正，可采用SQL語(yǔ)言操作模式或應(yīng)用軟件系統(tǒng)提供的專用工具，按確定的校正值對(duì)各校正單元逐一進(jìn)行計(jì)算。

(7)利用校正計(jì)算結(jié)果重新生成符號(hào)分級(jí)圖。

(8)觀察全圖，對(duì)部分校正結(jié)果不理想的單元，可通過(guò)上述步驟，對(duì)單元和校正值進(jìn)行調(diào)整，并重新計(jì)算，直到校正數(shù)據(jù)和成圖效果符合全局規(guī)律為止。

數(shù)據(jù)預(yù)處理景觀區(qū)及地質(zhì)構(gòu)造分區(qū)校正

由于不同的地理景觀和地質(zhì)因素的影響，在進(jìn)行大區(qū)域成圖或異常研究中，使得低背景或弱小異常被掩蓋，失去了許多有用成礦信息，這種校正的目的是突出弱小異常。

▲校正方法

通過(guò)圈定研究區(qū)內(nèi)景觀或構(gòu)造單元，建立單元與數(shù)據(jù)表的空間關(guān)系，利用標(biāo)準(zhǔn)單元單元層為基準(zhǔn)，對(duì)其他校正單元的數(shù)據(jù)進(jìn)行歸一化校正處理?！椒ㄔ?/p>

設(shè)有m個(gè)校正單元，第t個(gè)校正單元為標(biāo)準(zhǔn)單元，則：

其中，Cij’為第j個(gè)單元第i個(gè)樣品校正后的值，Cij為第j個(gè)單元第i個(gè)樣品原始值，nj為第j個(gè)單元的樣品數(shù)，n為全區(qū)樣品數(shù)J單元平均值總平均值或標(biāo)準(zhǔn)單元平均值數(shù)據(jù)處理與分析

景觀區(qū)及地