【发布时间】:2021-02-23 13:54:14
【问题描述】:
我有一个包含 3 列(Col 1:Long,col 2:Lat,col 3:Altitude)和 1095 行的 data.frame 列表。我想将此 data.frame 更改为 RasterLayer 公式。接下来我该怎么办?任何有帮助的都将不胜感激!
这是一个中国的网格单元(1 * 1度,经纬度),有高程值,感谢您的时间和帮助!
陶亮
【问题讨论】:
我有一个包含 3 列(Col 1:Long,col 2:Lat,col 3:Altitude)和 1095 行的 data.frame 列表。我想将此 data.frame 更改为 RasterLayer 公式。接下来我该怎么办?任何有帮助的都将不胜感激!
这是一个中国的网格单元(1 * 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")
【讨论】: