段宗恩，朱谷昌，張思穎，姬祥永，隗含濤，朱莉莉

(1．中南大學(xué)地球科學(xué)與信息物理學(xué)院，長(zhǎng)沙410083;2．有色金屬礦產(chǎn)地質(zhì)調(diào)查中心，北京100012)

哈密黑山口地區(qū)遙感蝕變信息提取特征波段優(yōu)選研究

段宗恩1，2，朱谷昌1，2，張思穎1，姬祥永1，隗含濤1，朱莉莉1

(1．中南大學(xué)地球科學(xué)與信息物理學(xué)院，長(zhǎng)沙410083;2．有色金屬礦產(chǎn)地質(zhì)調(diào)查中心，北京100012)

在遙感蝕變信息提取研究中，特征波段的選取是主要的技術(shù)關(guān)鍵。本文以新疆哈密市黑山口地區(qū)為研究區(qū)，利用ETM衛(wèi)星遙感數(shù)據(jù)，在光譜數(shù)據(jù)空間幾何結(jié)構(gòu)特征分析的基礎(chǔ)上，通過(guò)對(duì)七個(gè)多光譜波段的回歸偏度分析，評(píng)價(jià)蝕變信息的客觀存在性，從而選取出針對(duì)不同蝕變類(lèi)型的特征波段，然后采用主成分分析、最優(yōu)密度分割等方法，提取出礦化蝕變信息，并對(duì)結(jié)果信息進(jìn)行了評(píng)價(jià)。通過(guò)野外地質(zhì)調(diào)查和光譜測(cè)量的實(shí)地驗(yàn)證，提取的蝕變信息與實(shí)地情況吻合性好，取得了較理想的效果。結(jié)論認(rèn)為，次級(jí)斷裂構(gòu)造帶上含礦化蝕變的鐵鎂質(zhì)巖體找礦前景較好。

遙感蝕變信息二維散點(diǎn)圖回歸偏度特征波段黑山口

Duan Zong－en，Zhu Gu－chang，Zhang Si－ying，Ji Xiang－yang，Wei Han－tao，Zhu Li－li．Selection of characteristic bands in alteration information extraction from remote sensing images for the Heishankou area，Hami，Xinjiang［J］．Geology and Exploration，2012，48(4):0823－0828．

引言

遙感蝕變信息提取需要兩個(gè)條件:第一，研究區(qū)客觀上應(yīng)該存在著礦化蝕變的成礦地質(zhì)背景;第二，礦化蝕變信息在遙感圖像上至少應(yīng)該有弱的反應(yīng)(張遠(yuǎn)飛，2008)。遙感弱蝕變信息提取主要依賴于合適的算法模型，常用的方法有比值法和主成分分析方法。趙元洪等提出了波段比值的主成分復(fù)合法(趙元洪，1991)，劉慶生等采用“主成分分析+比值法+分類(lèi)+混合像元分解”成功地提取了內(nèi)蒙哈達(dá)門(mén)溝金礦區(qū)山前鉀化帶遙感信息(劉慶生，1999)。目前國(guó)內(nèi)外學(xué)者主要集中于蝕變信息提取算法研究(荊鳳等，2005)，缺乏對(duì)遙感信息結(jié)構(gòu)等內(nèi)在信息機(jī)理研究，本文從回歸偏度和二維散點(diǎn)圖等方面入手對(duì)蝕變信息進(jìn)行分析。

針對(duì)多光譜遙感蝕變信息提取技術(shù)，主要是以巖石礦物的電磁波特征反射譜帶作為理論基礎(chǔ)，運(yùn)用數(shù)學(xué)方法，從遙感圖像數(shù)據(jù)中提取出礦化蝕變信息(李紅，2010)。本文以新疆哈密市黑山口地區(qū)為研究區(qū)，利用ETM多光譜遙感數(shù)據(jù)，首先通過(guò)波段序偶的二維散點(diǎn)圖分析，對(duì)背景、干擾和蝕變信息三個(gè)主要對(duì)象的光譜空間結(jié)構(gòu)特征進(jìn)行研究，檢測(cè)蝕變信息的客觀存在性;再應(yīng)用回歸偏度理論，計(jì)算各波段的回歸偏度值，從而優(yōu)選出存在蝕變信息的特征波段;然后采用主成分分析、最優(yōu)密度分割等方法，提取出反映客觀實(shí)際的蝕變信息。

