Global Administrative Areas (GADM) is a spatial database of the location of the world's administrative areas. Files are provided as Google Earth KMZ, shapefiles (.shp), ESRI geodatabase and as I was exited to see: as .RData files.
Download or load directly from the website.
This will create an object called gadm:
[r]
library(sp)
# let's take a small country...
con <- url("http://gadm.org/data/rda/CHE_adm1.RData")
print(load(con))
close(con)
[/r]
Conveniently, the object turns out to be a SpatialPolygonsDataFrame.
[r]
str(gadm, max.level=2)
Formal class 'SpatialPolygonsDataFrame' [package "sp"] with 5 slots
..@ data :'data.frame': 27 obs. of 16 variables:
..@ polygons :List of 27
..@ plotOrder : int [1:27] 10 6 24 25 17 22 27 7 13 1 ...
..@ bbox : num [1:2, 1:2] 5.96 45.82 10.49 47.81
.. ..- attr(*, "dimnames")=List of 2
..@ proj4string:Formal class 'CRS' [package "sp"] with 1 slots
[/r]
Add some data (from http://www.all-about-switzerland.info/swiss-cantons.html)
[r]
language <- c("german", "german", "german", "german","german","german","french","french","german","german","french", "french", "german", "french","german","german","german","german","german","german","german","italian","german","french","french","german","german")
gadm$language <- as.factor(language)
col = rainbow(length(levels(gadm$language)))
spplot(gadm, "language", col.regions=col, main="Swiss Language Regions")
[/r]
nice!
Is there a way to set frame border to 0 o cut it from the plot? I didn't find any option in man page.
I suggest you check out ?par.