Кто-нибудь знает способ превратить вывод полигонов 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)