加载中…用ENVI和ARCGIS提取石漠化信息实验
(2019-03-09 15:46:18) 普定县的石漠化信息提取
根据习慧鹏等的文章《西南典型喀斯特石漠化时空演变特征》，2018，生态学报，提取石漠化信息步骤以及相关阈值如下：
表一：信息提取实验用到的文件
| 文件编号 | 文件名称 | 文件说明 | 获取方式 |
| F001 | PUDING.shp | 普定县的行政区划 | 下载 |
| F002 | LT51280422006100BKT00.tar.gz | 经纬度是：（105，26）的TM影像，影像资料来源于www.gscloud.cn | 下载 |
| F003 | ASTGTM_N26E105.img.gz | 经纬度是（105，26）的数字高程数据，来源于www.gscloud.cn | 下载 |
| F004 | PUDING-FILE-SLOPE.tif | 普定县所在的方形区域内的坡度信息 | 见步骤一 |
| F005 | PUDING-FILE-NDVI.tif | 普定县所在方形区域的NDVI信息 | 见步骤二 |
| F006 | PUDING-FILE-WATER.shp | 通过NDRI指数获取的水体信息矢量 | 见步骤三 |
| F007 | PUDING-FILE-NDRI.tif | NDRI修正以后的数据 | 见步骤四 |
| F008 | PUDING-FILE-WATER-OFF.shp | 去掉水体以后的普定行政矢量 | 见步骤五 |
表二：石漠化分类定级标准
一、根据DEM制作
1 用arcgis打开ASTGTM_N26E105，然后导出成为tif图像
2 在envi当中，通过subset工具，按照普定的行政区划进行剪切
3 通过arcgis打开剪切完成的栅格数据，然后通过3D analysis tool
-->Slope生成坡度栅格
3 导出成为tif文件备用：PUDING-FILE-SLOPE-PD.tif
二、生成NDVI
1 通过ENVI打开F002，F001
2 用F001的矢量对F002进行剪切
3 通过envi的NDVI工具箱对剪切完成的数据生成归一化植被指数
4 通过compute statistics工具统计该区域的情况
5 取1%和99%作为阈值，生成修正的NDVI指数
三、生成水体的矢量信息
1 执行(band5 - band3)/(band5 + band3)（这里的要点在于整数运算）
2 执行完成的结果可以把水体信息有效提取，然后通过(b1 eq 0)*1 + (b1 gt 0)*0得出水体信息
3 然后通过栅格转矢量生成evf文件
4 最后通过evf转shape生成文件备用F006
四、生成NDRI信息
1 在剪切的文件基础上，执行（float(b5) -
float(b3)）/(float(b5)+float(b3))生成ndri的数据栅格
2 用F006矢量文件制作排除掩模
3 通过compute statiscs在掩模的基础上进行统计
4 通过1%，99%进行信息矫正
5 生成并导出数据F007
五、去除水体信息以后的普定矢量
1 在arcgis当中打开普定行政区划F001和F006
2 通过arcgis的analysis tool -->Extract
-->Clip工具进行剪切，把F001作为输入项，F006作为剪切项
3 导出成为普定去掉水体矢量F008
六、生成石漠化信息
1 将F004，F005，F007在ENVI当中打开
2 通过ENVI的layer stacking工具将三个文件合并在一起
3 打开F008，生成掩模图层 （mask on）
4 通过subset工具将合并文件进行剪切
5 计算无石漠化信息：( bndri lt 0.2 ) * ( bndvi gt 0.7) *(bslope lt
5)，
问题：
这里基本上没有什么值，回头去观察ndri像元值的分布，发现低于0.2的只占很小的比例，而回头去看文章中的阈值，不知道20是从何而来？已经咨询作者
喜欢
0
赠金笔
前一篇：我的育儿导论
后一篇：关于简七理财中的基金定投的验证