ICA在甘肅黃草灘地區鐵染蝕變信息提取中的應用研究

付壘 趙鵬彬 于占超

摘要：本文以ICA（獨立成分分析）作為遙感蝕變信息提取的主要方法，結合傳統PCA（主成分分析）在光譜去相關性操作中的基本特性，通過互補式的ICA方法進行遙感蝕變信息的綜合提取。以甘肅黃草灘矽卡巖型磁鐵礦為研究區，使用ICA和PCA方法進行ETM+影像鐵染蝕變信息提取，通過對提取結果進行比較分析，結合野外實地查證，評估了兩種方法在特定地質背景下遙感蝕變信息提取的適應性和局限性，最后對遙感蝕變信息提取結果進行了歸納總結，結果表明：互補式的ICA能很好地反映研究區的蝕變異常分布狀況，可以將該方法應用到此類地區的地質找礦中，同時為進一步分析研究區成礦模式提供科學依據。

關鍵字：獨立成分分析；蝕變信息；矽卡巖型；地質找礦

Abstract： This paper using independent component analysis （ICA） as main method for extraction of alteration information from remote sensing data，combined the basic characteristics of traditional PCA （principal component analysis） in removing spectrum correlation，comprehensive extraction of alteration information from remote sensing data by complementary type ICA. Taking skarn magnetite in Huangcaotan of Gansu province as an example，using ICA and PCA method of ETM+ image for extraction of alteration information，through the results of a comparative analysis of the extraction，combined with field verification and evaluation of the two methods under the specific geological background in remote sensing alteration information from the adaptability and limitations. Finally，the remote sensing alteration information extraction results are summarized in the paper. The results showed that the complementary type of ICA can reflect the distribution of alteration anomaly primely，the method can be applied to such areas of geological prospecting，and provide a scientific basis for further analyzing and studying metallogenic model.

Key words：Independent Component Analysis；Alteration information extraction；Skarn；Geological prospecting

1. 引言

由于不同礦物的晶體結構各不相同，晶格振動所產生的光譜特征與其特有的晶體結構有關，因此不同的巖石具有不同的光譜特征。遙感技術可以識別這些光譜特征，從而判別出地物的礦物組成和巖性，這就為通過遙感技術識別蝕變礦物分布范圍提供了物理基礎[1]。目前遙感蝕變異常信息提取的方法有多種，較常用的有比值法、主成分分析法、光譜角法、歸一化植被指數法、混合像元分解法以及它們的混合法。每種方法都有其針對性和適用的范圍。其中主成分分析法提取蝕變信息是相對最為廣泛的。主成分分析（PCA：Principal Component Analysis）是基于信號二階統計特性的分析方法，由于所獲各主成分之間不相關，主成分之間信息沒有重復或冗余。PCA的這一基本性質在蝕變異常信息提取中被充分利用。多光譜遙感數據通過PCA所獲每一主成分常常代表一定的地質意義[2]，且互不重復，即各主成分的地質意義有其獨特性。但是由于蝕變礦物形成的影像特征在遙感圖像上往往表現得很微弱或不明顯，經常被主色調或“偽信息”所掩蓋，而且很多弱信息包含在高階統計量中。若僅僅采用PCA算法，對圖像進行處理可能會忽略很多關鍵信息。因此，本文在提取礦物蝕變信息中引進了互補式的獨立成分分析法（ICA：Independent Component Analysis）。ICA算法是一種基于高階統計量的去相關多元數據處理方法[3]，不僅能夠消除多變量數據中的二階相關關系，而且能夠消除數據間的高階相關。本文利用PCA算法在光譜去相關性操作中的基本特性，以ICA為遙感蝕變信息提取的主方法，綜合利用兩種不同方法提取的異常信息，同時結合研究區地質環境背景，提取了遙感蝕變異常信息，通過室內分析及野外驗證，取得了較好的成果。

2. 研究方法

2.1 ICA基本原理

ICA又稱盲源分離（Blind source separation，BSS）[4]，它假設觀察到的隨機信號x服從模型x=As，其中s為未知源信號，其分量相互獨立，A為一未知混合矩陣。ICA的目的是通過且僅通過觀察x來估計混合矩陣A以及源信號s。用數學語言描述，即，通過獨立成分分析，從一組觀察信號分離出一組統計獨立的信號源的估計值，也是統計獨立的。

假定第個觀測信號由個獨立成分線性混合而成：

（1）

用矢量表示觀測變量，矢量表示源變量，表示混合矩陣式（1）可表示成：

（2）

統計模型式（2）稱為獨立成分分析模型。該模型描述了觀測數據是如何由信息源混合生成的。源變量是隱藏變量，不能直接觀測到，而且混合矩陣也是未知。所有能觀測到的數據只有隨機變量，所以必須估計出混合矩陣和。

獨立成分分析模型的求解基于一個簡單的假設：源變量是統計獨立的，同時非高斯分布。在估計出混合矩陣后，可以進一步計算混合矩陣的逆，也即分離矩陣，從而得到獨立成分的估計。

（3）

ICA的目的就是估計出分離矩陣。目前獨立成分分析的實現方法有很多：信息最大化法、基于極大似然估計的ICA算法[5]、最小化互信息（負熵最大化）法[6]、FastICA算法[7]等。本文采用的是FastICA算法來做獨立成分分析。

遙感影像所表現的信息是不同地物對某個波段反射率的一種反映[8]，多光譜遙感影像的各波段可以認為是相互獨立的各種地物的光譜信息按一定的規律混合的結果。將ICA作用于多光譜影像，得到的每個獨立成分可集中體現某些地物的信息，使得不同地物間的可識別度增大。

2.2 技術方法

本文采用ICA方法提取波段信息中的高階信息和獨立成分，以獨立成分替代主成分提取鐵染蝕變信息。經ICA處理后的各成分相對于PCA的各成分來說是獨立的而且包含了一些基于高階統計的信息。從統計分析的角度看，簡單ICA和PCA一樣，同屬多變量數據分析的線性方法。因此可以按照PCA提取蝕變信息的準則來處理經ICA處理的成分。以TM數據為例，通過Crosta方法[9]說明其準則，通過TM1、TM3、TM4、TM5的波段組合進行主成分分析提取鐵染蝕變。由TM1、TM3、TM4、TM5作為輸入波段組合進行正交線性變換，變換后的某個新的組分圖像可能集中了鐵染蝕變信息[10]。對代表鐵染蝕變的主成分的判斷準則是：TM3的系數應與TM1、TM4的系數相反。鐵染蝕變存在于絕大多數成礦巖體中，提取這種蝕變信息基本可以確定研究區成礦巖石的分布情況[11]。

所以，本文的研究方法就是通過ETM+1、ETM+3、ETM+4、ETM+5波段組合進行ICA提取鐵染蝕變信息。為了充分發揮ICA提取結果的優越性，同時進行PCA處理，并對兩種方法提取的結果進行分析總結，結合研究區地質背景，形成研究區的蝕變信息分布圖。

基于ICA提取礦物蝕變信息的方法主要包括：圖像預處理、ICA分析、選擇有效成分、異常切割、蝕變信息處理等步驟[12]。原始數據經過一系列的圖像預處理，包括幾何糾正，大氣校正，掩膜去背景等；再選擇波段分別進行ICA和PCA處理，根據上述原則進行成分的選擇；將選擇的成分用相同的異常切割標準處理，突出蝕變信息；并與已知的礦化蝕變帶進行疊加分析，利用PCA提取的蝕變信息結果來強化和完善ICA提取的蝕變結果。

在做異常切割時，利用（X+kσ）確定異常下限和劃分異常強度等級[13]。X是某一成分的統計均值代表區域背景，σ是該成分的標準差。k取值一般是（1～3），有了這一標準，切割異常時可以減少主觀任意性，并使操作較為規范化。經ICA處理后得到的各成分是相互獨立的，每一成分可以看成是單一的地物。但是由于大氣，傳感器等影響，各成分并不只是單一的地物，所以在提取異常時還是利用了閾值切割法。

