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

问题描述

我创建了两个“ 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)

解决方法

暂无找到可以解决该程序问题的有效方法，小编努力寻找整理中！

如果你已经找到好的解决方法，欢迎将解决方案带上本链接一起发送给小编。

小编邮箱:dio#foxmail.com (将#修改为@）