Adds a number of utilities for handling and visualising geographical data - for example choropleth mapping with 'nice' legends.
1 2 3 4 5 6 7 8 9 10 11 | # Load up the libraries needed
library(maptools)
library(RColorBrewer)
# Read in map data and compute a rate for mapping
sids <- readShapePoly(system.file("shapes/sids.shp", package="maptools")[1],
proj4string=CRS("+proj=longlat +ellps=clrk66"))
sids.rate=10000*sids$SID74/sids$BIR74
# Create the shading scheme, plot a choropleth map and add a legend
shades = auto.shading(sids.rate)
choropleth(sids,sids.rate,shades)
choro.legend(-83.77,37.87,shades,fmt="%4.1f",title='Rate per 10,000')
|
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.