1 研究區(qū)概況

哈密黑山口地區(qū)位于哈密市東南約150 km，地理坐標(biāo)為東經(jīng)94°45'～95°;北緯41°10'～41°20'。區(qū)內(nèi)地形較緩，低山丘陵居多，無(wú)常年流水，干溝發(fā)育。

該區(qū)大地構(gòu)造位置處于塔里木板塊與哈薩克斯坦－準(zhǔn)噶爾板塊活動(dòng)大陸邊緣膠結(jié)帶，中亞造山帶東天山東段。近東西向的黃山斷裂是控制成礦巖體的主要斷裂構(gòu)造，侵入巖為華力西晚期的橄欖巖、蘇長(zhǎng)巖、輝石巖、閃長(zhǎng)巖等，含鎂鐵－超鎂鐵質(zhì)巖體侵位于下石炭統(tǒng)干洞組，干洞組主要為細(xì)碧巖、砂礫巖、砂礫灰?guī)r和生物碎屑灰?guī)r。該區(qū)的巖體鐵染和礦床圍巖碳酸巖化、綠泥石化強(qiáng)烈。本區(qū)賦礦巖體包含四套巖體組合，主要為含礦橄欖巖－輝長(zhǎng)巖組合(傅飄兒，2009)(見(jiàn)圖1)。

銅鎳硫化物礦床賦存于橄欖巖輝長(zhǎng)巖和閃長(zhǎng)巖組成的鎂鐵－超鎂鐵質(zhì)雜巖體中，含礦巖體巖相發(fā)育完全，橄欖巖相發(fā)育，是鐵鎂－超鐵鎂質(zhì)成礦巖漿作用的典型區(qū)域。

圖1 哈密黑山口地區(qū)區(qū)域地質(zhì)簡(jiǎn)圖Fig．1 Geological sketch map of the Heishankou area in Hami，Xinjiang

2 遙感波段數(shù)據(jù)分析

2．1 遙感數(shù)據(jù)特征

本文選用ETM數(shù)據(jù)。以ETM為代表的Landsat類(lèi)多光譜遙感數(shù)據(jù)，由于其地面覆蓋范圍寬，空間分辨率、光譜分辨率能滿足1∶50000區(qū)域地質(zhì)調(diào)查及地質(zhì)信息提取的要求，數(shù)據(jù)形式易于處理，多波段優(yōu)化組合的圖像信息豐富。ETM的各個(gè)波段各自針對(duì)著不同的探測(cè)對(duì)象和探測(cè)目的。對(duì)于地質(zhì)應(yīng)用而言各個(gè)波段都能提供地質(zhì)構(gòu)造、地形地貌信息。波段1～4能夠區(qū)分巖石中鐵、錳礦物和含鐵、錳礦物的相對(duì)含量，尤其是第4波段對(duì)于含三價(jià)鐵的礦物比較敏感，往往借此區(qū)分巖性;第5波段對(duì)綠簾石族礦物的特征譜帶敏感;第6波段識(shí)別地?zé)岙惓?、巖石和構(gòu)造的含水性，鑒別地質(zhì)構(gòu)造有一定的用途;第7波段對(duì)于羥基礦物等具明顯的光譜吸收性。

2．2 光譜空間結(jié)構(gòu)特征分析

光譜數(shù)據(jù)的二維散點(diǎn)圖所包含著豐富的信息:

(1)二維散點(diǎn)圖是兩個(gè)波段的聯(lián)合概率密度分布的幾何表達(dá)，表達(dá)了二維變量數(shù)據(jù)空間的聚類(lèi)結(jié)構(gòu);

