require(raster)
## Cargando paquete requerido: raster
## Cargando paquete requerido: sp
require(terra)
## Cargando paquete requerido: terra
## terra 1.8.21
require(sp)
comunas=shapefile("C:/Users/Claudia/Downloads/Comuna.shp")
comunas
## class       : SpatialPolygonsDataFrame 
## features    : 22 
## extent      : -76.59284, -76.46125, 3.331802, 3.505871  (xmin, xmax, ymin, ymax)
## crs         : +proj=longlat +datum=WGS84 +no_defs 
## variables   : 7
## names       : POLY_ID, casos, comuna,   nombre, zona_recol,         area, perimetro 
## min values  :       1,     1,      1, Comuna 1,         NA,  2329397.941,  7983.949 
## max values  :      22,    39,     22, Comuna 9,         NA, 12555929.024, 26480.361
plot(comunas)

tabla=comunas@data

comuna22=comunas[22,]

plot(comuna22)

plot(comunas)
plot(comuna22,add=TRUE,col="red")