一、MNF变换
4、点击ok,弹出对话框如下
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、在弹出的对话框中,进行如下设置
点击是。
四、
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文件
然后,进行如下操作:
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即可。
至此,端元提取结束。