(2)二維散點(diǎn)圖可以作為主成分分析的幾何解釋?zhuān)渲蠵c1為第一主成分軸，Pc2為第二主成分軸。

(3)橫坐標(biāo)表示具有吸收峰的波段，縱坐標(biāo)表示具有反射峰的波段。每個(gè)散點(diǎn)位置坐標(biāo)(bx，by)。坐標(biāo)點(diǎn)位于對(duì)角線上方，必定有，反之，對(duì)角線下方的點(diǎn)，必定有于對(duì)角線上方的點(diǎn)才表示存在蝕變信息。

基于波段序偶的二維散點(diǎn)圖，下面對(duì)研究區(qū)的遙感數(shù)據(jù)進(jìn)行光譜空間結(jié)構(gòu)幾何特征分析(圖2、圖3)。

圖中箭頭直線表示背景主軸，虛線表示對(duì)角線。大氣散射影響地物光譜反射曲線(主要為ETM1、ETM4波段)，由此導(dǎo)致(ETM1、ETM3波段)光譜空間結(jié)構(gòu)發(fā)生改變(背景主軸不經(jīng)過(guò)原點(diǎn))，大氣散射影響了地物在光譜特征空間中的分布位置，但并不影響典型地物間的空間分布關(guān)系。由于鐵化巖石在ETM1、ETM2或ETM4波段具有吸收帶，而在ETM3波段具有反射峰，那么它在ETM3波段的反射率必定大于其它具有吸收帶的波段(ETM1或ETM2或ETM4波段)，因此，在ETM3、ETM1二維光譜空間中，具有鐵化的巖石散點(diǎn)位于對(duì)角線的上方，則可以從圖像數(shù)據(jù)統(tǒng)計(jì)意義上判別該研究區(qū)有鐵化信息存在。由于泥化巖石在ETM5為強(qiáng)反射，在ETM7為強(qiáng)吸收，那么它在ETM5、7組合中必然會(huì)位于對(duì)角線的上方，從統(tǒng)計(jì)意義上可以判斷存在泥化蝕變信息。

黑山口地區(qū)地表主要為裸巖，圖2大的橢圓表示該區(qū)域的巖石背景。圖2中大部分散點(diǎn)位于對(duì)角線上方，該區(qū)存在鐵化異常信息的可信度較高。游離于背景橢圓之外的綠線橢圓位于對(duì)角線上方，為干擾或蝕變信息，滿足信息提取的基本條件。圖3中對(duì)角線與背景主軸幾乎重合，圖像左下角有散點(diǎn)分布在背景大橢圓周邊，遠(yuǎn)離背景橢圓中心，并且存在于對(duì)角線上方，可以推斷ETM5、ETM7波段含有泥化蝕變信息。

2．3 回歸偏度分析

對(duì)于遙感影像的參數(shù)與圖形分析，本文主要采用序結(jié)構(gòu)分析技術(shù)(張遠(yuǎn)飛，2009，2010)，序結(jié)構(gòu)分析是對(duì)序列數(shù)據(jù)的結(jié)構(gòu)形態(tài)特征與變化規(guī)律的研究。ETM七個(gè)波段數(shù)據(jù)可以構(gòu)成波段序列，波段之間的相似性大小主要取決于兩個(gè)波段的相近程度和對(duì)地物光譜響應(yīng)的差異程度。波長(zhǎng)越相近，則越相似，對(duì)地物光譜響應(yīng)差異程度越小則越相似，這種相似性在波段之間是漸變的。回歸偏度分析是序結(jié)構(gòu)分析技術(shù)的方法之一?；貧w偏度計(jì)算是基于兩波段，兩波段回歸散點(diǎn)圖的回歸軸兩側(cè)圖形是不對(duì)稱(chēng)的，回歸偏度值也是對(duì)這種不對(duì)稱(chēng)性的數(shù)學(xué)描述。從理論上講，這種兩側(cè)的不對(duì)稱(chēng)性是由于異常信息的存在引起的。由于異常信息存在而得到的回歸偏度值在理論上有正常范圍，且不同地域的回歸偏度值正常范圍不一，不同蝕變類(lèi)型的回歸偏度值正常范圍也有所不同。若兩個(gè)波段(比如ETM6和ETM7)不相關(guān)，計(jì)算出來(lái)的回歸偏度太大，這種大的回歸偏度值就不是因?yàn)楫惓Ｐ畔⒌拇嬖诙鴮?dǎo)致的。

