Plotear datos raster

Plotear datos usando usando “breaks”

Para mayor claridad y visibilidad del plot, podemos dividir nuestros datos. Es util primero conocer la distribicion de los datos mediante un grafico de barras utilizando las funciones mutate() y cut() combinadas. Para esto utilizaremos la libreria dplyr Los datos que utilizaremos datos de elevacion de un modelo digital del terreno.

summary(df_dem$dem)

   Min. 1st Qu.  Median    Mean 3rd Qu.    Max.    NA's 
    100    2509    3000    2887    3689    4500   25701 

Para conocer los valores de de corte de cada grupo. Un grupo se van encontrar todos los valores NA’s

unique(df_dem $ elevation)

[1] <NA>                (95.6,1.57e+03]     (1.57e+03,3.03e+03] (3.03e+03,4.5e+03] 
Levels: (95.6,1.57e+03] (1.57e+03,3.03e+03] (3.03e+03,4.5e+03]

Para contar los valores por cada grupo

df_dem%>%
  group_by(df_dem$elevation)%>%
  summarize(counts= n())

Podemos personalizar los valores de corte para cada grupo de datos. Como ehemplo vamos a realizar 3 cortes entre los valores menores a 1000, 1001-2000 y mayores a 2000

custom_bins <- c(100, 2000, 3000,4500 )
df_dem_cut <-   df_dem %>%
                mutate(elevation_2 = cut(df_dem$dem, breaks = custom_bins))
ggplot() +
    geom_bar(data = df_dem_cut, aes(elevation_2))

unique(df_dem_cut$elevation_2)

[1] <NA>            (100,2e+03]     (2e+03,3e+03]   (3e+03,4.5e+03]
Levels: (100,2e+03] (2e+03,3e+03] (3e+03,4.5e+03]

Podemos plotear el raster en donde cada grupo de datos tenga un color diferente

ggplot() +
 geom_raster(data =df_dem_cut  , aes(x = x, y = y, fill = elevation_2))

ggplot() +
 geom_raster(data =df_dem_cut , aes(x = x, y = y, fill = elevation_2)) + 
    scale_fill_manual(values = terrain.colors(3))

Si necesitamos crear múltiples gráficos utilizando la misma paleta de colores, podemos crear un objeto R (my_col) para el conjunto de colores que queremos usar.

y_col <- terrain.colors(3)
ggplot() +
 geom_raster(data =df_dem_cut , aes(x = x, y = y, fill = elevation_2)) + 
    scale_fill_manual(values = y_col)

Podemos crear un efecto de sombreado tridimensional. Vamos a generar un hillshade el cual mapea las sobras y las texturas que tiene el terreno. Utilizaremos una escala de gris

ggplot() +
 geom_raster(data = df_dem , aes(x = x, y = y, alpha = dem)) + 
    scale_alpha(range =  c(0.15, 0.90), guide = "none")