混合像元分解及其應(yīng)用(ENVI操作文本)

1、一 基于PPI的端元提取借助純凈像元指數(shù)（PPI）和n維可視化工具用于端元波譜收集 。第一步、獲取純凈像元這個(gè)步驟是在MNF變換的結(jié)果上計(jì)算純凈像元指數(shù)（PPI），之后選擇閾值范圍從PPI圖像上獲得感興趣區(qū)，感興趣區(qū)包含的像元就是比較純凈的像元。（1）打開(kāi)高光譜數(shù)據(jù)。 （2）在ENVI主菜單中，選擇Spectral ->MNF Rotation- > Forward MNF -> Estimate Noise Statistics From Dat a。在標(biāo)準(zhǔn)ENVI文件選擇對(duì)話框中，選擇高光譜圖像文件。 （3）打開(kāi)Forward MNF Transform Paramete

2、rs面板，選擇MNF輸出路徑及文件名，單擊OK執(zhí)行MNF變換。 （4）在波段列表中輸出MNF影像及特征曲線值。從圖中可以看出，大約20個(gè)波段以后的MNF的特征值很小 （5）MNF變換后，在ENVI主菜單中，選擇 Spectral-> Pixel Purity Index->FAST New Output Band。在打開(kāi)的Pixel Purity Index Input File對(duì)話框中，選擇MNF變換結(jié)果，單擊Spectral Subset按鈕，選擇前面10個(gè)波段（MNF后面波段基本為噪聲），單擊OK。 （6）在Display窗口中顯示PPI結(jié)果。選擇Overlay->Re

3、gion of Interest，在ROI Tool 面板中，選擇Options->Band Threshold to ROI，選擇PPI圖像作為輸入波段，單擊OK，打開(kāi)Band Threshold to ROI面板（圖14.19）。 Min Thresh Value：10，Max Thresh Value：空（PPI圖像最大值），其他默認(rèn)設(shè)置，單擊OK計(jì)算感興趣區(qū)，得到的感興趣區(qū)顯示在Display窗口中。 第二步、構(gòu)建n維可視化窗口（1）在ENVI主菜單中，選擇Spectral ->n-Dimensional Visualizer，在n-D Visualizer Input F

4、ile對(duì)話框中選擇MNF變換結(jié)果，單擊OK。（2）在n-D Controls面板中，選擇1、2、3、4、5波段，構(gòu)建5維的散點(diǎn)圖。第三步：選擇端元波譜（1）在n-D Controls面板中，設(shè)置適當(dāng)?shù)乃俣龋⊿peed），單擊Start按鈕，在n-D Visualizer窗口中的點(diǎn)云隨機(jī)旋轉(zhuǎn)，當(dāng)在n-D Visualizer窗口中的點(diǎn)云有部分聚集在一塊時(shí)，單擊Stop按鈕。（2）在n-D Visualizer窗口中，用鼠標(biāo)左鍵勾畫“白點(diǎn)”集中區(qū)域，選擇的點(diǎn)被標(biāo)示顏色。（3）在n-D Controls面板中，選擇Class->Items 1:20->White（用于刪除點(diǎn)），單擊Sta

5、rt按鈕，當(dāng)看到有部分選擇的點(diǎn)云分散時(shí)候，單擊Stop按鈕，在n-D Visualizer窗口中選擇分散的點(diǎn)，自動(dòng)會(huì)將選擇的點(diǎn)刪除。借助<-，->，New按鈕可以一幀幀從不同視角瀏覽以輔助刪除分散點(diǎn)。（4）在n-D Visualizer窗口中，單擊右鍵選擇New Class快捷菜單，重復(fù)（1）（3）選擇其他“白點(diǎn)”集中區(qū)域。 第四步、輸出端元波譜（1）在n-D Controls面板中，選擇Options->Mean All，在Input File Associated with n-D Scatter Plot對(duì)話框中選擇原圖像，單擊OK。（2）獲取的平均波譜曲線繪制在n_D

6、 Mean繪圖窗口中。（3） 識(shí)別每條波譜曲線對(duì)應(yīng)的地物類型。（4） 在n_D Mean繪圖窗口中，選擇File->Save Plot As->Spectral Library（或者ASCII），將端元波譜保存為波譜庫(kù)文件或者文本文件。二 豐度解混（以線性光譜解混為例）（1）Spectral ->Mapping Method-> spectral angle mapper 打開(kāi)要分類的圖像（2）在端元選擇器中打開(kāi)光譜庫(kù) Import->form Spectral Library File （3）選擇所有端元 （4）Apply后，填寫保存路徑即可三 利用端元進(jìn)行分類（以光譜角填圖（spectral angle mapper）為例）（1）Spectral ->Mapping Method->Liner spectral unmixing 打開(kāi)