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