回歸偏度值太高可能是不相關(guān)的體現(xiàn)。

設(shè)有線性回歸方程:

則實(shí)測(cè)數(shù)據(jù)中^yi的與回歸值^yi之間的偏差(yi－^yi)的平方和為

n為觀測(cè)點(diǎn)數(shù)，Q值改寫(xiě)成如下形式

Q值稱(chēng)之為回歸偏度，回歸偏度較大預(yù)示著蝕變信息存在。

2．4 蝕變信息客觀存在性評(píng)

對(duì)研究區(qū)波段數(shù)據(jù)作回歸偏度分析來(lái)選擇信息提取的特征波段，表1為六個(gè)波段的回歸偏度矩陣。ETM6波段屬于熱紅外(10．4～12．5μm)，與其他波段呈負(fù)相關(guān)，所以不對(duì)其進(jìn)行回歸偏度分析?；貧w偏度矩陣中，0值表示某波段和自身的回歸偏度。

表1 TM123457波段回歸偏度矩陣Table1 Bands(TM123457)partial regression matrix

ETM3、ETM5、ETM7三波段的回歸偏度較大，此三波段集中了主要光譜信息。由地物波譜的物理機(jī)制可知，ETM1和ETM3是兩個(gè)重要的鐵化蝕變特征波譜帶。ETM3的偏度較大，表明該地區(qū)的三價(jià)鐵離子在該波段存的反射特征很明顯，信息豐富，預(yù)示鐵化信息的存在。ETM3和ETM7的回歸偏度最高。ETM5和ETM7是羥基礦物蝕變特征波譜帶，從圖中可知，ETM5、ETM7這兩個(gè)波段偏度較大，說(shuō)明它們載荷了較多的泥化信息。本區(qū)的特征波段選擇方案如下:

(1)鐵化蝕變異常信息主要在ETM1、ETM3波段;

(2)泥化蝕變異常信息主要在ETM5、ETM7波段。

3 遙感蝕變信息提取

3．1 主成分分析

多光譜各波段是高度相關(guān)的，相鄰波段相關(guān)系數(shù)更高。因此，要去除一些相關(guān)信息。主成分分析(K－L變換)就是在一定程度上保證數(shù)據(jù)質(zhì)量，消除數(shù)據(jù)冗余，將相關(guān)的波段轉(zhuǎn)化成不相關(guān)的幾個(gè)主成份。將研究區(qū)數(shù)據(jù)作主成分分析得到:

對(duì)ETM1、ETM3、ETM4、ETM5做主成分分析(見(jiàn)表2)。一般而言，含鐵礦物在第1波段波譜具有吸收特性，在第3波段具有反射特性(張建國(guó)，2005)。表1中PC1和PC2載荷較多的背景信息。PC3、PC4主分量中ETM1和ETM3中符號(hào)相反。PC4中ETM3的絕對(duì)值大于PC3，因此可認(rèn)為PC4中載有較高的鐵的反射率，可以看出鐵化異常信息位于第四主分量，最后選擇PC4做最優(yōu)密度分割，提取鐵化信息。同樣選取ETM1、ETM4、ETM5、ETM7做主成分分析(見(jiàn)表3)，PC1和PC2載荷較多的背景信息。PC4中ETM5、ETM7波段符號(hào)相反，反射大于吸收，對(duì)PC4做最優(yōu)密度分割提取泥化蝕變信息。

表2 1、3、4、5波段主成分分析表Table2 Principal component analysis for bands 1，3，4，and 5

