问题描述
我有一个光栅网格,在网格中,我有一些感兴趣的值。但是,由于栅格单元内没有观测值,因此没有可使用的信息,因此可能会丢失某些值。
我想问什么是一些方法以及如何在 R 中实现它们以执行一些值的插值。
我知道当我有一些感兴趣的值和相关坐标时会执行基本的克里金法预测。然后,对栅格单元(网格)进行预测。
但是,是否有可能首先从栅格值开始,然后使用一些预测填充 NAN?
library(tidyverse)
library(sf)
library(raster)
# get data. using your example,we'll take spatial data from the raster package
Regions <- getData("GADM",country = "CZ",level = 0)
Regions <-
Regions %>%
st_as_sf()
grid_spacing <- 0.25
# Create a grid for the borber --------------------------------------------
polygony <-
st_make_grid(Regions,square = T,cellsize = c(grid_spacing,grid_spacing)) %>%
st_sf()
grid <- st_intersection(polygony,Regions)
# Add Values to the grid --------------------------------------------------
grid <-
grid %>%
mutate(Value = rnorm(n = nrow(grid),500,50))
# Creating NAN values to be predicted -------------------------------------
grid[1:23,'Value'] = NaN
grid[200:202,'Value'] = NaN
grid[80:100,'Value'] = NaN
grid[113:114,'Value'] = NaN
grid[133:144,'Value'] = NaN
# The NAN values here should be interpolated
grid %>%
ggplot(aes(fill = Value)) +
geom_sf()
解决方法
暂无找到可以解决该程序问题的有效方法,小编努力寻找整理中!
如果你已经找到好的解决方法,欢迎将解决方案带上本链接一起发送给小编。
小编邮箱:dio#foxmail.com (将#修改为@)