如何为克里格预测创建网格,如'meuse.grid'对象

如何为克里格预测创建网格,如'meuse.grid'对象,r,maps,polygon,raster,kriging,R,Maps,Polygon,Raster,Kriging,如何创建一个网格对象,然后能够在克里格预测模型中使用它。在R中的每一个克里金教程中,每个人都使用meuse.grid对象,但是我首先如何在我的数据上创建这个对象呢 我在下面准备了一个例子。这里,我有我的网格对象,我想在克里格预测中使用它。然而,当我运行它时,克里格模型的预测值是NAN。通过各种健全性检查,我得出结论,我的grid对象存在问题 如何将其转换为与meuse.grid对象完全相同的对象类型 library(raster) library(tidyverse) library(sf) l

如何创建一个网格对象,然后能够在克里格预测模型中使用它。在R中的每一个克里金教程中,每个人都使用
meuse.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

库(光栅)
图书馆(tidyverse)
图书馆(sf)
图书馆(星级)
图书馆(gstat)
图书馆(ggmap)

我认为,
meuse.grid
只是一个带有x和y坐标的data.frame来预测。在您的示例中,您没有执行任何克里格。