library(maptools)
nc.shp <- read.shape(system.file("shapes/sids.shp", package="maptools")[1])
data(ncshp)
nc1 <- as.SpatialRings.Shapes(nc.shp$Shapes, as.character(nc.shp$att.data$FIPS))
plotSpatialRings(nc1)
invisible(text(getSRSringsLabptSlots(nc1), labels=getSRSringsIDSlots(nc1), cex=0.6))
library(maps)
ncmap <- map("county", "north carolina", fill=TRUE, col="transparent", plot=FALSE)
data(ncmap)
IDs <- sapply(strsplit(ncmap$names, ","), function(x) x[2])
nc2 <- as.SpatialRings.map(ncmap, IDs)
plotSpatialRings(nc2)
invisible(text(getSRSringsLabptSlots(nc2), labels=getSRSringsIDSlots(nc2), cex=0.6))
library(RArcInfo)
fl <- "http://www.census.gov/geo/cob/bdy/co/co90e00/co37_d90_e00.zip"
download.file(fl, "co37_d90_e00.zip")
e00 <- zip.file.extract("co37_d90.e00", "co37_d90_e00.zip")
e00toavc(e00, "nc")
arc <- get.arcdata(".", "nc")
pal <- get.paldata(".", "nc")
pat <- get.tabledata("info", "NC.PAT")
data(co37_d90_arc)
data(co37_d90_pal)
sapply(pal[[2]], function(x) length(x[[1]]))
data(co37_d90_pat)
IDs <- paste(pat$ST[-1], pat$CO[-1], sep="")
nc3 <- as.SpatialRings.pal(arc, pal, IDs=IDs)
plotSpatialRings(nc3)
invisible(text(getSRSringsLabptSlots(nc3), labels=getSRSringsIDSlots(nc3), cex=0.6))
Run the code above in your browser using DataCamp Workspace