Read a spatial file from a local folder, optionally crop it by a bounding box, and return it as a sf object readLayer is used within TADA_addPolys and TADA_addPoints

Usage

readLayer(layerfilepath, gpkg, layer, bbox = NULL)

Arguments

layerfilepath: Local path to the data folder containing the .gpkg file
gpkg: name of the .gpkg file
layer: name of the layer within the .gpkg file
bbox: A bounding box from the sf function st_bbox; used to filter the query results. Optional; defaults to NULL.

Value

sf object containing the layer

Examples

if (FALSE) { # \dontrun{
# Load example dataset
utils::data(Data_TribalNations_Harmonized)

# Use the example dataset bounding box
bbox <- sf::st_bbox(
  c(
    xmin = min(Data_TribalNations_Harmonized$TADA.LongitudeMeasure),
    ymin = min(Data_TribalNations_Harmonized$TADA.LatitudeMeasure),
    xmax = max(Data_TribalNations_Harmonized$TADA.LongitudeMeasure),
    ymax = max(Data_TribalNations_Harmonized$TADA.LatitudeMeasure)
  ),
  crs = sf::st_crs(Data_TribalNations_Harmonized)
)

# Read the American Indian Reservations feature layer,
# filtered by the bounding box for the Data_TribalNations_Harmonized
# example dataset
layerfilepath <- "extdata"
gpkg <- "Tribal.gpkg"
layer <- "AmericanIndian"
readLayer(layerfilepath, gpkg, layer)
} # }