將提取的蝕變異常分布圖與已知的礦化蝕變帶進行疊加分析，計算出在一定區域內提取的蝕變信息量大小，用吻合率來表示。吻合率越大證明提取的有效信息量越大，偽異常信息量越小。

3. 應用實例

本文以甘肅黃草灘野馬井南磁鐵礦為研究區，采用ETM+作為遙感蝕變異常提取數據源，通過ICA方法提取蝕變信息，結合PCA的提取結果，運用互補式的ICA分析法，綜合圈定研究區遙感蝕變異常區帶。本次ETM+遙感數據采集的時相為2002/08/12，行帶號為137-31。

3.1 區域地質背景

研究區北部出露地層為寒武系雙鷹山組，中南部為平頭山組和少量震旦系洗腸井群，二者呈斷層接觸，平頭山組壓覆在洗腸井群之上，呈“飛來峰”。雙鷹山組巖性以深灰色硅質巖為主，局部夾灰黑色含炭質硅質巖、細砂巖、粉砂巖；平頭山組以中厚層灰白色白云巖、白云質大理巖為主，局部夾細砂巖、粉砂巖透鏡體。洗腸井群主要巖性為深灰色薄層狀硅質巖。侵入巖主要有泥盆紀雙峰山巖體，巖性主要有肉紅色中粗粒二長花崗巖、灰色細粒閃長巖、灰色～淺灰綠色中細粒閃長巖。上述地層受該巖體侵位影響，致使二者之間的接觸部位多呈港灣狀接觸，并多處形成矽卡巖帶或矽卡巖化帶，控制著矽卡巖型磁鐵礦產的分布。

研究區已知的礦床主要為鐵礦，次為錳礦、黃鐵礦、黃銅礦等，磁鐵礦點主要分布在野馬井之南的平頭山組、洗腸井群與晚泥盆世雙峰山石英閃長巖、閃長巖巖體和早泥盆世冰墩子山鉀長花崗巖熱液交代作用形成的矽卡巖化帶或矽卡巖帶。侵入巖比較發育，主要有早泥盆世冰敦子山鉀長花崗巖、花崗斑巖及晚泥盆世炭山子輝長巖等，見后期侵入的輝長輝綠巖脈（圖1）。

野馬井地區已知的礦床主要為鐵礦，次為錳礦、黃鐵礦、黃銅礦等，磁鐵礦點主要分布在野馬井之南的平頭山組、洗腸井群與晚泥盆世雙峰山石英閃長巖、閃長巖巖體和早泥盆世冰墩子山鉀長花崗巖熱液交代作用形成的矽卡巖化帶或矽卡巖帶。出露地層為青白口系園藻山群平頭山組和南華-震旦系洗腸井群，二者為逆沖斷層接觸，平頭山組壓覆在洗腸井群之上，呈“飛來峰”。平頭山組巖性主要為灰白色弱矽卡巖化中薄層白云巖、弱矽卡巖化白云質大理巖夾少量細碎屑巖；洗腸井群主要巖性為深灰色薄層狀硅質巖。侵入巖比較發育，主要有早泥盆世冰敦子山鉀長花崗巖、花崗斑巖及晚泥盆世炭山子輝長巖等，見后期侵入的輝長輝綠巖脈（圖1）。

根據研究區成礦地質背景，區內最主要的礦產是鐵礦和錳礦，鐵礦的富集使得遙感可以識別的鐵離子增多，研究區主要蝕變類型有磁鐵礦化、黃鐵礦化、黃銅礦化、綠簾石化、陽起石化、透閃石化等。因此，本文主要關注含鐵礦物蝕變的光譜特征[14]（圖2）。

3.2 鐵染信息提取

研究區鐵染蝕變較強，根據含鐵離子礦物的特征波譜，對ETM+1、3、4、5波段分別進行ICA和PCA處理，生成反映鐵離子變化的蝕變信息。由ICA提取的鐵染蝕變波段組合特征向量見表1。為了比較驗證ICA方法的準確性，并對ICA方法的完整性進行補充，輔以PCA進行特征向量的提取見表2。

