混合像元分解及其應用

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

2、及文件名，單擊OK執(zhí)行 MNF變換。(4)在波段列表中輸出MNF影像及特征曲線值。從圖中可以看出，大約20個波段以后的MNF 的 特征值很小F lie Edi i Opti c-itE Plct_Fnjnct i on HelpMNP rile: cur95effJntoln-gc%-3oln-gc%-310203040 5CEigenvalue Numb&r（5）MNF 變換后，在 ENVI （5）MNF 變換后，在 ENVI 主菜單中，選擇 Spectral- Pixel Purity Index-FAST New Output Bando 在打開的Pixel Purity Index I

3、nput File對話框中，選擇MNF變換結果，單擊Spectral Subset按鈕，選擇前面10個波段(MNF后面波段基本為噪聲)，單擊OK。Inpnit File : C： D de uni ent z aiM Settings A ilri inis ti_ at nrUutjiut File : C： i.Dlic uni erut s :eltl-I S etti ilee lA dm i m e tr at orC:eulc el1D00003DQ00II ill, Illi I.ll III, IlliI5Du 1DDD i 50Q 2000 2500 i-r-, ;-Eera

4、tion n|xUpti ons rlcit_FuxLc11.on HelpC:eulc el1D00003DQ00II ill, Illi I.ll III, IlliI5Du 1DDD i 50Q 2000 2500 i-r-, ;-Eeration n|xUpti ons rlcit_FuxLc11.on Helpuxd &口teration #2800; 31/121 11 I1 1 I11 1 1 I1 1 1 Ir（6）在 Display 窗口中顯示 PPI 結果。選擇 Overlay-Region of Interest,在 ROI Tool 面板中，選 擇 Options-Ba

5、nd Threshold to ROI，選擇 PPI 圖像作為輸入波段，單擊 OK，打開 Band Threshold to ROI 面板（圖 14.19）。Min Thresh Value： 10，Max Thresh Value：空（PPI 圖像最大值），其他默認設置， 單擊OK計算感興趣區(qū)，得到的感興趣區(qū)顯示在Display窗口中。馨 Band Threshold to RO-. -Min Thresh ValueM 邳 Thr esh V alueRClI Ham e Thresh (ppi_2. img)EOI Color Red.典匚第二步、構建n維可視化窗口（1）在 ENVI 主

6、菜單中，選擇 Spectral -n-Dimensional Visualizer，在 n-D Visualizer Input File 對話框中選擇MNF變換結果，單擊OK。（2）在n-D Controls面板中，選擇1、2、3、4、5波段，構建5維的散點圖。第三步：選擇端元波譜（1）在n-D Controls面板中，設置適當的速度（Speed），單擊Start按鈕，在n-D Visualizer 窗口中的點云隨機旋轉，當在n-D Visualizer窗口中的點云有部分聚集在一塊時，單擊Stop按鈕。（2）在n-D Visualizer窗口中，用鼠標左鍵勾畫“白點”集中區(qū)域，選擇的點被標示

7、顏色。（3）在 n-D Controls 面板中，選擇 Class-Items 1:20-White （用于刪除點），單擊 Start 按鈕， 當看到有部分選擇的點云分散時候，單擊Stop按鈕，在n-D Visualizer窗口中選擇分散的點，自動會 將選擇的點刪除。借助，New按鈕可以一幀幀從不同視角瀏覽以輔助刪除分散點。（4）在n-D Visualizer窗口中，單擊右鍵選擇New Class快捷菜單，重復（1）（3）選擇其他白 點”集中區(qū)域。萼 n-D isualizer萼 n-D isualizer第四步、輸出端元波譜（1）在 n-D Controls 面板中，選擇 Options-M

8、ean All,在 Input File Associated with n-D Scatter Plot對話框中選擇原圖像，單擊OK。（2）獲取的平均波譜曲線繪制在n_D Mean繪圖窗口中。（3）識別每條波譜曲線對應的地物類型。（4）在 n_D Mean繪圖窗口中，選擇 File-Save Plot As-Spectral Library（或者 ASCII），將 端元波譜保存為波譜庫文件或者文本文件。豐度解混（以線性光譜解混為例）(1) Spectral -Mapping Method- spectral angle mapper 打開要分類的圖像(2)在端元選擇器中打開光譜庫Import

9、-form Spectral Library FileLibrary Input FileCjaTlCcICjaTlCcISelect Input FillFile Informati ort：1 min class LibTJiTvFile! C : i.D dcuiti Hjit z and. Sett i ngs I Adm i m e trDims： 50 x 6 k 1 B5Q-Ei ze! Bouble Pi-eci sicn 2 400 bytes.File Tyi旭-ElP/I 舊ctral LibrryS e ns or Typ 已：LPriinBiV工lByte Uriie

10、r : Kost (Intel：Pr oj e c t i oil : ffoilsWwLtl言七k : 1.切。8 tci 2. 47J Mi crome+tirsUppTr Ln：tt Corner； 1.)1DE-:rip + i&rL： 1巨y spectr ：=J_ libry file Thu Apr 28 11:03:19 2011 一ITJLihrSVM_c1& 鼬甘鋁備 FFI.litraryMemoryZ4?？?或35晦1) cup95 eE. int（3）選擇所有端元NlUTitifcjr of 1 + HjfTiE EfcjlHC + eil：機NlUTitifcjr

11、of 1 + HjfTiE EfcjlHC + eil：機#2器靶瓠#6 5 5 E H 5 LI 5 5 E E 5 g H3.1M 如 H 叫 II II I 1- il II c c n- c c c D D D D c D - n n n n n n野 Input Spectral LibraryAv：iilabl SpectraSelect All I+ hjtis CLe：t- All I ternsWN5rLg+h Units Microiristt-rE |i Data Mult ipli er OUClOOU（4） Apply后，填寫保存路徑即可三利用端元進行分類（以光譜角填圖