（转）ENVI端元提取（MNF→PPI→n-D Visualizer）

仍以cai_jian1为例说明

    一、MNF变换

    1、选择Spectral→MNF Rotation→Forward MNF→Estimate Noise Statistics from Data,如图：

     2、在弹出的对话框MNF Transform Input File对话框中选择cai_jian1文件

     3、在弹出的对话框，如下图示，添加需要保存的mnf输出的路径，此处，保存文件名mnf

 点击OK，输出mnf文件。

4、并弹出MNF Eigenvalues示意图，如下：

 

二、PPI计算

1、如图：

2、在弹出的Fast Pixel Purity Index Input Data File中，选择步骤一中生成的mnf文件

3、在弹出的对话框Fast Pixel Purity Index Parameters中,选择保存路径，其他设置保存默认即可。

 4、点击OK之后进行迭代计算10000次，最终生成ppi图像。

三、

1、进行如下操作：

2、在弹出的Band Threshold to ROI Input Band对话框如下中，选择ppi图像，点击OK

 3、在弹出的对话框中,进行如下设置

 

4、点击ok，弹出对话框如下

点击是。

四、

1、选择Spectral→n-Dimensional Visualizer，如下图

2、在弹出的n-D Visualizer Input File中选择mnf文件

 

3点击Spectral Subset，选择mnf变换的前20个波段

4、弹出两个窗口如下：

 5、在n-D Controls窗口中的n-D Selected Bands中选择前7个波段，在左侧的n-D Visualizer中显示如下图示

6、通过Start按钮，旋转左侧n-D Visualizer中的点群，通过在右侧n-D Controls中的Class菜单下的Red、Green、Blue等颜色，来圈定左侧n-D Visualizer中的点密集的区域。

7、最终圈定的类别如下示：

 8、选择Options→Mean All.弹出对话框如下：

 

五、由于本文使用的图像与USGS波谱库中的sli中的光谱维分辨率不同，需要根据本文中的cai_jian1的图像对USGS波谱库中的sli波谱库文件进行重采样

1、选择Sprctral→Spectral Libraries→Spectral Libery Resampling

2、在弹出的Spectral Resampling Input File对话框中,选择USGS的usgs_min.sli文件

 

 注意： 如果该文件在Select Input File中没有出现，进行如下操作添加：

 然后，进行如下操作：

 

 

 选择usgs_min打开即可。

3、对步骤2中的对话框选择了usgs_min.sli波谱库后，点击OK

弹出如下对话框，并设置输出路径，点击OK，即可得到usgs_min.sli 重采样后的数据usgs_min_resampling

六、进行波段分析

1、

                                                                            

  2、在弹出的Spectral Analyst Input Spectral Library对话框中选择重采样的波谱库，如下：

 3、点击OK。在弹出的Edit Identify Methods Weighting对话框中选择默认方式,点击OK

4、弹出Spectral Analyst对话框

  5、 点击Apply按钮，弹出如下对话框

 6、选择n_D Class #4点击OK，即可得到结果如下：

可知，步骤四中的分步8中得到的光谱曲线图 n_D Means:cai_jian中的n_D Class #4 为goethit3(查询Google在线翻译知，此为针铁矿)

7、重复以上步骤5-6，可判定此影像中的其他三种矿物闪锌矿(Sphaler)、泻利盐(Epsonit)、赤铁矿(Hematit).

8、可以在步骤四中的分步8中得到的 光谱曲线图 n_D Means:cai_jian中，选择Edit→Data Parameters...弹出的如下对话框中

 通过点击n_D Class #1来修改其名字,点击Apply即可。

至此，端元提取结束。