Nothing
# Copyright 2003-5 (c) Roger Bivand
subset.polylist <- function(x, subset, fit.bbox=TRUE, ...) {
.Deprecated("", package="maptools", msg="objects other than Spatial objects defined in the sp package are deprecated")
if (!inherits(x, "polylist")) stop("x not a polylist object")
if (!is.logical(subset)) stop("subset not a logical vector")
if (length(x) < 1L) stop("zero length polylist")
if (length(x) != length(subset)) stop("x and subset different lengths")
res <- subset.default(x, subset)
attr(res, "region.id") <- subset.default(attr(x, "region.id"), subset)
old.ids <- 1:length(x)
new.ids <- match(old.ids, which(subset))
after <- new.ids[subset.default(attr(x, "after"), subset)]
area <- sapply(res, function(x) attr(x, "area"))
if (any(sapply(area, is.null)))
pO <- order(subset.default(attr(x, "plotOrder"), subset))
else pO <- order(area, decreasing=TRUE)
attr(res, "after") <- after
attr(res, "plotOrder") <- pO
class(res) <- "polylist"
attr(res, "maplim") <- attr(x, "maplim")
if (fit.bbox) attr(res, "maplim") <- maplimFromBbox(res)
res
}
maplimFromBbox <- function(plist) {
bboxes <- sapply(plist, function(x) attr(x, "bbox"))
mapxlim <- range(c(bboxes[c(1,3),]))
mapylim <- range(c(bboxes[c(2,4),]))
list(x=mapxlim, y=mapylim)
}
Any scripts or data that you put into this service are public.
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.