【ENVI基础】基于遥感的草原与沙漠化监测

又是夫也哥哥的投稿，大家有什么问题与疑问可以在文后留言。

---

以下是正文：

植被覆盖度是衡量地表植被状况和指示生态环境变化的一个重要指标，通过对干旱区草原植被覆盖度的定量反演，植被覆盖的变化监测，可以实现草原植被的高频率、大范围、实时的变化监测。

内容：

1、影像数据读取

2、辐射定标、图像配准、大气校正

3、植被反演、植被覆盖变化监测

步骤

一、影像数据读取

二、辐射定标、图像配准、大气校正

三、植被覆盖度反演

1、打开影像，波段选择432，影像中红色区域为植被

2、在Toolbox中搜索NDVI，双击，选择影像

波段默认3、4，选择保存路径，点OK

得到NDVI影像

3、植被覆盖度计算

公式：FC= [(NDVI- NDVIsoil) / (NDVIveg- NDVIsoil)]

实际运算表达式：

（b1gt 0.7）* 1 + (b1 lt 0.）*0 + (b1 ge 0 and b1 le 0.7) * ((b1 – 0.0) / (0.7-0.0))

表达式解释（上期教程中，有人询问，表达式具体含义是什么，这次详细解释一下。）

上期文章：【ENVI基础】如何进行水体提取？

B1 是指NDVI，数值0.7和0详见参考文献（其实就是一个经验值，但是为了教程的严谨性，大家可以阅读与自己实验地区相关的论文，得出数值）

（b1gt 0.7）* 1 指b1大于0.7就乘以1，就是认为NDVI大于0.7的植被覆盖度就是1

(b1 lt 0.）*0指b1小于0就乘以0，就是认为NDVI小于0的植被覆盖度就是0

(b1 ge 0 and b1 le 0.7) * ((b1 – 0.0) / (0.7-0.0)) 指NDVI在0到0.7之间，就采用(b1– 0.0) / (0.7-0.0)计算

将运算表达式输入BandMath，输出影像如下：

4、为了将草原进一步提取显示出来，我们将NDVI大于0.3的归为草原。（为什么不是大于0，因为存在混合像元的问题，裸地上也会出现植被）详见参考文献

相同的方法，换一个不同时期影像，推荐间隔3-5年，对比效果明显

5、使用Band math，对两个影像做差值，

6、右击差值后的影像，选择New raster color slice,得到草原植被变化图

图像分析：黄色为背景，绿色为植被增加区域，红色为植被减少区域。颜色可以自己修改

成果图：

参考文献： 

[1]马琳雅,崔霞,冯琦胜,梁天刚.2001-2011年甘南草地植被覆盖度动态变化分析[J].草业学报,2014,23(04):1-9. 

[2]阿多,赵文吉,宫兆宁,张敏,范云豹.1981—2013华北平原气候时空变化及其对植被覆盖度的影响[J].生态学报,2017,37(02):576-592.

[3]彭文甫,王广杰,周介铭,徐新良,罗怀良,赵景峰,杨存建.基于多时相Landsat5/8影像的岷江汶川-都江堰段植被覆盖动态监测[J].生态学报,2016,36(07):1975-1988.

那今天就到这里结束啦，欢迎留言讨论。文中的图片未经许可不要随便“引用”。

如果可以的话，希望能够转发分享，点个在看并且点个赞，给个赞赏~~也欢迎规范转载~

也希望大家和我多留言互动啊！（据说这样可以增加我的推送在你的订阅号里出现的概率）

文章推荐【ENVI基础】如何进行大气校正？【数据分享】如何获取全国矢量地铁站点与线路数据？【空间句法】20分钟学会空间句法高阶技能sDNA【城事数说】方向熵：上海与重庆的道路方向与混乱程度【数据整理】2019年行政区划调整的数据更新与分享【城事数说】时空圈—大城市之间的比较【城事数说】10年间建设用地如何变化？——全国、长三角、上海三个层面的比较【数据分享】扩张的世界，收缩的城市（分享世界人口密度数据）【数据分享·教程】20分钟做一个中国城市经济联系图【城事数说】中国的城市联系网络【GIS基础】谁都能学会的地形分析全流程（3）【数据分享】使用建筑轮廓数据估算现状建设指标【GIS基础】谁都能学会的地形分析全流程（2）【城事数说】大城市的人们都回哪去了？——人口迁徙视角（1）【GIS基础】谁都能学会的地形分析全流程（1）

END>

如需全文转载文章、投稿或者合作

可添加微信

(回复超慢！！！)

（不要添加我问各种问题，我大概率不会的==）

（入群请一定要备注入群）

（添加后会在晚上非工作时间通过，请稍安勿躁）

 公众号 微博

使用关键词搜索历史文章请点击【阅读原文】