表3 1、4、5、7波段主成分分析表Table3 Principal component analysis for bands 1，4，5，and 7

3．2 最優(yōu)密度分割

最優(yōu)密度分割方法是將數(shù)學(xué)地質(zhì)中有序地質(zhì)量的最優(yōu)分割法引入到遙感圖像的彩色密度分割中，并加以變化改進(jìn)(吳德文，2002)。密度分割是遙感信息分級(jí)的主要方法之一，該方法類(lèi)似于多級(jí)閾值法，即根據(jù)圖像的灰度值及其概率分布(直方圖)特征選擇幾個(gè)分割點(diǎn)，將灰度值分成幾個(gè)級(jí)別，分別代表不同的地物覆蓋類(lèi)型(張遠(yuǎn)飛，2010)。

對(duì)主成分分析的PC4做最優(yōu)密度分割。當(dāng)分割段數(shù)達(dá)到12后，曲線趨于平衡，因此選12為合理分割段數(shù)，最優(yōu)12段分割區(qū)間。分析對(duì)比認(rèn)為第1、2、3分割段為鐵化蝕變信息異常區(qū)，分別賦予紅綠藍(lán)三種顏色，該灰度段以外的灰度級(jí)作為背景處理，為消除噪聲影響對(duì)其進(jìn)行平滑濾波，而后將提取的鐵化信息與ETM1、ETM2、ETM3波段圖像進(jìn)行疊加，得到圖4、圖5鐵化、泥化蝕變結(jié)果。

4 野外實(shí)測(cè)光譜驗(yàn)證

設(shè)計(jì)13條穿越工作區(qū)地層的主干剖面線，采集野外巖石標(biāo)本，利用美國(guó)ASD(350－2500nm)光譜儀對(duì)標(biāo)本巖石進(jìn)行光譜測(cè)試。根據(jù)常見(jiàn)礦物的近紅外吸收特征，400～1100nm為Fe吸收譜帶，2300～2500nm為含羥基和碳酸鹽礦物的特征譜帶(童慶禧等，1999;甘甫平，2003)。

在圖4和圖5鐵化、泥化蝕變區(qū)里實(shí)測(cè)的光譜曲線中絕大部分都存在光譜特征吸收，樣品也有相應(yīng)的礦化蝕變。本區(qū)的鐵化蝕變主要在北部鉀質(zhì)花崗巖和東南方向的侵入巖體中心部位(肉紅色黑云母花崗巖)，巖體風(fēng)化強(qiáng)烈;干溝里也存在鐵化蝕變。結(jié)合野外光譜測(cè)試，提取的鐵化蝕變區(qū)主要發(fā)生鐵染，地表未見(jiàn)黃鐵礦化，干溝中的洪積物含有鐵離子。該區(qū)鐵化對(duì)成礦作用意義不大。該區(qū)泥化蝕變主要有綠泥石化、蛇紋石化、綠簾石化、碳酸巖化，偶見(jiàn)滑石化，并沿黃山斷裂帶分布，提取的泥化蝕變信息與實(shí)際情況相符。在南部的大巖體周邊，靠近斷裂帶處有綠泥石化蝕變分布，熱液活動(dòng)強(qiáng)烈;東部的花崗巖見(jiàn)有泥化，酸性巖體的泥化蝕變對(duì)銅鎳硫化物礦床成礦意義不大;北部的硅質(zhì)凝灰?guī)r泥化蝕變也比較強(qiáng)烈，成片狀分布;礦床附近的灰色閃長(zhǎng)巖和輝長(zhǎng)巖見(jiàn)有蛇紋石化、綠泥石化，對(duì)成礦有利。本次所提取的蝕變信息與實(shí)測(cè)蝕變區(qū)吻合性好。

5 結(jié)論

本文基于回歸偏度理論對(duì)蝕變信息客觀存在性進(jìn)行評(píng)估，準(zhǔn)確地鎖定含蝕變信息的波段，對(duì)含蝕變信息波段做主成分分析，最終利用最優(yōu)密度分割法提取出蝕變信息。

