R/check_atpol_square.R

In atpolR: ATPOL Grid Implementation

#' Reverse engineering of published ATPOL grids
#'
#' check_atpol_square() do a reverse engineering of published ATPOL grids species,
#' especially those published in \insertCite{zajacAtlasRozmieszczeniaRoslin2001;textual}{atpolR}
#' @importFrom sf st_buffer
#' @importFrom terra crop vect
#' @importFrom Rdpack reprompt
#' @param centroid Simple Feature point geometry for which the check is performed,
#' usually it corresponds to centroid of ATPOL 10km x 10km grid
#' @param raster geocoded raster, it has to be in EPSG:2180 projection
#' @param distance \code{st_buffer} distance from centroid point for which the check is done, default 1200 m
#' @returns "YES" or "?" for given SF point
#'
#' @usage check_atpol_square(centroid, raster, distance)
#' @export
#' @references
#'     \insertAllCited{}
#'
check_atpol_square <- function(centroid, raster, distance = 1200) {
  if (!grepl("ID\\[\"EPSG\",2180\\]\\]",terra::crs({{raster}}))) {
    warning("Raster has to be in EPSG:2180 projection")
    warning("Try terra::project(raster, \"EPSG:2180\")")
    stop()
  }
  c <- sf::st_buffer({{centroid}}, dist = {{distance}})
  mm <- terra::minmax(terra::crop({{raster}}, terra::vect(c)))
  y <- ifelse(mm[1] == 0 & mm[2] == 0, "YES", "?")
  return(y)
}