各獨立成分和主成分的特征向量是趨于一致的，這說明ICA和PCA在統計分析上是一致的，由此選擇了相同的切割標準對成分進行分割；同時，也表明ICA可以和PCA一樣用于提取蝕變信息。在特征向量表IC4/PC4中，由于1和3波段的貢獻系數最大，說明它們對IC4/PC4的貢獻最大，也就是說IC4/PC4的信息主要來源于ETM+1、3波段，并且ETM+3為負值，ETM+1、ETM+4為正值，具有ETM+3與ETM+1、4相反的貢獻標志，正符合鐵染信息在ETM+1波段具有強吸收，在ETM+3具有強反射的特點。因此該成分可以作為鐵染信息的指示分量。為了體現ICA方法在蝕變信息提取中的準確性和差異性，利用PCA提取蝕變結果對ICA進行驗證和補充，首先保證兩種方法提取的礦物蝕變異常信息的標準一致，將得到的IC4和PC4進行同一水平切割。根據研究區已有的礦化信息和蝕變帶來確定閾值。本次試驗采用（X+2.5σ）對IC4和PC4進行切割，然后對切割后的二值圖像進行窗口為3×3的中值濾波，最終將提取的蝕變結果與已知的礦化蝕變帶進行疊加分析（圖3），計算兩種方法分別在對應的蝕變帶中提取的蝕變總面積及與蝕變帶吻合的面積，最后計算出兩種方法的蝕變提取吻合率（表3）。

（紫色）和PCA（綠色）提取的鐵染蝕變信息與研究區的已知蝕變帶（橙色）重疊度較高。從圖中可以發現經ICA提取的鐵染蝕變信息大部分落在已知蝕變帶內及該蝕變帶附近，也有部分異常信息落在非蝕變帶區域附近，這部分零散的蝕變信息經實地驗證是偽蝕變信息。從圖3中也可以發現相對于PCA法，ICA提取的信息量和吻合率都有一定程度的提升，吻合率提升了近13%（表3）；通過分析圖3中各個已知蝕變帶的異常信息分布圖，發現經PCA提取的蝕變信息在4號蝕變帶附近只有部分異常，其中2號、6號蝕變帶周邊完全沒有PCA提取的信息，表3的結果正符合圖3所顯示的結果。從表3可以看出總體上經ICA提取的鐵染蝕變信息的吻合率大于PCA法，但是5號蝕變帶相反，ICA在5號蝕變帶附近提取的蝕變信息較弱，而PCA在該蝕變帶有較好的體現，這說明針對不同的區域、不同礦化蝕變類型，ICA法具有一定的局限性，為了更好地總結研究區礦化蝕變異常的提取方法，需要借鑒PCA的提取方法，綜合利用兩種方法的提取結果，采用一種互補式的ICA方法，全面準確地提取研究區含鐵礦物的蝕變信息。

3.3 野外查證及分析

本次工作在研究區共查證7處礦點，包括2個沉積型磁鐵礦點、5個矽卡巖型磁鐵礦礦點。

5個矽卡巖型磁鐵礦區的ICA異常都比較明顯，提取的蝕變帶呈團塊狀分布，特別是在巖體與地層的接觸帶上，ICA異常比較集中，例如在6號蝕變帶的矽卡巖化帶兩側均有ICA異常展布。而PCA異常在5個矽卡巖型磁鐵礦區的吻合率就相對較低，特別是在2、6號蝕變帶中基本沒有提取出相關蝕變（表3）。在2個沉積型磁鐵礦點的查證中，例如5號蝕變帶附近就少見ICA法提取的異常信息，而大部分異常信息為PCA法提取的（圖3），這說明了ICA法在不同成因類型的礦化信息提取中存在的局限性。在對研究區北部韌性剪切帶兩側的PCA異常進行查證中，并沒有發現蝕變礦物，說明在同一標準下，PCA提取的蝕變異常信息中含有更多的假異常信息，產生這些假異常信息主要是由于在干旱地區普遍存在“荒漠漆”現象造成假象褐鐵礦化蝕變異常所致，既由于氣候干燥，屬于構造破碎帶的巖石常年裸露，風化作用使巖石內部的Fe2+轉變成為Fe3+積聚在表面而形成，野外查證中較易發現這種鐵染異常假象。由此可以說明經ICA法提取的有效蝕變信息量要高于PCA法。