(1)在波段序結(jié)構(gòu)理論指導(dǎo)下，通過(guò)回歸偏度分析來(lái)理解光譜特征信息，這對(duì)遙感蝕變信息提取具有重要意義。通過(guò)回歸偏度分析確定蝕變信息潛在波段，對(duì)蝕變信息提取工作具有較高的實(shí)用價(jià)值。

(2)光譜曲線顯示的波譜反射吸收特征與主成分分析中的特征向量矩陣表示的蝕變信息機(jī)理存在一致性，將實(shí)測(cè)光譜數(shù)據(jù)與蝕變信息提取結(jié)果對(duì)比，吻合性好。

(3)超鐵鎂質(zhì)巖體的礦化蝕變與銅鎳硫化物礦床成礦關(guān)系密切，次級(jí)斷裂帶為熱液上升提供通道，位于次級(jí)斷裂構(gòu)造帶上的輝長(zhǎng)巖泥化蝕變對(duì)成礦有利。在以后工作中應(yīng)該注意次級(jí)斷裂構(gòu)造帶上鐵鎂質(zhì)巖體的蛇紋石化、綠泥石化、碳酸巖化信息。

Chen Song－ling，Lu Fu－h(huán)ong，Gao Guang－ming．2001．Inner Mongolia section gold wall－rock alteration of remote sensing recognition［J］．Remote Sensing for Land and Resources，48(02):13－18(in Chinese)

Fu bi－h(huán)ong，Chou Xiao－we．1995．Hydrocarbons of spectral reflectance characteristics and its significance in oil and gas exploration［J］．Remote Sensing Technology and Application，10(2):20－22(in Chinese)

Fu Piao－er．2009．Metallogenic magmatism of Huangshan Cu－Ni sulfide deposit in Xinjiang［J］GeoChimica，9:432－448(in Chinese)

Gan Fu－ping，Wang Run－sheng，Ma Ai－nai．2003．Spectral identification tree(SIT)for mineral extraction based on spectral characterristics of minerals［J］．Earth Science Frontier，10(2):445－454 (in Chinese with English abstract)

Hu Pei－qing．2010．Petrogenetic and ore－forming processes of Huanshangdong Cu－Ni sulfide depositin Hami，Xinjiang［J］．Mineral Deposits，2:0258－7106(in Chinese with English abstract)

Jing Feng，Chen Jian－ping．2005．Mineralized alteration of remote sensing information extraction method and reviewed in this paper［J］．remote sensing information，45(2):62－65(in Chinese with English abstract)

Li Hong，Zhu Gu－chang．2010．Pilot study on reflectance spectral features of objects in mineralization areas［J］．Remote Sensing for Land and Resources，83(1):89－95(in Chinese with English abstract)

Liu Qing－sheng，Yan Shou－xun，Ma Chao－fei．1999．Potassium remote sensing information extraction of Inner Mongolia Ha damen gold mining area［J］．Remote Sensing Technology and Application，14(3):7－11(in Chinese with English abstract)

Liu Yan－jun．1991．Remote sensing the principle and method of prospecting［M］．Beijing:Metallurgical Industry Press:3－7(in Chinese)

Meng Guang－lu．2008．Studies of Cu－Ni Sulfide Mineralization in the Huangshandong mafic－ultramafic intrusion，Hami，Xinjiang［D］．Lanzhou:Lanzhou Univesity:1－150(in Chinese)

Shen Li－xia，Liu Li－ping，Su Xin－xu．2008．Different vegetation coverage areas remote sensing mineralized alteration extraction research［J］．Geoscience，22(2):293－298(in Chinese with English abstract)

Tong Qing－xi，Tian Guo－liang．1990．China typical features and the spectrum analysis［M］．Beijing:Science Press:1－50(in Chinese)

