Cómo conseguir un mapa al estilo “Joy Division” con QGIS y R (II)

En esta ocasión voy a realizar una visualización al estilo «Joy Division» pero con los programas QGIS y R, totalmente gratuitos. Para ello vamos a seguir un flujo de trabajo de diferentes herramientas, voy a trabajar con datos de Tenerife, el MDT descargado de esta web. Lo voy a hacer rápidamente porque tenéis la información completa en el enlace que hay más abajo.

Así voy a crear un grid de tipo líneas , este sería el resultado:

 

Generaré los puntos a lo largo del grid:


 

También puedo realizar un clip si me interesa más sólo tener la isla:

 

Se utiliza a continuación la herramienta «Point Sampling Tool» si no la tienes en QGIS debes instalarla como complemento.

 

 

Y por último trabajamos con la herramienta «Add coordinates to point» de SAGA, y obtenemos un shape que posteriormente pasaremos a un formato tabla (.CSV) que será con el que trabajaremos en R.

 


 

Ya en R incluiremos este código obtenido de este tutorial 

 

library(ggplot2)
library(ggridges)
library(mapproj)
tenerife <- read.csv(file = "tenerife.csv", header = TRUE, sep = ",")
head(tenerife)
names(tenerife)[2:4] <- c("Elev", "Lon", "Lat")

head(tenerife)
tenerife_uno <- ggplot(tenerife, aes(x = Lon, y = Lat, group = Lat, height = Elev)) +
  geom_density_ridges(stat = "identity", scale=30)
tenerife_uno

tenerife_joy<- ggplot(tenerife, aes(x = Lon, y = Lat, group = Lat, height = Elev)) +
  geom_density_ridges(stat = "identity", scale=30,
                      fill = "transparent", color = "orange") +
  
  scale_x_continuous(name = "TENERIFE") +
  theme(panel.grid.major = element_blank(),
        panel.grid.minor = element_blank(),
        panel.background = element_rect(fill = "black"),
        axis.line = element_blank(),
        axis.text.x = element_blank(),
        plot.background = element_rect(fill = "black"),
        axis.ticks.x = element_blank(),
        axis.title.y = element_blank(),
        axis.text.y = element_blank(),
        axis.ticks.y = element_blank(),
        axis.title.x = element_text(colour = "white", size = 18)) +
  coord_map()
tenerife_joy
ggsave("tenerife2.png", dpi = 300)

Y este será el resultado, ojo también puedes hacer un clip y quedarte sólo con el área de Tenerife, os animo a que juguéis con el código:

Deja una respuesta

Tu dirección de correo electrónico no será publicada. Los campos obligatorios están marcados con *