In fdzul/rgeomex: The rgeomex provides the polygon data of INEGI 2019 and INE 2020

mapview::mapviewOptions(leafletWidth = 700, 
                        leafletHeight = 300)
knitr::opts_chunk$set(
  collapse = TRUE,
  comment = "#>"
)

library(rgeomex)

Para extraer una localidad de interes se usa la funcion extract_locality o la función extract_ageb, la primera extrae el polígono de la localidad y la segunda función extrae el polígono de la localidad y los AGEBs (Área Geostadística Básica) de la localidad. La primera regresa un objeto sf y la segunda una lista con dos objetos sf.

Para extraer una localidad con la función extract_locality, se usan dos argumentos, uno para definir la clave del estado (cve_edo) y el otro para definir el nombre de la localidad de acuerdo al INEGI 2020 (locality).

acapulco <- rgeomex::extract_locality(cve_edo = "12",
                                      locality = "Acapulco de Juárez")

mapview::mapview(acapulco, legend = FALSE)

Para extraer una localidad con la función extract_ageb, se usan dos argumentos, uno para definir la clave del estado (cve_edo) y el otro para definir el nombre de la localidad de acuerdo al INEGI 2020 (locality).

acapulco_ageb <- rgeomex::extract_ageb(cve_edo = "12",
                                      locality = "Acapulco de Juárez")

mapview::mapview(acapulco_ageb$locality, legend = FALSE)

Ambas funciones se pueden usar para todas las localidades urbanas de México, excepto la Ciudad de México para esta localidad se usa otro código para extraer toda la localidad.

rgeomex::loc_inegi19_mx|>
    dplyr::filter(CVE_ENT == "09")|>
mapview::mapview(legend = FALSE)

El argumento de la localidad acepta un vector con los nombres de las localidades que conforman el área metropolitana.

guadalajara <- rgeomex::extract_locality(cve_edo = "14",
                                      locality = c("Guadalajara",
                                                   "Zapopan",
                                                   "Tlaquepaque",
                                                   "Tonalá"))

mapview::mapview(guadalajara, legend = FALSE)

Tambien se puede usar base R para visualizar los bordes de la ciudad o el paquete ggpplot2, tmap o sfmap, etc.

ggplot2::ggplot()+
  ggplot2::geom_sf(data = acapulco_ageb$ageb,
                   fill = "gray80",
                   col = "white",
                   lwd = 0.1) +
  ggplot2::theme_void() +
  ggplot2::geom_sf(data = acapulco_ageb$locality,
                   fill = NA,
                   col = "black",
                   lwd = .3)