張楠楠，周可法，陳 曦，李 宏

(1.中國(guó)科學(xué)院新疆生態(tài)與地理研究所，烏魯木齊 830011;2.新疆維吾爾自治區(qū)遙感中心，烏魯木齊 830011;3.中國(guó)科學(xué)院新疆礦產(chǎn)資源研究中心，烏魯木齊 830011)

基于ETM+的遙感蝕變信息提取方法對(duì)比研究

張楠楠1，2，周可法1，3，陳 曦1，李 宏2

(1.中國(guó)科學(xué)院新疆生態(tài)與地理研究所，烏魯木齊 830011;2.新疆維吾爾自治區(qū)遙感中心，烏魯木齊 830011;3.中國(guó)科學(xué)院新疆礦產(chǎn)資源研究中心，烏魯木齊 830011)

為從微觀到宏觀了解、掌握和運(yùn)用圍巖蝕變的特征與原巖特征，尋找大型礦床或大型礦集區(qū)，在總結(jié)常用遙感蝕變信息提取模型的基礎(chǔ)上，結(jié)合主成分分析(PCA)和支持向量機(jī)(SVM)算法，建立了基于“掩模處理-主成分分析-支持向量機(jī)(MASK-PCA-SVM，MPS)”的遙感蝕變信息提取模型。在新疆哈圖地區(qū)進(jìn)行的基于ETM+的不同蝕變信息提取模型應(yīng)用對(duì)比研究表明，使用常用模型提取出的蝕變信息和礦點(diǎn)吻合效果較差，而基于MPS模型提取出的蝕變信息與已知礦點(diǎn)的吻合率可達(dá)86.51%;結(jié)合地質(zhì)專家知識(shí)、地質(zhì)圖和一個(gè)已知異常點(diǎn)的分析結(jié)果表明，基于MPS模型的蝕變信息提取精度較高，在研究區(qū)地質(zhì)應(yīng)用中效果較好，為今后該區(qū)的遙感蝕變信息提取工作提供了一種新思路。

ETM+;蝕變;信息提取;SVM

0 引言

從遙感圖像中深層次提取地學(xué)信息是地學(xué)定量分析的一項(xiàng)十分重要的工作。自遙感技術(shù)問世以來(lái)，遙感數(shù)據(jù)種類多、數(shù)量大、信息豐富，包含了地物的光譜信息、空間信息等，很好地利用這些信息必將有利于目標(biāo)信息的提?。?-4］;而如何從各類遙感數(shù)據(jù)源中有效、準(zhǔn)確地提取相關(guān)專題所需的目標(biāo)特征信息，一直是遙感界努力探索的目標(biāo)，國(guó)內(nèi)外遙感工作者都在不斷地設(shè)計(jì)、研制和總結(jié)對(duì)這種遙感信息的提取和識(shí)別技術(shù)。蝕變信息是找礦的一個(gè)重要標(biāo)志，尋找與成礦相關(guān)的蝕變信息在地質(zhì)找礦中有著舉足輕重的作用。絕大部分內(nèi)生礦床都伴隨有相應(yīng)的熱液圍巖蝕變;圍巖蝕變的特征與原巖特征從微觀到宏觀都有明顯的區(qū)別，了解、掌握和運(yùn)用這種特征對(duì)尋找大型礦床或大型礦集區(qū)具有特別重要的地質(zhì)意義［5］。因此，用數(shù)字處理技術(shù)從遙感圖像中提取遙感蝕變信息，可對(duì)縮小找礦范圍和提高找礦準(zhǔn)確度發(fā)揮重要作用。

