混合像元分解

注：由于计算量较大，所以将图像分成三块进行处理，另外此步骤在ENVI Classic里面完成

一：纯像元提取（PPI）

1、MNF变换

选择Spectral→MNF Rotation→Forward MNF→Estimate Noise Statistics from Data

在弹出的窗口中选择13-4-26a.img文件，点击ok

在弹出的窗口中设定保存路径，文件命名为13-4-26amnf，点击ok

2、PPI计算

选择Spectral→Pixel Purity Index→[FAST]New Output Band

在弹出的窗口中选择13-4-26amnf文件，点击ok

在弹出的窗口中设置输出文件，13-4-26appi，点击ok，这一步要等很久。。。

将获取的纯像元变为ROI：Basic Tools→Region of Interest→Band Threshold to ROI，选择13-4-26appi文件，点击ok

弹出的窗口中：

3、构建n维窗口，选择端元波谱

选择Spectral→n-Dimensional Visualizer→Visualize with New Data

在弹出窗口中选择13-4-26amnf文件，点击ok，选中1-5波段

（1）在n-D Controls面板中，设置适当的速度（Speed），单击Start按钮，在n-D Visualizer 窗口中的点云随机旋转，当在n-D Visualizer窗口中的点云有部分聚集在一块时，单击Stop 按钮。借助<-，->，New按钮可以一帧帧从不同视角浏览以辅助删除分散点。

（2）在n-D Controls面板中，选择Class->Items 1:20->red，可选择标志颜色

（3）在n-D Visualizer窗口中，用鼠标左键勾画“白点”集中区域，右键结束

4、输出端元波谱

在n-D Controls面板中，选择Options->Mean All，在Input File Associated with n-D Scatter Plot 对话框中选择原图像13-4-26a，单击OK

（2）获取的平均波谱曲线绘制在n_D Mean绘图窗口中。

（3）识别每条波谱曲线对应的地物类型。

（4）在n_D Mean绘图窗口中，选择File->Save Plot As->Spectral Library（或者ASCII），将端元波谱保存为波谱库文件或者文本文件。命名为13-4-26asp

二：丰度解混

Spectral ->Mapping Method-> spectral angle mapper 打开要分类的图像13-4-26a.img,点击ok 在弹出的窗口中，选择Import->form Spectral Library File

选择13-4-26asp文件

选择所有端元，点击apply，命名为13-4-26ax

三：利用端元分类

Spectral ->Mapping Method->Liner spectral unmixing 打开要分类的图像13-4-26a.img

在弹出的窗口中，选择Import->form Spectral Library File

选择13-4-26asp文件

选择所有端元

点击apply，命名为13-4-26aclass，注意将第一行选项改为yes，这样加权平均值才为1

AOI提取端元波谱

1、打开原图像（能目视解译的图像就行）右键原图像New Region of Interest新建一个roi 图层，修改roi name和color，在原图上画多边形勾出相应地类。点击添加地类

2、在tool窗口中optionS-send rois to n-D Visualizer ,点击ok

3、等一会儿，出现n-D Visualizer窗口，选中1-5，option-Mean all，选择要分类的图像

4、出现波普曲线，点击export-spectral library ,保存为波谱文件。提取端元波谱成功