compare_urban_metro.R

library(leaflet)
library(tigris)
library(sp)

compare_urban_metro <- function(metro_name) {
  
  cbsa <- core_based_statistical_areas(cb = TRUE, resolution = '20m')
  ua <- urban_areas(cb = TRUE)
  
  dfw_cbsa <- cbsa[cbsa$NAME == metro_name, ]
  o <- over(dfw_cbsa, ua, returnList = TRUE)
  
  dfw_ua <- ua[ua$NAME10 %in% o[[1]]$NAME10, ]
  
  leaflet() %>%
    addTiles() %>%
    addPolygons(data = dfw_cbsa) %>%
    addPolygons(data = dfw_ua, color = 'red', fillColor = 'red')
  
}

compare_urban_metro("Dallas-Fort Worth-Arlington, TX")