李娜 董新豐 甘甫平 閆柏琨

(中國自然資源航空物探遙感中心，北京 100083)

資源一號02D衛(wèi)星(又稱為5米光學(xué)業(yè)務(wù)衛(wèi)星)于2019年9月12日發(fā)射升空，該星作為中國自主建造并成功運行的首顆民用高光譜業(yè)務(wù)衛(wèi)星，是國家民用空間基礎(chǔ)設(shè)施中新型對地觀測衛(wèi)星發(fā)展的又一重要成果。其高光譜載荷可實現(xiàn)166個譜段的數(shù)據(jù)同時獲取和下傳，可應(yīng)用于自然資源調(diào)查監(jiān)測要素的定量化提取和反演。

礦物識別是高光譜遙感技術(shù)優(yōu)勢之一，已在地質(zhì)調(diào)查應(yīng)用中發(fā)揮了重要作用[1-3]，但以往受其數(shù)據(jù)源限制，制約了大面積業(yè)務(wù)化應(yīng)用推廣。同時高光譜遙感以其高維數(shù)據(jù)特點，在巖性-構(gòu)造填圖具有明顯技術(shù)優(yōu)勢，已在相關(guān)應(yīng)用中取得了較好的成果[4]，資源一號02D衛(wèi)星作為高光譜業(yè)務(wù)衛(wèi)星，其高光譜數(shù)據(jù)應(yīng)用效果備受關(guān)注。

為了更好的了解該衛(wèi)星高光譜數(shù)據(jù)在地質(zhì)調(diào)查工作中的可應(yīng)用性和適宜性，本文選取哈密遙感地質(zhì)試驗場為研究區(qū)，進行數(shù)據(jù)處理、巖性-構(gòu)造信息識別、礦物信息提取并評價其應(yīng)用效果，同時為其在地質(zhì)調(diào)查領(lǐng)域中應(yīng)用推廣提供參考。

1 衛(wèi)星數(shù)據(jù)及預(yù)處理

1.1 資源一號02D高光譜數(shù)據(jù)

資源一號02D配置了可見近紅外相機和高光譜相機，其高光譜相機共獲取可見近紅外短波范圍波段數(shù)166個，其中可見近紅外76個波段，光譜范圍395～1040 nm；短波紅外90個波段，光譜范圍1005～2501 nm，高光譜載荷的具體參數(shù)見表1。

表1 資源一號02D高光譜載荷基本參數(shù)Table 1 Basic characteristics of the ZY-1-02D satellite hyperspectral payload

1.2 數(shù)據(jù)預(yù)處理

1.2.1 輻射定標(biāo)和大氣校正

輻射定標(biāo)是根據(jù)每個波段的輻射定標(biāo)系數(shù)將影像像元亮度(DN)值轉(zhuǎn)換為輻亮度值，而地面目標(biāo)的輻亮度值并不是地表真實反射率的反映，而是包含了由大氣吸收，尤其是散射作用造成的輻射量誤差。需要進行大氣校正消除這些由大氣影響造成的輻射誤差，反演地物真實的地面反射率的過程。本文利用Modtran模型，對資源一號02D高光譜數(shù)據(jù)進行大氣校正[5]。根據(jù)影像獲取時間、地理信息以及影像質(zhì)量等信息，計算獲得地面反射率數(shù)據(jù)。為檢驗大氣校正質(zhì)量，在研究區(qū)選擇已知碳酸鹽巖地質(zhì)體影像光譜與美國地質(zhì)調(diào)查局(USGS)標(biāo)準(zhǔn)波譜庫對比，如圖1所示，兩者在碳酸鹽巖特征譜帶范圍內(nèi)譜形吻合度高，且特征吸收位置一致(2330 nm附近)，可滿足后續(xù)應(yīng)用要求。

圖1 碳酸鹽巖地質(zhì)體的影像光譜(紅色)與USGS標(biāo)準(zhǔn)波譜(黑色)對比Fig.1 Comparison of image spectrum (red) and USGS standard spectrum (black) of carbonate geological body

1.2.2 信息增強

信息增強可將資源一號02D高光譜數(shù)據(jù)166個波段中包含的重要地物信息集中到盡可能少的波段圖像中，有效地調(diào)節(jié)數(shù)據(jù)中的噪聲并降低波段間的相關(guān)性，可在巖性-構(gòu)造信息識別過程中突出不易辨別的細節(jié)信息。

最小噪聲分離(MNF)是高光譜數(shù)據(jù)降維與重建的主要手段，可達到數(shù)據(jù)增強的目的，其本質(zhì)上是兩次層疊的主成分變換[6]。第一次變換用于分離和重新調(diào)節(jié)數(shù)據(jù)中的噪聲，這一操作使得變換后的噪聲數(shù)據(jù)只有最小的方差且沒有波段間的相關(guān)。第二次變換是對噪聲白化數(shù)據(jù)的標(biāo)準(zhǔn)主成分變換，計算公式如下。

