基于Landsat 8 OLI的遙感蝕變提取

編者按：圍巖蝕變現象是熱液礦床成礦作用發生的重要標志之一，使用遙感數據進行蝕變信息提取，能降低找礦成本，提高找礦效率。文章通過將Usgs波譜庫的鐵染礦物與羥基礦物重采樣至Landsat 8 OLI波段，采用主成分分析法對研究區進行遙感蝕變提取。通過對分析研究區蝕變信息，圈定出兩個找礦靶區，以期減少找礦工作量以及提升找礦效率，為下一步的找礦活動提供更多支持。

引言

礦產資源是一種不可再生的資源，同時也是經濟建設和社會發展的重要物質基礎，利用遙感技術進行蝕變信息提取在礦產資源勘查中已得到廣泛應用，Landsat系列數據是最常被用到的數據，常用的遙感蝕變提取圖像處理方法主要有波段彩色合成、波段比值、主成分分析、光譜角﹑像元分解等。在這些方法中，主成分分析具有實現簡單、速度快﹑效果好﹑穩健等優點，本文利用Landsat 8 OLI數據，通過主成分分析法對研究區進行鐵染蝕變與羥基蝕變信息提取，通過對研究區蝕變異常結果分析，圈定出找礦靶區，為下一步的找礦活動提供更多支持。

研究概況

研究區自然概況

研究區范圍30°10′N-30°50′N;85°29′E-86°23′E，位于西藏自治區的措勒縣和昂仁縣交界處，措勤縣屬高原丘陵型和高原寬灘型地貌，境地處羌塘高原大湖盆地帶，是高原丘陵型，高原寬灘型地貌。昂仁縣地勢由東向西逐漸抬升，山脈占據全縣總面積的五分之三。

研究區地質概況

研究區主要出露地層有拉嘎組（P1l）、年波組（E2n）、昂杰組（P1a）、帕那組（E2p）、典中組（E1d）、永珠組（C2γ）、砸勒單元（K1Zl）、沙松南單元K1S、仲青單元（K1Z）、日貢拉組（E3r）、花崗閃長巖E2γσΠ、典中組（E1d）、花崗斑巖（E2γΠ）;黑云二長花崗巖（E2ηγ）、石英斑巖（E2λΠ）、第四系洪積物沖積物洪積物（Qhal），沖積物（Qhpl）。

數據處理

研究數據來源

本文數據源使用了Landsat 8 OLI 2018-01-16影像，波譜信息如表1：

研究數據預處理

為了消除大氣中各種成分吸收散射的影響，需要對西藏措勤縣尼雄地區的遙感數據進行預處理，對OLI數據進行輻射定標、大氣校正等手段，得到最終的研究區OLI（7、5、2波段組合）影像（見圖1）。

去干擾處理

Landsat 8 OLI中，植被的提取方法為歸一化植被指數（NDVI=Band5-Band4/ Band5+Band4）得到植被掩膜;水體提取使用的是倪衡等提出的新型水體指數ⅠNDw8 =Band1-（Band6+Band7）/Band1+（Band6+Band7）;在對第四系進行剔除，進而獲得Landsat 8 OLI綜合掩膜。

蝕變提取

波譜重采樣

研究區鐵染礦物主要為褐鐵礦與赤鐵礦，羥基對應的主要為黏土類礦物如高嶺石、綠泥石，將USGS的礦物波譜重采樣至landsat 8 OLI波段得到圖2。

主成分分析（PCA）

主成分分析法主要是利用各地物相應反射率差別較大的波段進行主成分分析，通常根據需要將影像數據按波段分組進行主成分變換﹐根據含鐵離子礦物的波譜特征吸收帶對應波段的主成分特征向量的載荷進行蝕變提取。

鐵染異常礦物在Band2、Band5有兩個吸收谷，Band4、Band6波段可見較高的反射峰;黏土礦物Band2可見反射峰，Band5處可見較小的吸收谷，Band6處的反射率較高，Band7處黏土礦物均出現吸收特征，選擇Band2、Band4、Band5、Band6波段與Band2、Band5、Band6、Band7波段分別作為鐵染異常和羥基異常的主成分分析波段，計算得鐵染異常（表2）與羥基異常（表3）信息，選擇第四主分量作為鐵染異常與羥基異常的主成分分量。

結果提取

本文采用主分量分析門限法劃分蝕變異常等級，將蝕變異常等級劃分為3級，即以σ、2.5σ、3σ為閾值從高到低劃分為一級﹑二級、三級蝕變異常，通過計算得到第四分量的標準差（見表4），得到研究區蝕變信息圖（見圖3）。

通過分析鐵染蝕變異常與羥基蝕變異常進而劃分出靶區Ⅰ、靶區Ⅱ，靶區位置位于影像的南東與南西部位。找礦靶區Ⅰ中蝕變呈現南西-北東分布，羥基蝕變主要位于靶區的南西部分，蝕變形狀呈團塊狀;鐵染異常在靶區分布較均勻，形狀呈現團塊狀、條帶狀;靶區Ⅱ中蝕變在靶區中均勻分布，鐵染蝕變主要集中于靶區的東南部，呈現為塊狀，且一級蝕變集中性較好，羥基蝕變整體分布與河流走向呈現一致，羥基蝕變主要位于地勢平坦的區域，綜上將其圈定為找礦靶區Ⅱ。

結論

本文利用主成分分析法（PAC）對Landsat8 OLI數據進行蝕變異常信息的提取，主要提取出了研究區內鐵染蝕變以及羥基蝕變的異常信息與分布情況;通過對于研究區蝕變異常信息的研究，分析研究區蝕變信息的分布情況，最終圈定出兩個找礦靶區以達到減少找礦工作量以及提升找礦效率的目的，以期可以為下一步的找礦活動提供更多支持。

參考文獻：

[1]江寧.基于遙感技術的山西五臺東腰莊金礦成礦前景分析[D].太原：太原理工大學，2011.

[2]唐超，周可法，張楠楠，等.基于Landsat-8 OLI和 ASTER數據集成和融合的礦化蝕變信息提取：以包古圖斑巖型銅礦為例[J].地質科技情報，2018，37（06）：211-217.

[3]汪子義，張廷斌，易桂花，等.Landsat8 OLI 數據斑巖銅礦遙感蝕變礦物組合提取研究[J].國土資源遙感，2018，30（03）：89-95.

[4]張玉君，姚佛軍.應用多光譜ASTER數據對ETM遙感異常的定性判別研究—以東昆侖五龍溝為例[J].巖石學報，2009，25（04）：963-970.

[5]賀金鑫，姜天，董永勝，等.基于Landsat 8的遼寧弓長嶺區遙感蝕變信息提取[J].吉林大學學報（地球科學版），2019，49（03）：894-902.

[6]張玉君，曾朝銘，陳薇.ETM+（TM）蝕變遙感異常提取方法研究與應用——方法選擇和技術流程[J].國土資源遙感，2003（02）：44-49，78.

[7]張遠飛，袁繼明，楊自安，等.基于物理意義的二維散點圖類型劃分與遙感蝕變信息提取[J].國土資源遙感，2013，25（02）：57-62.

[8]倪衡，劉曉穎.一種基于 Landsat 8遙感數據快速提取水體信息的新方法[J].科技資訊，2015，13（19）：7-9.

[9]王守志，邢立新，仲波，等.基于Landsat-8 OLI和GF-1 PMS數據融合的鐵染蝕變信息提取[J].遙感技術與應用，2016，31（05）：950-957.

作者簡介：王洪飛（1996—），四川遂寧人，碩士，現就讀于成都理工大學。