Áreas del Acuerdo de Destino para el Puente Vehicular Nichupté

Superficie total: 47005.08 m2

El Puente Vehicular Nichupté, el cual consiste en la construcción de un puente vehicular sobre el sistema lagunar Nichupté, con una longitud de 8.80 km, incluyendo la construcción de dos accesos. La sección transversal consta de una corona de 14.9 metros de ancho, diseñada para alojar tres carriles de 3.5 metros (uno por sentido más un tercer carril reversible) y un andador peatonal. Esta infraestructura permitirá la conexión entre el Blvd. Kukulcán, ubicado en la zona hotelera, y la Av. Bonampak, situada en el centro de Cancún, proporcionando así una vía alterna de acceso a la zona hotelera.

# Instalar tmap si no está instalado
if (!require(tmap)) install.packages("tmap")
if (!require(stars)) install.packages("stars")
if (!require(sf)) install.packages("sf")
if (!require(dplyr)) install.packages("dplyr")



library(sf)
library(tmap)
library(dplyr)
library(stars) 
library(tmaptools) 
library(dplyr)
#library (tidyverse)
#library (leaflet)

tmap_mode("view") #Para activar el modo interactivo del mapa 

shp_path <- "F:/R/PVN"

poligonos <- st_read(shp_path)

Reading layer `poligonos_4.7ha' from data source `F:\R\PVN' using driver `ESRI Shapefile'
Simple feature collection with 15 features and 2 fields
Geometry type: POLYGON
Dimension:     XYZ
Bounding box:  xmin: 517951.8 ymin: 2336434 xmax: 519161.8 ymax: 2337740
z_range:       zmin: 0 zmax: 0
Projected CRS: WGS 84 / UTM zone 16N

names(poligonos) #verificar en contenido del SHP

[1] "Layer"    "Poligono" "geometry"

str(poligonos) #visualizar el contenido del SHP

Classes 'sf' and 'data.frame':  15 obs. of  3 variables:
 $ Layer   : chr  "ZI" "ZI" "ZFMT" "ZFMT" ...
 $ Poligono: chr  "Poligono3_ZI" "Poligono2_ZI" "Poligono11_ZFMT" "Poligono10_ZFMT" ...
 $ geometry:sfc_POLYGON of length 15; first list element: List of 1
  ..$ : num [1:126, 1:3] 518378 518258 518252 518205 518205 ...
  ..- attr(*, "class")= chr [1:3] "XYZ" "POLYGON" "sfg"
 - attr(*, "sf_column")= chr "geometry"
 - attr(*, "agr")= Factor w/ 3 levels "constant","aggregate",..: NA NA
  ..- attr(*, "names")= chr [1:2] "Layer" "Poligono"

 #2. CALCULAR SUPERFICIE (ha)
#---------------------------------------------------------
poligonos <- poligonos %>%
  mutate(
    area_m2 = st_area(poligonos$geometry),
    area_ha = round (as.numeric(area_m2) / 10000,2)
  )

#---------------------------------------------------------
# 3. CARGAR ORTOFOTO (raster)
#---------------------------------------------------------
#ortho_path <- "C://PVN-Ortos/CANCUN PT1.tif"
#ortho <- read_stars(ortho_path)

#st_get_dimension_values(ortho, "band") #cuantas bandas tiene el Tif

#library(raster)

#r <- brick(ortho_path)
#r <- r[[1:3]]  # solo RGB
#ortho <- st_as_stars(r)

#print(ortho)
#st_get_dimension_values(ortho, "band") #ver la ortofoto y sus bandas, y dependiendo de eso se puede usar el codigo para leerlo

#---------------------------------------------------------
# 4. MAPA CON tmap
#---------------------------------------------------------
ls()

[1] "poligonos" "shp_path" 

#View(ortho) 
#names(poligonos) #ver contenido del SHP
#table(poligonos$Layer) #resumen de los caracteres que contiene el SHP
#table(poligonos$Layer)
#sf::st_read("F:/R/PVN", quiet = FALSE, options = "STRICT=YES") #muestra errores estructurales del SHP

poligonos$Layer <- as.factor(poligonos$Layer)
#-----CREAR MAPA-------------------------------------------------------------
#tm_shape(ortho) +
 #   tm_rgb(r = "red", g = "green", b = "blue") +
###leaflet() %>% 
  ###addProviderTiles(providers$OpenStreetMap) %>%
#------------# Crear etiqueta para la leyenda--------------------------
poligonos <- poligonos %>%
  mutate(leyenda = paste0(Layer, " - ", area_ha, " ha"))
#tmap_mode("plot")

tm_basemap("OpenStreetMap") +
tm_shape(poligonos) +
  tm_polygons(
    col = "Layer",
    alpha = 0.6,
      palette = "Set2", 
    border.col = "black",
    title = "Superficie (ha)"
  ) +
  tm_text("Layer", size = 0.4, col = "black",
          auto.placement = TRUE,
          bg.color = "WHITE", bg.alpha = 0.5) +
  tm_layout(
    legend.outside = TRUE,
    frame = TRUE,
    main.title = "Polígonos y superficies",
    main.title.size = 1.2
  )