Refs

Draft version, still lots to explore

library(maptools)
## Loading required package: sp
## Checking rgeos availability: TRUE
uk <- readShapeSpatial("map.shp", proj4string=CRS("+proj=longlat"))
plot(uk)

summary(uk)
## Object of class SpatialPolygonsDataFrame
## Coordinates:
##         min       max
## x -13.69139  1.764168
## y  49.86542 60.845822
## Is projected: FALSE 
## proj4string : [+proj=longlat +ellps=WGS84]
## Data attributes:
##      MM_UID            CODE                  NAME1               NAME2    
##  Min.   :  1.00   GB.ED  :  2   England         :112   Aberdeen     :  1  
##  1st Qu.: 48.75   GB.AB  :  1   Northern Ireland: 26   Aberdeenshire:  1  
##  Median : 96.50   GB.AD  :  1   Scotland        : 32   Anglesey     :  1  
##  Mean   : 96.50   GB.AG  :  1   Wales           : 22   Angus        :  1  
##  3rd Qu.:144.25   GB.AM  :  1                          Antrim       :  1  
##  Max.   :192.00   GB.AN  :  1                          Ards         :  1  
##                   (Other):185                          (Other)      :186  
##          DIVISION  
##  Kingdom     :144  
##  Principality: 22  
##  Province    : 26  
##                    
##                    
##                    
## 
library(maps)

france<-map(database="france")

names(france)
## [1] "x"     "y"     "range" "names"
head(france$names)
## [1] "Nord"           "Pas-de-Calais"  "Somme"          "Nord:1"        
## [5] "Ardennes"       "Seine-Maritime"
library(mapdata)

map(database="worldHires", "Portugal") 

library(RgoogleMaps)

lat <- c(37,42) #define our map's ylim
lon <- c(-9,-6) #define our map's xlim
center = c(mean(lat), mean(lon))  #tell what point to center on
zoom <- 7  #zoom: 1 = furthest out (entire globe), larger numbers = closer in

terrmap <- GetMap(center=center, zoom=zoom, maptype= "terrain", destfile = "terrain.png") 
# lots of visual options, just like google maps: 
#   maptype = c("roadmap", "mobile", "satellite", "terrain", "hybrid", "mapmaker-roadmap", "mapmaker-hybrid")
PlotOnStaticMap(terrmap)