Hice un gráfico de puntos utilizando datos de longitud y latitud con el siguiente código:
##Archivo csv
gto<- read_csv(sprintf("http://fs.planeacion.sep.gob.mx/cct/cct11.csv"))
##Vector de niveles educativos (para ordenar)
nivel_ed <- c("Primaria", "Secundaria", "Bachillerato", "Superior")
gtomap <- subset(gto, select=c(longitud, latitud, nnivel)) %>%
filter(nnivel%in%c("PRIMARIA", "SECUNDARIA", "BACHILLERATO",
"PEDAGOGICA, UNIVERSITARIA O TECNOLOGICA", "PROFESIONAL")) %>%
mutate(nnivel= case_when(nnivel== "PEDAGOGICA, UNIVERSITARIA O TECNOLOGICA" ~ "Superior",
nnivel=="PROFESIONAL" ~ "Superior",
nnivel=="SECUNDARIA" ~ "Secundaria",
nnivel=="PRIMARIA" ~ "Primaria",
nnivel=="BACHILLERATO" ~ "Bachillerato")) %>%
mutate(nnivel= factor(nnivel, levels=nivel_ed)) %>%
melt(id =c("longitud","latitud")) %>%
mutate(value= factor(value, levels=nivel_ed)) %>%
ggplot(aes(x = longitud, y = latitud, col= value)) +
geom_point() +
theme(legend.title=element_blank()) +
labs(x = "", y="")
gtomap
Posteriormente, descargué un mapa utilizando:
##Para importar geom_sf
devtools::install_github("tidyverse/ggplot2")
##Mapa
mapgto <- sf::read_sf("http://geoinfo.iplaneg.net/geoserver/wms/kml?layers=geonode%3Amgm_gto2010&mode=download", quiet=T)
map <- data.frame(mapgto)
##El mapa se grafica directamente con
ggplot() +
geom_sf(data=map)
Como podrán observar, los puntos y el mapa corresponden a las mismas coordenadas. Por tanto, quiero superponer el primer gráfico al segundo a fin de que el gráfico final contenga el mapa y los puntos que contienen la ubicación de mis variables de interés. Sin embargo, si intento algo como:
gtomap + geom_sf(data=map)
Error in FUN(X[[i]], ...) : objeto 'longitud' no encontrado
Se aprecia un error. He leído diversos tutoriales, pero no he encontrado alguno que me ofrezca una alternativa para mi problema específico. Agradezco cualquier orientación y comentario.
map
que sale después dedata.frame(mapgto)
. Es un data.frame anidado, tiene listas dentro con las coordenadas. Seguramentegeom_sf
tiene un método interno para desanidarlas y hacer el mapa. (No lo probé para no instalar una versión en desarrollo de ggplot). Pero al no asignar conaes()
las variables/columnas no coinciden con las degtomap
, por eso ese error. Sugerencia: para tu cartografía trata de usar un data.frame "puro" con columnaslatitud
ylongitud
a las que asignes ay
yx
. En ese caso te va a coincidir sin problemas con el otro df.tidymaps
para hacernos la vida más fácil. Si tengo tiempo en la tarde busco una respuesta propiamente dicha a tu pregunta.