(1)利用高通濾波器模板對整幅影像或具有同一性質(zhì)的影像數(shù)據(jù)進行濾波處理，得到噪聲協(xié)方差矩陣CN，將其對角化為矩陣

DN=UTCNU

(1)

式中：DN為CN的特征值按照降序排列的對角矩陣；U為由特征向量組成的正交矩陣，UT為U的逆矩陣。進一步變換式(1)可得

I=PTCNP

(2)

式中：I為單位矩陣；P為變換矩陣，P=UDN-1/2。當(dāng)P應(yīng)用于影像數(shù)據(jù)X時，通過Y=PX變換，將原始影像投影到新的空間，產(chǎn)生的變換數(shù)據(jù)中的噪聲具有單位方差，且波段間不相關(guān)。其中，X為原始影像，Y為經(jīng)過變換后投影到新的空間的影像。

(2)對噪聲數(shù)據(jù)進行標(biāo)準(zhǔn)主成分變換，公式為

CD-adj=PTCDP

(3)

式中：CD為影像X的協(xié)方差矩陣；CD-adj為經(jīng)過P變換后的矩陣，進一步將其對角化為矩陣

DD-adj=VTCD-adjV

(4)

式中：DD-adj為CD-adj的特征值按照降序排列的對角矩陣；V為由特征向量組成的正交矩陣。通過以上2個步驟得到MNF的變換矩陣TMNF=PV。

MNF變換克服了噪聲對影像質(zhì)量的影響，且波段間不相關(guān)，因此采用MNF變換后不同波段間組合，對資源一號02D高光譜數(shù)據(jù)進行信息增強處理，以突出影像中的細節(jié)信息。

2 研究區(qū)概況

哈密遙感地質(zhì)試驗場位于新疆東天山哈密市東南，距離哈密市區(qū)約140 km，海拔高程500～1500 m左右，屬丘陵、沙漠、礫石戈壁地貌。區(qū)內(nèi)基巖裸露好，出露地層可分為結(jié)晶基底和蓋層兩部分，結(jié)晶基底為薊縣系地層，在研究區(qū)東南部有少量出露，可分為兩個巖性段，下段巖性為大理巖、白云巖、白云質(zhì)大理巖，上段巖性為黑云母片巖，絹云母綠泥石英片巖。蓋層出露有古生界泥盆系、石炭系、二疊系、中生界侏羅系以及新生界，其中以石炭系出露最廣，呈近東西分布于康古爾塔格—黃山斷裂兩側(cè)，為一套海相火山沉積—正常沉積巖石建造，其余地層均為零星分布。區(qū)內(nèi)經(jīng)歷多期構(gòu)造變形，巖漿活動強烈，皺褶與斷裂系統(tǒng)極為發(fā)育。皺褶以區(qū)域性緊閉線狀皺褶為主；斷裂系統(tǒng)有東西向、北東向、北西向、北東向及近南北向五組，以近東西向及北西向的低序次斷裂構(gòu)造為主，為康古爾深大斷裂、苦水大斷裂及派生的次級斷裂。區(qū)內(nèi)侵入巖較發(fā)育，從基性-超基性到中酸性均有分布。區(qū)內(nèi)礦產(chǎn)資源相對豐富，主要為銅鎳硫化物礦床、釩鈦磁鐵礦和金礦。與成礦關(guān)系密切的蝕變礦物主要有褐鐵礦、絹云母、方解石、白云石、綠泥石/綠簾石等。

3 巖性-構(gòu)造信息識別

利用資源一號02D高光譜數(shù)據(jù)信息增強圖，結(jié)合研究區(qū)地質(zhì)資料，對研究區(qū)開展巖性-構(gòu)造信息識別工作。如圖2信息增強圖和1∶50 000地質(zhì)圖對比發(fā)現(xiàn)，不同巖性之間色調(diào)差異明顯，界限清晰，信息增強圖對于區(qū)域中第四系、石炭系地層以及黑云母花崗巖、閃長巖和輝長巖巖體都能進行很好的識別，而且地層和巖體界線準(zhǔn)確，并在區(qū)域右下方(圖中黑色圓圈處)發(fā)現(xiàn)有不同色調(diào)的新巖體信息。

圖2 信息增強圖和1∶50 000地質(zhì)圖地層巖體解譯Fig.2 Information enhancement map and interpretation of stratum and rock mass in 1∶50000 geological map

