（一）編圖試點(diǎn)研究簡(jiǎn)要史料。
 
  1973年3月，湖南地質(zhì)局物探大隊(duì)編制了《地球物理（化學(xué)）探礦圖件編制規(guī)范》。
 
  1974年4月，冶金部地質(zhì)司在廣西桂林舉行冶金地質(zhì)系統(tǒng)區(qū)域地球化學(xué)編圖及新方法、新技術(shù)交流會(huì)。
 
  1978年4月，物探所任天祥在《物探與化探》上發(fā)表了《關(guān)于1：20萬(wàn)區(qū)域化探資料整理利用方法的討論》，提出編制1：20萬(wàn)和1: 50萬(wàn)地區(qū)域化學(xué)圖的方法。
 
  1980年10月，謝學(xué)錦等發(fā)表了《皖浙贛邊區(qū)區(qū)域化探的資料整理與異常評(píng)價(jià)》，編制了皖浙贛邊區(qū)Cu、Pb、Zn等12個(gè)元素的地球化學(xué)圖。
 
  同年，桂林冶金地質(zhì)研究所與廣東地質(zhì)勘探公司940隊(duì)開展了“廣東大陸區(qū)域化探編圖及找礦研究”。
 
  1985年10月，地礦部物化探所、江西物探大隊(duì)和北京計(jì)算中心合作完成的《贛西北地球化學(xué)試驗(yàn)圖集》由地質(zhì)出版社出版，這是我國(guó)公開出版的第一份地球化學(xué)圖集。
 
  1987年8月，地礦部物探所研制的IBMPC/XT應(yīng)用開發(fā)及通用成圖系統(tǒng)，移植和開發(fā)了20多個(gè)常用化探數(shù)據(jù)處理程序，并研制了一套微機(jī)成圖系統(tǒng)。
 
  1988～1992年期間，北京計(jì)算中心承擔(dān)地質(zhì)行業(yè)科學(xué)技術(shù)發(fā)展基金資助的科研項(xiàng)目（劉如英、金諄浩、鄔寬廉、鄔可筠等），開發(fā)研制了“計(jì)算機(jī)地球化學(xué)圖編圖系統(tǒng)”（小比例尺）(GCMCS)，在我國(guó)率先比較系統(tǒng)地研究小比例尺計(jì)算機(jī)地球化學(xué)編圖系統(tǒng)。
 
  1990年11月，地礦部北京計(jì)算中心（劉如英、金諄浩、鄔可筠）與遼寧物探隊(duì)（白希杰等）合作，首次用1: 20萬(wàn)區(qū)域化探數(shù)據(jù)為遼寧編制全省1：50萬(wàn)地球化學(xué)圖，編制了37種元素的等值線圖、色塊圖及異常圖，以及8種元素的綜合異常圖，可以全面反映遼寧省有關(guān)元素的分布特征和規(guī)律，是一份全省性的基礎(chǔ)地球化學(xué)資料。
 
  1993年，受地礦部勘技司委托，北京計(jì)算中心劉如英、金諄浩、趙玉祥等在中國(guó)地質(zhì)大學(xué)（北京）舉辦培訓(xùn)班，推廣應(yīng)用“計(jì)算機(jī)地球化學(xué)編圖系統(tǒng)（小比例尺）”。
 
  1993年12月，成都地質(zhì)學(xué)院完成的地面放射性測(cè)量工作站通過專家鑒定，以微機(jī)或筆記本機(jī)為平臺(tái)，實(shí)現(xiàn)了地面放射性測(cè)量數(shù)據(jù)的采集、傳輸、處理、成圖成像和解釋全過程的計(jì)算機(jī)化，該系統(tǒng)為國(guó)內(nèi)首創(chuàng)。
 
  20世紀(jì)90年代后，大部分?。▍^(qū)、市）按省（區(qū)、市）域編制了省（區(qū)、市）地球化學(xué)圖，如浙江、新疆、湖北、貴州等（楊永忠，1997），同時(shí)按成礦區(qū)帶和大區(qū)域編制了區(qū)帶地球化學(xué)圖，如東昆侖、北祁連、內(nèi)蒙古北山、華南等區(qū)域性地球化學(xué)圖。主要采用色塊、等值線和數(shù)據(jù)加密插值的成圖方法。
 
  （二）編制彩色地球化學(xué)圖件中一些問題的處理。
 
  插值：由于輸出的彩色圖件為全面色，每個(gè)像元點(diǎn)都應(yīng)具有數(shù)據(jù)及其對(duì)應(yīng)的顏色，因此輸出彩圖之前必須求出每個(gè)像元點(diǎn)的插值。處理中常用的插值算法包括距離平方倒數(shù)、雙線性方法、雙三次樣條等。單一的插值算法往往存在缺點(diǎn)。距離平方倒數(shù)方法插值不連續(xù)，并且，當(dāng)插值點(diǎn)與數(shù)據(jù)點(diǎn)的距離太近時(shí)，可能出現(xiàn)誤差；而雙三次樣條則在插值中出現(xiàn)高于數(shù)據(jù)中的最大值和低于數(shù)據(jù)中的最小值、甚至出現(xiàn)負(fù)值的不合理情況。為了避免這些問題使用聯(lián)合算法插值，成功地解決了這一問題。
 
  色階劃分：對(duì)已插值數(shù)據(jù)進(jìn)行分級(jí)，并對(duì)各數(shù)量級(jí)指定顯示的顏色。分級(jí)方法包括以標(biāo)準(zhǔn)差的倍數(shù)作為相鄰級(jí)別的差值、用累積頻率（或分位數(shù)）分級(jí)、使用等比級(jí)數(shù)分級(jí)等。
 
  色彩的調(diào)配：抖動(dòng)算法。由于彩色墨水只有4種顏色(青C、品紅M、黃Y和黑K)，由CMYK的不同組合至多可以獲得8種顏色（包括白色）。為了顯示更多的顏色，必須用加大象元、犧牲分辨率的方法。例如使用4×4的抖動(dòng)矩陣，每個(gè)象元可以獲得大約5 000種顏色，但此時(shí)對(duì)于分辨率為300 dpi的Novajet繪圖儀，分辨率降到300/4=75 dpi。
 
  (三）全國(guó)范圍的地球化學(xué)編圖。
 
  進(jìn)入20世紀(jì)90年代，全國(guó)區(qū)域化探掃面計(jì)劃宜測(cè)地區(qū)的任務(wù)基本完成，積累了豐富而系統(tǒng)的高質(zhì)量數(shù)據(jù)，許多?。▍^(qū)、市）據(jù)之編制了?。▍^(qū)、市）域的地球化學(xué)圖和圖冊(cè)。于是，編制全國(guó)范圍地球化學(xué)圖的任務(wù)就提上日程。
 
  1993年，中國(guó)勘查技術(shù)院立設(shè)立“中國(guó)小比例尺地球化學(xué)圖編制”項(xiàng)目，由地科院物化探研究所組織實(shí)施，在謝學(xué)錦和任天祥主持下，向運(yùn)川負(fù)責(zé)編圖方法技術(shù)研究。這次編圖以1：20萬(wàn)和1：50萬(wàn)區(qū)域化探掃面數(shù)據(jù)為基礎(chǔ)，以1：5萬(wàn)和1: 2.5萬(wàn)圖幅為統(tǒng)計(jì)單元計(jì)算元素含量算術(shù)平均值，作為編圖的數(shù)據(jù)源。采用改進(jìn)的距離倒數(shù)加權(quán)法對(duì)數(shù)據(jù)進(jìn)行加密（插值），將數(shù)據(jù)劃分為16個(gè)色階，生成彩色像元，通過彩色噴墨繪圖儀繪制圖件。分別編制了39種元素的全國(guó)1：500萬(wàn)地球化學(xué)圖（1:2.5萬(wàn)圖幅數(shù)據(jù)集）和1：
 
  1 000萬(wàn)全國(guó)地球化學(xué)圖集（1:5萬(wàn)圖幅數(shù)據(jù)集）。編圖成果于1996年提交，這是我國(guó)首次編制的全國(guó)性39種元素的地球化學(xué)圖件。
 
  2001年，中國(guó)地質(zhì)調(diào)查局下達(dá)新的全國(guó)性編圖任務(wù)，即“全國(guó)區(qū)域地球化學(xué)系列圖編制”項(xiàng)目，由牟緒贊、奚小環(huán)負(fù)責(zé)，中國(guó)地質(zhì)調(diào)查局發(fā)展研究中心組織實(shí)施，謝學(xué)錦、任天祥等為顧問，向運(yùn)川為編圖組長(zhǎng)。這次編圖的特點(diǎn)是，直接以化探掃面的原始數(shù)據(jù)為編圖依據(jù)，首次匯集了全國(guó)28個(gè)省（區(qū)、市）的1: 20萬(wàn)和1: 50萬(wàn)區(qū)域化探資料，涉及1：20萬(wàn)圖幅1299個(gè)，1：50萬(wàn)圖幅18個(gè)，數(shù)據(jù)點(diǎn)約142萬(wàn)個(gè)，包括39種元素及氧化物的測(cè)試數(shù)據(jù)，覆蓋面積達(dá)650×104 km2。最終，分別編制出了39種元素與氧化物的1：500萬(wàn)全國(guó)地球化學(xué)圖和1：1 200萬(wàn)地球化學(xué)圖集。并且，在編圖同時(shí)建立全國(guó)區(qū)域地球化學(xué)數(shù)據(jù)庫(kù)。項(xiàng)目歷時(shí)4年，2005年6月提交編圖成果。
 
  由于采用原始數(shù)據(jù)編圖，數(shù)據(jù)量大，地球化學(xué)數(shù)據(jù)處理又有其特殊性，出現(xiàn)了一系列技術(shù)問題需要解決。編圖工作采用中國(guó)地質(zhì)調(diào)查局研制開發(fā)的“多元地學(xué)空間數(shù)據(jù)管理與分析系統(tǒng)GeoExpl”作為編圖平臺(tái)，使數(shù)據(jù)處理、編圖和輸出等過程一體化操作。第一，要分別對(duì)各元素進(jìn)行系統(tǒng)誤差的處理，以便得到很好的圖面效應(yīng)。為此，制定了4條系統(tǒng)誤差校正的原則：①分析誤差源，所展示的數(shù)據(jù)誤差與周邊數(shù)據(jù)值具有明顯的臺(tái)階狀。
 
  ②數(shù)據(jù)誤差在空間上具有區(qū)域性特點(diǎn)，省、大區(qū)域、圖幅或分析批次。③在數(shù)據(jù)值的分布上，掩蓋了地球化學(xué)特征和地質(zhì)特征展布的延續(xù)性和規(guī)律性。④在數(shù)據(jù)校正方法上，盡可能地選擇線性校正，通過簡(jiǎn)單的計(jì)算可以復(fù)原數(shù)據(jù)。第二，數(shù)據(jù)網(wǎng)格化處理。大區(qū)域編圖難以利用原始點(diǎn)位直接構(gòu)圖，需要對(duì)原始數(shù)據(jù)進(jìn)行網(wǎng)格化處理，通過對(duì)不同網(wǎng)格化參數(shù)的選擇、對(duì)比和研究，確定了效果較好的網(wǎng)格化處理方法，即：網(wǎng)格距4 kmX4 km的網(wǎng)度；數(shù)距搜索方式采用搜索半徑為12 km的圓域；數(shù)據(jù)計(jì)算方法采用改進(jìn)的距離倒數(shù)加權(quán)法，即指數(shù)距離冪倒數(shù)加權(quán)法。第三，地球化學(xué)專題圖的生成。應(yīng)用GeoExpl提供的專題圖生成模塊，制作各元素等量線、等值區(qū)圖，以及相關(guān)的特征圖等。采用累計(jì)頻率與對(duì)數(shù)分級(jí)相結(jié)合，同時(shí)參考數(shù)據(jù)的特征值如平均值、中位數(shù)等，將數(shù)據(jù)分為19級(jí)，繼而劃分級(jí)次和色階，分級(jí)統(tǒng)計(jì)計(jì)算特征值，利用系統(tǒng)提供的專題圖創(chuàng)建模塊，生成地球化學(xué)等值線和等值區(qū)，分別以線圖層和區(qū)圖層的形式，作為系統(tǒng)顯示窗的圖形顯示數(shù)據(jù)。用這套方法制成的圖件，應(yīng)用了GIS技術(shù)和現(xiàn)代數(shù)據(jù)處理技術(shù)，是近十幾年來計(jì)算機(jī)技術(shù)在地球化學(xué)數(shù)據(jù)管理、處理、成圖方法研究的結(jié)晶。