微信公众号搜"智元新知"关注
微信扫一扫可直接关注哦!

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

如何解决如何为克里金法预测创建网格,如“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 举报,一经查实,本站将立刻删除。