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: