这可能是一个愚蠢的问题,但到目前为止我找不到答案。 我有一个这样的数据框:
df <- data.frame(lat = (-10: 10),
lon = (-10: 10),
A = runif(21, 1 ,10),
B = runif(21, 20 ,30))
其中lat和lon是地理坐标,A和B是两个连续变量。 我想做的是创建一个5乘5网格并计算网格每个单元格中A和B的平均值。这是我尝试的方法,但没有成功......
library(raster)
library(maptools)
coordinates(df) <- ~ lat + lon
rast <- raster(ncol = 5, nrow = 5)
extent(rast) <- extent(df)
rasterize(df, rast, df$A, fun = mean)