Masks

Collection of different masks used in my simulations.

Forest Masks

my-picture1

Forest masks were made by Matthieu Baron from the forest dataset of IGN BD Foret V2. The code and a note [2] to reproduce the masks can be found in his git depot [3].

To summarize, the BD foret dataset present polygones of vegetation cover of different types. Different vegetation types are selected according to characteristics choosen by Matthieu B. In the current version, the vegetation formations called “landes” and “formation herbacées” are removed, all the others are used to build the mask. The different polygones are then rasterized to the simulation grid.

RGI Glacier Masks

my-picture1

The glaciers mask is simply from a rasterisation of the RGI glacier inventory to the simulation grid.

Lakes Rivers Cities Masks

my-picture1

Lake, rivers and cities masks are hand picked from IGN shapefile BD TOPO. Lakes and rivers are taken form the hydrographic shapefiles, cities from the “commune” shapefile.

Geomorphons Masks

Geomorphons detection are methods for classification and mapping of landform elements from a DEM based on the principle of pattern recognition rather than differential geometry. [1]

I used a dedicated software to perform this analysis, Whitebox Tools, producing the pixel classification . The geomorphons concept is based on line-of-sight analysis for the eight topographic profiles in the cardinal directions surrounding each grid cell in the input DEM.

my-picture1

I used the command line interface of the whitebox_tools software (see website to install) with the -r=Geomorphons option. This command take as input a DEM (need to be on the same grid as the simulation), the DEM zone must be larger than the zone of interest because we are using a search radius to detect the different landform feature. I chose a search radius of 2.5 km (100 pixel at 250m resolution so option search=100). The larger the search radius is, the less dependant to resolution the classification is. I also choose a flat threshold of 1 degree with option threshold=1. Finally the output is presented for the 10 most common landforms as in [1] with the option forms.

Example of command:

./whitebox_tools -r=Geomorphons -v --wd=. --dem=DEM_FRANCE_L93_250m_bilinear.tif -o=output_100_r.tif --residuals --search=100 --threshold=1 --forms

Then to convert the geomorphon tif to netcdf (to be used as input by snowtools), you can use gdal_translate

gdal_translate -of NetCDF <input_filename> <output_filename>

Footnotes