我有一个sf对象，其中包含对多个单位的变量的年度观测值。我想对此进行栅格化处理，以便从1946-2005年（观察的起始年不同）每年有一个RasterLayer。数据如下：
> head(data)
Simple feature collection with 6 features and 4 fields
geometry type: POINT
dimension: XY
bbox: xmin: 0.5 ymin: 26.63 xmax: 0.5 ymax: 26.63
epsg (SRID): 4326
proj4string: +proj=longlat +datum=WGS84 +no_defs
# A tibble: 6 x 5
id year var1 geometry
<dbl> <int> <dbl> <POINT [°]>
1 615 1957 1 (0.5 26.63)
2 615 1958 1 (0.5 26.63)
3 615 1959 1 (0.5 26.63)
4 615 1960 1 (0.5 26.63)
5 615 1961 1 (0.5 26.63)
6 615 1962 1 (0.5 26.63)
每个RasterLayer应该包含所有几何点不存在var1的信息。
我尝试使用以下代码，但这将创建两层的RasterBrick而不是我想要的59层的RasterStack。
r <- raster::rasterize(data, raster, field = c("var1", "year"), fun = "first")