knitr::opts_chunk$set( echo = FALSE, warning = FALSE, message = FALSE ) m

GEOVISUALIZADOR DEL ÁREA NATURAL PROTEGIDA NEVADO DE TOLUCA

1. LIBRERÍAS

library(sf) library(sp) library(leaflet) library(leafem) library(ggplot2) library(dplyr) library(terra) library(htmlwidgets)

2. CARGA DE SHAPEFILES

ANP <- st_read(“C:/FINAL SIG II/ANP Nevado de Toluca.shp”)

Municipios <- st_read(“C:/FINAL SIG II/Municipios.shp”)

Localidades <- st_read(“C:/FINAL SIG II/localidades.shp”)

UsoSuelo <- st_read(“C:/FINAL SIG II/Uso de suelo.shp”)

Vegetacion <- st_read(“C:/FINAL SIG II/Vegetación.shp”)

Cobertura <- st_read(“C:/FINAL SIG II/Cobertura arborea.shp”)

Hidrologia <- st_read(“C:/FINAL SIG II/Hidrología.shp”)

Agua <- st_read(“C:/FINAL SIG II/Cuerpos de agua.shp”)

Incendios <- st_read(“C:/FINAL SIG II/Incendios Forestales.shp”)

Clima <- st_read(“C:/FINAL SIG II/Clima.shp”)

3. TRANSFORMAR A WGS84 (EPSG 4326)

ANP <- st_transform(ANP,4326) Municipios <- st_transform(Municipios,4326) Localidades <- st_transform(Localidades,4326) UsoSuelo <- st_transform(UsoSuelo,4326) Vegetacion <- st_transform(Vegetacion,4326) Cobertura <- st_transform(Cobertura,4326) Hidrologia <- st_transform(Hidrologia,4326) Agua <- st_transform(Agua,4326) Incendios <- st_transform(Incendios,4326) Clima <- st_transform(Clima,4326)

4. EXTENSIÓN DEL MAPA

caja <- st_bbox(ANP) %>% as.vector()

5. LOGOS

logo_uaemex <- “C:/FINAL SIG II/log_uni.png”

6. MAPA LEAFLET

m <- leaflet() %>%

addProviderTiles(“OpenStreetMap”) %>%

# ANP

addPolygons( data = ANP, color = “red”, weight = 4, fill = FALSE, label = ~NOMBRE, group = “ANP Nevado de Toluca” ) %>%

# Municipios

addPolygons( data = Municipios, color = “black”, weight = 1, fillOpacity = 0.15, label = ~NOMGEO, group = “Municipios” ) %>%

# Localidades

addPolygons( data = Localidades, color = “orange”, weight = 1, fillOpacity = 0.5, label = ~NOMGEO, group = “Localidades” ) %>%

# Uso de suelo

addPolygons( data = UsoSuelo, color = “gold”, weight = 1, fillOpacity = 0.4, label = ~DESCRIPCIO, group = “Uso de suelo” ) %>%

# Vegetación

addPolygons( data = Vegetacion, color = “darkgreen”, weight = 1, fillOpacity = 0.4, label = ~DESCRIPCIO, group = “Vegetación” ) %>%

# Cobertura arbórea

addPolygons( data = Cobertura, color = “green”, weight = 1, fillOpacity = 0.4, label = ~DESCRIPCIO, group = “Cobertura arbórea” ) %>%

# Hidrología

addPolylines( data = Hidrologia, color = “blue”, weight = 2, label = ~DESCRIPCIO, group = “Hidrología” ) %>%

# Cuerpos de agua

addPolygons( data = Agua, color = “cyan”, fillColor = “cyan”, fillOpacity = 0.6, label = ~nombre, group = “Cuerpos de agua” ) %>%

# Incendios forestales

addPolygons( data = Incendios, color = “hotpink”, fillColor = “hotpink”, fillOpacity = 0.6, label = ~TIPOINC, group = “Incendios forestales” ) %>%

# Clima

addPolygons( data = Clima, color = “purple”, weight = 1, fillOpacity = 0.3, label = ~DESC_CLIMA, group = “Clima” ) %>%

7. CONTROL DE CAPAS

addLayersControl( overlayGroups = c( “ANP Nevado de Toluca”, “Municipios”, “Localidades”, “Uso de suelo”, “Vegetación”, “Cobertura arbórea”, “Hidrología”, “Cuerpos de agua”, “Incendios forestales”, “Clima” ), options = layersControlOptions(collapsed = FALSE) ) %>%

8. LEYENDA

addLegend( position = “bottomright”, colors = c( “red”, “black”, “orange”, “gold”, “darkgreen”, “green”, “blue”, “cyan”, “hotpink”, “purple” ), labels = c( “ANP”, “Municipios”, “Localidades”, “Uso de suelo”, “Vegetación”, “Cobertura arbórea”, “Hidrología”, “Cuerpos de agua”, “Incendios forestales”, “Clima” ), title = “Leyenda”

) %>%

9. ZOOM

fitBounds( lng1 = caja[1], lat1 = caja[2], lng2 = caja[3], lat2 = caja[4] )

10. LOGOS

file.exists(“C:/FINAL SIG II/logo_fgeo.jpg”)

file.exists(“C:/FINAL SIG II/logo_uaemex.jpg”)

logo_fgeo <- “C:/FINAL SIG II/logo_fgeo.jpg” logo_uaemex <- “C:/FINAL SIG II/logo_uaemex.jpg”

m <- m %>% leafem::addLogo( img = logo_uaemex, width = 120, position = “topleft” )

m <- m %>% leafem::addLogo( img = “C:/FINAL SIG II/logo_fgeo.jpg”, width = 100, position = “bottomleft” ) m

12. EXPORTAR A HTML

saveWidget( m, file = “C:/FINAL SIG II/Geovisualizador_Nevado_Toluca.html”, selfcontained = TRUE )

file.exists(“C:/FINAL SIG II/Geovisualizador_Nevado_Toluca.html”)

```