library(sp)
library(ggmap)
demo(meuse, ask=FALSE, echo = FALSE)
merc = CRS("+init=epsg:3857")
WGS84 = CRS("+init=epsg:4326")
meuse = spTransform(meuse, WGS84)
bgMap = get_map(as.vector(bbox(meuse)), source = "google", zoom = 13) # useless without zoom level
# plot with ggmap-google bg:
plot(spTransform(meuse, merc), bgMap = bgMap, pch = 16, cex = .5)
# spplot with ggmap-google bg:
spplot(spTransform(meuse, merc), c("zinc", "lead"), colorkey = TRUE,
sp.layout = list(panel.ggmap, bgMap, first = TRUE))
# plot with ggmap-osm bg:
bb = t(apply(bbox(meuse), 1, bbexpand, .04))
bgMap = get_map(as.vector(bb), source = "osm") # WGS84 for background map
plot(spTransform(meuse, merc), bgMap = bgMap, pch = 16, cex = .5)
# RgoogleMaps:
center = apply(coordinates(meuse), 2, mean)[2:1]
library(RgoogleMaps)
g = GetMap(center=center, zoom=13) # google
par(mar = rep(0,4)) # fill full device
plot(spTransform(meuse, merc), bgMap = g, pch = 16, cex = .5)
spplot(spTransform(meuse, merc), c("zinc", "lead"), colorkey = TRUE,
sp.layout = list(panel.RgoogleMaps, g, first = TRUE),
scales = list(draw = TRUE))
# Norway boundary example:
library(cshapes)
cshp = cshp(as.Date("2000-01-1"))
norway = cshp[cshp$ISO1AL2 == "NO",]
bgMap = get_map(as.vector(bbox(norway))) # no is already in WGS84
plot(spTransform(norway, merc), bgMap = bgMap, border = 'red')
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.