计算栅格时序数据的斜率并提取特定斜率范围的栅格

• 前言
  本文所说的时序数据，本质上是多波段数据，理论上可以有无数个波段，但是需要合成为一个栅格文件。这个步骤可以在ENVI中进行。

  本文以DMSP数据为例，计算某地多年DMSP数据的斜率，得到一个值为斜率的栅格数据。

1 数据和软件

本文使用多年DMSP数据，需要掩膜提取成统一范围，建议读者设定投影坐标。
使用ENVI和QGIS，用ENVI计算斜率，用QGIS提取特定的斜率范围数据。

2 计算时序栅格斜率

2.1 在ENVI中安装斜率计算器插件

在ENVI的app store中搜索并安装时序数据变化趋势（斜率计算器）

安装后在extensions中可以看到Slope K Calculator工具即视为安装成功

2.2 波段合成

将掩膜提取后的所有栅格数据进行波段合成，即ENVI中的Build Layer Stack，在波段合成时需要注意，各个年份的DMSP的放置顺序必须沿着时间线，这样计算出的斜率才是有意义的。

• 波段合成后的数据：
  

2.3 计算斜率

输入刚刚波段合成后的数据，设定输出数据位置与名称

• 得到计算结果，这个计算结果是单波段数据，只有Slope K波段，栅格值就是这个像元位置时序数据的斜率

2.4 导出斜率栅格为tif

点击SlopeK数据，选择左上角的File—Save As—Save As TIFF

输出时候一定选择TIFF格式

3 提取特定斜率范围的栅格

在QGIS中选择刚刚添加进来的斜率栅格，右键导出—另存为


在栅格图层另存为窗口，勾选无数据值，填写非目标值域，比如我想保留-0.1到0.1的斜率像元，那么我就把非该值域范围的值填入其中

---

感谢阅读