library(tidyverse)
library(dplyr)
library(ggplot2)
library(sf)
library(leaflet)
library(sp)
library(rgdal)
mun_antioquia <- st_read("../GEOMATICA/05_ANTIOQUIA/ADMINISTRATIVO/MGN_MPIO_POLITICO.shp")
## Reading layer `MGN_MPIO_POLITICO' from data source `C:\Users\User\Desktop\GEOMATICA\05_ANTIOQUIA\ADMINISTRATIVO\MGN_MPIO_POLITICO.shp' using driver `ESRI Shapefile'
## Simple feature collection with 125 features and 9 fields
## Geometry type: POLYGON
## Dimension:     XY
## Bounding box:  xmin: -77.12783 ymin: 5.418558 xmax: -73.88128 ymax: 8.873974
## Geodetic CRS:  WGS 84
mun_antioquia$KM2 <- st_area(st_transform(mun_antioquia, 3116))/1E6
mapa