R/covGIS.R

In gcostaneto/envirotype: GIS platforms to support envirorytping in breeding programs

#'@title  Collect covariates from raster files
#'
#'
#' @description Imports environmental information from GIS spatial data
#' @author Germano Martins F Costa Neto <germano.cneto@usp.br>
#' @param .covraster vector containg site identification
#' @param .lat vector (numeric) containing latitude values
#' @param .lon vector (numeric) containing longitude values
#' @param .env vector(chr) containing environmental id
#' @param .path output directorie
#' @importFrom raster extract
#'

covGIS <-function(.covraster=NULL,.reference=NULL, .lon=NULL, .lat=NULL, .env=NULL, .crs=NULL,.path=NUL){
  
  leg  = .reference
  loc  = data.frame(x=leg[,.lon],y=leg[,.lat])
  site = leg[,.env]
  coordinates(loc)= ~x+y
  proj4string(loc) = CRS("+proj=longlat +datum=WGS84")
  
  for(i in 1:length(names(.covraster))){site = cbind(site,data.frame(extract(.covraster[[i]], loc)))}
  names(site)[-1] = names(.covraster)
  return(data.frame(unique(site)))
}