The code enables users to visualize spatial data in the cryosphere in 3D animations, including a model of the glacier as well as graphical data visualization. The code is based on the original Rayshader (https://github.com/tylermorganwall/rayshader). Necessary input data is a DEM of the area, spatial data of the glacier surface changing in time (e.g. data retrieved from satellites or models) and data to visualize on top including glacier velocities, surface features like lakes etc.
Note that distributed glacier thickness data exists for all glaciers globally (consult: https://www.glims.org/RGI/rgi60_dl.html). Velocity datasets also become more and more accesible (see for example https://nsidc.org/data/golive/)
Prerequisites are a working installation of R Studio (https://rstudio.com/) as well as working packages as detailed below.
The code input_cryoRayshader.R needs to be edited and run by the user, cryoRayshader.R includes the model core.
We use data from our research studies to show examples of how the model can be used below.
Please get in touch via firstname.lastname@example.org if you have other kinds of datasets you'd like to see work in this or there are bugs we can resolve.
Visualize modelled mass loss over a clean ice and a debris-covered glacier
Model outputs from Wijngaard et al. (2019) (data available on request) allows us to visualize change in space. We do that below for a clean ice (Hintereis, European Alps) and debris-covered glacier (Langtang, Himalaya).
Simple mass loss data is visualized with orthophotos draped on top and the extent of the respective ice cover enhanced.
Visualize modelled mass loss and associated data
In combination with the 3D visualization the data can be visualized alongside as a graph.
Additionally spatial data like mass loss or velocities can be draped over the glacier surface.
Visualize changing pond cover on a debris-covered glacier
as well as the total and average area covered by the features ...