R/Rios.R

In klauswiese/Oikos: Oikos

#Rios Cercanos a Proyecto
#K. Wiese Agosto 2018
#########################

Rios <- function(lat, lon, Rios, radio, proyecto, fondo = NULL){
  
  df_xy <- data.frame(proyecto=proyecto,lon=lon,lat=lat)
  xy <- df_xy[,c(2,3)]
  
  spdf <- SpatialPointsDataFrame(coords = xy, data=df_xy,
                                 proj4string = CRS("+proj=longlat +datum=WGS84 +ellps=WGS84 +towgs84=0,0,0"))
  buffer <- gBuffer(spdfp, width = radio + 0.2*radio)
  spdfp <- spTransform(spdf, CRSobj = crs(Rios))  
  Distancias <- as.vector(gDistance(spdfp, Rios, byid= TRUE))
  
  Sel <- Rios[as.vector(gDistance(spdfp, Rios, byid= TRUE)) <= radio,]
  Dist <- Distancias[as.vector(gDistance(spdfp, Rios, byid= TRUE)) <= radio]
  #Grafico
  plot(extent(buffer), col=NULL, xlab="X", ylab="Y", main=paste("Rios en la Cercania del Proyecto", proyecto, sep=" "))
  plot(fondo, border="#aaaaaa" ,col="grey",add=TRUE)
  plot(Sel, col="darkblue",add=TRUE, lwd=2)
  pointLabel(coordinates(fondo), labels = fondo$NEIGH_UPPE, col= "white", cex = 0.7)
  points(spdfp, pch=16, col="black")
  legend("bottomright", 
         legend = paste("Ubicación", proyecto, sep=" "),
         pch = 16,
         col = "black",
         cex = 1)
  axis(1)
  axis(2)
  axis(4)
  box()
  
  #Resultado
  return(data.frame(Rio = Sel$NOMBRE, Distancia = Dist))
}