DMEs

1 Paquetes a usar:

library(raster)
library(sf)
library(rgdal)
library(tidyverse)
library(leaflet)

1.1 Carga de datos del DME 253 + 500 A:

points <- read.csv("DME 253+500 A.csv", header = F)
names(points) <- c("id","x","y")
point <- st_as_sf(points,coords = c("x","y"),crs=32719)
polig <- st_as_sf(points,coords = c("x","y"),crs=32719) %>% 
  summarise(geometry=st_combine(geometry)) %>% 
  st_cast("POLYGON")

1.2 Carga de datos del DME 253 + 500 B:

points2 <- read.csv("DME 253+500 B.csv", header = F)
names(points2) <- c("id","x","y")
point2 <- st_as_sf(points2,coords = c("x","y"),crs=32719)
polig2 <- st_as_sf(points2,coords = c("x","y"),crs=32719) %>% 
  summarise(geometry=st_combine(geometry)) %>% 
  st_cast("POLYGON")

1.3 Aplicando Web mapping

1.3.1 Usando Basemap de Esri:

leaflet() %>% addTiles() %>% addProviderTiles(providers$Esri.WorldImagery) %>% 
  addPolygons(data=st_transform(polig,crs = 4326),
              label = "235+500 A",
                popup = "DME 235 + 500 A") %>% 
  addPolygons(data=st_transform(polig2,crs = 4326),
              label = "235+500 B",
              popup = "DME 235 + 500 B") %>% 
  setView(zoom = 16, lat = -13.674,ln=-70.47)

1.3.2 Usando Basemap de OpenStreetMap:

leaflet() %>% addTiles() %>% 
  addProviderTiles(providers$OpenStreetMap) %>% 
  addPolygons(data=st_transform(polig,crs = 4326),
              label = "235+500 A",
                popup = "DME 235 + 500 A") %>% 
  addPolygons(data=st_transform(polig2,crs = 4326),
              label = "235+500 B",
              popup = "DME 235 + 500 B")

1.4 Graficando con ggplot2

ggplot()+
  geom_sf(data=polig, color="deepskyblue", fill="transparent")+
  geom_sf(data=point,color="red")+
  geom_sf(data=polig2, color="darkgreen", fill="transparent")+
  geom_sf(data=point2,color="green")+
  theme_bw()+
  coord_sf(datum = 32719)+
  labs(x="Este",y="Norte", title = "DME 253 + 500 A y B")