Tired of waiting the perfect image with no cloud to show the snow cover? This monthly snow report script is for you.
This code will find where the snow is persistent for the last 30 days (from the chosen date). In order to well represent the land-cover, the script will store each pertinent date in a list and will represent the median value.
As the aim of the script is to represent the snow cover, all the other land-cover are saturated in green in order to easily see the snow.
The limitations are essentially:
- a few white rooftops
- when no uncloudy pixel is available in the previous 30 days
Thus, it is possible to make the script working with 90 days in order to have a trimonthly synthesis by changing settings to: “numberOfMonthsToUse = 3”.
Persistent snow cover in February 2019 in Corsica
Persistent snow cover in Pyrénées in March 2019
Mer de Glace / Alpes persistent snow cover in the previous 30 days of 24 april 2019.
See the supplementary material for more examples.
- Harel Dan temporal script
- NDSI/red thresold was adapted from here
- Snow representation by Simon Gascoin