if (FALSE) {
dhs <- rings(fdhs, N = c(100,200,300,400,500))
radius.N300 <- krige('r.radius', dhs, N = 300, nb.cells = 50)
prev.krige <- krige(r.wprev ~ 1, dhs, N = c(100, 300, 500))
plot(prev.krige, lty = 0)
library(ggplot2)
ggplot(prev.krige) +
aes(fill = r.wprev.N300.RInf) +
geom_sf(colour = "transparent") +
scale_fill_gradientn(colors = prevR.colors.red()) +
theme_prevR_light()
# Export r.wprev.N300.RInf surface in ASCII Grid
r <- terra::rast(stars::st_rasterize(prev.krige))
# writeRaster(r[[2]], "wprev.N300.asc")
}
Run the code above in your browser using DataLab