ENVI端元提取（MNF→PPI→n-D Visualizer）_文玲de学习时间

http://blog.sina.com.cn/u/2107714272

首页博文目录关于我

个人资料

微博

加好友发纸条

写留言加关注

博客等级：
博客积分：

博客访问：
关注人气：
获赠金笔：0支
赠出金笔：0支
荣誉徽章：

正文字体大小：大中小

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

(2014-02-21 09:02:47)

分类：遥感分类

一、MNF变换

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

http://img2.ph.126.net/JLdtw4nH5XcaJmxAKslG7g==/6598169278285256747.jpgVisualizer）" />

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

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

http://img0.ph.126.net/eiIyCgD0CvKraU-gxdTeaw==/6597709682425051482.jpgVisualizer）" />

点击OK，输出mnf文件。

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

http://img1.ph.126.net/wuNi3BMd_KBcKmn_NVKacQ==/6597278673866466280.jpgVisualizer）" />

二、PPI计算

1、如图：

http://img2.ph.126.net/Z-P48gfPmYpHpvIYNAEKnw==/6597875708680859876.jpgVisualizer）" />

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

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

http://img1.ph.126.net/duLHFI-boSCU8QKlKCNjlw==/6597782250192489889.jpgVisualizer）" />

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

三、

1、进行如下操作：

http://img1.ph.126.net/3ExzDx6WGzdufw-fQUdRaw==/6597963669611084347.jpgVisualizer）" />

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

http://img2.ph.126.net/CbpoTF1AAqtFgWIWXXwlLQ==/6598293523099203644.jpgVisualizer）" />

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

http://img2.ph.126.net/TMr676Up13j4UBodzqWjgQ==/6597722876564588456.jpgVisualizer）" />

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

http://img1.ph.126.net/CcASEPj01opfvGSGIObTFA==/6598149487075951474.jpgVisualizer）" />

点击是。

四、

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

http://img0.ph.126.net/T6Xr6BDsYRuHvF7Ok4kSqg==/6598179173889907128.jpgVisualizer）" />

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

http://img0.ph.126.net/79SB_R8UzkVpm8ulIUCDOg==/6598184671448045439.jpgVisualizer）" />

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

http://img2.ph.126.net/NxifJWniRs6niWKq2x3_Cw==/6598114302703860221.jpgVisualizer）" />

4、弹出两个窗口如下：

http://img1.ph.126.net/-j9WS1HnTsiN2dHBrvzF6w==/6597859216006437865.jpgVisualizer）" />

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

http://img2.ph.126.net/Tlk3qu0uVokD9qzjOen-UA==/6597783349704118022.jpgVisualizer）" />

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

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

http://img0.ph.126.net/hblfM8CXSb18p2Cr5Phghg==/6597540357634375857.jpgVisualizer）" />

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

http://img0.ph.126.net/i_LWL7PXwYmJNWc6ZcB_lg==/6598196766075948246.jpgVisualizer）" />

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

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

http://img1.ph.126.net/YvrY9-MxLtWIQWCTr8w_VA==/6597711881448307624.jpgVisualizer）" />

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

http://img2.ph.126.net/-yoaeYIfVzjdZpk-Wxt9bA==/6597926286215735291.jpgVisualizer）" />

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

http://img0.ph.126.net/knpFNDid3p10amM_2JJkaQ==/6597807538959937155.jpgVisualizer）" />

然后，进行如下操作：

http://img0.ph.126.net/xFwpqizsPertLTFBLtsRLw==/6597951574983176082.jpgVisualizer）" />

http://img0.ph.126.net/CvKagtkEBJfjx54yh8DJsw==/6597979062773865042.jpgVisualizer）" />

选择usgs_min打开即可。

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

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

http://img1.ph.126.net/E07lmwX4vifla76fIAxGPg==/6598231950448035253.jpgVisualizer）" />

六、进行波段分析

1、

http://img1.ph.126.net/8jwJGpGL6n_LjsId8gACaQ==/6597926286215735383.jpgVisualizer）" />

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

http://img0.ph.126.net/_TphcFyXhdqWmqogp8XKvw==/6597707483401796745.jpgVisualizer）" />

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

http://img1.ph.126.net/JVpM_ESCYNOtifxIX0oSeg==/6597330350912974705.jpgVisualizer）" />

4、弹出Spectral Analyst对话框

http://img1.ph.126.net/a2X_3XGogf6cO6ZN8Ko2KA==/6598227552401524279.jpgVisualizer）" />

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

http://img1.ph.126.net/pGOjsE9ULRL24d3WE4wr0w==/6597692090239012387.jpgVisualizer）" />

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

http://img1.ph.126.net/tOPzVVr9zvcanOdwUN55ug==/6597854817959927248.jpgVisualizer）" />

可知，步骤四中的分步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...弹出的如下对话框中

http://img0.ph.126.net/9fbpzbqH_NpUm4bumDA8qg==/6598090113448052055.jpgVisualizer）" />

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

至此，端元提取结束。

http://wlb321xy.blog.163.com/blog/#m=0&t=1&c=fks_084070082087083065087084085095082081087068083081081

阅读┊ 收藏 ┊ 喜欢 ▼ ┊打印┊举报/Report

前一篇：[转载]WGS84坐标与北京54坐标转换（转）

后一篇：arcgis中计算多边形面积

新浪BLOG意见反馈留言板　欢迎批评指正