Oikos: Oikos

#K. Wiese Agosto 2018
#####################

BarriosTGU <- function(lat, lon, Barrios, radio, proyecto){

  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"))
  
  spdfp <- spTransform(spdf, CRSobj = crs(Barrios))
  buffer <- gBuffer(spdfp, width = radio)
  
  
  
  #Barrios <- spTransform(Barrios, 
  #                       CRSobj = "+proj=utm +zone=16 +datum=WGS84 +units=m +no_defs")
  
  intersects <-  gIntersects (buffer, Barrios, byid=TRUE)
  
  # subset
  sel <- Barrios[as.vector(intersects),]
  
  #Grafico
  plot(extent(sel), col=NULL, xlab="X", ylab="Y", main=paste("Barrios en la Cercania del Proyecto", proyecto, sep=" "))
  plot(Barrios, border="#aaaaaa" ,col="grey",add=TRUE)
  plot(sel, col="red", add=TRUE)
  pointLabel(coordinates(sel), labels = sel$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)# border
  axis(1)
  axis(2)
  #axis(3)
  axis(4)
  box()
  
  #Resultado
  return(data.frame(Colonia = sel$NEIGH_UPPE, Hogares = sel$Hogares, Personas = sel$Poblacion))
}