## Not run:
# library(rgdal)
#
# # get a map using names
# g = gmap('Australia')
# plot(g, inter=TRUE)
# gs = gmap('Sydney, New South Wales, Australia', type='satellite')
# plot(gs, inter=TRUE)
# gs = gmap('Sydney, Australia', type='satellite', exp=3)
# plot(gs, inter=TRUE)
# gs = gmap('Sydney, Australia', type='hybrid', zoom=10, scale=2)
# plot(gs, inter=TRUE)
#
# # from a maxtrix with lon/lat points
# x = runif(30)*10 + 40
# y = runif(30)*10 - 20
# xy = cbind(x, y)
# g = gmap(xy, type='hybrid')
# plot(g, inter=TRUE)
# points(Mercator(xy) , col='red', pch=20)
#
# # or from an Extent object
# e = extent( -121.9531 , -120.3897 , 35.36 , 36.61956 )
# # you can also get an Extent object by clicking on the map twice after using:
# # drawExtent()
# r = gmap(e)
# plot(r, interpolate=TRUE)
#
# # transform points to Mercator for plotting on top of map:
# pt <- matrix(c(-121, 36), ncol=2)
# ptm <- Mercator(pt)
# points(ptm, cex=3, pch=20, col='blue')
# Mercator(ptm, inverse=TRUE)
#
# # transform Spatial objects to Mercator for plotting on top of map
# # here for points, but particularly relevant for lines and polygons
# pt <- data.frame(pt)
# coordinates(pt) <- ~X1 + X2
# proj4string(pt) <-"+proj=longlat +datum=WGS84 +ellps=WGS84"
# ptm2 <- spTransform(pt, CRS("+proj=merc +a=6378137 +b=6378137 +lat_ts=0.0
# +lon_0=0.0 +x_0=0.0 +y_0=0 +k=1.0 +units=m +nadgrids=@null +no_defs"))
# points(ptm, col='red', pch='x', cex=3)
#
# # styles:
# g <- gmap("Brooklyn", style="feature:road.local|element:geometry|hue:0x00ff00|saturation:100
# &style=feature:landscape|element:geometry|lightness:-100", type='roadmap')
# plot(g)
# ## End(Not run)
Run the code above in your browser using DataLab