data-raw/icao_firs.R

In euctrl-pru/pruatlas: The Atlas and Utilities for PRU Maps

# remotes::install_github("yonghah/esri2sf")
library(esri2sf)
library(dplyr)
url <- "https://gis.icao.int/ArcGIS/rest/services/FIRMSD/MapServer/18"
df <- esri2sf(url,
              outFields = c("FID", "Shape", "FIRname", "REGION",
                            "KIND", "UPPER", "LOWER", "ULC",
                            "ICAOCODE", "RESP", "NOM_COMP", "HISTORIC",
                            "REMARKS", "REMARKS2", "REMARKS3",
                            "centlong", "centlat",
                            "AREAsqkm", "PERIMEkm",
                            "SUPP_REGIO", "ET_ID")) %>%
  rename(fir_id = FID, fir_name = FIRname,
         icao_code = ICAOCODE,
         area_sqkm = AREAsqkm,  perimeter_km = PERIMEkm,
         centroid_lon = centlong, centroid_lat = centlat,
         region = REGION, supporting_region = SUPP_REGIO,
         responsible = RESP) %>%
  janitor::clean_names() %>%
  # make " " NA's
  mutate(across(where(is.character), str_trim),
         across(where(is.character), dplyr::na_if, ""))

df %>%
  write_sf("inst/extdata/icao_firs.geojson",
           delete_dsn = TRUE,
           delete_layer = TRUE)

# df %>%
#   filter(REGION == "EUR") %>%
#   filter(ULC == "LC", str_starts(ICAOCODE, "E") | str_starts(ICAOCODE, "L")) %>%
#   ggplot() +
#   geom_sf()