混合像元分解及其应用

1、一基于PPI的端元提取借助纯净像元指数（PPI）和n维可视化工具用于端元波谱收集第一步、获取纯净像元这个步骤是在MNF变换的结果上计算纯净像元指数（PPI），之后选择阈值范围从PPI图像上获得感兴 趣区，感兴趣区包含的像元就是比较纯净的像元。（1）打开高光谱数据。（2）在 ENVI 主菜单中，选择 Spectral -MNF Rotation- Forward MNF - Estimate Noise Statistics From Dat a。在标准ENVI文件选择对话框中，选择高光谱图像文件。(3)打开Forward MNF Transform Parameters面板，选择MNF输出路径

2、及文件名，单击OK执行 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计算感兴趣区，得到的感兴趣区显示在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窗口中，用鼠标左键勾画“白点”集中区域，选择的点被标示

7、颜色。（3）在 n-D Controls 面板中，选择 Class-Items 1:20-White （用于删除点），单击 Start 按钮， 当看到有部分选择的点云分散时候，单击Stop按钮，在n-D Visualizer窗口中选择分散的点，自动会 将选择的点删除。借助，New按钮可以一帧帧从不同视角浏览以辅助删除分散点。（4）在n-D Visualizer窗口中，单击右键选择New Class快捷菜单，重复（1）（3）选择其他白 点”集中区域。萼 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后，填写保存路径即可三利用端元进行分类（以光谱角填图