r - 如何为克里金预测创建一个网格，如“meuse.grid”对象

Question

如何创建一个可以用于克里金预测模型的网格对象。在 R 中的每个克里金教程中，每个人都使用该meuse.grid对象但是我如何首先在我的数据上创建这个对象？

我在下面准备了一个例子。在这里，我有我grid想在克里金预测中使用的对象。但是，当我运行它时，克里金模型的预测值是 NAN。通过各种健全性检查，我得出结论，我的grid对象存在问题。

如何将其转换为与对象完全相同类型的meuse.grid对象？

library(raster)
library(tidyverse)
library(sf)
library(stars)
library(gstat)
library(ggmap)

Regions <- getData("GADM", country= "CZE", level = 1)

Regions <- Regions[Regions@data$NAME_1 == 'Zlínský',]


# Create a Grid for the kriging predictions -------------------------------

Regions <-
  Regions %>% st_as_sf()

grid_spacing <- 0.02

polygony <- st_make_grid(Regions, square = T,
                         cellsize = c(grid_spacing, grid_spacing)) %>%
  st_sf()

plot(polygony, col = 'white')
plot(st_geometry(Regions), add = T)

grid = st_intersection(polygony, Regions)
plot(grid)

grid


data("meuse.grid")

meuse.grid

标签： rmapspolygonrasterkriging