R:从具有经纬度和高度值的数据帧创建RasterLayer对象

Question

我有一个包含3列(第1列:长，第2列:纬度，第3列:高度)和1095行的data.frame列表。我想将这个公式更改为RasterLayer data.frame。我接下来该怎么办？任何有帮助的东西都将非常感谢！！

这是一个具有高程值的中国网格单元(1 *1度，经度和纬度)，感谢您的时间和帮助！

陶亮

Answer 1

您可能可以使用rasterFromXYZ

示例数据

library(raster)
r <- raster(nrow=5, ncol=5, xmn=0, xmx=10, ymn=0, ymx=10, crs="")
set.seed(1)
values(r) <- sample(1:25)
r[r < 15] <- NA
xyz <- rasterToPoints(r)

解决方案

rst <- rasterFromXYZ(xyz)

或使用terra

library(terra)
x <- rast(xyz, type="xyz")