R语言实现M-K空间趋势分析并检验其显著性

1、安装需要的packages，并library他们

2、将需要分析的栅格按时间顺序叠加到一起。必须是按照时间顺序，叠加之前查看数据名称在R环境下是否与时间顺序一致，因为有的人命名不规则与windows下的顺序可能不一致。下面的例子是以39年的降水数据为例，先读第一年的栅格，后边的年份利用循环叠加的第一层上。

3、定义创建返回多年均值、slope和pValue的函数。

4、利用calc函数利用fun函数对叠加好的栅格逐像元运算。

5、输出栅格结果。