通過分析總結：對于矽卡巖型的磁鐵礦，熱液交代作用形成的圍巖蝕變信息較弱，空間展布不規則，干擾因素也較多，遙感光譜反射形成的可識別信號往往包含于高度相關的信息中，而ICA通過去相關操作，理論上保留了一些弱化的高階信息，從而使得在此類地質背景下的弱化信息得到了加強，提取效果就要優于PCA法。而對于沉積型的磁鐵礦，其含礦層位、礦體寬度、厚度及延伸都較為穩定，主要受地層和構造控制，圍巖蝕變相對集中，因而傳統PCA法提取的效果就相對較好。

ICA法作為本次遙感蝕變信息提取的主要方法，采用互補式的ICA法進行遙感蝕變帶的圈定，既能提高ICA法的準確性，發揮其它方法（PCA法）的特殊優勢，又能提高ICA法在蝕變信息提取中的全面性，為進一步把握研究區的找礦方向奠定了理論基礎。

4. 結語

研究區地處半干旱戈壁荒漠區，蝕變礦物主要為巖漿熱液交代形成的矽卡巖型磁鐵礦。結合研究區地質背景，通過遙感技術方法，在獲取特征波段之間的獨立成分前提下，利用ICA方法提取隱藏在高階統計信息中的有用信號（特定的礦物蝕變信息），更加全面的揭示了多光譜數據間的本質結構。

在借鑒了傳統PCA法的提取結果后，采用互補式的ICA方法提取礦物蝕變信息，經實地查證取得了較好的效果，在統一標準下比傳統的PCA方法效果要好，經ICA提取蝕變信息的吻合率相比傳統PCA法要高近13%，這也充分說明了ICA法在研究區提取蝕變的準確性和適用性。這一方法在類似區域開展遙感蝕變信息提取工作具有一定的指導意義。

將互補式的ICA方法用于提取遙感蝕變信息，以下幾點需要不斷研究與探索：

（1）在對選擇的成分進行分割時，閾值的選擇具有一定的局限性，根據已知的地質資料并不能完全滿足要求，而且閾值的確定對提取的結果有很大的影響，如何實現快速確定閾值有待進一步的研究。

（2）表3中的5號蝕變帶PCA法提取信息的吻合率要高于ICA法，說明了ICA方法也有一定的適用范圍，并不是所有經ICA提取的蝕變信息就一定比PCA法提取的有效信息量大，這可能與閾值的選擇及研究區地質背景有關。

（3）蝕變信息的提取，要盡可能地采用多種不同的方法進行對比研究，找到一種適合該區域的方法或幾種方法的組合來提高蝕變信息提取結果的可信度。將ICA方法與其它方法組合可能會取得更好的應用效果。

參考文獻：（References）：

[1] Zhang Yujun，Yang，Jianmin，Chenwei.A Study of the method for Extraction of Alteration Anomalies from the ETM+（TM） Data and Its Application： Geologic Basis and Spectral Preconditon[J].Remote Sensing for Land & Resources，2002，54（4）：30-36.張玉君，楊建民，陳薇.ETM+（TM）蝕變遙感異常提取方法研究與應用—地質依據和波譜前提[J].國土資源遙感 2002，54（4）：30-36.

[2] Zhang Nannan，Zhang Xiaofan，Zhou Kefa，et al.Extraction method of alteration information in arid area from TM image[J].Arid Land Geography，2006，29（3）：398-402.張楠楠，張曉帆，周可法等.干旱區TM 圖像蝕變信息提取方法研究[J].干旱區地理.2006，29（3）：398-402.

