PlotOnMapTiles function

plots on map tiles by "stitching" them together

Counterpart to PlotOnStaticMap for map tiles

PlotOnMapTiles(mt, lat, lon, center, size = c(768, 
    768), add = FALSE, FUN = points, mar = c(0, 0, 
    0, 0), verbose = 0, ...)

Arguments

• mt: list returned by GetMapTiles
• lat: latitude values to be overlaid, if any
• lon: longitude values to be overlaid, if any
• center: optional center
• size: size (in pixels) of "stitched" map
• add: start a new plot or add to an existing
• FUN: plotting function to use for overlay; typical choices would be points and lines
• mar: outer margin in plot; if you want to see axes, change the default
• verbose: level of verbosity
• ...: further arguments to be passed to FUN

Returns

nothing returned

Author(s)

Markus Loecher

Examples

if (0){
  lat = c(40.702147,40.718217,40.711614);
  lon = c(-74.012318,-74.015794,-73.998284);
  center = c(mean(lat), mean(lon));
  zoom <- min(MaxZoom(range(lat), range(lon)));
  bb=qbbox(lat,lon)
  
  manhattan_osm = GetMapTiles(latR =bb$latR , lonR=bb$lonR,zoom=zoom,verbose=1)
  PlotOnMapTiles(manhattan_osm,lat=lat,lon=lon,pch=20,col=c('red', 'blue', 'green'),cex=2)
  
  manhattan_goo = GetMapTiles(latR =bb$latR , lonR=bb$lonR,zoom=zoom,
                   tileDir= TRUE, type="google" )
  PlotOnMapTiles(manhattan_goo,lat=lat,lon=lon,pch=20,col=c('red', 'blue', 'green'),cex=2)
  
}