Statistical features from NDVI mosaics

Statistics NDVI mosaics

These mosaics are produced by further processing S2ind-mosaics [https://ckan.ymparisto.fi/dataset/sentinel-2-image-index-mosaics-s2ind-sentinel-2-kuvamosaiikit-s2ind](https://ckan.ymparisto.fi/dataset/sentinel-2-image-index-mosaics-s2ind-sentinel-2-kuvamosaiikit-s2ind) for the full months. The following mosaics are produced:

* Yearly maximum value
* Yearly minimum value
* Yearly mean
* Yearly median
* Yearly 10-quantile
* Yearly 25-quantile
* Yearly amplitude (pixelwise maximum - 25-quantile from current year and two previous years)
* Yearly sum
* Monthly mosaics for April, May, June, July, August, September and October

These mosaics are not produced at once, but S2ind mosaics are first tiled into a 12x8 grid, which are again combined after the processing. During processing, all nodata values are ignored. The processing steps are the following:

1. Snowy areas are masked from the data using NDSI threshold of 140 (approximately 0.4)
2. Two base mosaics are constructed, one for spring and one for autumn. The values for these mosaics are median values of all available data (2016 to 2023). Spring mosaic is constructed from dates between 1.4. and 31.5., while the autumn mosaic is constructed from dates between 15.9. and 31.10.
3. All gaps in mosaics are filled with the maximum value of the same month and pixel from two previous years
4. All remaining nodata for April and May mosaics are filled with corresponding values from the base spring mosaic
5. All remaining nodata for October are filled with corresponding values from the base autumn mosaic
6. All remaining nodata for May, June, July, August and September are filled with the mean value of the adjacent months of the same year
7. Yearly statistics are collated
8. All sea, river and lake areas are masked from the data
9. All produced mosaics are clipped to Finnish borders

Because amplitude requires two previous years to derive, it is available only from 2020 onwards. Nodatavalues are 0 for other statistic aside from sum and amplitude, which have Nodatavalue of -999.

The mosaics are availabe from the interface of NSDC using URL, like https://pta.data.lit.fmi.fi/sen2/phenology/s2_phenology_statistics_year_feature.tif, where year is one of [2018, 2019, 2020, 2021, 2022, 2023] and feature is [amp, max, mean, median, min, quantile_10, quantile_25, sum].

The pixel size for the dataset is 10m. Datatype is unsigned 8 bit integer for all mosaics except yearly sum and amplitude, whose datatype is unsigned 16 bit integer. These datasets can be used according to open data license (CC BY 4.0)

Dataset updates yearly.

---

Vuosittaiset Sentinel-2 indeksimosaiikit

Mosaiikit lasketaan S2ind-mosaiikeista [https://ckan.ymparisto.fi/dataset/sentinel-2-image-index-mosaics-s2ind-sentinel-2-kuvamosaiikit-s2ind](https://ckan.ymparisto.fi/dataset/sentinel-2-image-index-mosaics-s2ind-sentinel-2-kuvamosaiikit-s2ind) täysiltä kuukausilta. Aineistoista tuotetaan seuraavat mosaiikit:

* Vuosittainen maksimi
* Vuosittainen minimi
* Vuosittainen keskiarvo
* Vuosittainen mediaani
* Vuosittainen 10-kvantiili
* Vuosittainen 25-kvantiili
* Vuosittainen amplitudi (Vuosittainen maksimi - vuoden ja kahden edellisen vuoden 25-kvantiili)
* Vuosittainen summa
* Kuukausittaiset mosaiikit huhti-, touko-, kesä-, heinä-, elo-, syys- ja lokakuulle

Kokonaisia mosaiikkeja ei prosessoida kerralla, vaan S2ind-mosaiikit jaetaan 12x8 hilaan, jotka yhdistetään toisiinsa lopuksi. Prosessoinnin aikana nodata-arvoja ei huomioida tunnuslukujen laskennassa. Prosessoinnissa suoritetaan seuraavat vaiheet:

1. Aineistoista poistetaan lumiset alueet käyttäen vastaavaa NDSI-mosaiikkia ja kynnysarvoa 140 (noin 0.4)
2. Luodaan kaksi perusmosaiikkia, toinen keväälle ja toinen syksylle. Näiden mosaiikkien pikseliarvot ovat mediaaniarvot kaikesta käytettävissä olevasta aineistosta (2016-2023). Kevätmosaiikki kootaan aikavälin 1.4. ja 31.5. kuvista, ja syksymosaiikki aikavälin 15.9. ja 31.10. kuvista.
3. Kaikki aukot mosaiikeissa täytetään kahden edellisen vuoden saman pikselin ja kuukauden maksimiarvolla
4. Kaikki jäljelläoleva nodata huhti- ja toukokuulle täytetään kevätmosaiikin vastaavien pikselien arvoilla
5. Kaikki jäljelläoleva nodata lokakuulle täytetään syksymosaiikin vastaavien pikselien arvoilla
6. Kaikki jäljelläoleva nodata touko-, kesä-, heinä-, elo-, ja syyskuulle täytetään saman vuoden vierekkäisten kuukausien vastaavien pikselien keskiarvoilla
7. Lasketaan tilastolliset tunnusluvut
8. Vesi-, meri ja järvialueet leikataan pois aineistosta
9. Kaikki tuotetut mosaiikit leikataan Suomen rajohin

Koska amplitudille vaaditaan kaksi aiempaa vuotta, se on saatavilla vasta vuodesta 2020 eteenpäin. Aineistojen nodata-arvo on nolla lukuunottamatta summa- ja amplitudimosaiikkeja, joille se on -999.

Aineisto on saatavilla Kansallisen satelliittidatakeskuksen käyttöliittymästä URL-osoitteen avulla, esim. https://pta.data.lit.fmi.fi/sen2/phenology/s2_phenology_statistics_vuosi_piirre.tif, missä vuosi on yksi joukosta [2018, 2019, 2020, 2021, 2022, 2023] ja piirre on yksi joukosta [amp, max, mean, median, min, quantile_10, quantile_25, sum].

Aineiston pikselikoko on 10m. Kaikkien mosaiikkien datatyyppi on 8-bittinen etumerkitön kokonaisluku lukuunottamatta vuosittaisia summa- ja amplitudimosaiikkeja, joiden datatyyppi on 16-bittinen etumerkitön kokonaisluku. 

Aineisto päivittyy vuosittain.

Description of the index computation process is described in S2ind metadata.

Organization: Suomen ympäristökeskus (Syke)