r - 空间插值（克里金法），多边形而不是栅格输出

Question

有没有一种（简单的）方法可以将 SpatialPixelsDataFrame（来自 krige）转换为例如 SpatialPolygonsDataFrame（矢量图形而不是像素图）。

设置值范围并将栅格插入多边形或使用另一种生成 SpatialPolygonsDataFrame 的 krige 方法会很好。我期待一个简单的例子。

krige 示例：即书中https://oscarperpinan.github.io/spacetime-vis/ 8.1.5 Spatial Interpolation，完整源码https://github.com/oscarperpinan/spacetime-vis/blob/master/bubble.R

library(gstat)

airGrid <- spsample(NO2sp, type="regular", n=1e5)
gridded(airGrid) <- TRUE
airKrige <- krige(mean ~ 1, NO2sp, airGrid)

spplot(airKrige["var1.pred"],
       col.regions=colorRampPalette(airPal)) + ...

Answer 1

像这样的东西可能会起作用：

library(raster)
x <- raster(airKrige["var1.pred"])
y <- cut(x, c(10,20,30,40,50,60,70))
z <- rasterToPolygons(y, dissolve=TRUE)
spplot(z)