Maps

Defining a bounding box

bbox_2 <- osmdata::getbb("Jayanagar, Bangalore, India")
bbox_2
       min      max
x 77.58034 77.58507
y 12.93683 12.94356
locations <-
  osmdata::opq(bbox = bbox_2) %>%
  osmdata::add_osm_feature(
    key = "amenity",
    value = c("restaurant", "atm", "college")
  ) %>%
  osmdata_sf() %>% # Convert to Simple Features format
  purrr::pluck("osm_points") # Pull out the data frame of interest
###
dat_buildings <-
  osmdata::opq(bbox = bbox_2) %>%
  osmdata::add_osm_feature(key = "building") %>%
  osmdata_sf() %>%
  purrr::pluck("osm_polygons")

Inspecting OSM data

str(locations)
Classes 'sf' and 'data.frame':  37 obs. of  22 variables:
 $ osm_id            : chr  "493493804" "493493806" "493493810" "1217571540" ...
 $ name              : chr  NA NA NA NA ...
 $ amenity           : chr  NA NA NA NA ...
 $ brand             : chr  NA NA NA NA ...
 $ brand:en          : chr  NA NA NA NA ...
 $ brand:hi          : chr  NA NA NA NA ...
 $ brand:kn          : chr  NA NA NA NA ...
 $ brand:pa          : chr  NA NA NA NA ...
 $ brand:pnb         : chr  NA NA NA NA ...
 $ brand:ur          : chr  NA NA NA NA ...
 $ brand:wikidata    : chr  NA NA NA NA ...
 $ brand:wikipedia   : chr  NA NA NA NA ...
 $ brand:wikipedia:pa: chr  NA NA NA NA ...
 $ name:en           : chr  NA NA NA NA ...
 $ name:hi           : chr  NA NA NA NA ...
 $ name:kn           : chr  NA NA NA NA ...
 $ name:pa           : chr  NA NA NA NA ...
 $ name:pnb          : chr  NA NA NA NA ...
 $ name:ur           : chr  NA NA NA NA ...
 $ operator          : chr  NA NA NA NA ...
 $ operator:wikidata : chr  NA NA NA NA ...
 $ geometry          :sfc_POINT of length 37; first list element:  'XY' num  77.6 12.9
 - attr(*, "sf_column")= chr "geometry"
 - attr(*, "agr")= Factor w/ 3 levels "constant","aggregate",..: NA NA NA NA NA NA NA NA NA NA ...
  ..- attr(*, "names")= chr [1:21] "osm_id" "name" "amenity" "brand" ...
str(dat_buildings)
Classes 'sf' and 'data.frame':  733 obs. of  26 variables:
 $ osm_id          : chr  "105734855" "105734857" "105734859" "130863019" ...
 $ name            : chr  NA NA NA "BBMP Complex (Jayanagar 2nd Block)" ...
 $ addr:city       : chr  NA NA NA NA ...
 $ addr:housenumber: chr  NA NA NA NA ...
 $ addr:postcode   : chr  NA NA NA NA ...
 $ addr:street     : chr  NA NA NA NA ...
 $ alt_name        : chr  NA NA NA NA ...
 $ amenity         : chr  NA NA NA "public_building" ...
 $ atm             : chr  NA NA NA NA ...
 $ brand           : chr  NA NA NA NA ...
 $ brand:wikidata  : chr  NA NA NA NA ...
 $ brand:wikipedia : chr  NA NA NA NA ...
 $ building        : chr  NA "yes" "yes" "yes" ...
 $ building:levels : chr  NA NA NA NA ...
 $ level           : chr  NA NA NA NA ...
 $ lit             : chr  NA NA NA NA ...
 $ name:en         : chr  NA NA NA NA ...
 $ name:kn         : chr  NA NA NA "ಬಿಬಿಎಂಪಿ ಸಂಕೀರ್ಣ" ...
 $ office          : chr  NA NA NA NA ...
 $ opening_hours   : chr  NA NA NA NA ...
 $ operator        : chr  NA NA NA NA ...
 $ phone           : chr  NA NA NA NA ...
 $ shop            : chr  NA NA NA NA ...
 $ smoking         : chr  NA NA NA NA ...
 $ website         : chr  NA NA NA NA ...
 $ geometry        :sfc_POLYGON of length 733; first list element: List of 1
  ..$ : num [1:14, 1:2] 77.6 77.6 77.6 77.6 77.6 ...
  .. ..- attr(*, "dimnames")=List of 2
  .. .. ..$ : chr [1:14] "1217571544" "1217571540" "493493806" "1217571541" ...
  .. .. ..$ : chr [1:2] "lon" "lat"
  ..- attr(*, "class")= chr [1:3] "XY" "POLYGON" "sfg"
 - attr(*, "sf_column")= chr "geometry"
 - attr(*, "agr")= Factor w/ 3 levels "constant","aggregate",..: NA NA NA NA NA NA NA NA NA NA ...
  ..- attr(*, "names")= chr [1:25] "osm_id" "name" "addr:city" "addr:housenumber" ...

Plotting Base Maps with ggplot/ggformula

ggplot() +
  geom_sf(data = locations, fill = "gold", color = "grey", linewidth = 0.025) +
  geom_sf(data = dat_buildings, fill = "purple") +
  coord_sf()

gf_sf(geometry = ~geometry, data = locations) %>%
  gf_sf(geometry = ~geometry, dat = dat_buildings)

Plotting with tmap

tm_shape(locations) +
  tm_dots(size = 0.5) +
  tm_shape(dat_buildings) +
  tm_polygons() +
  tm_compass(position = c("right", "top"))

Adding User Data to a Base Map

Back to top