近30 a來(lái)，已有大量學(xué)者利用 ETM+，ASTER或SPOT數(shù)據(jù)進(jìn)行過遙感蝕變信息提取研究。在國(guó)外，早在1976年，戈茨利用短波紅外波段1.6 μm和2.2 μm波長(zhǎng)的兩譜帶反射比的比值定量劃分了蝕變巖和未蝕變巖;Abrams［6］利用 TM圖像波段比值+主成分變換方法排除了植被覆蓋因素干擾，提取了蝕變信息;Loughlin［7］利用Landsat-5 TM圖像數(shù)據(jù)的4個(gè)波段(TM1，TM3，TM5和 TM7或 TMl，TM4，TM5和TM7)，設(shè)計(jì)了主成分變換+特定主因子求反的方法，填制了巴西熱帶地區(qū)殘積土壤中的三價(jià)鐵和羥基蝕變巖信息圖;Rokos等［8］利用TM，SPOT數(shù)據(jù)，用B5/B7，B3/B1等比值和主成分分析等數(shù)學(xué)運(yùn)算進(jìn)行了蝕變帶的識(shí)別;Kusky等［9］利用B5/B7(R)，B5/B1(G)，B5/B4(B)的假彩色合成及分析，在阿拉伯—努比亞地盾的干旱氣候下提取了金礦化蝕變信息;Tommasoa等［10］在阿根廷利用ASTER數(shù)據(jù)的光譜信息及比值方法進(jìn)行了蝕變信息提取及斑巖銅礦的識(shí)別。在國(guó)內(nèi)，劉燕君［11］在MSS圖像上成功地把礦化點(diǎn)含礦圍巖蝕變強(qiáng)烈部位增強(qiáng)出來(lái);趙元洪等［12］針對(duì)浙江新昌地區(qū)的植被覆蓋問題采用TM4/TM3比值+主成分變換技術(shù)壓抑植被干擾達(dá)到了提取蝕變信息的目的;馬建文［13］提出TM掩模+主成分變換+分類的識(shí)別、提取礦產(chǎn)弱信息方法，在遼寧二道溝金礦提取了三價(jià)鐵蝕變火山碎屑巖信息，在河北華北地臺(tái)北緣提取了含金鉀化帶信息;劉素紅等［14］利用Gramschmidt投影方法從高山區(qū)提取了TM數(shù)據(jù)中的含礦蝕變帶信息;張玉君等［15-16］總結(jié)了ETM+(TM)蝕變遙感異常提取方法研究與應(yīng)用的地質(zhì)依據(jù)、波譜前提、方法選擇和技術(shù)流程;洪金益等［17］開始嘗試?yán)肧VM方法對(duì)TM數(shù)據(jù)進(jìn)行蝕變信息提取;耿新霞等［18］利用ASTER數(shù)據(jù)和主成分分析方法在新疆阿巴宮提取了鉛鋅礦蝕變信息;刁海等［19］利用ASTER數(shù)據(jù)和主成分分析及分形模型計(jì)算異?；叶鹊拈撝?，提取了蝕變異常信息;鄧吉秋等［20］利用TM數(shù)據(jù)提取了錳礦化蝕變信息并利用該信息進(jìn)行了成礦預(yù)測(cè);另外，平仲良［21］、張滿郎［22］、劉慶生等［23］、李培軍等［24］、呂鳳軍等［25］、張遠(yuǎn)飛等［26］、程知言等［27］在金礦蝕變信息的識(shí)別和增強(qiáng)提取方面也做了不少工作，獲得了不少新認(rèn)識(shí)。

1995年，Vapnik［28］提出了支持向量機(jī)(support vector machines，SVM)算法，我國(guó)最早見到的 SVM方面的研究是在2000年以后［29-30］。隨后，SVM在各領(lǐng)域的應(yīng)用受到高度重視，并被廣泛地應(yīng)用于許多方面［31-34］。但用SVM方法提取遙感蝕變信息的研究在國(guó)內(nèi)外較少，尚處于起步階段;把SVM和傳統(tǒng)方法相結(jié)合進(jìn)行有效的蝕變信息提取的研究則更少，只是在近兩年陸續(xù)有人開始嘗試?yán)肧VM進(jìn)行蝕變礦物信息提取。本文旨在利用比較容易獲取到的ETM+數(shù)據(jù)，結(jié)合SVM和傳統(tǒng)的主成分分析，建立基于“掩模處理 -主成分分析 -支持向量機(jī)(MASK-PCA-SVM，MPS)”的蝕變礦物提取模型，并在新疆哈圖地區(qū)將本文提出的MPS方法與傳統(tǒng)的方法進(jìn)行了應(yīng)用對(duì)比研究，探索出一種適合本研究區(qū)的遙感蝕變信息提取和識(shí)別方法。

1 研究區(qū)地質(zhì)概況

作為研究區(qū)的新疆哈圖地區(qū)位于準(zhǔn)噶爾盆地西北緣，屬于西準(zhǔn)噶爾華力西期褶皺帶扎依爾一達(dá)拉布特復(fù)向斜東段的北翼［35］。區(qū)域地層以石炭系為主，包括下石炭統(tǒng)太勒古拉組(C1t)、包古圖組(C1b)和希貝庫(kù)拉斯組(C1x);泥盆系地層出露較少，為中泥盆統(tǒng)巴爾雷克組(D2b)及喀依爾巴依組(D2k);第三系、第四系地層出現(xiàn)于山前及溝谷中［36］(圖1)。區(qū)內(nèi)自北向南分布著哈圖、安齊、達(dá)拉布特3條走向北東、大致平行的區(qū)域性斷裂。區(qū)域內(nèi)巖漿活動(dòng)強(qiáng)烈，侵入巖發(fā)育，主要有鉀長(zhǎng)花崗巖(如哈圖巖體、阿克巴斯套巖體及克拉瑪依北巖體等)、斜長(zhǎng)花崗巖、閃長(zhǎng)巖及超基性巖等［37］。

圖1 哈圖地區(qū)地質(zhì)簡(jiǎn)圖(根據(jù)1∶50萬(wàn)地質(zhì)圖修編)Fig.1 Geological sketch-map of Hatu area(based on the geological map at 1∶500 000 scale)

同大多數(shù)金礦一樣，礦區(qū)內(nèi)圍巖蝕變較發(fā)育，主要為黃鐵礦化及毒砂化，其次為碳酸鹽化、硅化、絹云母化和綠泥石化等［38］。黃鐵礦化及毒砂化較普遍，分布于脈巖的本身及其兩側(cè)，寬約1～4 m，沿破碎帶也有黃鐵礦化分布。凡見有黃鐵礦化和毒砂化的地方，均含有金;此兩種蝕變程度愈強(qiáng)，含金亦愈富。在地表氧化帶中，因受風(fēng)化作用影響，黃鐵礦和毒砂為褐鐵礦所代替，并使巖石變成褐—紅色，此種情況通稱為“鐵化”。硅化、絹云母化、綠泥石化和碳酸鹽化與成礦關(guān)系也極為密切。絹云母一類的礦物與其他金礦的重要蝕變巖——黃鐵絹英巖中的絹云母相比，哈圖金礦的絹云母頗似淺灰黃色臘狀的粘土礦物，多出現(xiàn)在高品位的含金石英脈中，分布不均勻。本區(qū)蝕變帶中的綠泥石因受到成礦后的強(qiáng)烈動(dòng)力變質(zhì)作用而嚴(yán)重片理化。

2 遙感蝕變信息提取模型綜述

2.1 比值模型

比值分析(ratioing)是一種多光譜圖像處理方法，由2個(gè)波段對(duì)應(yīng)像元的灰度值之比或幾個(gè)波段組合的對(duì)應(yīng)像元灰度值之比獲得［39-40］。通常根據(jù)特定目標(biāo)的反射率特征，選取反射率差值最大的2個(gè)波段來(lái)生成比值圖像。利用蝕變礦物的吸收及反射光譜特征，有針對(duì)性地?cái)U(kuò)大要提取信息的對(duì)應(yīng)波段的差異，突出目標(biāo)信息，壓制、消除或減弱地形陰影和植被的干擾。

根據(jù)蝕變礦物在ETM+相關(guān)波段的光譜特征可總結(jié)出若干比值模型:含鐵離子的蝕變礦物可用ETM+3/ETM+1提取;黃鐵鉀礬和赤鐵礦可用ETM+5/ETM+4或ETM+3/ETM+4提取;含錳離子的蝕變礦物、含羥基的蝕變礦物和含碳酸根蝕變礦物可用ETM+5/ETM+7提取。

2.2 主成分模型

主成分分析(PCA)在數(shù)學(xué)上常被稱為KL變換，是在統(tǒng)計(jì)上的多維正交線性變換［41］。PCA主要通過一次或多次坐標(biāo)軸旋轉(zhuǎn)變換后，使得原來(lái)在單波段圖像上難以識(shí)別的目標(biāo)信息在主成分圖像上得到增強(qiáng)和分離。PCA實(shí)際上能對(duì)圖像數(shù)據(jù)進(jìn)行集中和壓縮，可將多光譜圖像各波段中高度相關(guān)的信息集中到少數(shù)幾個(gè)波段并盡可能保證這些波段的信息互不相關(guān)，即用幾個(gè)綜合性波段代表原圖像的多個(gè)波段，使處理的數(shù)據(jù)量大幅度減少。目前，在遙感蝕變信息提取中常用的PCA是Crosta等［42］提出的Crosta特征PCA 方法，即利用ETM+1，3，4，5 波段進(jìn)行PCA提取鐵化蝕變信息，利用 ETM+1，4，5，7波段進(jìn)行PCA提取粘土化蝕變信息。

