R中corLocal获取Kendall/Pearson相关的P值与Sen斜率栅格

你通过corLocal生成的对象p是一个RasterStack（栅格栈）：当使用method = "kendall"时，它包含三个栅格层（相关系数、P值、Sen斜率）；若使用pearson方法，则包含相关系数和P值两个层。

具体操作步骤： #

• 首先确认栅格栈的层名称，明确各指标对应位置：

names(p)

执行后会输出类似[1] "correlation" "p.value" "sen.slope"的结果。

• 提取并保存单个栅格层：
  • 保存P值栅格文件：

  writeRaster(p[["p.value"]], filename = "kendall_pvalue.tif", format = "GTiff", overwrite = TRUE)
  

  • 保存Sen斜率栅格文件：

  writeRaster(p[["sen.slope"]], filename = "kendall_sen_slope.tif", format = "GTiff", overwrite = TRUE)
  

  若使用pearson方法，仅需提取保存p.value层即可。

补充提示： #

• 可修改filename参数自定义输出文件的路径与名称；
• overwrite = TRUE用于覆盖同名已有文件，可根据需求调整；
• 若需提取相关系数层，使用p[["correlation"]]即可。

内容的提问来源于stack exchange，提问作者Tenzing Ingty