经纬度栅格

Question

给定 raster对象 r ，如何创建具有相同范围和分辨率的新栅格，单元格值等于 r 中相应单元格的纬度(或经度) ?

例如，r可能看起来像:

r <- raster(matrix(runif(100), ncol=10))

Answer 1

您可以使用 init功能

library(raster)
r <- raster(nrow=10, ncol=10)
lon <- init(r, 'x')
lat <- init(r, 'y')