2.3 光譜角模型

光譜角制圖(spectral angle mapper，SAM)方法是一種光譜匹配技術(shù)，這種技術(shù)基于估計(jì)像元光譜與樣本光譜或是混合像元中端元組分(endmember)光譜的相似性來(lái)區(qū)分各像元點(diǎn)的光譜曲線［43］。SAM方法的原理是把光譜作為矢量投影到N維空間上(其維數(shù)N為試驗(yàn)時(shí)選取的所有波段數(shù)，各光譜曲線被看作有方向有長(zhǎng)度的矢量，而各光譜之間形成的夾角叫作光譜角)，該方法考慮的是光譜矢量的方向而非光譜矢量的長(zhǎng)度。由于光譜矢量的長(zhǎng)度關(guān)聯(lián)到影像亮度而SAM方法對(duì)亮度值的影響并不敏感，所以當(dāng)計(jì)算影像光譜角時(shí)，SAM方法僅僅需要知道光譜矢量的方向而不關(guān)注影像本身的亮度。

2.4 MPS模型

本文在常用模型和前人研究的基礎(chǔ)上，提出了基于MPS的遙感礦化蝕變信息提取模型。MPS的含義是:M代表MASK(即掩模處理)，P代表PCA(即主成分分析)，S代表SVM(即支持向量機(jī)處理)。SVM是從線性可分情形下的最優(yōu)分類超平面發(fā)展起來(lái)的，其基本思想是要使訓(xùn)練集的間隔最大化［44］(圖2)。

圖2 SVM的基本思想示意圖Fig.2 Sketch map of the basic idea of SVM

圖2中的黑色圓點(diǎn)和白色圓點(diǎn)分別代表2類樣本;H是分類面，H1和H2分別是穿過距離2類分類面最近的樣本且平行于分類面的直線;2條直線之間的距離叫做分類間隔(margin)。所謂最優(yōu)分類面就是要求分類面不但能將2類樣本正確分開(訓(xùn)練錯(cuò)誤率為0)，而且可使分類間隔最大。

基于MPS模型提取遙感礦化蝕變信息的整體思路(圖3)是:在通過掩模處理去除干擾信息的基礎(chǔ)上進(jìn)行PCA;利用PCA的結(jié)果與已知礦點(diǎn)進(jìn)行比對(duì)，在圖像與已知礦點(diǎn)相吻合的區(qū)域采集SVM所需的訓(xùn)練樣本，對(duì)樣本進(jìn)行訓(xùn)練后對(duì)整個(gè)研究區(qū)進(jìn)行礦化蝕變異常信息提取(具體過程將在另文中詳述)。

圖3 基于MPS模型的遙感蝕變信息提取流程圖Fig.3 Flow chart of extracting alteration information from remote sensing image based on the MPS model

3 模型應(yīng)用對(duì)比

為了進(jìn)行不同模型提取遙感礦化蝕變信息的效果評(píng)價(jià)，同樣用其他模型對(duì)哈圖金礦地區(qū)進(jìn)行了蝕變信息提取，幾種模型的提取結(jié)果對(duì)比如下:

1)分別用比值ETM+3/ETM+1提取“鐵化”信息，用比值ETM+5/ETM+7提取“泥化”信息，提取結(jié)果見圖4(a);

2)分別用常規(guī)的 PCA(ETM+1，3，4，5)和 PCA(ETM+1，4，5，7)提取“鐵化”信息和“泥化”信息，提取結(jié)果見圖4(b);

3)在已知礦點(diǎn)和蝕變帶處采集樣本光譜，用SAM提取“鐵化”信息和“泥化”信息，光譜角閥值設(shè)為0.25(該閥值最大設(shè)置為1，為盡可能減少假異常，本文設(shè)置的閥值相對(duì)較低)，提取結(jié)果見圖4(c);

4)在選取核函數(shù)和最優(yōu)參數(shù)后，用訓(xùn)練模型對(duì)整個(gè)圖像進(jìn)行識(shí)別，最終得到基于MPS模型的研究區(qū)蝕變信息提取圖(圖4(d))。

圖4(a)—(d)中，紅色圖斑代表“鐵化”信息，藍(lán)色圖斑代表“泥化”信息，黃色的點(diǎn)代表已知礦點(diǎn)。

