data-raw/imap.R

In elipousson/mapmaryland: Easy Access to Maryland Spatial Data

library(getdata)
library(sfext)
library(dplyr)

state_fips <- 24
pkg_crs <- 3857
pkg_state_abb <- "MD"
pkg_state_name <- "Maryland"

# "https://data.imap.maryland.gov/datasets/maryland-waterbodies-rivers-and-streams-detailed"

url <-
  get_imap_url("rivers_and_streams_detailed")

md_streams_detailed <-
  read_sf_url(url = url)

md_streams_detailed <-
  md_streams_detailed %>%
  dplyr::transmute(
    type = layer
  ) %>%
  sf::st_transform(pkg_crs) %>%
  sf::st_join(dplyr::select(md_counties, countyfp, county = name), largest = TRUE) %>%
  rename_sf_col()

# NOTE: md_streams_detailed is no longer being exported

url <-
  get_imap_url("rivers_and_streams_generalized")

md_streams <-
  read_sf_url(url = url) %>%
  sf::st_transform(pkg_crs) %>%
  sf::st_join(dplyr::select(md_counties, countyfp, county = name), largest = TRUE) %>%
  rename_sf_col()

# NOTE: md_streams is no longer being exported

url <- get_imap_url("county_boundaries_detailed")

md_counties_detailed <-
  read_sf_url(url = url) %>%
  sf::st_transform(pkg_crs) %>%
  rename_sf_col()

md_counties_detailed <-
  dplyr::transmute(
    md_counties_detailed,
    statefp = "24",
    countyfp = stringr::str_pad(county_fip, width = 3, pad = "0"),
    name = county
  ) %>%
  dplyr::left_join(
    dplyr::select(
      sf::st_drop_geometry(md_counties),
      countyfp, geoid, namelsad
    ),
    by = "countyfp"
  )

usethis::use_data(
  md_counties_detailed,
  overwrite = TRUE
)

md_mpos <-
  get_imap_data(
    location = NULL,
    nm = "metropolitan_planning_organizations_boundaries",
    crs = 3857
  )

md_mpos <-
  md_mpos %>%
  getdata::rename_with_xwalk(
    xwalk = list(
      "id" = "mpo_id",
      "name" = "mpo_name",
      "abb" = "acronym",
      "url" = "mpo_url"
    )
  ) %>%
  dplyr::rowwise() %>%
  dplyr::mutate(
    states = paste0(c(state, state_2, state_3), collapse = ";"),
    states = stringr::str_remove_all(states, ";N/A")
  ) %>%
  dplyr::select(
    name, abb,
    geoid = id, url, states, geometry
  ) %>%
  dplyr::ungroup()

md_mpos <- sf::st_transform(as_sf(dplyr::ungroup(md_mpos)), crs = 3857)

usethis::use_data(md_mpos, overwrite = TRUE)

url <- "https://geodata.md.gov/imap/rest/services/BusinessEconomy/MD_IncentiveZones/FeatureServer/2"

md_arts_districts <- esri2sf::esri2sf(url, crs = 3857)

md_arts_districts <- md_arts_districts |>
  dplyr::select(
    name = sitename,
    desc = sitedesc,
    org_name = orgname,
    org_url = website,
    org_phone = phonenumbe,
    org_address = streetaddr,
    org_city = city,
    org_county = county,
    zipcode = zip,
    lon = longitude,
    lat = latitude,
  ) |>
  sfext::rename_sf_col()

usethis::use_data(md_arts_districts, overwrite = TRUE)