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")