圖4 不同模型提取遙感礦化蝕變信息結(jié)果Fig.4 Mineralization alteration information extracted by different models using remote sensing methods

4 結(jié)果驗(yàn)證

本文用兩種方法對(duì)上述提取結(jié)果進(jìn)行驗(yàn)證。一種方法是利用地質(zhì)圖上的已知礦點(diǎn)及蝕變帶作為查證依據(jù)，主要通過礦點(diǎn)及已知蝕變帶與所提結(jié)果的吻合程度來(lái)進(jìn)行對(duì)比驗(yàn)證;第二種方法是利用專家知識(shí)及已有的野外典型異常查證資料來(lái)進(jìn)行驗(yàn)證。

4.1 與地質(zhì)圖對(duì)比驗(yàn)證

研究區(qū)地質(zhì)圖(圖1)上標(biāo)出的蝕變有黃鐵礦化、綠泥石化、褐鐵礦化、混合巖化等。對(duì)比地質(zhì)圖和圖4(d)可以看出，除了阿克巴斯套巖體邊緣的混合巖化體現(xiàn)較弱外，其他已知蝕變的地區(qū)在圖4(d)上基本都有所反映。本文收集的該區(qū)已知礦點(diǎn)共126個(gè)，疊加到圖4(d)上后，除了17個(gè)礦點(diǎn)吻合效果不好，其他礦點(diǎn)均與所提取的蝕變信息吻合較好，吻合率達(dá)86.51%。

比較圖4(a)—(c)可以看出，用比值模型、PCA模型和SAM模型提取出的遙感蝕變異常中假異常都較多。其中，用比值模型提取出的假異常信息最多，提取出的“鐵化”信息中大部分為植被信息，而“泥化”信息中大部分為高亮度巖體信息，與礦點(diǎn)吻合效果不好(只有24個(gè)礦點(diǎn)周邊有蝕變信息，吻合率只有3.37%);用PCA模型提取出的遙感蝕變信息與比值模型相比，減少了部分假異常信息，但仍然有植被和高亮度巖體等信息存在，蝕變信息與礦點(diǎn)吻合率為69.05%;用SAM模型提取出的遙感蝕變信息與PCA模型相比，減少了植被的信息，但仍然存在一些高亮度巖體信息，蝕變信息與礦點(diǎn)吻合率為65.08%。MPS模型圖(4(d))與其他3種模型相比，所提取出的遙感蝕變信息中幾乎沒有植被等假異常信息，而且和礦點(diǎn)吻合效果最好(與已知礦點(diǎn)的吻合率達(dá)86.51%)。

4.2 用專家知識(shí)驗(yàn)證

根據(jù)專家知識(shí)可知，圖4(d)提取結(jié)果中的1—5號(hào)黃框內(nèi)區(qū)域?yàn)橐巴獾刭|(zhì)專家確定的蝕變區(qū)，其中1號(hào)區(qū)域在哈圖斷裂附近，有著很好的容礦和導(dǎo)礦條件，因此發(fā)現(xiàn)礦點(diǎn)數(shù)量較多，也是蝕變數(shù)量多且蝕變強(qiáng)烈的地段，并且和地質(zhì)圖上的褐鐵礦礦化及黃鐵礦化蝕變的位置吻合較好;2號(hào)區(qū)域是齊Ⅰ金礦床的所在區(qū)域，結(jié)合已有的地質(zhì)資料，在齊Ⅰ金礦西南約5 km的安齊斷裂中部，存在Ap49異常點(diǎn)(其地表照片見圖5)，地表可見褐鐵礦化。該區(qū)位于安齊斷裂的上盤，處于構(gòu)造破碎帶，蝕變類型為硅化、黃鐵礦化及碳酸鹽化。

圖5 Ap49異常區(qū)地表照片F(xiàn)ig.5 Field photo of Ap49 abnormal area

從圖4(d)中截取的圖6可以看出，所提信息與地質(zhì)資料描述相符。

圖6 Ap49異常區(qū)提取的蝕變信息Fig.6 Alteration information extracted in Ap49 abnormal area

3號(hào)和4號(hào)區(qū)域?yàn)樾r體周邊，發(fā)現(xiàn)的礦床也較多，并且該區(qū)的小巖體周邊已成為找礦的重點(diǎn)地段;5號(hào)區(qū)域的東側(cè)也發(fā)現(xiàn)很多礦床蝕變比較強(qiáng)烈。6號(hào)區(qū)域則是從地質(zhì)圖上已知的有該區(qū)典型蝕變礦物綠泥石分布的區(qū)域。

