最近看到一些博主分享关于 Sen+MK 检验的代码，对于新手来说可能有点复杂。我们编写了一段 MATLAB 代码，能够一次性解决这些问题，简化操作流程。我们还准备了几个关于趋势检验的空间分布图，供大家参考。

一、Sen's Slope和Mann-Kendall检验

在研究中，我们常结合 Theil-Sen 斜率估计（Sen's Slope）和 Mann-Kendall (MK) 检验对栅格数据进行趋势分析。以下是这些方法的具体公式：

Liu, Z.; Wang, H.; Li, N.; Zhu, J.; Pan, Z.; Qin, F. Spatial and Temporal Characteristics and Driving Forces of Vegetation Changes in the Huaihe River Basin from 2003 to 2018. Sustainability 2020, 12, 2198. https://doi.org/10.3390/su12062198

二、MATLAB代码

只需要更改栅格数据的输入输出路径、年份和栅格的有效值。

 

 

 

三、ArcGIS制图

ABC是MATLAB输出的栅格，我们需要把D做出来。以下是具体的步骤：

1. 将得到三个栅格数据，图中分别是MK检验、SEN趋势和MK显著区域。

2. 对significant_mk栅格转面。

3. 更改矢量的样式。