空间插值（克里格），多边形而不是光栅输出

有没有一种（简单）的方法来将SpatialPixelsDataFrame（从krige）转换为一个SpatialPolygonsDataFrame（vectorgraphic而不是pixelgraphic）。空间插值（克里格），多边形而不是光栅输出

设置值范围并将栅格内插到多边形或使用另一个生成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)) + ...

enter image description here

来源

2015-04-03 ckluss

我可能不要知道答案，但我试图运行代码，并最终卡住试图找到所需的zip文件从http：//www.madrid.or g/nomecalles/Callejero_madrid.icm或http://www.madrid.org/nomecalles/DescargaBDTCorte.icm – 2015-04-04 03:24:23

像这样的东西可能会奏效：

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)

来源

2015-04-04 04:17:24 RobertH

空间插值（克里格），多边形而不是光栅输出

回答

相关问题