R：创建电子表格和 netcdf 数据集之间的关联图答案

【问题标题】：R: Create correlation map between spreadsheet and netcdf datasetR：创建电子表格和 netcdf 数据集之间的关联图
【发布时间】：2023-04-03 12:35:02
【问题描述】：

我被卡住了，无法在任何地方找到解决方案。我想创建一个简单的相关图，例如this one在我的数据之间，它基本上是一个来自某个位置的每日降雨量值的列表/ Excel 电子表格，用于2010 年 1 月的月份和 netcdf 数据集，其中包含同一时期西非的每日降雨量值。 netcdf 数据集 TRMM 可用here 并可以通过单击蓝色链接“netCDF”下载。我的电子表格数据来自加纳的阿克拉（纬度：5.6N，经度：-0.22W），如下所示：

- Date, Rainfall
- 2010-01-01, 5.5
- 2010-01-02, 3
- 2010-01-03, 0
- 2010-01-04, 7
- 2010-01-05, NA
- 2010-01-06, 0
- ...
- 2010-01-31, 4.5

我将 netcdf 文件加载为raster::brick，但我不确定如何处理来自加纳的数据以计算两个数据集之间的相关性。

【问题讨论】：

标签： r correlation netcdf

【解决方案1】：

解决此类问题的第一步是用简单的数据创建一个小示例，以便您了解原理。例如：

# example data 
library(raster)
set.seed(424)
b <- brick(nrow=10, ncol=10, nl=100)
values(b) <- matrix(runif(100*100, max=100), 100, 100)   
rain <- runif(100, max=100)

# function to compute correlation between rain and RasterBrick values
fun <- function(x) cor(as.vector(x), rain)

# use the function with calc
r <- calc(b, fun)
plot(r)

【讨论】：