## Not run:
# # Load data
# data("com")
#
# # Get isochones with lon/lat coordinates, default breaks
# iso <- osrmIsochrone(loc = c(5.936036, 49.24882))
# plot(iso)
# points(5.936036, 49.24882, pch = 20, col = "red")
#
# # Map
# if(require("cartography")){
# osm <- getTiles(spdf = iso, crop = TRUE, type = "osmgrayscale")
# tilesLayer(osm)
# breaks <- sort(c(unique(iso$min), max(iso$max)))
# cartography::choroLayer(spdf = iso, df = iso@data,
# var = "center", breaks = breaks,
# border = NA,
# legend.pos = "topleft",legend.frame = TRUE,
# legend.title.txt = "Isochrones\n(min)",
# add = TRUE)
# }
#
# # Get isochones with a SpatialPointsDataFrame, custom breaks
# iso2 <- osrmIsochrone(loc = src[7,], breaks = seq(from = 0,to = 30, by = 5))
#
# # Map
# if(require("cartography")){
# osm2 <- getTiles(spdf = iso2, crop = TRUE, type = "osmgrayscale")
# tilesLayer(osm2)
# breaks2 <- sort(c(unique(iso2$min), max(iso2$max)))
# cartography::choroLayer(spdf = iso2, df = iso2@data,
# var = "center", breaks = breaks2,
# border = NA,
# legend.pos = "topleft",legend.frame = TRUE,
# legend.title.txt = "Isochrones\n(min)",
# add = TRUE)
# }
# ## End(Not run)
Run the code above in your browser using DataLab