r - 如何将轮廓线变成填充轮廓？

Question

有谁知道像contourLines那样转动filled.contours多边形的输出以绘制为填充轮廓的方法。是否有顺序必须绘制多边形才能查看所有可用层？这是无效的示例代码片段:

#typical plot
filled.contour(volcano, color.palette = terrain.colors)

#try
cont <- contourLines(volcano)
fun <- function(x) x$level
LEVS <- sort(unique(unlist(lapply(cont, fun))))
COLS <- terrain.colors(length(LEVS))
contour(volcano)
for(i in seq(cont)){
    COLNUM <- match(cont[[i]]$level, LEVS)
    polygon(cont[[i]], col=COLS[COLNUM], border="NA")
}
contour(volcano, add=TRUE)

Answer 1

使用raster包(调用rgeos和sp)的解决方案。输出是一个SpatialPolygonsDataFrame，它将覆盖网格中的每个值:

library('raster')
rr <- raster(t(volcano))
rc <- cut(rr, breaks= 10)
pols <- rasterToPolygons(rc, dissolve=T)
spplot(pols)

Here's a discussion，它将向您展示如何简化(“美化”)结果多边形。