¿Cómo superponer las coordenadas de los puntos en el shapefile de tmap?

Question

Hasta ahora, he conseguido presentar visualmente los datos de diferentes factores socioeconómicos a nivel de distrito en Pakistán con la ayuda del tmap paquete en RStudio para crear la siguiente salida.

Al mismo tiempo, he utilizado el ggmap para trazar las coordenadas geográficas de varias torres para diferentes operadores de redes móviles en Pakistán (aproximadamente 13000 torres). El resultado es el siguiente.

Mi siguiente requisito es superponer el segundo gráfico sobre el primero, mostrando esencialmente las torres sobre las puntuaciones de electricidad por distrito.

¿Alguna idea sobre qué paquetes podría utilizar (así como código de ejemplo) para conseguirlo?

Answer 1

Esto también se puede hacer con tmap . Echa un vistazo a la tmap viñeta especialmente el Sección de ploteo con elementos tmap .

En tu caso, tu código debería ser básicamente así:

tm_shape(shapefile_with_electricity_data) + 
  tm_fill("Electricity_in_Schools_variable") + 
tm_shape(shapefile_with_MNO_data) + 
  tm_dots("MNO_variable")

Answer 2

Puede probar el st_as_sf de la función sf paquete.

Lo utilicé en uno de mis proyectos:

boston <- st_read("../data/boston.geojson")
proj4 <- st_crs(boston)$proj4string

student_points <- student_info %>% 
   st_as_sf(coords = c('lng', 'lat'), crs = proj4)

tm_shape(boston) +
  tm_polygons('#f0f0f0f0', border.alpha = 0.2) +
tm_shape(student_points_rastered) +
  tm_raster(alpha = 0.7, title = '# of students')