#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))
}
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.