R/parse_location.R

In neotoma2: Working with the Neotoma Paleoecology Database

#' @title parse_location
#' @author Socorro Dominguez \email{s.dominguez@ht-data.com}
#' @import gtools
#' @import lubridate
#' @import geojsonsf
#' @returns `sf` object to parse as location.
#' @importFrom methods new
#' @description
#' Retrieve location from datasets WTK, GeoJson, bounding box
#' @keywords internal
#' @noRd
#' @param x location object
parse_location <- function(x) {
  if (is.numeric(x)) {
    # We're getting a numeric vector of coordinates:
    assertthat::assert_that(length(na.omit(x)) == 4,
                            msg = "Numeric coordinates need to be an array of 4 units, c(xmin, xmax, ymax, ymin)")
    my_bbox <- sf::st_bbox(c(xmin = x[1], xmax = x[3],
                             ymax = x[4], ymin = x[2]),
                           crs = st_crs(4326))
    my_bbox <- sf::st_as_sfc(my_bbox)
    new_geojson <- geojsonsf::sfc_geojson(my_bbox)
    loc <- new_geojson[1]
    return(loc[1])
  } else if (class(x)[1] == "sf") {
    geom <- sf::st_geometry(x)
    if ("sfc" %in% class(geom)){
        loc <- geojsonsf::sfc_geojson(geom)
    } else {
      loc <- geojsonsf::sf_geojson(x)
    }
    return(loc)
  } else if ("sfg" %in% class(x)) {
    geo <- sf::st_sfc(x)
    loc <- geojsonsf::sfc_geojson(geo)
    return(loc)
  } else if (is.character(x)) {
    loc1 <- try(geojsonsf::geojson_sf(x), silent = TRUE)
    if ("try-error" %in% class(loc1)) {
      loc <- sf::st_as_sf(wk::new_wk_wkt(x))
      loc <- parse_location(loc)
    } else {
      loc <- parse_location(loc1)
    }
    return(loc)
  }
}