多语言展示
当前在线:261今日阅读:113今日分享:31

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

R语言实现栅格空间Mann-Kendall分析,计算逐像元Sen's slope,并实现MK显著性检验(mk test)。
工具/原料

电脑,R或者Rstudio

方法/步骤
1

安装需要的packages,并library他们

2

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

3

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

4

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

5

输出栅格结果。

注意事项

有问题留言,没问题点赞。

推荐信息