ENVI扩展工具:植被覆盖指数计算器

0
分享 2016-06-19
更新记录:
20151224:修复BUG一枚,当矢量属性与统计结果名称有冲突时报错的问题,错误信息如下:

新版下载地址:http://pan.baidu.com/s/1pKoexBH


在上一个博文中(http://blog.sina.com.cn/s/blog_764b1e9d0102wet7.html)介绍了如何利用ENVI进行植被覆盖指数的计算工作,建议先浏览此博文了解具体操作流程。
从预处理到统计结果,中间流程较多,且重复操作较多,为了提高工作效率,特地编写专门实现植被覆盖指数的计算工具。
注:此工具与上述博文中的流程完全一致,如有变动,需自行修改源代码。

使用方法
下载地址:请下载新版补丁,链接在文章开头。
压缩包内包含补丁文件(.sav)和IDL源代码(.pro)。

首先将vegetation_cover_index_calculator.sav补丁拷贝到如下安装路径(以ENVI 5.2 SP1为例):
  • ENVI 5.2 SP1:C:\Program Files\Exelis\ENVI52\extensions\
  • ENVI Classic:C:\Program Files\Exelis\ENVI52\classic\save_add\


注1:由于在统计植被覆盖指数时,调用ENVI 5.2 SP1版本提供的新API,如果想基于Shapefile统计得到植被覆盖指数,需要ENVI最低版本为ENVI 5.2 SP1。
注2:低于5.2 SP1版本的ENVI,可以进行植被覆盖栅格图像的计算,然后在ArcMap中进行分区统计。
注3:本次提供的sav文件适用于ENVI 5.2及更高版本,如果需要在较低版本中使用,请自行将工程导入IDL工作空间,右键构建工程(或构建项目)生成对应版本的sav即可。

操作步骤
1) 通过如下菜单或工具启动补丁:
Ÿ ENVI Classic:Spectral > Vegetation Analysis > Vegetation Cover Index Calculator
Ÿ ENVI 5.x:Toolbox/Extensions/Vegetation Cover Index Calculator
2) 在弹出的Vegetation Cover Index Calculator Parameters面板中设置如下参数:
Ÿ Select HDF Files:通过多选输入20景MOD13/MYD13数据,可多次添加;
Ÿ Select HDF Directory:通过输入路径自动搜索HDF文件;
Ÿ Full Path:切换列表中显示完整文件路径;
Ÿ File Name:切换列表中只显示文件名;
Ÿ Delete:删除列表中选中的文件,可多选;
Ÿ Clear All:清空列表;
Ÿ Number of Files:显示目前已加入列表的文件个数;
Ÿ Resampling:重投影到Albers投影时的重采样方法;
Ÿ Warp Method:重投影到Albers投影时的几何校正方法;
Ÿ Grid Spacing:采样间隔;
Ÿ Select Shapefile:选择矢量文件,用于分区统计植被覆盖指数;
Ÿ Zone Field:选择矢量文件中用于统计的分区字段;
Ÿ Enter Output Raster Filename:单击Choose按钮,设置输出栅格文件名。
图:Vegetation Cover Index Calculator Parameters面板
3) 单击OK按钮进行处理,会进行如下几个处理步骤:
Ÿ 依次打开20景HDF文件;
Ÿ 进行波段运算,计算植被覆盖指数栅格结果;
Ÿ 利用Shapefile进行分区统计,得到文件名为“输出栅格文件名+_stats.csv”的表格统计结果,并自动在Excel中打开,如下图所示,MEAN一列即为每个县市的植被覆盖指数。
注:分区统计时,需要ENVI 5.2 SP1版本支持。


图:植被覆盖指数统计结果

文章来源:http://blog.sina.com.cn/s/blog_764b1e9d0102wfc8.html

0 个评论

要回复文章请先登录注册