基于arcgis栅格计算器进行连续时间序列NDVI的Slope趋势分析

一、趋势分析法简介

趋势分析法是一种通过对随时间变化的变量进行线性回归分析，从而来预测其变化趋势的方法。其计算方法如下：

式中，Slope为像元回归方程的斜率，NDVI i为第i年的NDVI的平均值，n为研究的时间长度,视自身情况而定。当Slope>0时，表示该像元NDVI为增加趋势；当Slope=0，表示该像元NDVI基本不变；当Slope<0时，表示该像元NDVI为减少趋势。

二、基于ArcGIS进行趋势分析

1、趋势分析前提

趋势分析的前提是要有NDVI年均值图像，可以直接通过arcmap处理得到，过程简单，不再赘述。
步骤：ArcToolBox——>Spatial Analyst Tools——>Local——>Cell Statistics
该像元统计工具还可以做最大值合成，其本质都是对两景或两景以上影像求均值。实际上还是要根据自己的需求进行前期处理。

2、趋势分析

工具：栅格计算器
位置：ArcToolBox——>Spatial Analyst Tools——>Map Algebra——>Raster Calculator
步骤：
（1）加载年均值图像。我的是2016-2020五年时间，所以公式中n取5。

（2）打开栅格计算器，输入趋势分析公式。

由于公式累加较长，故先计算分子左边：

计算分子右边（注：与i有关的累加直接口算即可，如果多的话也可以通过excel等进行计算）

计算Slope，分母结果为常数，直接算出来即可：

所得结果：

（3）重分类，划分等级，出趋势分析图。
重分类工具位置：ArcToolBox——>Spatial Analyst Tools——>Reclass——>Reclassify
（注意：数据管理工具——>栅格——>栅格处理中也有一个重分类工具，不要用这个）

在重分类中设置分的类别数以及断点：

OK后，按照个人喜好调节色带：右击图层——>属性——>符号，结果图如下：

3、出专题图
打开布局，按照需要设置图例，比例尺，指南针等，就可以得到自己所需的NDVI趋势分析图。

（以上用栅格计算器进行趋势分析的思路来自于B站up主TR繁鑫~
个人认为这种方法只适合于少量年份的分析处理，如果多的话单纯用ArcGIS是很费时间的，所以还是要学好matlab/IDL/python编程，才能事半功倍！）