R/find_nearest_point.R

In aws.wrfsmn: Data Processing of SMN Hi-Res Weather Forecast from 'AWS'

#' Temporal series of closest location
#'
#' @description
#' Location of nearest point to lon/lat and temporal serie of location
#'
#' @usage find.nearest.point(data.spat.raster = data.spat.raster, lon = lon, lat = lat)
#'
#' @param data.spat.raster Spat Raster of WRF SMN (only one or several)
#' @param lon Longitude location of nearest point to find
#' @param lat Latitude location of nearest point to find
#'
#' @importFrom terra extract
#' @importFrom terra vect
#' @return a vector with the nearest location (lon/lat) and time serie of that location
#'
#'
#' @export
#'
find.nearest.point <- function(data.spat.raster = data.spat.raster,
                               lon = lon, lat = lat){

  aux <- extract(data.spat.raster, vect(cbind(lon, lat)), xy = TRUE, ID = FALSE)

  nearest.point <- data.frame(closest.lon = aux$x, closest.lat = aux$y)

  nearest.serie <- aux[, -c(length(aux)-1, length(aux))]

  return(cbind(nearest.point, nearest.serie))
}