Question

我陷入困境，无法在任何地方找到解决方案。我想在我的数据之间创建一个简单相关性地图，例如this one ，这基本上是来自特定位置的每日降雨量值的列表/ excel电子表格2010年1月的和netcdf数据集，其中包含同期西非的每日降雨量值。 netcdf数据集TRMM可用here，可以通过点击蓝色链接下载＆＃34; netCDF＆＃34;来下载。我的spreedsheet数据来自加纳阿克拉（纬度：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，但我不知道如何处理来自加纳的数据以计算两个数据集之间的相关性。

Answer 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)

R：在电子表格和netcdf数据集之间创建相关性映射

1 个答案: