R/clean_dup_by_year.R

In luismurao/hsi: "Habitat suitability index"

#' Function to clean duplicated coordinates data
#' @description Clean duplicated longitude and latitude data by year using a threshold distance.
#' @param this_species Species Temporal Data object see \code{\link[hsi]{sp_temporal_data}}.
#' @param threshold A numerc value representig the distance (in grades) between coordinates
#'            to be considered as a duplicate.
#' @seealso For the current version of this functions that work with data frames go to \code{\link[hsi]{clean_dup_by_year_df}}
#' @export
#'

clean_dup_by_year <- function(this_species,threshold){
  stopifnot(inherits(this_species, "sp.temporal.modeling"))
  df_occs_year <- this_species$occs_data
  df_occs_year$this_species <- as.numeric(as.character(df_occs_year$years_in_occs))
  clean_by_year <- df_occs_year %>%
    split(.$years_in_occs,drop=T) %>%
    purrr::map_df(~clean_dup(data = .x,
                             longitude = this_species$lon_lat_vars[1],
                             latitude = this_species$lon_lat_vars[2],
                             threshold = threshold))

  sp.temp.data.clean <- list(sp_coords = clean_by_year[,this_species$lon_lat_vars],
                             sp_occs_year = clean_by_year$years_in_occs,
                             oocs_data = clean_by_year,
                             lon_lat_vars = this_species$lon_lat_vars ,
                             layers_path_by_year = this_species$layers_path_by_year)
  class(sp.temp.data.clean) <- c("list", "sp.temporal.modeling")


  return(sp.temp.data.clean)
}
. <- NULL