Wu De－wen，Zhang Yuan－fei，Zhu Gu－chang．2002．Remote sensing rock information exstraction basing on optimal density segmentation［J］．Remote Sensing For Land＆Resources，(4):51－54(in Chinese)

Zhang Jian－guo．2005．Remote sensing spectrum information extraction and metallogenic metallogenic prediction［D］．Changsha:Centrel South University:1－50(in Chinese)

Zhang Yuan－fei，Wu De－wen，Zhu Gu－chang．2008．Remote－sensing detection of alteration information background and interference studies［J］．Remote Sensing For Land＆Resources，(2):22－26 (in Chinese)

Zhang Yuan－fei，Wu Jian－sheng．1999．Based on remote sensing image extraction Alteration Information［J］．Non－ferrous metals mining and exploration，8(6):604－606(in Chinese)

Zhang Yuan－fei，Yang Zi－an，Zhu Gu－chang．2009．Remote－sensing image information on detection of altered spatial structure of spectral data analysis［J］．Remote Sensing Information，(1):3－9(in Chinese)

Zhang Yuan－fei，Zhu Gu－chang，Wu De－wen．2007．Geological survey of remote sensing alteration mineral multi－level separation and extractionof information technology and application［J］．China Association of Remote Sensing Applications，Annual Meeting:(in Chinese)

Zhang Yuan－fei，Wu De－wen，Zhang Yin－in，Zhu Gu－chang，Li Hong．2010．High－band spectral data analysis and application of coherent structures［J］．Remote Sensing For Land＆Resources，21 (1):30－38(in Chinese with English abstract)

Zhang Zong－gui，Wang Yun－sheng，Guo Xiao－fang．2003．Based on the spectral characteristics of spectral imaging remote sensing mineral recognition［J］．Earth Science Frontiers，10(2):437－443(in Chinese with English abstract)

Zhao Peng－da，Cheng Yong－qin．1999．Theory and practice of geoanomaly in mineral exploration［M］．Wuhan:China University of Geosciences Press:1－150(in Chinese)

Zhao Yuan－h(huán)ong．1991．The application of principal component intergration of band ration to extracting hydrothermal alteration information［J］．Remote Sensing for Land＆Resources，(3):83－87(in Chinese)

［附中文參考文獻(xiàn)］

陳松嶺，盧福宏，高光明．2001．華北地臺(tái)北緣內(nèi)蒙古段金礦圍巖蝕變的遙感識(shí)別［J］，國(guó)土資源遙感，48(02):13－18

傅碧宏，丑曉偉．1995．烴類(lèi)物質(zhì)的反射光譜特征研究及其油氣勘探意義［J］．遙感技術(shù)與應(yīng)用，10(2):20－22

傅飄兒．2009．新疆黃山銅鎳硫化物礦床成礦巖漿作用過(guò)程［J］．地球化學(xué)，9:432－448

甘甫平，王潤(rùn)生，馬藹乃．2003．基于特征譜帶的高光譜遙感礦物譜系識(shí)別［J］．地學(xué)前緣，10(2):445－454

胡沛青．2010．新疆哈密黃山東銅鎳硫化物礦床成巖成礦作用［J］．礦床地質(zhì)，2:0258－7106

荊鳳，陳建平．2005．礦化蝕變信息的遙感提取方法綜述［J］．遙感信息，45(2):62－65

李紅，朱谷昌，張遠(yuǎn)飛．2010．礦化蝕變區(qū)典型地物光譜特征分析－以內(nèi)蒙古突泉縣扎魯特旗成礦帶為例［J］．國(guó)土資源遙感，83(1):89－95

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

劉燕君．1991．遙感找礦的原理和方法［M］．北京:冶金工業(yè)出版社:3－7

孟廣路．2008．新疆哈密黃山東銅鎳硫化物礦床成巖成礦作用研究［D］．蘭州:蘭州大學(xué):1－50

沈利霞，劉麗萍，蘇新旭．2008．不同植被覆蓋率地區(qū)遙感礦化蝕變提取研究［J］．現(xiàn)代地質(zhì)，22(2):293－298

