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

【问题标题】：Spatial interpolation (Kriging), polygon instead of raster output空间插值（克里金法），多边形而不是栅格输出
【发布时间】：2015-06-08 23:21:50
【问题描述】：

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

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

克里格示例：即在书中https://oscarperpinan.github.io/spacetime-vis/ 8.1.5 空间插值，完整源代码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)) + ...

【问题讨论】：

我可能不知道答案，但我尝试运行代码并最终卡在试图从madrid.org/nomecalles/Callejero_madrid.icm 或madrid.org/nomecalles/DescargaBDTCorte.icm 找到所需的zip 文件

标签： r polygon geospatial lattice kriging

【解决方案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)

【讨论】：