[3] Zhong Jiaqiang，Wang Runsheng.Multitemporal Remote Sensing Images Change Detection Based on ICA[J].Journal of Electronics & Information Technology，2006，28（6）：994-998.鐘家強，王潤生.基于獨立成分分析的多時相遙感圖像變化檢測[J].電子與信息學報，2006，28（6）：994-998.

[4] COMON P.Independent component analysis，a new concept？[J].Signal Processing，1994，36（3）：287-314.

[5] Feng Yan，He Mingyi，Song Jianghong，et a1.ICA-Based Dimensionality Reduction and Compression of Hyperspectral Images[J].Journal of Electronics & Information Technology，2007，29（12）：2871-2875.馮燕，何明一，宋江紅，等.基于獨立成分分析的高光譜圖像數據降維及壓縮[J].電子與信息學報，2007，29（12）：2871-2875.

[6] Zeng Shenggen，Wang Xiaomin，Fan Ruibin，et a1.Remote Image Classification Based on Independent Component Analysis[J].Journal of Remote Sensing，2004，8（2）：150-157.曾生根，王小敏，范瑞彬，等.基于獨立分量分析的遙感圖像分類技術[J].遙感學報，2004，8（2）：150-157.

[7] Wu wei，Peng hua，Zhang fan.Performance analysis of FastICA and RobustICA on blind sources sepatation[J].Application Research of Computer，2014，31（1）：95-98.吳微，彭華，張帆.FastICA和RobustICA算法在盲源分離中的性能分析[J].計算機應用研究，2014，31（1）：95-98.

[8] Hyvrinen A.Fast and robust fixed-point algorithms for independent component analysis[J].IEEE Trans.on Neural Networks，1999，10（3）：626-634：

[9] Chen Jiangping，Wang Qian，Dong Qingji，et al.Extraction of Remote Sensing Alteration Information in Tuotuohe，Qinghai Province[J].Earth Science Journal of China University of Geosciences，2009，34（2）：314-318.陳建平，王倩，董慶吉，等.青海沱沱河地區遙感蝕變提取[J].地球科學—中國地質大學學報，2009，34（2）：314-318.

[10] Li Jianguo，Mao Debao.The Methods for Extraction Anomalies based on the ETM+ and Aster Data-A Case Study of the Mandula Area[J].Geological Survey and Research，2007，30（3）：234-240.李建國，毛德寶.基于ETM+與ASTER數據的礦化蝕變信息提取-以滿都拉地區為例[J].地質調查與研究，2007，30（3）：234-240.

[11] Liu Cheng，Wang Danli，Li Xiaomei.Extracting Clay Alteration Information of Medium Vegetation Covered Areas based on Liner Model of Spectral Mixture Analysis[J].Remote Sensing Technology and Application，2003，18（2）：95-98.劉成，王丹麗，李笑梅.用混合像元線性模型提取中等植被覆蓋區的粘土蝕變信息[J].遙感技術與應用，2003，18（2）：95-98.

[12] Su Zhicheng，Nv Hongwei.Remote Image Classification Based on Independent Component Analysis[J].Science Technology and Engineering，2007，7（23）：6244-6247.蘇志成，呂宏偉.基于獨立分量分析的遙感影像分類方法[J].科學技術與工程.2007，7（23）：6244-6247.

[13] Wang Qian，Chen Jianping.Extraction and Grading of Remote Sensing Alteration Anomaly based on the Fractal Theory[J].Geological Bulletin of China，2009，28（2-3） ：285-288.王倩，陳建平.基于分形理論的遙感蝕變異常提取與分級[J].地質通報，2009，28（2-3）：285-288.

[14] Jin Feng，Cheng Jianping.MODIS Remote Sensing Data Processing and Acquisition of Atmosphere Products[J].Remote Sensing Information，2005（2）：62-65，57.荊鳳，陳建平.礦化蝕變信息的遙感提取方法綜述[J].遙感信息，2005（2）：62-65，57.