data-raw/transit_routes.R

# prepare transit routes

requireNamespace("readxl", quietly = TRUE)
requireNamespace("fs", quietly = TRUE)
requireNamespace("tigris", quietly = TRUE)
requireNamespace("janitor", quietly = TRUE)

library(dplyr)
library(fs)
library(sf)
library(tigris)
library(janitor)


# ## transit routes ---------------
# temp <- tempfile()
# download.file("ftp://ftp.gisdata.mn.gov/pub/gdrs/data/pub/us_mn_state_metc/trans_transit_routes/gpkg_trans_transit_routes.zip",
#   destfile = temp
# )
# trans_routes <- sf::read_sf(unzip(temp, "trans_transit_routes.gpkg")) %>%
#   st_transform(4326)
# 
# fs::file_delete("trans_transit_routes.gpkg")
# 
# usethis::use_data(trans_routes, overwrite = TRUE)
# 


## transit stops ---------------
temp <- tempfile()
temp2 <- tempfile()
download.file("ftp://ftp.gisdata.mn.gov/pub/gdrs/data/pub/us_mn_state_metc/trans_transit_stops/shp_trans_transit_stops.zip",
  destfile = temp
)

unzip(zipfile = temp, exdir = temp2) # list.files(temp2)
trans_stops <- sf::read_sf(paste0(temp2, pattern = "/TransitStops.shp")) %>%
  filter(
    busstop_yn == "Y",
    board_flag == 1
  ) %>%
  st_transform(4326)

# trans_stops %>% ggplot() + geom_sf(aes(col = busstop_yn))

usethis::use_data(trans_stops, overwrite = TRUE)
Metropolitan-Council/economic-values-atlas documentation built on Dec. 17, 2021, 4:13 a.m.