如何为克里格预测创建网格,如'meuse.grid'对象
如何创建一个网格对象,然后能够在克里格预测模型中使用它。在R中的每一个克里金教程中,每个人都使用如何为克里格预测创建网格,如'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
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来预测。在您的示例中,您没有执行任何克里格。