EVOLUTION-MANAGER

Edit File: webmap.R

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')