童慶禧，田國(guó)良．1990．中國(guó)典型地物波譜及其特征分析［M］．北京:科學(xué)出版社:1－50

吳德文，張遠(yuǎn)飛，朱谷昌．2002．遙感圖像巖石信息提取的最優(yōu)密度分割方法［J］．國(guó)土資源遙感，(4):51－54

張建國(guó)．2005．遙感多光譜成礦信息提取與成礦預(yù)測(cè)［D］．長(zhǎng)沙:中南大學(xué):1－50

張遠(yuǎn)飛，吳德文，張艮中，朱谷昌，李紅．2010．高光譜數(shù)據(jù)的波段序結(jié)構(gòu)分析與應(yīng)用研究［J］．國(guó)土資源遙感，21(1):30－38

張遠(yuǎn)飛，吳德文，朱谷昌．2008．遙感蝕變信息檢測(cè)中背景與干擾問(wèn)題研究［J］．國(guó)土資源遙感，(2):22－26

張遠(yuǎn)飛，吳健生．1999．基于遙感圖像提取礦化蝕變信息［J］．有色金屬礦產(chǎn)與勘查，8(6):604－606

張遠(yuǎn)飛，楊自安，朱谷昌．2009．遙感圖像蝕變信息檢測(cè)中的光譜數(shù)據(jù)空間結(jié)構(gòu)分析［J］．遙感信息，(1):3－9

張遠(yuǎn)飛，朱谷昌，吳德文．2007．地質(zhì)礦產(chǎn)調(diào)查的遙感蝕變信息多層次分離提取技術(shù)與應(yīng)用［J］．中國(guó)遙感應(yīng)用協(xié)會(huì)2007年學(xué)術(shù)年會(huì)

張宗貴，王潤(rùn)生，郭小方．2003．基于地物光譜特征的成像光譜遙感礦物識(shí)別方法［J］．地學(xué)前緣，10(2):437－443

趙鵬大，陳永清．1999．地質(zhì)異常成礦預(yù)測(cè)理論與實(shí)踐［M］．北京:地質(zhì)大學(xué)出版社:1－150

趙元洪．1991．波段比值的主成分復(fù)合在熱液蝕變信息提取中的應(yīng)用［J］．國(guó)土資源遙感，(3):83－87

Selection of Characteristic Bands in Alteration Information Extraction from Remote Sensing Images for the Heishankouyong Area，Hami，Xinjiang

DUAN Zong－en1，2，ZHU Gu－chang1，2，ZHANG Si－ying1，JI Xiang－yong1，WEI Han－tao1，ZHU Li－li1
(1．Centre South University，Changsha，Hunan410083; 2．China Non－ferrous Metals Resource Geological Survey，Beijing100012)

Selection of characteristic bands is a critical process for extracting mineral alteration information from remote sensing images．Taking the Heishankou area in Hami of Xinjiang as an example，we make geometric structural analysis of spectral data space and use ETM data to study the partial regression of the 7 ETM bands，and to evaluate mineral alteration information from these data．Then we select the characteristic bands for varied alteration types to extract alteration information，and use the principal component analysis method and the optimal density separation method to extract mineral alteration information．Field geological survey and spectral measurements prove that the extracted alteration information is consistent with the real situation，indicating the ideal result of this method．And this work suggests that the mineralized alteration–bearing ferruginous rock bodies on the secondary fault zone shows a good prospecting potential．

remote sensing alteration information，planar scatter plot，partial regression，characteristic band，Heishankou

book=7,ebook=228

P627

A

0495－5331(2012)04－0823－6

2011－12－30;

2012－03－20;［責(zé)任編輯］郝情情。

國(guó)家高科技研究發(fā)展計(jì)劃(863)項(xiàng)目(編號(hào):2006AA06Z112)資助。

段宗恩(1988年—)，男，黃岡市人，中南大學(xué)，在讀碩士研究生，主要從事GIS和RS的應(yīng)用研究。E－mail:dze2008@163．com。