经纬度栅格

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

例如，r 可能如下所示:

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

最佳答案

如果您的问题是关于创建一个与另一个光栅对象具有相同范围和分辨率的新光栅对象，您可以使用命令 templatetemplate 是用于设置范围的 Raster* 或 Extent 对象(如果是 Raster* 对象，则为 CRS)。如果不是 NULL，则忽略参数 xmn、xmx、ymn、ymx 和 crs(除非模板是 Extent 对象)

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

关于经纬度栅格，我们在Stack Overflow上找到一个类似的问题：https://stackoverflow.com/questions/35580730/