5 結(jié)論

通過對(duì)3種傳統(tǒng)模型和本文提出的MPS模型在新疆哈圖地區(qū)進(jìn)行蝕變遙感信息提取的對(duì)比研究可以看出:

1)比值模型和主成分模型都容易造成假異常過多現(xiàn)象，光譜角模型受光譜分辨和空間分辨率的限制對(duì)數(shù)據(jù)要求較高。

2)本文提出的基于“掩模處理-主成分分析-支持向量機(jī)(MASK-PCA-SVM，MPS)”的遙感蝕變信息提取模型的應(yīng)用效果較好，與已知礦點(diǎn)吻合效果較好，避免了明顯的假異常，突出了目標(biāo)信息即礦化蝕變信息，在研究區(qū)與已知礦床吻合率達(dá)86.51%，并與地質(zhì)圖上的已知蝕變信息相吻合。另外，利用專家知識(shí)發(fā)現(xiàn)，所提的礦化蝕變信息與野外調(diào)查情況相符，與已有的野外典型異常的查證資料相符;所提取的遙感異常與野外Ap49點(diǎn)的異常特征描述相符、位置吻合。因此，本文的MPS蝕變信息提取模型，在很大程度上提高了蝕變異常信息提取的精度，更有效地去除假異常信息，為蝕變信息提取提供了一種新思路。這種模型在其他地區(qū)也可以推廣使用。

3)MPS模型也存在一些問題，例如參數(shù)的選擇需要在不同地區(qū)進(jìn)行調(diào)整，處理整景圖像時(shí)運(yùn)行速度較慢，這些需要在今后的研究中進(jìn)一步改進(jìn)。

［1］趙英時(shí).遙感應(yīng)用分析原理與方法［M］.北京:科學(xué)出版社，2003:5-35.

［2］周成虎.遙感影像地學(xué)理解與分析［M］.北京:科學(xué)出版社，1999:11-20.

［3］陳述彭，童慶禧，郭華東.遙感信息機(jī)理研究［M］.北京:科學(xué)出版社，1998:18-20.

［4］趙鵬大.定量地學(xué)方法及應(yīng)用［M］.北京:高等教育出版社，2004:8-22.

［5］陳淳福，梁有彬.內(nèi)生金礦床圍巖蝕變及其找礦意義［J］.地質(zhì)找礦論叢，1999，5(3):63-70.

［6］Abrams M J.Landsat-4 Thermatic Mapper and Thematic Mapper Simulator Data for a Porphyry Copper Deposit［J］.Photogrammetric Engineering and Remote Sensing，1984，50(8):1171-1173.

［7］Loughlin W P.Principal Component Analysis for Alteration Mapping［J］.Photogrammetric Engineering and Remote Sensing，1991，57(9):163-1169.

［8］Rokos D，Argialas R，Mavrantza，et al.Structural Analysis for Gold Mineralization Using Remote Sensing and Geochemical Techniques in a GIS Environment:Island of Lesvos，Hellas［J］.Natural Resources Research，2000，9(4):277-293.

［9］Kusky T M，Talaat M R.Structural Controls on Neoproterozoic Mineralization in the South Eastern Desert，Egypt:an Integrated Field，Landsat TM，and SIR-C/X SAR Approach［J］.Journal of African Earth Sciences，2002，35(1):107-121.

［10］Tommasoa I D，Rubinstein N.Hydrothermal Alteration Mapping Using ASTER Data in the Infiernillo Porphyry Deposit，Argentina［J］.Ore Geology Reviews，2007，32(1-2):275-290.

［11］劉燕君.礦產(chǎn)信息的遙感地面模式［M］.北京:地質(zhì)出版社，1993:25-50.

［12］趙元洪，張福祥，陳南峰.波段比值的主成份復(fù)合在熱液蝕變信息提取中的應(yīng)用［J］.國(guó)土資源遙感，1991(3):12-18.

［13］馬建文.利用TM數(shù)據(jù)快速提取含礦蝕變帶方法研究［J］.遙感學(xué)報(bào)，1997，1(3):208-213.

［14］劉素紅，馬建文，藺啟忠.通過Gram-Schmidt投影方法在高山區(qū)提取TM數(shù)據(jù)中含礦蝕變帶信息［J］.地質(zhì)與勘探，2000，36(5):62-65.

