data-raw/sf_green.R

In EcoDynIZW/d6berlin: A collection of template maps for visualizing ecological studies in Berlin

## code to prepare `sf_green` dataset goes here

d6berlin::download_data_berlin()

## FILE PATHS ----------------------------------------------------------------
shp_path  <- "./data-raw/geo-raw/berlin_shapes"

## PREPARE DATA --------------------------------------------------------------

sf_green <-
  suppressMessages(
    ## Berlin Landuse Categories (WGS 84)
    sf::read_sf(dsn = paste0(shp_path, "/gis_osm_landuse_a_free_1.shp"),
                layer = "gis_osm_landuse_a_free_1") |>
    ## Berlin Natural Areas (WGS 84)
    rbind(sf::read_sf(dsn = paste0(shp_path, "/gis_osm_natural_a_free_1.shp"),
                      layer = "gis_osm_natural_a_free_1")) |>
    dplyr::filter(fclass %in% c("forest", "grass", "meadow", "nature_reserve",
                                "scrub", "heath", "beach", "cliff")) |>
      dplyr::mutate(osm_id = as.factor(osm_id)) |>
      sf::st_intersection(d6berlin::sf_districts)
  )

usethis::use_data(sf_green, overwrite = TRUE)