Ⅰ 高光譜礦物信息提取
8.6.1 方法與流程
8.6.1.1 產品生成業務化流程
常用的高光譜礦物填圖標准方法：先運用最大雜訊分離（MNF）變換對反射率數據進行波譜降維，用像元純度指數（PPI）分析進行空間降維，然後用N維可視化方法來確定圖像端元，進行端元識別，確定礦物波譜，最後用合適的方法進行礦物填圖，並對結果進行優化操作。礦物填圖流程圖如圖8.40所示。
1）數據降維。利用MNF變換進行光譜數據減維，分離數據中的雜訊，減小運算量。觀察最終特徵值和MNF圖像，確定數據的固有維數，選擇合適的MNF波段進行PPI指數計算。
2）計算圖像的PPI指數。計算MNF圖像的PPI指數，最終產生PPI圖像，圖像中像元的DN值代表像元被記錄為極值的次數，從直方圖中選擇閾值，僅選擇最純的像元以保證被分析的像元數最小，這些像元被輸入到分離特定光譜端元的互動式可視化演算法中。
3）N維可視化進一步提純純凈像元。在N維可視化中，光譜可視為n維散點圖中的一個點，n是波段數。通過選擇n維散點圖的頂點和拐角，進一步提純純凈像元。將最終確定的端元輸入到圖像中的感興趣區（ROI），從圖像中提取每個感興趣區平均反射率光譜曲線作為成像光譜礦物填圖的候選端元。
圖8.40 常見高光譜礦物填圖流程
4）建立真實波譜庫。將採集端元的波譜曲線輸出，經過波譜分析（用波譜特徵擬合演算法與USGS標准波譜庫進行匹配識別）得到的波譜曲線圖，確定端元礦物並建立該礦物的真實波譜庫。
5）選擇合適的填圖方法，輔以自行建立的真實波譜庫，提取出各類礦物。並對提取結果做最小圖斑去除、類合並等操作，優化礦物提取結果。
8.6.1.2 礦物信息提取的高光譜數據預處理
先將成像高光譜數據進行輻射定標、大氣校正等過程得到光譜反射率數據。無論是高光譜成像儀還是傳統的多光譜感測器，它們所記錄的數據都是地面觀測目標的反射或輻射能量的光譜輻射絕對值，與地物目標的光譜反射率或光譜輻亮度值是不一致的。因此，輻射定標和光譜重建是地物識別和定量分析不可缺少的環節。只有經過輻射標定、輻射校正和大氣校正，剔除由於大氣散射、吸收、地形起伏及感測器本身不穩定帶來的各種失真，將記錄的圖像值轉換為地面的反照率值，重建像元地面光譜，才能根據光譜特徵，有效地識別地物，反演地物成分。
8.6.1.3 典型礦物光譜吸收特徵分析
首先打開波譜資料庫系統圖8.41調出所研究的典型礦物的光譜進行分析。
根據主要譜帶的位置將礦物劃分為含Fe2＋礦物、Fe3＋礦物、Mn2＋礦物、碳酸鹽礦物、含Al-OH鍵礦物、含Mg-OH鍵礦物等大類。
Al-OH：診斷譜帶一般位於2165～2215 nm附近；
Mg-OH：診斷譜帶一般位於2315～2335 nm附近；
Fe2＋：診斷譜帶一般位於1000～1100nm附近；
圖8.41 波譜庫系統主界面
Fe3＋：診斷譜帶一般位於600～900nm附近；
Mn2＋：診斷譜帶一般位於450～600nm附近。
以主吸收譜帶和次要吸收譜帶的組合特徵，並考慮到吸收譜帶在不同岩石中的變異，識別礦物族。如對Al-OH鍵礦物可分為明礬石族、蒙脫石族、白雲母族、高嶺石族；含Mg-OH礦物可分為綠泥石族、黑雲母族、綠簾石族、蛇紋石族、滑石-透閃石族等；碳酸鹽礦物包括方解石族、白雲石族等。
1）對於Al-OH鍵礦物：
明礬石族：兩個吸收特徵分別位於2175nm（2165～2175nm，± 10nm）和2325nm。
蒙脫石族：位於2215nm（2205～2215 nm，± 10nm）的主要吸收特徵伴隨有一個2440nm的次級吸收。
白雲母-伊利石族：位於2205nm（2195～2215 nm，± 20nm）的主要特徵吸收伴隨有三個次級吸收峰，分別位於2355nm，2440nm和2115nm。
高嶺石族：2200nm附近的OH-譜帶，在該譜帶左翼存在次一級的譜帶（2160nm）。
圖8.42為典型Al-OH鍵礦物波譜特徵圖。
圖8.42 幾種典型Al-OH鍵礦物波譜特徵圖
2）對於碳酸鹽礦物：
方解石族：位於2330 nm～2340 nm的主要譜帶。
菱鐵礦族：同時具有Fe2＋在1000nm與
圖8.43為典型碳酸鹽礦物波譜特徵圖。
圖8.43 幾種典型碳酸鹽礦物波譜特徵圖
3）對於Mg-OH鍵礦物：
綠泥石族：位於2330 nm主要吸收特徵伴隨有一個2386 nm及2224 nm兩個次級吸收。
綠簾石族：位於2330 nm主要吸收特徵伴隨有一個2225～2226 nm的次級吸收。
蛇紋石族：位於2320 nm附近主要吸收特徵伴隨有2110 nm為中心的較寬的弱譜帶。
滑石-透閃石族：位於2290～2310 nm附近的較寬的雙峰譜帶和2385 nm附近的次級吸收。
圖8.44為典型Mg-OH鍵礦物波譜特徵圖。
圖8.44 幾種典型Mg-OH鍵礦物波譜特徵圖
8.6.1.4 礦物識別技術
目前，基於成像光譜數據的礦物識別的方法主要有：基於光譜波形參數、基於光譜相似性測度、基於混合光譜模型、基於地質統計規律、基於光譜知識的智能識別等。
（1）光譜波形參數提取與蝕變識別的技術方法
岩石礦物單個診斷性吸收特徵峰可以用吸收波段位置（λ）、吸收深度（H）、吸收寬度（W）、吸收面積（A）、吸收對稱性（d）、吸收的數目（n）和排序參數作完整地表徵（陳述彭等，1998）。根據端元礦物的單個診斷性吸收波形，從成像光譜數據中提取並增強這些參數信息，可直接用於識別岩礦類型。如HSI編碼與吸收波段圖（Kruse，1988）是利用連續統去除後的光譜圖像，定義出波段吸收中心位置圖像，波段深度圖像及波段半極值寬度（FWHM）圖像，並分別賦予HSI空間的明度（H）、強度（I）和飽和度（S），然後逆變換到RGB色度空間。代表性的方法還有光譜特徵擬合（SFF）等。
光譜特徵擬合（SFF，Spectral feature fitting）：光譜特徵擬合是一種基於光譜吸收特徵的方法，使用最小二乘法技術將圖像中每一個像元的光譜曲線與所選擇的參考光譜曲線的吸收特徵進行擬合。是選擇包含目標礦物特定吸收譜帶的光譜區間，利用最小二乘擬合方法，比較像元光譜與目標光譜吸收特徵的整體形態和吸收深度。該方法對噪音和地形有較好地抑製作用。
（2）基於相似性測度的識別技術方法
成像光譜最大的優勢在於利用有限細分的光譜波段，去再現像元對應物的光譜曲線。這樣，利用整個光譜曲線進行礦物匹配識別，可以在一定程度上改善單個波形的不確定性影響（如光譜漂移、變異等），提高識別的精度（甘甫平等，2000）。基於整個波形的識別技術方法是在參考光譜與像元光譜組成的二維空間中，合理地選擇測度函數度量標准光譜或實測光譜與圖像光譜的相似程度。這類方法主要有光譜角制圖（SAM）技術、光譜二值編碼（Binary encoding）技術等。
（3）基於光譜知識模型識別的技術方法
基於光譜知識模型識別的技術方法是建立在一定的光學、光譜學、岩礦結晶學和數學理論之上的信號處理技術方法。它不僅能夠克服利用單一譜形識別所存在的缺陷，而且從地物光譜學原理入手、從本質上理解認識岩礦光譜的物理機制與物理過程，建立光譜數學物理模型，識別並定量提取岩礦信息。這在一定程度上能精確地量化地表物質的組成及深入地描述地物組成的物理特性，以進一步探測地物所蘊藏的成生環境本質。例如，建立在Hpake（1981）光譜雙向反射理論基礎之上的線性混合光譜分解模型（SMA/SUM），可以根據不同地物或者不同像元光譜反射率響應的差異，構造光譜線性分解模型，從而識別地物，量化地物成分，挖掘地物成生環境信息。
（4）基於地質統計特徵的分類識別方法
該類方法是基於地物在圖像上的統計分布規律，建立地質模型進行圖像分類識別。比較典型的有（Ieohku，etal，1996），概率模型（Porbabilistic Model），幾何光學模型（Geomertic-optical Model），隨機幾何模型（Stochastic Geometric Model）及非參數地質統計模型（Non-paramertic GeostatisticTechnique）。
（5）基於光譜知識的智能識別方法
傳統的及上述的成像光譜識別方法利弊共存。對於高維與超大容量的成像光譜數據及大量的實驗室光譜研究結果等迫切要求新的高效的遙感定量分析技術。因此，專家系統、人工神經網路、模糊識別等基於光譜知識的智能識別應運而生。神經網路技術是應用最廣泛的識別技術（Gong，1996；Jimenez，et al.，1998；Benediktsson，et al.，1995；Garcia Haro，et al.，1998）。Goetting 和 Lyon（1986）建立一個試驗性專家系統；Kruse等（1993）建立了基於知識的成像光譜礦物自動匹配製圖系統。以Dempster-Shafer證據理論為基礎格架的證據推理方法也得到了一定程度的應用。基於光譜知識的智能識別技術方法與系統將是成像光譜遙感地物信息識別、提取與量化及實用化的最佳選擇與發展方向，極具潛力。
8.6.2 結果與分析
8.6.2.1 蛇紋石礦物的填圖及效果分析
對Hyperion反射率數據進行MNF變換，根據變換結果人工確定變換維數，利用像元純度指數（PPI）進行分析，根據n維光譜空間進行特徵端元採集。運用光譜特徵擬合的方法對端元波譜曲線進行波譜匹配，建立真實的波譜庫，如圖8.45所示。
圖8.45 蛇紋石波譜匹配圖
圖8.46 五種填圖方法結果對比
用前文簡述的常用五種方法提取蛇紋石的具體位置，並進行對比分析，填圖結果如圖8.46所示。
用SAM礦物填圖方法，可以快速提取端元採集出的所有地物，並用彩色分類影像來顯示，使人們很直觀地看出該礦物的分布范圍。
光譜特徵擬合將為每個參照波譜輸出一幅比例圖像和 RMS 圖像或一幅合成的「擬合」圖像（Scale/RMS）。較高的擬合值表明該像元與參照波譜匹配較好。但是波譜特徵擬合效果並不是很好，因為如果輸入了錯誤的參照端元或使用了錯誤的波長范圍，也會出現一個遠遠大於1的比例值。
MTMF（混合調制匹配濾波）結果將以一系列灰階圖像的形式出現，兩幅圖像對應一個被選波譜。浮點型結果提供了像元與參照波譜相對匹配程度的估計方法（1.0表示完全匹配），以及亞像元的權重。不可行性（Infeasibility）值以sigma雜訊為單位，顯示了匹配濾波結果的可行性。用於准確制圖的像元有一個大於背景分布值的較高的匹配濾波值和一個較低的不可行性值。對得出的浮點圖進行密度分割等操作，也能很好的得出蛇紋石的分布范圍。
MF（匹配濾波結果）將以一系列灰階圖像的形式出現，每幅圖像對應一個選擇的端元。浮點型結果提供了像元與參照波譜相對匹配程度的估計方法（1.0表示完全匹配），以及亞像元的權重。對MF的結果圖做密度分割等操作，也可明顯地看出蛇紋石的分布范圍。
線性光譜分解的結果將以一系列灰階圖像的形式出現，每幅灰階圖像加上一個RMS誤差圖像對應於一個端元。經過發現，此方法的提取效果並不是很好，雜訊較多。因此在實際應用中，並不建議使用此方法。
通過以上對比分析，發現SAM，MTMF，MF三種方法在本地區都比較適合蛇紋石這種礦物的填圖，只是MTMF，MF兩種方法還需對結果圖進行密度分割處理，稍顯繁復，但填圖效果良好。
8.6.2.2 白雲母礦物的填圖及效果分析
對Hyperion反射率數據進行MNF變換，根據變換結果人工確定變換維數，利用像元純度指數（PPI）進行分析，根據n維光譜空間進行特徵端元採集。運用光譜特徵擬合的方法對端元波譜曲線進行波譜匹配，建立真實的波譜庫（圖8.47）。因為隨著白雲母中Al含量的減少，Al-OH的譜帶位置逐漸向長波方向移動，所以可以分辨出貧鋁白雲母和富鋁白雲母的波譜曲線。
圖8.47 白雲母波譜匹配圖
用上述五種方法提取蛇紋石的具體位置，填圖結果如圖8.48所示。
經過以上五種填圖方法的對比分析發現，SAM填圖方法效果最佳，不僅明確填出白雲母礦物的范圍，而且避免了誤差引起的填圖錯誤。此外，對光譜特徵擬合後的結果圖做密度分割等操作，所提取的礦物范圍也較好。因此，兩種方法均適合進行白雲母的礦物填圖。
8.6.2.3 綠簾石礦物的填圖及效果分析
對Hyperion反射率數據進行MNF變換，根據變換結果人工確定變換維數，利用像元純度指數（PPI）進行分析，根據n維光譜空間進行特徵端元採集。運用光譜特徵擬合的方法對端元波譜曲線進行波譜匹配，建立真實的波譜庫，如圖8.49所示。
用上述五種方法提取蛇紋石的具體位置，填圖結果如圖8.50所示。
經過以上五種填圖方法的對比分析，發現SAM填圖方法效果最好。但光譜特徵擬合、MTMF和MF三種方法對綠簾石也有較好的顯示，對他們的結果圖做密度分割等操作，也能較好的提取出蛇紋石的范圍，因此，以上四種方法都適合蛇紋石礦物的提取填圖。
圖8.48 五種填圖方法結果對比
圖8.49 綠簾石波譜匹配圖
8.6.2.4 綠泥石礦物的填圖及效果分析
選取HyMap數據，選取USGS標准波譜庫中綠泥石的光譜曲線作為參考光譜。利用上面的方法提取綠泥石的具體位置，填圖結果如圖8.51所示。
8.6.2.5 多種礦物填圖情況
（1）HyMap數據1情況分析
依據地面區域地質資料，在HyMap反射率數據上選取典型的訓練樣區，利用這些訓練樣區提取五種礦物的分布位置。
SAM方法填圖：該方法需要設置閾值max angle threshod，通過不同的嘗試，綠泥石與富鋁白雲母分布較廣，針對兩種閾值的設定情況對比如下（圖8.52 ，圖5.3）。
Binary Encoding方法填圖：該方法需要設置一個閾值minimum encoding threshod，通過不同的嘗試，綠泥石與富鋁白雲母分布較廣，針對兩種閾值的設定情況對比（圖8.54 ，圖8.55）。實驗發現富鋁白雲母（黃顏色圖例）在閾值設置為0.97和0.96時，分類效果變化非常大，由於軟體對該閾值的設定只支持小數點後2位，所以不能進一步的精確，對比分析兩種閾值的分類效果，得出二進制編碼方法對富鋁白雲母的提取效果不是很好。
Minimun Distance（MD）方法填圖：該方法需要設置兩個閾值：Max stdev from Mean；Max Distance Error。通過不同的嘗試，綠泥石與富鋁白雲母分布較廣，針對兩種閾值的設定情況對比如圖8.56和圖8.57所示。
Spectral Information Divergence（SID）方法填圖：該方法需要設置閾值：Maximun Divergence Threshod。通過不同的嘗試，參數設置和填圖效果如下（圖8.58至圖8.60）。
（2）HyMap數據2情況分析
依據地面區域地質資料，在HyMap反射率數據上選取典型的訓練樣區，利用這些訓練樣區提取五種礦物的分布位置。
SAM方法填圖：該方法需要設置閾值max angle threshod，通過不同的嘗試，綠泥石與富鋁白雲母分布較廣，針對兩種閾值的設定情況對比如圖8.61所示。
最小距離填圖：最小距離填圖的參數設置及填圖效果差異如圖8.62和圖8.63所示。
（3）模擬數據情況分析
依據地面區域地質資料，在模擬的反射率數據上選取典型的訓練樣區，利用這些訓練樣區提取五種礦物的分布位置。
圖8.50 五種填圖方法結果對比
圖8.51 三種填圖方法結果對比
圖8.52 局部區域填圖效果
圖8.53 整體區域填圖效果
圖8.54 填圖效果1
圖8.55 填圖效果2
圖8.56 局部區域填圖效果截圖
圖8.57 整體區域填圖效果
圖8.58 不同閾值局部地區填圖效果
圖8.59 填圖效果（Threshold＝0.001）
圖8.60 填圖效果（Threshold＝0.005）
圖8.61 SAM 方法填圖
圖8.62 參數1下效果圖
圖8.63 參數2下效果圖
SAM方法填圖：該方法需要設置閾值max angle threshod，通過不同的嘗試，閾值的設定情況和填圖效果如圖8.64所示。
圖8.64 SAM 填圖效果
Binary Encoding方法填圖：該方法需要設置一個閾值minimum encoding threshod，通過不同的嘗試，閾值的設定和填圖效果情況如圖8.65所示。多次試驗發現該方法對提取chlorite和chlorite＋serpentine的效果不是很理想。
圖8.65 填圖效果
Minimun Distance（MD）方法填圖：該方法需要設置兩個閾值：Max stdev from Mean；Max Distance Error。通過不同的嘗試，閾值的設定和填圖效果情況如圖（其中Max Distance Error＝2000 ，圖8.66）。
圖8.66 填圖效果
Spectral Information Divergence（SID）方法填圖：該方法需要設置閾值：Maximun Divergence Threshod。通過不同的嘗試，參數設置和填圖效果如圖8.67所示。
圖8.67 填圖效果（Threshold＝0.005）
8.6.2.6 礦物信息填圖結果分析
綜合前面的分析結果，開展了礦物信息提取結果圖，為了分析礦物信息填圖效果，搜集了前人對此試驗區（圖8.68）進行的礦物填圖結果（圖8.69）。
圖8.68 礦物信息研究的區域圖
圖8.69 前人的礦物提取結果填圖結果
為開展蝕變礦物識別精度分析，可對比研究區本次提取的蝕變礦物分布圖（圖8.70）與前人提取的蝕變礦物信息分布圖結果，逐像素進行對比（由於解析度不同需做像素變換），如果研究區總的點數為N，本次與前人礦物分布圖的結論相同就認為該點取值1，最後統計結果中1的個數n，這樣蝕變礦物識別精度可定義為
高光譜遙感技術原理及礦產與能源勘查應用
對比分析本次礦物信息填圖結果和前人礦物信息的填圖結果，發現其結果並非完全一樣。其中蛇紋石的填圖效果一致性最高；綠簾石的分布范圍一致，但本次數據的填圖結果顯示的范圍較大；白雲母的分布范圍基本一致，但貧鋁白雲母和富鋁白雲母的分布范圍有些混淆，經過對比分析，本次數據和前人數據填圖范圍的一致性達到89%，基本滿足高光譜數據的填圖要求。對出現誤差的主要原因分析如下：
1）前人數據是機載高光譜數據，空間解析度可達到3～12m，此處的Hymap數據的空間解析度為5m，而本次數據是模擬星載高光譜數據，空間解析度為30m。隨著解析度的降低，單個像元所對應的地面面積將增大，導致每個像元中包含更多的礦物類型，礦物間的影響性也會增加。一方面，由於混合像元的平均效應，目標礦物在像元中的等效豐度會下降，光譜信息減弱。當像元中目標礦物的等效豐度下降到檢出限以下時，礦物將不能被識別，而造成礦物分布區的外圍含量較低的地段和含量較低分布區填繪面積的減少，點狀集群分布區的漏識別，線狀分布區的斷續分布。這種效應相當於檢出限的下降。另一方面，當像元中目標礦物的等效豐度在檢出限以上時，會使面狀礦物分布區的范圍擴大、空洞的充填、相鄰小區的連接，點狀集群分布區的成片，線狀區域的斑點效應。這兩種效應的綜合作用結果，空間解析度的減小則會使強異常區（包括高豐度區和光譜反襯度較高的礦物分布區）范圍的擴大、小區的相連、點群的結合，而使異常更加醒目，但其細節特徵會因此消失；弱異常區（包括低豐度區和光譜反襯度較低的那些礦物分布區）和小異常區面積會縮小或漏檢；線狀異常可能會形成斷續分布的小斑塊而使線狀特徵和其走向變得不清晰。因此，在填圖的細致程度上，解析度越高的數據填圖效果越好，這是造成填圖效果差異的最主要的原因之一。
圖8.70 本次試驗礦物信息提取填圖結果
2）輻射校正的精確程度。無論是高光譜成像儀還是傳統的多光譜感測器，它們所記錄的數據都是地面觀測目標的反射或輻射能量的光譜輻射絕對值，與地物目標的光譜反射率或光譜輻亮度值是不一致的。因此，輻射校正和光譜重建是地物識別不可缺少的環節。由於Hymap數據是機載高光譜數據，Hyperion數據是星載高光譜數據，對兩種數據進行輻射校正的參數也是不一致的，因此，得出的反射率數據也是有差異的，這也是造成填圖差異的原因之一。
總之，用本次數據進行礦物填圖是可行的，它可以在一定程度上對礦物的種類及分布進行識別。
Ⅱ 成像光譜方法技術
一方面，高光譜解析度的成像光譜遙感技術是對多光譜遙感技術的繼承、發展和創新，因此，絕大部分多光譜遙感數據處理分析方法，仍然可用於高光譜數據；另一方面，成像光譜技術具有與多光譜技術不一樣的技術特點，即高光譜解析度、超多波段（波段＜1000，通常為100～200個左右）和甚高光譜（Ultra Spectral）解析度（波段＞1000，主要用於探測大氣化學成分）的海量數據。因此，常規多光譜數據處理方法不適合於成像光譜數據的定量分析，於是成像光譜數據處理和分析技術應運而生。在成像光譜數據處理和分析方法中，關鍵性的技術問題是地物光譜重建，光譜特徵的量化及提取，混合像元的分解和定量分析及模型識別。
2.3.2.1 光譜重建技術
按照不同的模型及演算法，從成像光譜數據中把地物的光譜特性反演出來的過程就是地物光譜重建技術。根據不同的工作情況及條件，採取不同反演模型來重建地物光譜，是實現成像光譜數據遙感定量化分析的第一步。若對其不進行反演，則沒有一個統一物理量進行對比。目前，光譜反演模型大體可分為三大類型：基於大氣傳輸理論的大氣傳輸模型，基於統計分析的統計模型以及基於地面地物同步觀測的經驗回歸模型。
2.3.2.1.1 基於大氣傳輸理論的模型
該模型實質上就是用理論模型消除大氣中氣體分子、水蒸氣、氣溶膠及塵埃等分子顆粒對地反射輻射能量吸收與散射以及大氣程輻射效應，並將其還原成地物的反射輻射光譜。這是一種比較復雜的同時必須進行地物光譜及大氣參量測量的絕對反射率生成方法，也就是對成像光譜數據進行絕對輻射標定的再反演。在這一反演過程中，關鍵是建立大氣傳輸的模型。自1960年，Chndrasekhar提出了輻射傳輸理論以來，相繼發展了許多方法，如Ordinate方法和Variational方法等來解決輻射傳輸問題。目前，常見的大氣傳輸模型有5 S（the Simulation of the Satellite Signal in the Solar Spectrum）、6S（the Second Simulation of the Satellite Signal in the Solar Spectrum）、LOWTRAN 7及MODTRAN（Teillet，P.M.，1989；Vermote，E.，Tanṙen，D.，Deuże，J.L.et al.，1994，1996；Bo⁃Cai Gao，K.B.Heidebrecht and A.F.H.Goetz，1997；Z.QIN，A.Karnieli and P.Berliner，2001）。其中，6S模型是由法國Tanré等人研究開發的，是目前世界上發展比較完善的大氣輻射校正模型演算法之一。該演算法既能合理地處理大氣散射、吸收，又能產生連續光譜，避免在光譜反演中較大的定量誤差。它還充分利用了分析表達式和預選大氣模式，使計算時間大大縮短。許多遙感專家使用此模式進行地物光譜反演後認為，該模型較其他模型計算精度高。不足之處是必須開展試驗區典型地物光譜反射率觀測以及大氣環境參量實測，如：大氣光學厚度、溫度、氣壓、水蒸氣含量、大氣分布狀況等。相對來說，盡管LOWT⁃RAN 7和MODTRAN模型計算精度低一些，但它不需要地面實測典型地物的反射率。這些模型一般用於對感測器選定標定場，開展數據絕對輻射標定。
2.3.2.1.2 基於統計分析的模型
該模型的建立是在分析不同地物光譜遙感信息在不同光譜波段的傳輸特點基礎上，利用計算機對典型地物的光譜特性進行統計分析後，得到的地物光譜特性反演模型。對成像光譜數據進行地物光譜反演常用模型有平滑域反射率模型 FFR（Flat Field Reflectance）（Goetz and Srivastava，1985；Conel，1985；Crowley at al.，1988；Rast et al.，1991），內在平均相對反射率模型IARR（Internal Average Relative Reflectance）（Kruse et al.，1985；Kruse，1988；Mackin and Munday，1988；Zamoudio and Atkinson 1990），對數剩餘模型LRC（Logarithm Resial Correction）（Green，R.O.etal.，1985；Gower，J.F.R.etal.，1992）。在這3種模型中，FFR模型是在圖像上選取光譜和地貌特徵都均勻平滑（平滑性是指地物無光譜吸收譜帶，光譜曲線平直）的地物平均值，來消除大氣輻射衰減和儀器的零響應；White模型是根據整幅圖像的平均光譜曲線平均值對圖像歸一化處理，然後計算每個像元光譜曲線與平均光譜曲線的比值，也就是地物光譜特性；LRC模型是經Lyon和Lanze修正後，對太陽輻射衰減、大氣效應及地形影響都有所消除。Green和Graige提出的對數剩餘糾正公式如下：
lg（Rij）=lg（DNij）-lg（aveDNi）+lg（DNi）+lg（DNg）（2-3-1）
這里Rij是第i波段、第j個像元的剩餘值；DNij是第i波段、第j個像元的亮度值；DNi是第i波段所有像元數據的平均值；DNj是第j像元在所有波段上數據的均值；DNg是所有波段及所有像元的均值。該方法完全基於圖像本身特徵，不需要野外地物光譜測量。在前兩種模型中，FFR模型優於IARR模型，它克服了IARR模型因受像福強吸收特徵的影響而出現的假反射峰的弱點，而且計算量較小。
2.3.2.1.3 經驗線性回歸模型
利用該方法重建地物光譜技術實質就是通過開展典型地物的同步反射率觀測，根據成像光譜數據DN值與地面實測地物反射率值，經最小二乘法求出回歸方程Rij=Aj·DNij+Bj（這里Aj，Bj是感測器第j波段的線性回歸系數），然後，根據此方程反演地物的反射光譜。這種模型的數學和物理意義明確，方法簡便，運算量少，應用廣泛（Roberts et al.，1985；Conel et al.，1988；Elvidge，1988；Green et al.，1988；Kruse et al.，1990；Zamoudio and Atkinson 1990）。例如：美國JPL的Abrams利用該模型在美國Nevada州的Cuprite礦區進行礦物學填圖；美國科羅拉多大學的Zamudio等人，利用該模型在美國 Nevada 州東部進行礦物識別和岩相分析；日本Pasco公司的Mochizuki，利用該模型在美國Navada州進行蝕變礦物的反射光譜研究等。該模型的不足之處是要開展野外地物光譜觀測，成本比較高，回歸精度的高低依賴於對野外概實測的精度。
除上述這些典型光譜重建模型之外，還有 UA RT Code，JPL Code，連續內插波段比演算法（CIBR），背景法等模型（De Jong，1998）。
2.3.2.2 岩礦光譜特徵的量化、提取，定量分析及識別模型
成像光譜數據經過光譜重建模型處理後，獲得了地物的光譜特徵譜線。不同地物光譜具有不同的診斷特徵譜帶，如吸收譜帶，特徵譜線的微積分變化，波形變化等等。如何有效地開展地物特徵定量分析和識別地物，首先要弄清楚如何去量化及提取地物的光譜特徵。因此，開展基於地物特徵譜的量化提取是十分必要的。
2.3.2.2.1 地物光譜特徵度量、提取與匹配識別模型
（1）就地物光譜特徵（這里指地物反射輻射光譜）而言，不外乎兩大類型：吸收譜帶（或反射谷）和光譜曲線的斜率變化（含波形變化）。針對這兩類光譜特徵的形態、結構，分別採取不同的度量方式。目前，對吸收譜帶的分析度量方法是外殼系數法，它通過把光譜曲線歸一化後去測量吸收譜帶的波長位置（position）、吸收深度（depth）、吸收寬度（FWHM）和對稱性（sym⁃metry）（Lyon et al.，1985；F.A.Kruse，A.B.Lefkoff，1993）。這種外殼系數法可以由外殼凸形曲線與光譜之比來表示，也可以由外殼值去減相應波長上的光譜反射率值來求得。由於吸收峰的非對稱性，採用RBD方法難以准確描述其特徵。連續插值小組段演算法（CIBR，Continuun interpolated band algorithm）（De Jong，1998）和光譜吸收指數圖像（SAI，spectral absorption index image）（王晉年等，1996）與相對吸收深度圖方法類似，但引入了對稱度因子，使其對吸收特徵的描述更為合理。除了這些測量參數外，對植被光譜有多種度量參量，如植被、綠度指數等等。對於光譜曲線斜率變化的特徵，表徵和提取的方法有基於地物光譜的總體波形特徵度量，如傅里葉變換的波形分析方法是利用有限級次的諧波振幅和初位相度量地物波譜特徵；基於切比雪夫多項函數的波形分析是利用多項式函數對地物光譜曲線進行擬合，提取有限項的系數來表示或組合特徵，或用其比值來表示地物光譜波形特徵的參量；基於波形相似性（總體或分段）分析的光譜角度量；還有光譜曲線特徵的微分度量、積分度量及二值度量等等。當然，對吸收譜帶及斜率特徵度量還可以用統計特徵量去表徵度量，如均值、方差、協方差矩陣、特徵值、特徵向量、特徵因子及組內離差等。
（2）光譜匹配識別模型不同於多光譜的模式識別，它是根據光譜特徵度量參數來進行匹配識別的，是成像光譜數據處理分析的特色之一。這種特色模型在處理過程中往往是採用可視化的互動式的圖像與光譜、光譜與標准光譜形式進行的。目前，光譜匹配識別模型有：編碼匹配識別法（均值編碼匹配法，坡向編碼匹配法，比值編碼匹配法L吸收峰編碼匹配法、波形匹配法和光譜角匹配法等等。在這些匹配識別模型中，編碼匹配基本上都是按二值（0和1）進行編碼、匹配與識別的；吸收特徵編碼匹配是根據外殼系數法歸一化後，對每個特徵的吸收深度與波長位置進行編碼；波形匹配包括相似度、傅里葉變換參量，切比雪夫等方法。近幾年發展起來的小波變換分析在成像光譜數據分析處理中應用相當多，尤其是對原始信號按不同小波尺度，分解成不同的小波進行波形分析，突出低頻弱信息，有利於信息增強，比如用小波變換進行圖形圖像插值、融合及混合像元分解等。
2.3.2.2.2 成像光譜數據的定量分析及識別模型
定量化分析及識別模型化是當今遙感技術的發展方向之一，應用於成像光譜數據處理。定量化分析與識別模型，除了不斷完善和改進已有基於統計分析的定量化及識別模型（如：改進的主成分分析、最佳波段組合、改模型最大似然法、基於決策邊界特徵矩陣的變換和正交子空間投影），其他學科的新思想、新方法也在不斷地引人遙感數據分析和理解之中，如人工智慧的專家系統，模糊邏輯映射，證據推理、神經網路、分形和分維等。
人工智慧專家系統技術是目前比較流行的信息處理技術，尤其對比較復雜問題的解決有獨到之處，Gotting和Lyon在1986年就已建立的光譜信息專家識別系統用於分析實驗室和野外光譜，它是結合已有的地物光譜特徵知識，由專家確定判別規則決策樹（Decision Tree）以達到識別地物或地物類別目的。決策樹這一基於知識的判別准則層次是建立專家系統成敗的關鍵。通過這個系統進行編碼匹配，他們成功地從大量的實驗室光譜中識別出11種礦物。1993年，美國地球空間研究中心（CSES）和美國環境科學研究所的F.A.Kruse和A.B.Lefkoff研製了基於知識的成像光譜地質制圖專家系統。選定各種特徵在識別過程中的作用並賦予相應的權值，或根據專家對判別知識和經驗建立判別准則進行識別。
目前，神經網路模型在遙感地物分析和識別方面備受青睞，有著廣泛的應用（Golen Giser，1996；Giles，M.F.et al.，1995；郭小方，1998；王潤生等，2000）。由於神經網路分類規則對訓練樣本的數量及分布特徵沒有特定要求，因而可以在特徵空間形成非線性判別邊界，並且還有一定的抗雜訊、抗干擾和自適應能力，適用於大數據量的分類研究，最為常用的分類准則是後向傳播（BP）網路模型。
目前，從成像光譜遙感數據分析與識別的各種新理論、新方法的引入來看，大多數模型的研究和應用還是一種嘗試，在如何將模型與成像光譜數據相結合的研究方面，分析不夠深入。
2.3.2.3 混合像元分解模型
由於空間解析度不高的原因，在圖像像元內會出現不同成分（end member）的地物，即混合像元。不同的地物具有不同光譜特徵，因此需要通過混合光譜分解技術來提高識別精度。混合像元問題是遙感技術的研究難點和熱點。由於成像光譜技術的光譜解析度已從微米（μm）提高到納米級（nm），因此，其混合像元分析、分解及其模型研究就顯得更為重要。
目前，開展高光譜遙感混合像元研究的方法技術，首先從實驗著手，進行地物混合光譜的測試、分析、數字模擬、分解模型開發研究，然後將其外推到遙感圖像上，進行典型地物混合像元分析，主要包括空-地同步觀測獲取典型地物（或可通過人工布標）數據，經模型分析後，對混合像元的地物進行分解，或混合光譜模擬合成。在實驗室里，通過對不同礦物光譜混合含量測試發現，不透明礦物或暗色礦物，其光譜按比例混合到其他礦物中，混合的反射率急劇下降，而不是逐漸下降，說明其混合光譜與其混合的端員礦物光譜是非線性關系（磁鐵礦和橄欖石）。當兩種礦物的色調相近時，實驗測試的混合光譜與線性模型合成的混合光譜都呈線性逐漸變化，說明混合光譜可以按線性模型分解端員礦物光譜，如橄欖石和紫蘇輝石，且吸收譜帶的波長位置也是逐漸從一個波長位置逐漸過渡到另一波長位置。不僅如此，還發現在可見光、近紅外這一波長上，低成分端員混合時呈線性趨勢，當成分增加時，線性關系劇烈變成非線性關系。在這三種情況中，第一種非線性關系是由於組合混合光譜的端員成分之間互相作用、互相影響後光譜被光譜儀檢測到；第二種線性關系是由於各端員成分之間無互相影響作用，各自獨立地反射電磁波能量貢獻於混合光譜；第三種情況是兩種關系都存在，二者之間存在臨界條件（邊界條件）。目前，有關此方面的研究極少。根據這些分析，混合像元分解模型大體分為線性模型和非線性模型。在遙感混合像元中，絕大多數反射率相似的地物，可以用線性模型來分解端員成分，如：土壤與植被、不同含水量的耕地、岩石露頭與草地、荒地等等。在一幅圖像中，事先知道有N種端員（地物種類），並且也知道各種端員的光譜反射率，那麼就可以用線性模型：
成像光譜岩礦識別方法技術研究和影響因素分析
這里DNc是波段C上混合像元的DN值或反射率；Fi是第i種端員在混合像元中所佔比例（或權系數）；DNi，c是C波段上第i種端員的DN值（或反射率）；Ec是C波段上擬合誤差。對每個像元都按照最小二乘法解方程，進行分解。在圖像中，端員的DN值（或反射率值）要麼可以從訓練區取值，要麼在地面實測。端員成分的確定過程實質上是一個迭代過程，迭代結果使M個波段上總誤差ε最小（且N≤M）。
成像光譜岩礦識別方法技術研究和影響因素分析
求得版中各種端員成分之後，就可以定量或半定量地對端員豐度製作豐度等專題圖件。
用非線性模型開展混合像元分解不多見，但已有這方面的研究，如模糊分割模型（Jin Ⅱkim，1996），概率鵬模型，幾何光學模型（Charles Ichoku，1996）及基於神經網路模型的混合像元分解（王喜鵬，張養貞等，1998）等等。
目前開發的模型有：
——光譜吸收指數模型SAI（王晉年，童慶禧等，1996）：
SAI=∑fiSAIi，∑fi=1，fi＞0 （2-3-4）
——高斯模型法MGM：該模型是基於礦物和岩石的反射、吸收光譜性質模擬反射光譜的各種模擬方法。它是一種確定性的而不是統計性的方法。高斯改進模型MGM 是近幾年在分析反射光譜的基礎上發展起來的分析技術（Cloutis，1989，Veverka，J.et al.，2000）。
m（x）=S*exp（-（xn-μn）2/2σ2）， （2-3-5）
通常取n=-1。
光譜識別與分類技術（Spectral Classification）：主要是利用地物高光譜特徵的量化參數，結合其在圖像空間上分布進行提取有利的信息，達到分類的目的。主要的分類方法有：
——最大似然法MLC：
g（x）=-ln｜∑i｜+（x-mi）t∑i-1（x-mi），（2-3-6）
——人工神經網路技術ANN：一般採用前饋網路模型，即第一隱層的節點輸入等於輸入層諸節點輸出的加權和。迭代的次數以系統的平均誤差為最小時為准。
成像光譜岩礦識別方法技術研究和影響因素分析
——光譜角制圖法SAM（Spectral Angle Mapper）：該方法是通過計算測試樣本光譜矢量（像元光譜）與參考光譜矢量（訓練的端員樣品光譜，或標准光譜庫的光譜），在n維空間（n波段）上的角度來確定它們兩者的相似度。一般兩矢量之間的角度越小，兩光譜向量越相似，進而可識別兩種地物為同類，否則視為異類。數學模型是：
成像光譜岩礦識別方法技術研究和影響因素分析
這里i=1，2，3，……，n，n為波段數。
——光譜維特徵提取法（Spectral Dimension Feature Extraction）：在高光譜遙感分類中，使用該方法對多波段、高相關、數據冗餘度高的數據進行降維處理。相關的有統計方法，如主成分、典型變數及改進的PCA法等。
——光學模型（Optical Modeling）：除了前述的數據分析及模型外，植被因其特有反射性質，還有獨特分析模型（光學模型）。該模型主要利用高光譜遙感數據預測或估計植被的多種生物物理、化學參量，如葉面指數LAI、總生物量、覆蓋度等；葉綠素、水分、N、P、K含量等。該模型也屬於經驗性的統計模型方法。一般性通用模型為：
S=f（λ；θs，Φs；θv，Φv；С）， （2-3-9）
這里S為預測的生物物理、化學參數；λ是波長；θs，Φs，θv，Φv是入射光和感測器探測幾何位置參數，C是描述植被冠層為特性參數。依靠法的模型有葉子光學性質光譜模型PROSPECT，葉子的任意斜散射模型SAIL，即生化參量反演的LIBERTY模型等。
高光譜在植被應用中除了生物、化學參量的反演分析外，還注重利用植被光譜特性譜線的藍邊、反射峰、黃邊、紅光吸收谷、紅邊、近紅外反射高原區等變化及數據的歸一化、對數、微分等變換，來監測植被的長勢及病蟲害，進行森林識別、分類、制圖（Clark，R.N.，Roush.T.L.，1984）。
2.3.2.4 光譜數據應用處理分析軟體
通過開展岩礦高光譜特性測試分析和成像光譜方法技術及應用分析研究，已發展並開發了如下數據處與分析軟體：
2.3.2.4.1 光譜資料庫及分析軟體（400～2500φ）
國外：美國地質調查所USGS和JPL的標准礦物光譜庫（含機載光譜）及光譜分析管理軟體SPAM，IRIS，日本地調所的岩石礦物光譜庫等（http://speclib.jpl.nasa.gov；http://speclab.cr.usgs.gov；Kruse F A et al.1993）。
國內：中國科學院安徽光學研究所、中國科學院遙感應用研究所、原地質礦產部航空物探遙感中心等科研單位都已建自己的光譜庫（王潤生等，2000）。
2.3.2.4.2 圖像處理分析軟體
目前國內外常用的光譜圖像處理分析軟體有：Erdas、PCI、ENVI等。其中PCI和ENVI都有高光譜分析處理功能（ENVI User』s Guide.，2000）。此外，還有像Tetricorder（Clark，R.N.，G.A.Swayze，K.E.Livo，2003）。國內通過高光譜遙感方法技術及示範應用研究，中國科學院遙感應用研究所、國土資源部航空物探遙感中心相繼建立了成像光譜數據分析處理系統，如：HIPAS，ISDPS等。
Ⅲ 遙感光譜數據的獲取
遙感技術從航空攝影測量逐步演變發展起來，大致經歷了3個發展階段：
1.航空攝影測量發展階段
目前仍保存著的最早一幀航空相片是1860年J.W.布萊克從氣球上拍攝的波士頓市的相片。在地質上的應用則始於1913年，有人在飛機上用攝影機對著非洲利比亞的本格遜油田攝影成像，並用這套骯空相片編制了本格遜油田地質圖。航空攝影遙感主要以飛機或者氣球為運載工具，用航空攝影機對目標獲取信息，然後再經過負片和正片過程得到最終的航空相片。航空攝影利用的是電磁波可見光全色波段，用感光膠片接受所攝目標物反射來的太陽光線感光、成像，一般感光片的感光范圍是0.3～0.9μm。航空攝影大多數情況下是垂直攝影，即航空攝影機主軸保持沿鉛垂方向進行拍照；在特殊情況下，利用專門相機進行斜傾攝影。航空攝影按所利用的電磁波波段、相應的感光片及所成圖像的特點，分成4種，即：航空可見光全色黑白圖像；航空可見光真彩色圖像：航空紅外假彩色圖像：航空紅外黑白圖像。其中，航空可見光全色黑白圖像和航空紅外假彩色圖像最為常用，它們主要利用地物波譜的寬波段反射強度特性。
2.多光譜衛星遙感階段
數字衛星成像首先是從氣象衛星開始的，在1960年TIROS-1氣象衛星提供了非常粗糙的衛星圖像，主要用來展示雲的樣式。隨後，在1970年代，美國國家海洋和大氣管理局（NOAA）發射了甚高解析度輻射感測器（AVHRR）進行氣象預報，它的地面解析度是1.1km，我們在電視氣象預報節目中看到它所獲得的雲圖。同時，從1970年代開始，相繼發射了一些搭載更高解析度感測器的衛星。如：1972年7月23日，美國國家航空和宇宙航行局（NASA）發射了第一顆專門用來進行地球表面監測和填圖的地球資源技術衛星（ERTS-U），1975年被更名為陸地衛星（Landsat）。在Landsatl-3上都裝有多光譜掃描儀（MSS），該掃描儀有4個波段，即綠、紅和兩個紅外波段，地面解析度約為80m。1982年，Landsat4搭載了專題制圖儀（TM），它有7個波段，比MSS覆蓋波譜范圍更寬，波段寬度劃分得更細些，更能反映地物反射光譜特性的變化規律，其地面解析度除第6波段為120m外，均為30m。多光譜遙感的最典型特徵是能夠利用多個波段同時獲取同一目標的多個波譜特徵。這樣就大大提高了遙感識別地物的能力。隨後各國紛紛效仿，感測器的光譜范圍從可見光、紅外直至微波波段，應用范圍也不斷擴大。
3.成像光譜遙感技術發展階段
成像光譜遙感技術是多光譜技術發展的一次跨越。Hunt的研究結果表明特徵礦物的吸收寬度大約在20～40nm，而多光譜遙感數據（例如，MSS和TM）的光譜解析度僅為100nm左右，因此遙感科學家們開始研究高光譜解析度和空間解析度的遙感感測器。1981年，一台太空梭多光譜紅外輻射計（SMIRR）隨著美國太空梭「哥倫比亞」號對地球表面進行了一次有限航帶的觀測，第一次實現了從空間通過高光譜解析度遙感直鑒別碳酸鹽岩以及粘土高嶺土礦物，由此拉開了成像光譜遙感岩性識別的新篇章。繼JPL的AIS-1和AIS-2以及AVIRIS航空成像光譜儀研製成功之後，加拿大也先後研製成功了FIL/PML，CAS1及SFSI等幾種成像光譜儀（童慶禧等，1993）。其他的還有：HIRIS（high resolution imaging spectrometer）成像光譜儀，在0.4～2.5μm范圍內有192 個光譜波段，地面解析度30m，在0.4～1.0μm波長范圍光譜解析度為9.4nm，1.0～2.5μm范圍內為11.7nm（Goetz＆ Herring 1989；Kerekes ＆ Landgrebe，1991）。美國地球物理環境研究公司（Geophysical and Environ-mental Research Corporation）的63通道成像光譜儀（GER）是專門為地質遙感研究設計的，被多次用於岩性填圖（鄭蘭芬等，1992；Bamaby W rockwell，1997）。除航空成像光譜儀外，美國和歐洲空間局（ESA）已制定了發展航天成像光譜儀的計劃，其中美國的中解析度成像光譜儀（MODIS）已經加入地球觀測系統（EOS）發射入軌，對地球實現周期性的高光譜解析度遙感觀測。歐空局的中解析度成像光譜儀（MERIS）也將於同時發射（童慶禧等，1993）。
從1990～1995年，Roger N.Clark等人先後利用AVIRIS數據在美國內華達州，卡普來特試驗場進行了礦物和岩性的識別和填圖，他們發現成像光譜儀不僅能區分地表發射光譜中總體亮度和坡度差異（多光譜技術MSS，TM和SPOT區分地物的基礎），而且能得出用於識別特殊地物的光譜吸收波段，成像光譜數據的光譜分析可以對任何在測量光譜范圍內有獨特吸收特徵的物質（礦物、植被、人T物體、水體、雪等）進行識別和填圖（Clark，R.N.et al.，1996）。
中國科學院上海技術物理研究所是我國成像光譜儀的主要研製機構。1983年研製成功了第一台工作於短波紅外光譜區（2.05～2.5μm）的6通道紅外細分光譜掃描儀，其光譜解析度在30～50nm之間。1987年，在國家和中國科學院黃金找礦任務的驅動下，該儀器發展到12個通道，其波段位置更趨於與地面粘土礦物、碳酸鹽岩礦物的吸收波段相一致，因而在地質岩性識別方面具有更大的能力（童慶禧等，1993）。另外還有熱紅外多光譜掃描儀（TIMS），19 波段多光譜掃描儀（AMSS）以及71波段多光譜機載成像光譜儀（MATS）等。這些光譜儀的數據主要用於油氣資源遙感（朱振海，1993）和礦物制圖（王晉年等，1996）等方面，數據的處理技術和礦物識別的理論研究都取得了不同程度的進展（李天宏，1997）。
綜觀遙感光譜數據的獲取，具有幾個新的發展：
①擴展了應用光譜范圍，增加了光譜波段；②提高了光譜和空間解析度；③具有獲得立體像對的功能，打破了只有航空相片才能有立體像對的能力（如SPOT圖像）；④改進了探測器性能或探測器器件，即線、面陣CCD器件；⑤提高了圖像數據精度；⑥應用領域縱向發展，如用TM圖像數據直接可以識別赤鐵礦、針鐵礦等礦物。
在20世紀末和21世紀初，空間高光譜成像衛星已成為遙感對地觀測中的一項重要前沿技術，在研究地球資源、監測地球環境中發揮越來越重要的作用。
高光譜解析度遙感技術的發展是20世紀末的最後兩個10年中人類在對地觀測方面所取得的重大技術突破之一，是當前乃至21世紀初的遙感前沿技術、通過高光譜成像所獲取的地球表面的圖像包含了豐富的空間、輻射和光譜三重信息。進入20世紀90年代後期，伴隨著高光譜遙感應用的一系列基本問題，如高光譜成像信息的定標和定量化、成像光譜圖像信息可視化及多維表達、圖像-光譜變換、大數據量信息處理等的解決、高光譜遙感已由實驗研究階段逐步轉向實際應用階段，而作為高光譜遙感應用這一熱點中的重點就是高光譜數據信息挖掘技術的提高和與之緊密相連的應用領域的擴展。
高光譜遙感數據最主要的特點是：將傳統的圖像維與光譜維信息融合為一體，在獲取地表空間圖像的同時，得到每個地物的連續光譜信息，從而實現依據地物光譜特徵的地物成分信息反演與地物識別。它由以下3部分組成：
（1）空間圖像維
在空間圖像維，高光譜數據與一般的圖像相似。一般的遙感圖像模式識別演算法是適用的信息挖掘技術。
（2）光譜維
從高光譜圖像的每一個象元可以獲得一個「連續」的光譜曲線，基於光譜資料庫的「光譜匹配」技術可以實現識別地物的目的。同時大多數地物具有典型的光譜波形特徵，尤其是光譜吸收特徵與地物化學成分密切相關，對光譜吸收特徵參數（吸收波長位置、吸收深度、吸收寬度）的提取將成為高光譜信息挖掘的主要方面。
（3）特徵空間維
高光譜圖像提供一個超維特徵空間，對高光譜信息挖掘需要深切了解地物在高光譜數據形成的二維特徵空間中分布的特點與行為，研究發現：高光譜的高維空間是相當空的，數據分布不均勻，且趨向於集中在超維立方體空間的角端，典型數據的差異性，可以映射到一系列低維的子空間，因此迫切需要發展有效的特徵提取演算法去發現保持重要差異性的低維子空間，從而有效地實現信息挖掘。
Ⅳ 光譜特徵分析與提取
6.1.1 基於光譜重排的光譜特徵提取方法
首先，針對光譜吸收特徵受雜訊影響較大的問題，對數據進行最小雜訊分量正變換，消除雜訊後，再將最小雜訊分量特徵空間的數據變換回原數據空間，即最小雜訊分量反變換；然後針對單個吸收不穩定、光照等對光譜幅值影響較大等問題，提出在連續去除的基礎上，利用所有吸收特徵並將光譜吸收特徵按吸收深度由強至弱重排，從而實現穩定、可靠的光譜特徵提取。
（1）最小雜訊分量變換
在實際應用中，地物光譜吸收特徵對雜訊敏感，因此，在進行特徵提取之前，研究中引入了最小雜訊分量變換（Minimum Noise Fraction，MNF），去除雜訊對特徵提取影響的同時去除數據相關性。
MNF變換是Green等人在主成分分析理論的基礎上改進得到的。通常被用來去除數據中的雜訊成分，同時確定高光譜數據的本徵維數，從而減少後續處理的運算量。
該方法以雜訊協方差的估計矩陣為基礎，調整雜訊的取值並去除其波段間的相關性。在結果數據中雜訊的方差為1，並且在波段間無相關性。假設高光譜數據X＝［x1，x2，…，xm］T可以表示為
X ＝ Z ＋ N （6.1）
式中：矩陣Z，N分別是理想信號和雜訊矩陣，且彼此不相關；第i 波段的雜訊分量定義為NFi ＝
∑X，∑Z和∑N分別為可觀測信號、理想信號及雜訊的協方差矩陣，並且有
∑X ＝∑XZ ＋∑N （6.2）
假設F為∑N的白化矩陣，∑N的特徵值矩陣為
高光譜遙感技術原理及礦產與能源勘查應用
高光譜遙感技術原理及礦產與能源勘查應用
式中：I為單位矩陣，矩陣
假設∑w＝F T∑X F為雜訊白化之後的觀測數據的協方差矩陣，∑w矩陣特徵值組成的對角矩陣為
高光譜遙感技術原理及礦產與能源勘查應用
於是得到最小雜訊分量變換矩陣：
高光譜遙感技術原理及礦產與能源勘查應用
由式（6.6）得觀測信號最小雜訊分量變換後的矩陣為
T ＝ MTX （6.7）
經過式（6.7）變換之後，可觀測信號各個波段間彼此不相關，且各個波段按信噪比由大到小排列
（2）光譜重排
不同地物的光譜信息是不相同的，因此，高光譜遙感提供的地物精細的光譜信息可以直接作為特徵提取與目標識別的依據，比如利用紅邊、綠峰、NDVI等特徵可以提取植被。但當不同地物之間的光譜在形狀、幅值、變化趨勢等指標大致相同的時候（即光譜特徵相似），提取區分不同地物顯著特徵是非常困難的，即地物之間的不相關性均勻地分布在各個波段；此外，由於單個光譜吸收特徵容易受到光照條件、大氣等影響使得提取的光譜特徵參量不穩定。因此，針對以上問題，研究中提出了基於光譜重排的特徵提取方法，根據光譜吸收深度的由強到弱排列，剩餘的沒有吸收特徵波段則按波長由小到大排列。
光譜重排的實現過程如下：
1）通過不同階數的微分值確定的光譜彎曲點、最大最小光譜反射率及其波長位置，計算連續統去除後目標光譜的吸收位置λM 及其反射率值ρM、吸收深度H、吸收左右肩（ρL，ρR）及其反射值（λL，λR），並且吸收深度H的計算公式如下：
H ＝ d × ρL ＋（1-d）× ρR-ρM （6.8）
d ＝（λL-λM）/（λR-λL） （6.9）
2）將目標光譜按照吸收深度H由強至弱進行排列，若無吸收特徵，則按波長由小到大進行排列；
3）以目標光譜為基譜，將圖像數據光譜按照目標光譜重排後的波長進行排序。
該方法有效地利用了高光譜遙感數據提供的地物所有吸收特徵，增加了特徵提取的穩定性和可靠性；並且通過大量的實驗發現，任何兩種不同地物的光譜通過光譜重排之後，區分不同地物的顯著特徵更加明顯，增加了類別間的可分性。
（3）演算法實現
基於光譜重排的抗雜訊光譜特徵提取方法的實現流程如圖6.1所示。該方法中為了消除雜訊對光譜吸收特徵參數提取的影響，引入了MNF變換；為了有效抑制由於光照條件、感測器等因素產生的光譜幅值變化對光譜特徵提取的影響，引入了連續統去除操作；為了克服單一特徵不穩定、不同地物光譜特徵相似等問題，提出了光譜重排的方法。
（4）實驗分析
為了驗證上述研究中方法的有效性和可行性，採用AVIRIS航空高光譜數據進行實驗分析，並利用光譜之間的光譜角進行可分性的定量化分析。
實驗數據為1995年7月在美國內華達州Cuprite礦區AVIRIS航空高光譜數據，並且使用ATREM方法校正得到了地表反射率，波段范圍為1990～2480nm，空間解析度20m，光譜解析度10nm，數據大小為255 × 350 × 50。
圖6.1 光譜特徵提取方法實現流程
該研究區域的礦物分布圖如圖6.2（a）所示，從數據中提取高嶺石光譜曲線如圖6.2（b）所示，光譜重排後的光譜如圖6.2（c）所示。高嶺石、明礬石、布丁石及熱液硅石特徵提取前的光譜比較如圖6.3（a）所示，以高嶺石光譜為基譜，光譜重排後四種礦物的光譜特徵如圖6.3（b）（圖中的光譜曲線縱坐標做了平移處理）所示。利用光譜角的方法進行四種礦物光譜重排前後可分性的比較，結果如表6.1和表6.2所示。
圖6.2 高嶺石礦物光譜比較
圖6.3 四種礦物光譜比較
表6.1 原始光譜數據四種礦物的可分性
表6.2 重排後光譜數據四種礦物的可分性
由圖6.2和圖6.3可以看出，經光譜重排後，高嶺石礦物光譜吸收特徵按吸收深度的強弱進行了重新排列，較好的顯現了高光譜所有吸收特徵及主次吸收特徵的變化；並且明礬石與高嶺石礦物在2200 nm的光譜特徵由於吸收寬度等不同而能將二者較好的區分。由圖6.3與表6.2可以看出，經過光譜重排後，高嶺石與其他三種礦物的可分性均存在不同程度的增大，特別是，高嶺石與明礬石的可分性從0.1978增加為0.225；為後續礦物識別與分類等處理奠定了良好的基礎。
圖6.4 SAM方法礦物識別結果
為了進一步驗證該方法的性能，進行了利用該方法以及基於SAM方法的礦物識別結果比對分析。利用原始光譜進行光譜角匹配識別的結果如圖6.4所示。利用基於光譜重排的抗雜訊特徵提取方法得到的數據進行礦物識別，結果如圖6.5 所示。可以看出，兩種方法均能實現四種主要蝕變礦物的識別，但是，採用原始光譜進行識別的結果中存在著一定程度的礦物混淆，並且布丁石的識別結果混淆尤其明顯；而在研究方法中進行特徵提取基礎上得到的礦物識別結果礦物混淆明顯降低，取得了較好的識別結果，證明了上述研究中提出的方法的優越性能。
圖6.5 基於光譜重排特徵提取方法礦物識別結果
6.1.2 吸收波長加權匹配方法
光譜曲線往往包含了許多由雜訊引入的無效特徵，利用同類地物光譜特徵求交，實現了有效吸收波長、吸收深度的提取；常用的SAFP匹配方法中，只有參考光譜和測試光譜的特徵在相同的波長位置時，兩條光譜才被判為相同，匹配准則比較苛刻，導致由於雜訊等因素影響光譜特徵而無法匹配，吸收波長加權匹配法利用偏移加權矩陣實現了吸收波長的容偏匹配，大大增加了匹配的准確性，降低了外界因素對吸收參量特徵的影響。
對同類地物光譜曲線特徵求交，得出識別地物的有效特徵；地物光譜的診斷吸收特徵總是出現在特定的波段上，在某些情況下會有局部的偏移；對吸收特徵的中心波長進行匹配，並容許一定程度的波段偏移，容許程度用偏移加權矩陣來度量，能夠對地物光譜實現精確的識別。考慮到實際應用雜訊及系統誤差引入的干擾，用吸收深度對單個中心波長進行加權，吸收深度小的吸收特徵對整體相似度的貢獻小，吸收深度大的吸收特徵對整體相似度的貢獻大，這樣一定程度上抑制了無法去除的非有效特徵的影響。
（1）吸收波長加權匹配的實現
有效吸收特徵的精確提取和容偏匹配實現流程如圖6.6所示，具體包含以下幾個步驟：
1）對參考光譜連續統去除。利用導數法確定各吸收特徵的中心位置和左右肩對應的波長後，利用下列公式提取吸收特徵中心波長和吸收深度：
高光譜遙感技術原理及礦產與能源勘查應用
式中：
沒有標准參考光譜時，參考光譜通過訓練樣本得到。通過上述方法提取各條參考光譜的吸收中心波長和吸收深度後，對所有訓練樣本的吸收特徵參數求交，方法如下：
光譜A和B的所有吸收特徵為feature_a，feature_b，A的第i個波段上存在特徵，對feature_b計算：
judge ＝ Weight·feature_b（［i-BandOffset：i ＋ BandOffset］） （6.11）
如果，judge＞0 ，則光譜A的第i個波段上的特徵為有效特徵。
得到參考光譜共有的有效特徵，此處需要記錄的是有效特徵的位置和吸收深度的大小，保存在向量EffFeatureIndex和Depth中。
2）提取未知光譜所有吸收位置和對應的吸收深度特徵，記錄在FeaturePos和FeatureDepth；
圖6.6 中心波長加權匹配流程圖
未知光譜特徵與參考光譜有效特徵按位匹配，匹配方法包含兩個參數，容許波段偏移數BandOffset和偏移加權矩陣Weight。
3）找到參考光譜第i個特徵位置，生成特徵檢驗區間：
TestIndex ＝（i-BandOffset）：（i ＋ BandOffset） （6.12）
計算特徵檢驗值：
TestValue ＝ Weight·FeatureDepth（TestIndex） （6.13）
TestValue不為0 ，則說明未知光譜對應位置存在有效特徵，反之則不存在，未知光譜中的識別特徵所在波段記錄在向量EffIndex中。
4）重復3）的過程，直到對未知光譜的所有有效特徵進行了檢測，未知光譜中識別特徵存在的波段記錄在向量EffIndex中。
5）對吸收位置用吸收深度加權匹配，匹配度的計算公式如下：
高光譜遙感技術原理及礦產與能源勘查應用
6）根據匹配度degree的值判斷未知光譜與參考光譜的近似程度，閾值Thresh手動選擇，根據經驗，在用吸收深度加權的匹配方法中，Thresh＝0.8就能獲得較高的識別率。
用吸收深度加權對吸收特徵中心波長進行容偏匹配的關鍵在於：有效吸收特徵的准確提取和偏移加權矩陣Weight或容許波段偏移數目BandOffset的選擇，反射率曲線所有吸收特徵的精確提取是前提，偏移加權矩陣的確定需要根據對像光譜的采樣間隔來確定，Weight的分量的個數為2 × BandOffset＋1；並且有效特徵提取和特徵識別過程使用的偏移加權矩陣Weight可以不同，光譜采樣間隔較大時，可以選擇Weight的各個分量服從高斯分布。
（2）基於USGS光譜庫數據的實驗與結果分析
圖6.7（a）為USGS礦物光譜庫中六條綠泥石連續統去除後反射率曲線；波段偏移參數BandOffset＝1，對應的容偏矩陣Weight＝［1，1，1］；即兩條光譜的特徵相差一個波段以下認為該特徵為有效特徵；綠泥石的有效特徵見圖6.7（b），用方框標記出了吸收谷的波長位置；圖6.7（c）給出了利用吸收波長加權匹配方法得到的綠泥石有效特徵；圖6.7（d）給出了綠泥石和陽起石反射率光譜。
圖6.7 有效特徵提取
匹配加權矩陣Weight＝［1，1，1］表示容許兩端偏移，Weight＝［0.1，1，0.1］表示不容許偏移；兩情況對應的相似度見表6.3和表6.4。對比表6.3和表6.4的相似度值可以看出，容許波段偏移後，綠泥石光譜間的相似度明顯變大。利用圖6.7（c）的有效特徵對圖6.7（d）所示的陽起石和綠泥石光譜進行Weight＝［1，1，1］匹配，近似度見表6.5，用綠泥石的有效光譜能有效的識別出綠泥石光譜與陽起石光譜的差異。
表6.3 綠泥石光譜識別Weight＝［1，1，1］
表6.4 綠泥石光譜識別Weight＝［0.1，1，0.1］
表6.5 陽起石和綠泥石識別Weight＝［1，1，1］
（3）基於AVIRIS數據的實驗與結果分析
利用內華達州Cuprite礦區的AVIRIS數據進行基於吸收波長加權提取方法實現礦物匹配識別研究。利用的礦物端元光譜如圖6.8所示，識別結果如圖6.9所示。
從地質圖6.2（a）與結果圖6.9比較可以看出，該方法對具有明顯光譜吸收特徵的明礬石和高嶺石礦物具有較高精度的識別效果，但是對於吸收特徵較寬、較淺的白雲母和布丁石的識別效果則較差。
圖6.8 演算法中用到的端元光譜
圖6.9 基於吸收波長加權特徵提取的礦物匹配識別結果
Ⅳ 端元光譜選取與信息提取
3.1.3.1 端元數目的確定
一個混合像元可能包含了幾種地物類別，所以確定端元的數目是進行光譜線性分解的前提，也是整個光譜線性解混技術中必不可少的一個環節。對於多光譜數據常用的方法是根據主成分分析（PCA）中協方差的大小來判定端元數，但是其分析方法比較粗糙，而對於含有上百個窄波段的高光譜遙感影像來說，主成分分析方法很容易把細微的光譜信息歸到雜訊部分（Chang C，2007）。因此，目前常用的針對高光譜遙感影像確定端元數目的方法是基於Neyman－Pearson探測理論的特徵閾值分析方法（Harsanyi et al.，1994），簡稱HFC（Harsanyi，Farrand，Chang）。本章也採用此方法進行端元數目的確定。
HFC原理是通過對影像進行矩陣的相關計算，分別得到其相關矩陣Rm×n和協方差矩陣Km×n及其特徵值，並把特徵值分別記作 。
如果影像的信號能量為正，則有
高光譜遙感影像信息提取技術
式中：n即為端元數目；m為高光譜影像的波段數。
即便是沒有採用主成分分析的方法，也不可避免的有部分弱信號被認作是雜訊信號被排除掉，因此，為了盡可能地減少這種現象的發生，在計算端元數目之前最好進行白化處理（Gruninger et al.，2004）。
3.1.3.2 端元的提取
高光譜影像端元提取的方法目前研究的較為深入，研究者從不同的角度提出了很多實用性的提取方法，其中較為常用的有純像元指數（PPI）、內部最大體積法（N－FINDR）、頂點成分分析（VCA）、單形體投影方法（SPM）、順序最大凸錐（SMACC）、迭代誤差分析（IEA）、外包單形體收縮（SSWA）、最小體積單形體分析（MVSA）、凸錐分析（CCA）、光學實時自適應光譜識別系統（ORASIS）、自動形態學（AMEE）、最大距離法（MaxD）、最大體積法（MaxV）、最大零空間投影距離法（NSP）、定量化獨立成分分析法（ICA）等（張兵等，2011）。本章採用的端元提取的方法是順序最大凸錐（SMACC），它提供了更快、更自動化的方法來獲取端元波譜，但是它的結果近似程度較高，精度較低。由於本章的研究目的是一種改進型的線性分解方式，不是集中在端元選擇問題上的研究，因此使用SMACC雖然不是最好的演算法，但是完全可以達到本章實驗的目的。
SMACC演算法（Gruninger et al.，2004）在提取端元的同時可以獲得豐度反演的結果，它的基本原理是通過迭代的方法來獲取端元，經過數次迭代，每次都不斷地計算和調整各個端元在混合像元中所佔的比例，並且利用投影變換消除端元之間的相互影響。其中，最關鍵的步驟就是判斷該像元中是否有該端元，並且是否需要進行斜交投影（或正交投影）。具體演算法如下：
設原始像元集表示為 ，其第j次迭代前的像元集合表示為 ，第j次迭代前的端元集合表示為 ，wj表示每次迭代時的投影方向，Xj－1為最長的光譜向量，則Xj－1在wj方向上的投影系數為
高光譜遙感影像信息提取技術
ej在xi中的比例系數為
高光譜遙感影像信息提取技術
式中：βij為調整系數，當βij＝1時為正交投影，否則為斜交投影。經過多次迭代，最終可以得到的ej在xi中的比例系數 。
像元的投影結果為
高光譜遙感影像信息提取技術
其中，系數βij調整的原則為：當Oij≤0時，βij＝0，表示沒有該端元；否則，根據
高光譜遙感影像信息提取技術
通過計算vk得到其中的最小值，記為vmin。當vmin＞1，則βij＝1，為正射投影；否則為斜交投影，βij＝vmin。
3.1.3.3 混合像元地物信息提取及其分類
經過對混合像元光譜曲線構建的矩陣進行分解，可以得到每一類端元光譜在混合像元中的豐度值fj（j＝1，2，…，n），但是對高光譜影像的分類或地物信息提取是按照以像元為基本單位進行劃分的，也就是說在分類或地物提取中，一個像元不可能被塗上不同的顏色，因此，為了便於分類，選取fj中值最大的一個對應的端元作為該混合像元的地物種類，即Max（fj）（j＝1，2，…，n）對應的端元。
Ⅵ 如何根據高光譜圖像數據提取其空間特徵
特徵提取和特徵選擇是模式識別的關鍵問題之一，它影響到分類器的設計及其性能．高光譜圖像數據
是超高維多特徵數據集，如何實現高維特徵空間的特徵壓縮和特徵提取是一個重要課題．基於高光譜圖像譜圖合一、數據維度高的數據結構特點，該文從光譜和圖像兩個層面分別綜述了主成分分析、最小雜訊分離、獨立成分分
析等光譜特徵提取方法以及基於顏色、紋理、形狀等圖像特徵提取方法．還詳細介紹了核主成分分析和投影尋蹤方
法這兩種高光譜特徵提取新方法，並給出了以上方法的應用實例．特徵提取和特徵選擇的研究將為後續的高光譜圖
Ⅶ 光譜維特徵提取方法
特徵是對象所表現出來的各種屬性與特點。在遙感圖像分析中特徵提取可以從兩個意義上來實施：一種是按照一定的准則直接從原始空間中選出一個子集（即子空間），實踐中的波段選擇即屬於此類；另一類是在原始特徵空間和新特徵空間之間找到某種映射關系P，P：x→y，將原始特徵空間x=｛x1，x2…，xn｝ 映射到維數降低了的特徵空間y中去，y=｛y1，y2…，ym｝，m＜n。對於用於分類目的的特徵提取，好的特徵提取方法能使同類物質樣本的分布具有密集性，而不同類物質的樣本在特徵空間中能夠隔離分布，為進一步分類打下良好基礎。因為高光譜數據具有波段多、波段間相關性高及數據冗餘度高等特點，所以對高光譜遙感數據的特徵提取具有特殊意義。遙感圖像特徵提取包含的內容非常廣泛，提取方法也很多，光譜維特徵提取和空間維特徵提取是表現圖像特徵提取的兩種主要方法。這里主要介紹適用高光譜數據的一些光譜維特徵提取方法，主要涉及主成分分析法，典範變數分析法及改進的CA方法。
主成分分析是一種把原來多個指標化為少數幾個相互獨立的綜合指標的一種分析技術。對波段間高度相關的數據非常有效（Cloutis，1996）。PCA技術已被用在不同的地質遙感項目，包括寬波段和高光譜數據（Lee等，1990；Resmini等，1997，Fujimura ＆ Kiyasu，1994）。由於高光譜數據波段間的相關性、高冗餘度，直接利用所有的原始波段作分類或特徵提取顯得很不經濟。因此先對原始數據作PCA變換，然後對少數幾個綜合指標（成分）分析將會收到事半功倍的效果。在高光譜數據分析中，PCA技術可將總體大部分方差集中在前面少數幾個主成分中。於是，人們利用這少數幾個主成分做一些地質分析，如利用前3個主成分的假彩色合成圖判讀地質礦物信息，進而成圖。但在主成分合成圖上的彩色在不同的圖像上是變化的，並不代表一定的地質礦物成分，除非有相似的地質露頭和覆蓋，更困難的是，我們不能根據岩石、土壤和礦物等反射光譜作指示來判讀主成分合成圖上的彩色。另外，確定每個主成分的物理意義也相當困難。再者，礦物的顯著變異可能僅引起光譜的細微差異，這種細微差異常被淹沒在高方差的主成分中而被忽略。因此PCA可能較適宜用來粗略地識別光譜差異顯著的礦物和岩性類別，而不是定量的礦物識別和制圖（Coutis，1996）。
Jia＆Richards（1999）發展的分塊主成分分析法用於特徵提取，取得了一定的分類和顯示效果。類似3.3，4中建立的SMLDF判別函數的思想，他們將全部波段的相關陣按照相鄰波段的相關性分成若干塊。一般高度相關的塊沿對角線分布，而相關性低的塊遠離對角線。塊矩陣本身集合了相鄰波段間相關性高的波段。因此沿相關矩陣對角線可分成若干塊（波段組），並對每組進行主成分變換，最後將每組的重要特徵（主成分）再重新組合在一起作為進一步主成分分析與特徵選擇之用。
典範分析也是將較多的變數化為少數幾個典範變數，通過這較少的典範變數之間的相關性來綜合地描述兩個多元隨機變數之間關系的一種數學方法（唐守正，1986）。假如我們有兩個多元隨機變數（設x為p維隨機變數，y是q維隨機變數），如何描述這兩個多元隨機變數之間關系的緊密程度呢？直接的方法就是逐一計算兩個多元隨機變數各分量之間的相關系數或其他相似系數，可計算出p×q個相關系數。但這樣做既繁瑣，也不能本質地說明這兩個隨機變數總體相關水平。類似主成分分析，從每個多元隨機變數中造就數個「綜合變數」——典範變數。在求算兩隨機變數各自典範變數過程中得到的特徵根即為對應典範變數對的典範相關系數。如果將非零特徵根按從大到小排列，則最大的特徵根即為第一對典範變數（分別對應x和y）的相關系數，如果典範相關系數越大，則說明這一對典範變數關系越緊密。一般在實踐中只取前面k（k＜p，q）個典範相關系數和典範變數進行分析，捨去後面的典範變數已無關緊要，這達到了特徵提取的目的。在利用高光譜數據進行岩性識別分類時，首先可以將高光譜圖像數據分成若干未知岩性的類別，然後在每個類別中抽取一定的樣本（象元）數，同時抽取已知岩性一定的樣本數，組成兩個類似前述的多元隨機變數（一個為已知岩性，另一個為未知岩性）的樣本並計算它們前面數對典範變數。如前面數對典範變數（通常為3對）關系緊密（反映在它們相應的特徵根上），這個未知岩性的類別就很有可能與已知岩性為同類岩性，反之就可能不是。實踐中前面數對典範變數關系緊密程度靠經驗裁定。以此類推，可以將研究區內所有未知岩性的類別與已知岩性的類別求算兩兩典範變數對，並根據它們各自的相關緊密程度和判據決定未知岩性類別的歸屬。
PCA主要想最大限度地將不同類別分開，而典範分析則是在低維變數間尋找能代表高維變數的相關性，以達到分類、識別目標物的目的。典範分析在高光譜地質應用中潛在的功用與PCA技術大部分是一致的。
劉建貴（1999）在分析K-L變換性質的基礎上，根據高光譜數據用於城市目標物識別提取的特點，提出了面向分類的特徵提取的CA改進方法。選擇適當的變換矩陣，同時考慮類內與類間距離的CA方法，設法使原特徵空間的各類的樣本點在光譜維上的投影能使類間距離與類內距離的比值達到最大。這種比值稱為廣義瑞利商。根據這一原則來決定變換矩陣的選擇。劉建貴（1999）用這種改進的方法實施對北京市沙河鎮城市地物特徵的提取。具體處理過程：①對原始高光譜圖像進行預處理，得到相對反射率圖像；②選擇最終成分光譜，即需要分類的類別數，最終成分光譜根據圖像及地面調查的情況進行，共找出11個類別；③對每一類統計出均值向量和協方差陣；④求出每兩個類別對之間的類間和類內距離，利用CA變換方法求出變換特徵；⑤求出該兩個類別對應於每個原始波段的巴氏距、載入系數（劉建貴，1999）以及變換域每個特徵的巴氏距、每個特徵度量維上的類間類內距離比，載入系數可通過CA變換成分與波段間的相關性求算；⑥找出最優特徵。實驗結果表明這種特徵提取方法非常有效。採用這種方法，提取的特徵能夠增加樣本的類內凝聚度和擴大類間距離，同時消除波段問的相關性，因而能改善分類性能。
Ⅷ 基於光譜特徵的信息提取與分類過程
地物光譜特徵研究是現代遙感技術的重要組成部分。它既是感測器波段選擇和設計的依據，又是遙感數據分析解譯的基礎。遙感探測是成像空間地理實體的電磁波譜和輻射能特徵。具有明確的物理意義，而基於光譜特徵的信息提取與分類，是通過遙感光譜數據的變化規律來識別和研究地物類型。因此，研究不同地理實體的光譜表達模型，是有效地提取專題信息的關鍵。
1.金屬硫化物礦床近礦圍岩岩石光譜特徵研究
圖3.3.1是東秦嶺地區比較有代表性的幾種近礦蝕變岩及含礦體的反射波譜曲線。從波譜曲線形態可以看出，金屬硫化物蝕變帶在藍綠波段（0.4～0.6μm）和近紅外（0.85～1.1μm,2.2～2.4 μm）波段呈強吸收特徵；在紅光（0.6～0.85μm）波段和近紅（1.28～1.46μm）波段出現強反射峰。硅化蝕變岩及含羥基的蝕變岩類反射波譜曲線形態比較接近，其吸收帶仍然位於0.4～0.55μm、0.85～1.1μm、1.9～2.3μm波段內，而在0.6～0.85 μm、1.48～1.88 μm波段內出現反射肩。
圖3.3.1東秦嶺地區近礦蝕變岩與礦化體反射波譜曲線
①高嶺土化蝕變岩；②硅化蝕變岩；③金屬硫化物蝕變岩
圖3.3.2是東秦嶺地區近礦蝕變帶與近礦圍岩的反射波譜對比曲線，成礦圍岩與金屬硫化物帶波譜曲線相比，波譜響應趨向平緩，不出現大幅度跳躍的波峰或波谷。從波譜曲線形態可以看出，隨金屬硫化物帶（礦體）的遠離（礦化蝕變的減弱），波譜曲線走向平緩的趨勢更加明顯。其中安山岩類和大理岩類波譜形態相似，除在可見光藍綠光段（0.4～0.6μm）出現一個較強吸收帶外，紅光到近紅外光段基本是光滑的弧頂向上的曲線，雖然在1.4 μm、1.9μm處有弱吸收顯示，但不出現明顯的吸收谷和反射肩；片麻岩和花崗岩類為一條平滑的直線，基本無異常顯示。從不同蝕變類型與不同近礦圍岩的波譜對比中還可以看到，波譜曲線在0.4～1.4μm和1.9～2.5μm兩個波段區間集約，僅在1.4～1.9 fμm波段區間呈離散狀態，顯示出較大的反射差。
圖3.3.2東秦嶺地區近礦蝕變岩及近礦圍岩反射波譜曲線
①千枚岩；②片麻岩；③大理岩；④安山岩；⑤硅化蝕變岩
上述反射波譜曲線特徵表明，蝕變岩與非蝕變岩類的反射波譜有較明顯的差異，其中以1.4～1.9 μm波段離散程度最好，即Landsat-TM5波段應為提取礦化蝕變信息的最基本波段。1.9～2.4μm區段亦有一定的離散傾向，故Landsat-TM7可選作輔助波段；0.4～0.6μm波段相對1.4～2.5μm波段為強吸收帶，Landsat-TM1、2可作為理想的匹配波段。
圖3.3.3是闊葉林（接骨木）在不同狀態下的反射波曲線。從曲線形態可以看出，隨著植物病害程度的加重在0.6～0.7μm、1.4～1.6μm和1.9～2.4μm波段的吸收逐漸躍起，相反從0.7～1.4 μm段的反射峰跌落。從曲線的離散、集合特徵來看，0.6～1.8μm、2.0～2.5μm波段離散程度最好。因此，處於其間的Landsat-TM4、7波段可作為基本應用波段，0.4～0.5μm波段相對為強吸收帶，處於其中的Landst-TM1可作為基本匹配波段。
根據代數運算的原理，當波段間反射率差值相近而曲線斜率不同時，反射波段與吸收波段的比值處理，可在一定程度上擴展地物波譜的差異性，顯示出動態范圍。表3.3.1是根據地物反射波譜數據計算的不同組合比值數值表，從表中可以看出，作為提取蝕變岩帶信息的基本應用波段、輔助波段和匹配波段Landsat-TM 5、7、1可以最大限度的顯示出以鐵帽、硅化、高嶺土化、絹雲母化為特徵的蝕變岩與圍岩背景的波譜差。如需要進一步區
圖3.3.3河南上宮金礦區植被（接骨木）不同狀態下的反射波譜曲線
1—正常；2—弱毒化期；3—強毒化期
分蝕變帶的類型，則以7/1、7/2、5/1、5/2的比值合成處理可在一定程度上突出以Fe3＋為主體的蝕變岩信息。據實驗研究，混合比值處理［(TM3×TM4)-K]/TM7則極大地壓縮非蝕變背景信息，比值合成TM5/1(R)＋TM7/1(B)＋［(TM3×TM4)－K]/TM7(G）可突出蝕變信息，特別是以羥基（OH）礦物為代表的蝕變岩信息。在比值合成圖像上，金屬硫化物帶應為暗紅色，高嶺土、絹雲母化帶應為亮黃色，非蝕變背景接近青色（見彩圖）。
表3.3.1河南豫西地區遙感成礦特徵波段比值一覽表
表3.3.1是根據植被在不同狀態下的反射波譜數據計算出的不同組合比值數值。從組合比值數值中可以看出，作為提取植被受重金屬離子毒害信息在Landsat-TM4、7波段得到了比較好的波譜差。如果突出受害嚴重的區域，需要有混合波段比值處理，如COSTM4×TM7-K等。經實驗研究，比值合成［COSTM4×TM3-K］(R)＋TM4/TM1(B)＋TM7/TM1(G）圖像反映重害區為橘黃區，輕害區為接近白色，正常區接近青色；也可以用TM4(R)+TM7(B)+[COSTM4×TM7-K](G）圖像，正常區為紅色，輕害區接近粉紅色，重害區接近白色。
2.反射波譜特徵模式的應用原則
由於岩石反射波譜測試大多是在室內條件下進行的，而遙感感測器所記錄的岩石反射率則為自然條件下的反映。岩石的反射波譜特徵在自然環境中受植被、土壤、水分等因素的干擾，往往會產生較強的畸變。對於近礦蝕變帶來說，遙感衛星所提供的礦化蝕變信息，常常與礦化蝕變帶的水、土及植被等信息有關。當礦化蝕變帶中含有較高的As、Hg、Pb等有毒元素和重金屬元素時，植被因此可以出現較大范圍的毒化反映；當蝕變帶遭受強烈的風化剝蝕時，殘留在地表多是硅質和含羥基的粘土質礦物集合體，而這種風化殘積物的規模可遠遠超過蝕變帶的分布范圍。因此，在蝕變信息的提取過程中必須考慮干擾因素的存在，根據當地礦化蝕變岩的出露特徵確定提取的目標物（植被毒化、Fe3＋、OH－），根據反射波譜的特徵模式採取不同的運算公式，加以補救和篩選。
如秦嶺山地植被蓋率在70%左右，林木繁盛地區可達90%以上。該區主要植被類型為櫟類與油松、華山松、白樺等混交林及灌草叢。據區域地球化學研究，崤山、伏牛山地區有色金屬、貴金屬金礦床多伴生有以 As為主的有害元素和以Cu、Pb、Zn、Mo為主的重金屬元素，這些元素會在金礦蝕變帶的分散暈圈半徑內產生植被的有限毒害暈。根據這一特徵，選擇在風化殘積物較薄、植被受毒害相對明顯的伏牛山北坡為試驗區，以植被生長狀態反應比較敏感的Landsat-TM3、4、5、7波段為基礎進行混合比值處理，工作程序見圖3.3.4。
圖3.3.4河南省伏牛山地區植被毒化暈的提取工作程序圖
TM5/TM4、TM4/TM3比值是最佳的植被指數。它們對植被毒害程度的反映是一個反演序列，即植被從正常發育到毒害變態反射率分別在0.36～0.61～1.20(TM5/TM4）、6.25～3.3～1.3(TM4/TM3）降低（或稱藍移現象）。TM5/TM1對鐵帽反應比較敏感，對植被中度以上毒害反應突出，它們的比值系數為3.38、2.90，分別高出背景值一倍以上。植被毒化區及蝕變區為高頻域，通過濾波可壓縮背景低頻信息，高頻域均可以鮮明的色調給以增高突出。
如圖版Ⅰ.1所示，TM4/TM3(B）高值區位於植被毒害和蝕變岩區，R、G近等量合成後呈亮黃色調的局部斑塊，即代表與礦化有關的信息。應指出的是，這種黃色斑並不確切指出礦化蝕變的現存位置，因分散暈的遷移造成的位置差，色斑往往位於礦化蝕變岩的下游。
熊耳山－外方山區，礦化蝕變均沿構造破碎帶發育。構造破碎帶在多數情況下為負地形，除硅化體礦化蝕變帶有斷續出露外，大部分地段被褐鐵礦化及粘土質物質所充填。其主要成礦圍岩為太古界太華群花崗－綠岩系和元古界熊耳群安山岩類。因此，需要增強的信息是發育在構造帶中的富含羥基的粘土礦物集合體。
根據地物反射波譜資料，在Landsat-TM3、5、7波段目標物和背景的離散比較好。如進一步對蝕變帶進行劃分，還需要對以Landsat-TM3、5、7波段為基礎進行混合像元處理，工作程序如下圖3.3.5所示：
圖3.3.5羥基粘土礦物暈的提取工作程序圖
從岩石反射波譜特徵模式中可以了解到，TM7/TM1、TM5/TM4的比值結果突出了蝕變暈帶中的核心——金屬硫化物氧化鐵帽，TM3/TM4的比值結果突出了植被受重金屬毒害信息。傅氏變換將比值數據轉成頻率域，然後通過高濾波壓抑低頻背景，使蝕變信息得到很高的增強；經反變換（IFT）將信息還原到空間域，用已知蝕變區作為樣本確定彩色分割的閥值，以此閥值為標准進行假彩色漫遊或假彩色密度分割，即可得到反映圖版Ⅰ.2。
在反射波譜特徵模式中提出的壓縮背景混合像元處理方法，即［(TM3×TM4)－K］和［COS-TM4×TM7]-K，這里的K值代表礦化蝕變帶圍岩的反射率（背景值）。壓縮背景的方法在圖像處理中叫「分段線性擴展」，即將圖像亮度值的整個動態范圍分成若干區段，按區段進行不同程度的擴展（圖3.3.6）。
圖中的L1為原圖像的亮度值變數，L2為變換後影像的亮度值變數。a1、a2、a3分別為所選擇的分段斷點。斷點之間的斜率，控制區段內亮度值的變換。
圖中k1、k2、k3分別為對應區段內變換曲線的斜率。適當選擇斷點和斜率，可以獲得特定亮度值區內目標圖像的對比度增強，或者壓縮某些目標的對比度。
圖3.3.6分段線性擴展示意圖
圖版Ⅰ.3是為提取導礦斷裂信息而設計的圖像處理程序得到的圖像。為了突出斷裂構造帶中不同地段的植被富水性及蝕變等特徵信息，分別應用TM4/TM3、TM5/TM1、TM5/TM2進行比值處理。對比值圖像數據進行例拉伸（SCALE）得到灰階為0～255的灰度圖像。然後分別在比值圖像中找出斷裂特徵信息的亮度值區間（192、128、115）為斷點，將非目標區壓縮為0，將目標區給予較大的擴展斜率，將擴展後的圖像合成具斷裂構造意義的線性體，以鮮明色調給以確切的顯示。如將蝕變暈斑疊加其上，該斷裂的控礦意義更加明確。
Ⅸ 高光譜岩性信息提取
8.5.1 方法與流程
8.5.1.1 岩性信息產品生成業務化流程
採用的高光譜岩性信息分類填圖方法：先將成像高光譜數據進行輻射定標、大氣校正、光譜重建等過程得到光譜反射率數據，然後根據礦物特徵吸收峰分布情況對高光譜反射率數據進行波譜降維，對照已有地質圖選取其中的典型岩性並結合像元純凈度指數圖在高光譜數據區域中建立感興趣區，確定岩性分類的先驗樣區，最後用合適的方法進行岩性信息分類填圖，並對結果進行優化操作。岩性信息分類填圖流程圖如圖8.20。
圖8.20 高光譜岩性信息分類填圖流程
8.5.1.2 岩性分類信息提取的高光譜數據預處理
無論是高光譜成像儀還是傳統的多光譜感測器，它們所記錄的數據都是地面觀測目標的反射或輻射能量的光譜輻射絕對值，與地物目標的光譜反射率或光譜輻亮度值是不一致的。因此，輻射定標和光譜重建是地物識別和定量分析不可缺少的環節。只有經過輻射標定、輻射校正和大氣校正，剔除由於大氣散射、吸收、地形起伏及感測器本身不穩定帶來的各種失真，將記錄的圖像值轉換為地面的反照率值，重建像元地面光譜，才能根據光譜特徵，有效地識別地物，反演地物成分。
8.5.1.3 端元選擇
對於一個地區的未知岩性分類信息提取，一般要進行岩性端元信息的提取。所謂端元，指的是成分單一的岩性像元。只有提取出端元才能進一步對研究區域進行岩性信息的分類分析。目前端元選擇的方式概括起來主要有兩種：①根據野外波譜測量或從已有的地物波譜信息庫中選擇端元。通過這種途徑選擇的端元稱為「參考端元」。②直接從待分類岩性的圖像上選擇端元，然後不斷對其修改、調整，確定端元，這種圖像上選擇的端元區域稱為先驗樣區。
野外獲取岩性端元信息一般要經過實地勘察，先選擇好樣區，然後選擇合適的時間進行量測。一般情況下，要求在獲取影像的同時進行地物波譜量測，但這種難度很大，實際應用中很少能做到。通過野外測量方式獲取的參考端元理論上比較精確，但遙感圖像上地物的波譜曲線受到大氣、地形和感測器等的影響，這些與野外實地測量的地物波譜曲線存在很大差別，即使對影像進行各種糾正進一步消除這些因素的影響，也不會與野外測量的地物波譜曲線很好地吻合。對於大部分研究區域而言，有與研究區對應的實際地物波譜庫的區域很少，也很少投入大量人力、物力進行野外波譜測量。因此在缺乏野外波譜測量數據情況下，從影像本身獲取端元是目前獲取端元的主要方式。目前，在端元從影像本身的像元獲取方式中，除了對遙感影像像元的色調、波譜特徵進行目視解譯直接分析外，還可藉助以下分析方法進行。
基於圖像的端元選擇是假定圖像中存在有基本上僅反映一種岩性或礦物光譜的「純像元」或「非混合像元」，用數學方法自動或交互地從圖像中提取這些「純像元」作為端元，用同一類「純像元」的典型光譜或平均像元光譜作為端元光譜。PPI方法由於推出較早，計算簡便直觀，並且處理高光譜數據中使用最廣泛的遙感圖像處理系統ENVI中有該功能，因而得到較廣泛的應用。但是在進行PPI之前，需先對高光譜數據進行降維處理。我們常用的降維處理手段主要是最大雜訊分離（MNF）變換。
（1）最大雜訊分離（MNF）變換
成像光譜的光譜解析度很高，波段多，數據海量，且波段之間的相關性很強，數據冗餘度高，需要在端元選取和礦物識別之前對數據作減維，並弱化雜訊。最常用的方法是「最大雜訊組分變換」。
MNF變換是利用圖像的雜訊組分矩陣（∑N∑-1）的特徵向量對圖像進行變換，使按特徵值由大到小排序的變換分量所包含的雜訊成分逐漸減小，圖像質量順次提高。∑為圖像的總協方差矩陣；∑N為圖像雜訊的協方差矩陣。MNF有兩個重要的性質，一是對圖像的任何波段作比例擴展，變換結果不變；二是使圖像矢量、信息分量和加性雜訊分量互相垂直，達到信息分離的目的。乘性雜訊可通過對數變換轉換為加性雜訊，變換後可針對性地對各分量圖像進行去噪，或舍棄雜訊占優勢的分量，以達到成像光譜數據減維和去噪的目的。
（2）像元純度指數（PPI）分析
基於上述MNF變換，排序低的MNF波段被暫時忽略，僅選擇高序次波段進一步處理。PPI設計指定光譜極值像元，對應為混合光譜端元。通過反復投影n維散點圖到隨機單元矢量來計算PPI。記錄每次投影的極值像元，注記每個像元被標定為極值的總次數。PPI圖像產生，其中，每個像元的DN值對應像元被記錄為極值的次數。這些圖像的直方圖顯示被PPI「擊中」（hit）的分布。從直方圖中選擇閾值，用於選擇最純的像元以保證被分析的像元數最小。這些像元被輸入到分離特定光譜端元的互動式可視化演算法中。
（3）N維可視化（n-Dimensional Visualization）
由於PPI演算法本身確定的並不是最終的端元，而是從圖像眾多像元中選出包含所有端元像元的較小像元子集。所以，一般要將PPI的處理結果輸入到N維可視化（n-Dimensional Visualization）工具中，選擇出最終的端元像元。
在N維可視化中，光譜可視為n維散點圖中的一個點，n是波段數。對給定的像元，n維空間中組成n值的點的坐標是每個波段的光譜反射率。這些點在n維空間的分布可用於估計光譜端元數和它們的純光譜特徵。在兩維空間，如果只有兩個端元混合，混合光譜將落入直方圖的線中。純端元將落入混合線的兩端；如果三個端元混合，混合像元將落入四面體中；余類推。混合的端元落在純端元之間，處在純端元勾畫的多面體中。這種混合光譜的凸面幾何特徵，可用於確定端元光譜數並估算它們的光譜特徵，經過PPI選擇出的潛在端元光譜輸入n維散點圖中進行反復旋轉以識別出純端元。根據前面的分析，較好的端元通常會出現在n維散點圖的頂點和拐角處，當一系列的端元點被確定後，就可以將其輸入到圖像中的感興趣區（ROI），從圖像中提取每個感興趣區平均反射率光譜曲線作為成像光譜礦物填圖的候選端元。
鑒於PPI演算法得出的結果是從圖像眾多像元中選出包含所有端元的較小像元子集（圖8.21），在已有的地質圖中結合PPI結果圖可選取盡量純凈的端元區域作為後期監督分類的感興趣區（ROI，Region of Interest）。
圖8.21 選擇出的純凈像元的波譜曲線
8.5.2 結果與分析
目前，從是否需要先驗樣區可將岩性分類技術分為兩大類：非監督分類和監督分類。非監督分類是指人們事先對分類過程不施加任何的先驗知識，而僅憑數據遙感影像地物的光譜特徵的分布規律，即自然聚類的特性，進行「盲目」的分類；其分類的結果只是對不同類別達到了區分，但並不能確定類別的屬性；其類別的屬性是通過分類結束後目視判讀或實地調查確定的。非監督分類也稱聚類分析。一般的聚類演算法是先選擇若干個模式點作為聚類的中心。每一中心代表一個類別，按照某種相似性度量方法（如最小距離方法）將各模式歸於各聚類中心所代表的類別，形成初始分類。然後由聚類准則判斷初始分類是否合理，如果不合理就修改分類，如此反復迭代運算，直到合理為止。監督分類又稱訓練場地法，是以建立統計識別函數為理論基礎，依據典型樣本訓練方法進行分類的技術。即根據已知訓練區提供的樣本，通過選擇特徵參數，求出特徵參數作為決策規則，建立判別函數以對各待分類影像進行的圖像分類，是模式識別的一種方法。要求訓練區域具有典型性和代表性。判別准則若滿足分類精度要求，則此准則成立；反之，需重新建立分類的決策規則，直至滿足分類精度要求為止（圖8.22）。
圖8.22 地質圖
對原始數據的211個波段進行目視解譯，去掉其中明顯的雜訊波段得到剩下的178個波段，對這178個波段的數據信息用上述方法進行分類得到所需的結果，然後結合各種岩性地質特徵，得知大多數礦物特徵吸收峰均位於1300 nm以後的波段中，故在這178個波段中繼續裁剪得93個波段，並對裁剪結果用相同的方法及參數重新分類，將得到的結果與第一次分類結果進行對比。
8.5.2.1 非監督分類填圖結果
（1）K-Means
預選待分類數據可分為八類，最大迭代次數為2，且設變化閾值為5%，即當每一類像元數變化小於閾值時結束迭代過程，最大允許標准差和最大允許距離誤差不輸入，即所有像元都參與分類。從而，當達到閾值5% 或迭代達到兩次時則分類結束（圖8.23 ，圖8.24）。
圖8.23 所用數據：dts_178，K-Means填圖結果
圖8.24 所用數據：resize_dts_93，K-Means填圖結果
（2）ISODATA
預選待分類數據可分為5～10類，最大迭代次數為2 ，變化閾值為5%，每一類最少含有像元數為500 ，最大允許標准差為10 ，即如果一類的標准差大於10 ，則該類被拆分為兩類。類均值間允許最小距離為5，能夠被合並成對的最大對數為5，即當類均值間距離小於5 時，這一類就會被合並，而合並後的成對類的最大值為5（圖8.25 ，圖8.26）。
圖8.25 所用數據：dts_178，ISODATA 填圖結果
圖8.26 所用數據：resize_dts_93，ISODATA填圖結果
8.5.2.2 監督分類填圖結果
首先，對照地質圖和PPI圖像選取兩種岩性（Cgammabeta：石炭紀黑雲母花崗岩，Cdelta：石炭紀閃長岩）的感興趣區。
（1）平行六面體
設置最大允許標准差（Max stdev from Mean）為1.4（圖8.27，圖8.28）。
圖8.27 所用數據：dts_178，平行六面體填圖結果
圖8.28 所用數據：resize_dts_93，平行六面體填圖結果
（2）最小距離
設置最大標准差為10 ，最大允許距離誤差為2500 ，則分類過程中由兩者中較小的一個判定像元是否參與分類，若一旦大於任何一個值則該像元不參與分類，歸屬為無類別（圖8.29 ，圖8.30）。
圖8.29 所用數據dts_178，最小距離填圖結果
圖8.30 所用數據resize_dts_93，最小距離填圖結果
（3）光譜角制圖
設置最大允許角度為0.05°，即當像元波譜與終端端元波譜間夾角大於0.05°時，則不參與分類（圖8.31 ，圖8.32）。
圖8.31 所用數據：dts_178，光譜角制圖填圖結果
圖8.32 所用數據：resize_dts_93，光譜角制圖填圖結果
（4）光譜信息散度
設置最大散度閾值為0.002（圖8.33）。
圖8.33 所用數據：dts_178，光譜信息散度填圖結果
（5）二值編碼
設置最小二進制閾值為0.95 ，決定了哪些像元參與分類（圖8.34 ，圖8.35）。
圖8.34 所用數據：dts_178，二進制編碼填圖結果
圖8.35 所用數據：resize_dts_93，二進制編碼填圖結果
（6）最小距離（Hymap數據）
該方法需要設置兩個閾值：Max stdev from Mean；Max Distance Error。通過不同的嘗試，針對兩種閾值的設定如圖8.36所示。
圖8.36 局部填圖效果
（7）光譜信息散度（Hymap數據）
該方法需要設置閾值：Maximun Divergence Threshod。通過不同的嘗試，參數設置和填圖效果如圖8.37所示。
圖8.37 Threshod＝0.005 時的譜信息散度填圖結果
8.5.2.3 岩性信息填圖結果分析
為了分析高光譜星載模擬數據的岩性填圖效果（圖8.38），可對比研究區的岩性信息分布圖與地質圖逐像素進行對比，如果研究區總的點數為N，地質圖中與岩性信息分布圖的結論相同就認為該點取值1，最後統計結果中1的個數n，這樣岩性信息識別精度可定義為
高光譜遙感技術原理及礦產與能源勘查應用
圖8.38 岩性提取結果圖
筆者發現岩性分類結果與地質圖吻合度達到了89%，對出現誤差的主要原因分析如下（圖8.39）：
1）實驗中所用數據是高光譜星載模擬數據，其空間解析度為30m，隨著解析度的降低，單個像元所對應的地面面積將增大，導致每個像元中包含更多的岩性類型，岩性間的影響性也會增加。一方面，由於混合像元的平均效應，目標岩性在像元中的等效豐度會下降，光譜信息減弱。當像元中目標岩性的等效豐度下降到檢出限以下時，岩性將不能被識別，而造成岩性分布區的外圍含量較低的地段和含量較低分布區填繪面積的減少，點狀集群分布區的漏識別，線狀分布區的斷續分布。這種效應相當於檢出限的下降。另一方面，當像元中目標岩性的等效豐度在檢出限以上時，會使面狀岩性分布區的范圍擴大、空洞的充填、相鄰小區的連接，點狀集群分布區的成片，線狀區域的斑點效應。這兩種效應的綜合作用結果，空間解析度的減小則會使強異常區（包括高豐度區和光譜反襯度較高的礦物分布區）范圍的擴大、小區的相連、點群的結合，而使異常更加醒目，但其細節特徵會因此消失；弱異常區（包括低豐度區和光譜反襯度較低的那些岩性分布區）和小異常區面積會縮小或漏檢；線狀異常可能會形成斷續分布的小斑塊而使線狀特徵和其走向變得不清晰。因此，在填圖的細致程度上，解析度越高的數據填圖效果越好，這是造成填圖效果差異的最主要原因之一。
2）輻射校正的精確程度。無論是高光譜成像儀還是傳統的多光譜感測器，它們所記錄的數據都是地面觀測目標的反射或輻射能量的光譜輻射絕對值，與地物目標的光譜反射率或光譜輻亮度值是不一致的。因此，輻射校正和光譜重建是地物識別不可缺少的環節。但校正後數據必然丟失一定的信息，故再進行後續處理也會有一定的誤差。
3）感興趣區的選取是進行監督分類的一個重要環節，但由於一般情況下並不能得到純度足夠高的感興趣區，造成選取的樣區含有多種岩性，從而對後期處理造成不可避免的誤差。
4）在數據一定的情況下，各種分類方法中參數的選取決定了分類效果，但由於無法遍取各種參數進行嘗試，實驗過程中僅是進行有限嘗試後選擇了具有較好的結果的參數，但並不能保證所用參數是最適合的參數。
總之，用高光譜星載模擬數據進行礦物填圖是可行的，它可以在一定程度上對岩性的種類及分布進行識別。
圖8.39 對比地質圖（上圖是下圖黑框中部分）