［15］張玉君，楊建民，陳 微.ETM+(TM)蝕變遙感異常提取方法研究與應(yīng)用——地質(zhì)依據(jù)和波譜前提［J］.國(guó)土資源遙感，2002(4):30-37.

［16］張玉君，曾朝銘，陳 微.ETM+(TM)蝕變遙感異常提取方法研究與應(yīng)用——方法選擇和技術(shù)流程［J］.國(guó)土資源遙感，2003(2):44-49.

［17］洪金益，姚學(xué)恒，潘 東.基于SVM遙感圖像礦化信息提取試驗(yàn)［J］.礦業(yè)研究與開發(fā)，2004，24(5):63-65.

［18］耿新霞，楊建明，姚佛軍，等.新疆阿巴宮鉛鋅礦遙感找礦信息提取及地化多元信息綜合分析［J］.地質(zhì)與勘探，2011，46(5):942-952.

［19］刁 海，張 達(dá)，狄永軍，等.基于主成分分析和分形模型的ASTER蝕變異常信息提取［J］.國(guó)土資源遙感，2011(2):75-80.

［20］鄧吉秋，謝 楊，張寶一，等.ETM+圖像錳礦化蝕變信息提取與找礦預(yù)測(cè)［J］.國(guó)土資源遙感，2011(1):102-105.

［21］平仲良.從陸地衛(wèi)星TM數(shù)據(jù)提取膠東某地區(qū)某種類型金礦的圍巖蝕變信息［J］.國(guó)土資源遙感，1993(4):42-45

［22］張滿郎.金礦化蝕變信息提取中的主成分分析［J］.遙感技術(shù)與應(yīng)用，1996，11(3):1-6.

［23］劉慶生，燕守勛，馬超飛，等.內(nèi)蒙哈達(dá)門溝金礦區(qū)山前鉀化帶遙感信息提?。跩］.遙感技術(shù)與應(yīng)用，1999，14(3):7-11.

［24］李培軍，龍玄耀，劉 立.基于ASTER數(shù)據(jù)的蛇綠巖組分識(shí)別:以德爾尼礦區(qū)為例［J］.巖石學(xué)報(bào)，2007，23(5):1175-1180.

［25］呂鳳軍，郝躍生，石 靜，等.ASTER遙感數(shù)據(jù)蝕變遙感異常提取研究［J］.地球?qū)W報(bào)，2009，30(2):271-276.

［26］張遠(yuǎn)飛，袁繼明，朱谷昌，等.基于遙感數(shù)據(jù)隨機(jī)模型的空間結(jié)構(gòu)分與蝕變信息提?。跩］.國(guó)土資源遙感，2010(4):34-39.

［27］程知言，趙培松.云南播卡金礦遙感蝕變信息提取及成礦預(yù)測(cè)［J］.國(guó)土資源遙感，2011(2):135-140.

［28］Vapnik V N.The Nature of Statistical Learning Theory［M］.New York:Springer-Verlag，1995:18-78.

［29］Vapnik V N.統(tǒng)計(jì)學(xué)習(xí)理論的本質(zhì)［M］.張學(xué)工，譯.北京:清華大學(xué)出版社，2000:5-35.

［30］邊肇琪.模式識(shí)別［M］.北京:清華大學(xué)出版社，2001:10-45.

［31］Joachims T.Text Categorization with Support Vector Machines:Learning with Many Relevant Features［C］//Proceedings of ECML-98，10th European Conference on Machine Learning.1998:1225-1235.

［32］Zien A，R?tsch G，Mika S，et al.Engineering Support Vector Machine Kernels that Recognize Translation Initiation Sites［C］//Proceedings of the 24th Annual International Conference on Machine Learning.Oxford:Oxford University Press，2000:799-807.

［33］盧增祥，李衍達(dá).交互支持向量機(jī)學(xué)習(xí)算法及其應(yīng)用［J］.清華大學(xué)學(xué)報(bào):自然科學(xué)版，1999，39(7):93-97.

［34］凌旭峰，楊 杰，葉晨洲.基于支撐向量機(jī)的人臉識(shí)別技術(shù)［J］.紅外與激光工程，2002，30(5):318-327.

［35］劉鵬鸚.新疆哈圖金礦地質(zhì)地球化學(xué)特征及其成礦作用［J］.黃金，1995，16(9):2-6.

