R terra -将数组分配给不包含所有单元格值的SpatRaster对象的最有效方法

Question

假设我们有一个数组modeldata(数据来自地面模型)，其维度为:

> dim(modeldata)
[1] 67420   518

第一个维度包括网格单元，第二个维度包括来自1500:2017的时间序列
第一个维度的异常长度是由于单独存在地面单元以节省空间。

在raster包中，我以以下方式处理它:

> coords
[,1]   [,2]
[1,] -179.75 -16.25
[2,] -179.75  65.25
[3,] -179.75  65.75
[4,] -179.75  66.25
[...,]  ...     ...
[67420,] 179.75  71.25
> wgs84 <- sp::CRS("+proj=longlat +ellps=WGS84 +datum=WGS84 +no_defs")
> modeldata_spdf <- sp::SpatialPixelsDataFrame(coords,
data = data.frame(modeldata),
proj4string = wgs84)
> modeldata_brick <- raster::brick(modeldata_spdf)

请不要以这种方式评判我，
我对使用terra包的类似(性能)方法更感兴趣。
另一种方法也是很好的，将使用SpatRaster蒙版层代替坐标。

谢谢:-)

Answer 1

下面的方法适用于在已知的规则栅格上计算位置子集(感兴趣的单元格)的值的情况。其他情况参见values<-、rast(, xyz=TRUE)和rasterize。

示例数据

library(terra)
r <- rast(res=10)
set.seed(0)
cells <- sample(ncell(r), 5)
xy <- xyFromCell(r, cells)
xy
#        x   y
#[1,] -165 -25
#[2,]   25  55
#[3,] -135 -55
#[4,] -155 -45
#[5,]  -75   5

在这种情况下，我们有5个位置，下面的模型数据md有2个时间步(层)，每个时间步(层)有5个值

md <- cbind(A=1:5, B=5:1)
# compute cell numbers (not needed here as we have them already)
cells <- cellFromXY(r, xy)
# create empty raster
x <- rast(r, nlyr=2)
# assign values
x[cells] <- md

在您的情况下，保留单元格编号而不是坐标可能更有效。

以上应该也适用于raster，只需稍加修改。