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

问题描述

如何创建一个可用于克里金预测模型的网格对象。在 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

解决方法

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

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

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