［36］魏少妮，朱永峰.新疆包古圖中酸性巖漿侵位的P-T-fO2條件及巖體地球化學(xué)研究［J］.地質(zhì)學(xué)報(bào)，2010，84(7):1017-1029.

［37］李華芹，陳富文，蔡 紅.新疆西準(zhǔn)噶爾地區(qū)不同類型金礦床Rb-Sr同位素年代研究［J］.地質(zhì)學(xué)報(bào)，2000，74(2):181-191.

［38］王冶華，張慧玉，王科強(qiáng)，等.ETM圖像數(shù)據(jù)中礦化蝕變信息的提取——以新疆哈圖地區(qū)為例［J］.黃金，2009，30(2):6-9.

［39］朱亮璞.遙感地質(zhì)學(xué)［M］.北京:地質(zhì)出版社，1999:25-78.

［40］劉志杰，韓先菊，劉蔭椿.比值-特征主成分混合分析提取金礦蝕變信息——以贛南遂川地區(qū)為例［J］.黃金地質(zhì)，1999，4(1):74-77.

［41］李玉珍，王宜懷.主成分分析及算法［J］.蘇州大學(xué)學(xué)報(bào):自然科學(xué)版，2005，21(1):32-36.

［42］Crósta A P，Moore J M.Enhancement of Landsat Thematic Mapper Imagery for Residual Soil Mapping in SW Minas Gerais［C］//Proceedings of the 7th(ERIM)Thematic conference:Remote Sensing for Exploration Geology.Calgar，1989:225-255.

［43］安 斌，陳書海，嚴(yán)衛(wèi)東.SAM法在多光譜圖像分類中的應(yīng)用［J］.中國(guó)體視學(xué)與圖像分析，2005，10(1):55-60.

［44］Vapnik V N.Statistical Learning Theory［M］.New York:Wiley，1998:3-55.

A Comparative Study of Extraction Methods for Alteration Information Based on ETM+

ZHANG Nan-nan1，2，ZHOU Ke-fa1，3，CHEN Xi1，LI Hong2
(1.Xinjiang Institute of Ecology and Geography，Chinese Academy of Sciences，Urumqi 830011，China;2.Remote Sensing Center of Xinjiang Uyghur Autonomous Region，Urumqi 830011，China;3.Xinjiang Research Center for Mineral Resources，Chinese Academy of Sciences，Urumqi 830011，China)

In order to understand，grasp and apply the characteristics of altered wall rock and the original rock from microcosmic and macroscopic angles and look for large deposits or large ore districts，the authors summarized the common extraction models of alteration information，combined the principal component analysis(PCA)with support vector machines(SVM)to establish the MASK process-PCA-SVM(MPS)model，applied and comparatively studied these models in Hatu area of Xinjiang，and made use of the two methods to verify the extraction results.The results show that the extracting results of the common models are rather poorly consistent with the known deposits，but the alteration information extracted by MPS model is better concordant with the known mineral deposits and altered belts，with the coincidence rate arriving at 86.51%.The combination of geological experts'knowledge with geological maps and a known anomaly spot reveals that the result of MPS model is in accord with the geological fact.Therefore，the extracting precision based on MPS model is higher than that of the common models，and the effect of application is fairly good in the study area.The model proposed in this paper provides a new idea for extracting remote sensing alteration information in the future.

ETM+;alteration;information extraction;SVM

TP 75

A

1001-070X(2012)02-0034-07

張楠楠(1980-)，女，博士，副研究員，畢業(yè)于中國(guó)科學(xué)院新疆生態(tài)與地理研究所，主要從事遙感地質(zhì)及成礦預(yù)測(cè)方面的研究。E-mail:znn_0802@163.com。

(責(zé)任編輯:劉心季)

10.6046/gtzyyg.2012.02.07

2011-05-26;

2011-07-02

中國(guó)科學(xué)院數(shù)字地球重點(diǎn)實(shí)驗(yàn)室開放基金項(xiàng)目(編號(hào):2011LDE005)、國(guó)家重點(diǎn)基礎(chǔ)研究發(fā)展計(jì)劃(編號(hào):2007CB411308)、國(guó)家科技支撐計(jì)劃(編號(hào):2006BAB07B07-04)、新疆自然科學(xué)基金項(xiàng)目(編號(hào):200821185)和新疆科技廳項(xiàng)目(編號(hào):200815116)共同資助。