R 空间分析 - 使用 for 循环创建每日克里格地图的问题

Question

我一直在尝试使用 R 通过数据集创建 24 小时 PM2.5 值的每日克里格地图。它应该提供 61 张地图作为输出。通过下面的代码，我能够获得 61 张地图的数据点。但是，它只显示了一张传单格式的地图。我不确定我哪里出错了。非常感谢任何帮助。谢谢。

enddate <- as.Date("2020-10-31")

dates <- seq(enddate , startdate , by = -1)

maps <- list()

for (i in seq_along(dates)) {
  currentd <- subset(lab3eleutm, as.Date(date) == dates[i])
  
`  krigelab <- fitme(pm25 ~ 1 + Matern(1|X+Y),
                    data = lab31utm,
                    fixed=list(nu=0.5))`
`  
  predslab <- predict(krigelab,
                     ` newdata=pred3utm %>%
            `            st_drop_geometry(), 
                      variances=list(linPred=TRUE,predVar=TRUE))
  `
  predele <- tibble(pred3utm,
                    first_max_value_pred = predslab[,1],
                    first_max_value_var = attributes(predslab)$predVar) %>%
    st_as_sf(coords=c("X","Y"),
             crs=st_crs(lab3eleutm)) 
  
  predrast3 <- predele %>%
    st_transform(st_crs(lab3eleutm)) %>%
    st_rasterize(template = lab3gridraster)
  
  finallab3  <- tm_shape(predrast3) +
 `   tm_raster() +
   ` tm_shape(lab3eleutm) +
    tm_dots(col = "pm25") +`
    tm_layout(frame = FALSE)`
  `
  maps[[i]] <- finallab3
I have attached the code above. I can share the images if needed.
`
[enter image description here](https://i.stack.imgur.com/NbnIs.png)
The above is the image generated.