研究區(qū)內(nèi)斷裂構(gòu)造多以東西和北東產(chǎn)狀發(fā)育，基本為不同地層巖體的分界面呈線狀延伸，在信息增強圖上呈明顯的線狀特征，構(gòu)造兩側(cè)色調(diào)差異明顯，紋理清晰，容易識別，并且發(fā)現(xiàn)有褶皺地層發(fā)育(如圖3(a)中黑色折線所示)，比1∶50 000地質(zhì)圖信息更加豐富。

圖3 信息增強圖和1∶50 000地質(zhì)圖斷裂構(gòu)造解譯Fig.3 Information enhancement map and 1∶50000 geological map for fault structure interpretation

基于以上解譯原理，對研究區(qū)巖性-構(gòu)造信息進行綜合識別，如圖4所示，雖然資源一號02D高光譜數(shù)據(jù)空間分辨率為30 m，但經(jīng)信息增強后，其光譜維優(yōu)勢能夠反應(yīng)詳細的地層、巖體及構(gòu)造信息，可為今后1∶50 000地質(zhì)填圖提供重要技術(shù)支撐。

圖4 巖性-構(gòu)造綜合解譯圖Fig.4 Comprehensive interpretation of lithology and structure

4 礦物信息提取

高光譜礦物信息提取方法主要有光譜匹配、光譜特征參量以及光譜解混等三大類方法。其中光譜匹配方法常見的有距離法[7]、光譜角法[8]、光譜信息散度[9]、匹配濾波、混合調(diào)制匹配濾波[10]等，光譜特征參量代表性的方法有光譜特征擬合[11]、光譜吸收指數(shù)[12]、吸收譜帶定位分析等。光譜解混主要是利用線性或非線性解混方法，對礦物含量進行反演[13]。但在實踐應(yīng)用中，上述方法均存在不能滿足礦物精細識別以及準(zhǔn)確度較低的情況。為此文獻[14]在總結(jié)分析了已有礦物信息提取方法優(yōu)缺性的基礎(chǔ)上，提出了光譜特征增強匹配度算法。文獻[15]依據(jù)分層譜系思路，提出了綜合光譜特征增強匹配度和特征參量的高光譜礦物信息提取方法，取得了較好的礦物識別效果。

本文同樣基于礦物分層譜系識別思想，采用綜合光譜特征增強匹配度和特征參量信息提取方法，對研究區(qū)礦物信息進行提取，其主要步驟如下。

(1)根據(jù)礦物的譜形特征，分析其主要特征吸收譜帶范圍和特征吸收位置；

(2)計算不同礦物特征譜帶與其標(biāo)準(zhǔn)光譜相應(yīng)譜帶范圍的光譜匹配度；

(3)計算特征譜帶中最小吸收深度和位置等參量；

(4)對上述參量設(shè)定閾值獲得相應(yīng)礦物信息，形成礦物分布圖。

根據(jù)上述礦物識別方法，最終在研究區(qū)內(nèi)提取了褐鐵礦、方解石、白云石、(高、中、低)鋁絹云母、鎂綠泥石、鎂鐵綠泥石/綠簾石、鐵綠泥石、角閃石等11種礦物信息，如圖5所示。為進一步驗證礦物提取信息的準(zhǔn)確性，在研究區(qū)內(nèi)結(jié)合已有地質(zhì)資料，選擇具有一定分布面積，且易于對比的地物，進行光譜對比驗證。本文選取研究區(qū)內(nèi)方解石、白云石、中鋁絹云母、鐵綠泥石礦物信息，對其識別結(jié)果進行隨機選點驗證，對比其影像光譜與USGS標(biāo)準(zhǔn)波譜庫光譜。對比結(jié)果顯示(見圖6)，兩者在上述礦物的特征譜帶譜形一致性較好，特征吸收位置準(zhǔn)確，可見信息提取結(jié)果可靠。

圖5 礦物信息分布圖Fig.5 Distribution of mineral information

圖6 礦物光譜對比圖(黑色為USGS標(biāo)準(zhǔn)光譜，紅色為影像光譜)Fig.6 Comparison of mineral spectra (USGS standard spectrum in black and image spectrum in red)

5 結(jié)束語

本文利用資源一號02D高光譜數(shù)據(jù)，開展了巖性-構(gòu)造信息識別和礦物信息提取應(yīng)用評價。結(jié)果表明：資源一號02D高光譜數(shù)據(jù)經(jīng)信息增強后，可對不同地層、巖體以及構(gòu)造信息反映清晰，也可識別較小的構(gòu)造或脈體，可為今后1∶50 000地質(zhì)填圖提供重要技術(shù)支撐；在礦物填圖方面，通過典型礦物光譜對比分析，認(rèn)為其譜形特征明顯，吸收位置準(zhǔn)確，礦物信息提取準(zhǔn)確率高，可為以后找礦預(yù)測、資源潛力評價等提供重要參考信息。