Creo el mapa con los polígonos que se usan en la EPH

Preprocesamiento
Dejo esto para que cualquiera pueda reproducirlo
#cargo los datos de la EPH resumen
load("base_ep_resumen.RData")

#cargo los aglomerados de la EPH
radios.eph <- read_sf("radios_eph.json") 
radios.eph <- st_transform(radios.eph,'+proj=longlat +datum=WGS84') #4326) #pongo en codigo leaflet

#los junto
radios.eph$AGLOMERADO <- toupper(radios.eph$eph_aglome)

#cree objectos con centroides
radios.eph$centroid <- st_centroid(radios.eph$geometry)
radios.eph <- cbind(radios.eph,st_coordinates(radios.eph$centroid ))
#reome thingswithout location
radios.eph <- radios.eph[!is.na(radios.eph$X),]
radios.eph <- distinct(radios.eph)
radios.eph$aglomerado_y_loc <- paste0(radios.eph$AGLOMERADO,"_",radios.eph$localidade)

#combinar geometrías por aglomerado 
radios.eph.aglom <- radios.eph %>%
  group_by(AGLOMERADO) %>% 
  summarize(geometry = st_union(geometry))

unique(individual_03.hoy$AGLOMERADO[individual_03.hoy$AGLOMERADO %!in% radios.eph.aglom$AGLOMERADO ])

grep("JU",unique(radios.eph.aglom$AGLOMERADO),value = T)

radios.eph.aglom$AGLOMERADO[radios.eph.aglom$AGLOMERADO=="PARTIDOS DEL GBA" ] <- "PARTIDOS GBA"
radios.eph.aglom$AGLOMERADO[radios.eph.aglom$AGLOMERADO=="GRAN PARANÁ"] <- "GRAN PARANA"
radios.eph.aglom$AGLOMERADO[radios.eph.aglom$AGLOMERADO=="GRAN CÓRDOBA"] <- "GRAN CORDOBA"
radios.eph.aglom$AGLOMERADO[radios.eph.aglom$AGLOMERADO=="GRAN TUCUMÁN - TAFI VIEJO"] <- "GRAN TUCUMAN - TAFI V."
radios.eph.aglom$AGLOMERADO[radios.eph.aglom$AGLOMERADO=="SAN LUIS - EL CHORRILLO"] <- "GRAN SAN LUIS"
radios.eph.aglom$AGLOMERADO[radios.eph.aglom$AGLOMERADO=="SALTA"] <- "GRAN SALTA"
radios.eph.aglom$AGLOMERADO[radios.eph.aglom$AGLOMERADO=="MAR DEL PLATA - BATÁN" ] <- "MAR DEL PLATA"
radios.eph.aglom$AGLOMERADO[radios.eph.aglom$AGLOMERADO=="SANTIAGO DEL ESTERO - LA BANDA" ] <-"S. DEL ESTERO - LA BANDA"
radios.eph.aglom$AGLOMERADO[radios.eph.aglom$AGLOMERADO== "SAN NICOLAS - VILLA CONSTITUCIÓN" ] <- "SAN NICOLAS - V. CONSTITUCION"
radios.eph.aglom$AGLOMERADO[radios.eph.aglom$AGLOMERADO== "SAN NICOLAS - VILLA CONSTITIUCIÓN" ] <- "SAN NICOLAS - V. CONSTITUCION"
individual_03.hoy$AGLOMERADO[individual_03.hoy$AGLOMERADO== "USUAHIA - RIO GRANDE" ] <- "USHUAIA - RIO GRANDE" #error en eph 
radios.eph.aglom$AGLOMERADO[radios.eph.aglom$AGLOMERADO== "VIEDMA - CARMEN DE PATAGONES" ] <- "VIEDMA - C. DE PATAGONES"
radios.eph.aglom$AGLOMERADO[radios.eph.aglom$AGLOMERADO=="COMODORO RIVADAVIA - RADA TILLY"  ] <- "C. RIVADAVIA - R. TILLY"
radios.eph.aglom$AGLOMERADO[radios.eph.aglom$AGLOMERADO=="NEUQUÉN - PLOTTIER"  ] <-"NEUQUEN - PLOTTIER"
radios.eph.aglom$AGLOMERADO[radios.eph.aglom$AGLOMERADO=="JUJUY - PALPALÁ"  ] <- "JUJUY - PALPALA"

#Algunos datos no tienen aglomerado... los voy a dejar por fuer por ahora 
individual_03.hoy[individual_03.hoy$AGLOMERADO == "N/A",]

save.image("mapa.RData")

radios.eph.aglom.map <- left_join(individual_03.hoy,radios.eph.aglom)
radios.eph.aglom.map <- st_as_sf(radios.eph.aglom.map,sf_column_name ="geometry",crs='+proj=longlat +datum=WGS84')

radios.eph.aglom.map <- radios.eph.aglom.map[radios.eph.aglom.map$AGLOMERADO != "N/A",]

radios.eph.aglom.map <- st_transform(radios.eph.aglom.map,'+proj=longlat +datum=WGS84')
mapa_versión1

2023-07-11

Creo el mapa con los polígonos que se usan en la EPH

Así se ve el mapa por ahora
