R和SF中两个栅格之间的空间相关性

Question

我创建了两个“ sf”克里金法对象（点向量），一个用于温度，另一个用于农业产量。为了制作网格并进行点插值，我保留在“ sf”包中。

我现在想在这两个变量之间创建空间局部相关性“栅格”，如本网页https://statnmap.com/2018-01-27-spatial-correlation-between-rasters/ https://statnmap.com/2018-01-27-spatial-correlation-between-rasters/所示。但是，在该示例中，他们使用了“光栅”包和“焦点”功能。我想知道是否可以在“ sf”内实现此目的，即无需更改班级？如果没有，将这些对象转换为栅格类的最佳方法是什么？

以下是我的克里金代码摘录，以供参考：

library(sf)
sf_data <- st_as_sf(x = data, coords = c("longitude", "latitude"), crs = 4326)
library(gstat)
vgm_utci <- variogram(UTCI~1, sf_data)
utci_fit <- fit.variogram(vgm_utci, vgm("Gau"), fit.kappa = TRUE)
plot(vgm_utci, utci_fit)
istria <- read_sf(“./Istria_Boundary.shp")
istria <- istria$geometry
istria.grid <- istria %>% 
 st_make_grid(cellsize = 0.05, what = "centers") %>% 
 st_intersection(istria) 
library(ggplot2)
ggplot() + geom_sf(data = istria) + geom_sf(data = istria.grid)
library(stars)
utci_krig <- krige(formula = sf_data$UTCI ~ 1, locations = sf_data, 
                  newdata = istria.grid, model = utci_fit)