7.3 Lectura de shapefile
## leer shape del pais
ShapeSAE <- read_sf("01 Modelo de area/CHL/2017/ShapeDeptoCHL/comunas.shp")
ShapeSAE %<>% mutate(dam2 = str_pad(
string = cod_comuna,
width = 5,
pad = "0"
),
dam = str_pad(
string = codregion,
width = 2,
pad = "0"
)
) %>%
filter(dam != "00", !dam2 %in% c("05201", "13603"))