Harmonized Landsat satellite image mosaic timeseries
The purpose of this data is to provide harmonized dataset of yearly satellite image mosaics based on Landsat imagery and starting from year 1984. Both yearly reflectance mosaics of visible, near and shortwave infrared bands, as well as several image index mosaics like Normalized Difference Vegetation Index are provided. So far, years 1984 - 2024 are processed; years 1984 – 2014 are based on Landsat-5 Thematic Mapper (launched 1.3.1984, decommissioned 5.6.2013) and Landsat-7 Enhanced Thematic Mapper Plus (launched 15.4.1999, scan line corrector failure 31.5.2003, imaging suspended 19.1.2024), and years 2013 – 2024 based on Landsat-8 (launched 11.2.2013) and Landsat-9 (launched 27.9.2021) Operational Land Imager imagery.
The USGS provides several differently processed Landsat datasets; Landsat Collection-2 Level-2 imagery has been used in this timeseries. These images provide surface reflectance estimates combined with the best globally available geometric orthocorrection for Landsat imagery. Mosaics have been done using visible B1-B4, near B5 and shortwave infrared B6 & B7 bands (numbering of bands is based Landsat-9 Operational Land Imager bands, wavelengths in micrometers):
- B1: 0.433-0.453
- B2: 0.450-0.515
- B3: 0.525-0.600
- B4: 0.630-0.680
- B5: 0.845-0.885
- B6: 1.560-1.660
- B7: 2.100-2.300
The number of bands and their numbering varies between instruments; Landsat-8 OLI and -9 OLI-2 has additional blue band B1 that Landsat-7 ETM+ and Landsat-5 TM do not have. Also, the widths of bands are a bit different in different instruments, the details of band numbering and widths can be found from following WWW-pages:
Individual images have been downloaded to National Satellite Data Centre (NSDC: https://nsdc.fmi.fi/). The mosaicking was performed using CalFin-processing cluster at NSDC. The cloud masks provided by images (based on CFmask v3.3.1) and some other checks were used to select valid values for mosaicking. In case of Landsat-5 and -7, pixel has been selected for processing if it was outside masks cloud, cloud shadow, fill and cloud and cloud shadow confidences are low, distance to cloud was greater than 0.6 km, NDVI was between -1 and 1, and surface reflectance of bands more than 0.005. In case of Landsat-8 and -9 there were some small changes; pixel was selected for processing if it was outside masks cloud, cloud_shadow, cirrus, dillated_cloud and designated_fill, if distance to cloud was greater than 1 km and surface reflectance of bands more than 0.005.
The time window used to select observations for mosaicking varied in different parts of Finland, because growing season length is different. The time windows and areas were:
- South (lat. 59-64 deg.): 1.5.-30.9.
- Central (lat. 64-68 deg.): 15.5.-15.9.
- North (lat. 68-70 deg.): 1.6.-31.8.
The mosaicking was based on 95% percentile of NDVI-values. This was used as an attempt to avoid shadow areas in final mosaic.
Postprocessing GDAL/Python-script merged parts together and made coordinate transformation to TM35Fin-coordinate system (EPSG 3067) with 20 meter pixel size. Area of mosaics is (TM35Fin, EPSG 3067):
- ULE: 50000
- ULN: 7800000
- LRE: 842000
- LRN: 6600000
Following image indices have been computed:
- Normalized Difference Moisture Index NDMI = ( TM_B4 – TM_B5 ) / ( TM_B4 + TM_B5 ) = ( OLI_B5 – OLI_B6 ) / ( OLI_B5 + OLI_B6 )
- Normalized Difference Tillage Index NDTI = ( TM_B5 – TM_B7 ) / ( TM_B5 + TM_B7 ) = ( OLI_B6 – OLI_B7 ) / ( OLI_B6 + OLI_B7 )
- Normalized Difference Vegetation Index NDVI = ( TM_B4 – TM_B3 ) / ( TM_B4 + TM_B3 ) = ( OLI_B5 – OLI_B4 ) / ( OLI_B5 + OLI_B4 )
- Normalized Difference Water Index NDWI = ( TM_B2 – TM_B4 ) / ( TM_B2 + TM_B4 ) = ( OLI_B3 – OLI_B5 ) / ( OLI_B3 + OLI_B5 )
- TM_Bx: Landsat-5 TM and -7 ETM+ band x
- OLI_Bx: Landsat-8 OLI and -9 OLI-2 band x
Other indices can be computed by downloading individual bands and making the computations by him/herself.
Other published imagery consists of
- META: The pixel value is the number of day of selected observation.
- NUMOBS: The pixel value is the number of available valid observations. NOTE: The common area of succeeding images from same day is counted twice.
The mosaics are available by downloading from the archive of the National Satellite Data Centre (https://nsdc.fmi.fi/):
Reflectance band mosaics of Landsat-5 and -7: https://pta.data.lit.fmi.fi/lans/landsat57/bands/landsat57_refl_bn_year.tif
- bn: b1, b2, b3, b4, b5 or b7
- year: between 1984 - 2014
Index mosaics of Landsat-5 and -7: https://pta.data.lit.fmi.fi/lans/landsat57/index/landsat57_ind_year.tif
- ind: ndmi, ndti, ndvi, ndwi
- year: between 1984 - 2014
Metadata mosaics of Landsat-5 and -7:
Reflectance band mosaics of Landsat-8 and -9: https://pta.data.lit.fmi.fi/lans/landsat89/bands/landsat89_refl_bn_year.tif
- bn: b1, b2, b3, b4, b5, b6 or b7
- year: between 2013 - 2024
Index mosaics of Landsat-8 and -9: https://pta.data.lit.fmi.fi/lans/landsat89/bands/landsat89_ind_year.tif
- ind: ndmi, ndti, ndvi, ndwi
- year: between 2013 - 2024
Metadata mosaics of Landsat-8 and -9:
The drawbacks of this timeseries are that there are gaps of data due to lack of images, and there are atmospheric effects still visible in made mosaics like clouds that the cloud masking has not been detected, or cloud shadow areas. It was noticed that the NDVI of cloud shadow areas can be higher than neighboring shadowless areas, leading to situations that shadow areas are selected instead of shadowless areas into final mosaic. It seems that the atmospheric correction corrects visible bands, especially red, too much leading to too low reflectance of red band and therefore increased NDVI. Also, water areas can have negative reflectance, and these have been treated as outliers and removed from mosaicking process.
Reference: Landsat Collection 2 Level-2 Science Products, https://www.usgs.gov/landsat-missions/landsat-collection-2-level-2-science-products
Original Landsat satellite images by USGS, processing by SYKE. This SYKE’s dataset can be used according to open data license (CC BY 4.0).
Suomeksi:
Harmonisoitu Landsat satelliittikuvamosaiikkiaikasarja
Tämän datasetti käsittää Landsat-satelliittisarjan kuviin perustuvia vuosittaisia kuvamosaiikkeja alkaen vuodesta 1984. Vuosittainen data koostuu näkyvän valon, lähi- ja keski-infrapuna-alueen kanavien maanpinnan reflektanssin mosaiikeista, sekä muutamasta kuvaindeksimosaiikista kuten NDVI (Normalized Difference Vegetation Index). Tällä hetkellä (24.10.2024), vuodet 1984-2024 on prosessoitu; vuosien 1984-2014 mosaiikit perustuvat Landsat-5 Thematic Mapper (laukaistu kiertoradalle 1.3.1984, käyttö lopetettu 5.6.2013) ja Landsat-7 Thematic Mapper Plus (laukaistu kiertoradalle 15.4.1999, keilauslinjan korjaaja rikkoutui 31.5.2003, käyttö lopetettu 19.1.2024) kuviin, ja vuosien 2013-2024 mosaiikit perustuvat Landsat-8 Operational Land Imager (laukaistu 11.2.2013) ja Landsat-9 Operational Land Imager-2 (laukaistu 27.9.2021) kuviin.
Yhdysvaltojen Geologian tutkimuslaitos USGS hallinnoi Landsat-ohjelmaa ja hoitaa kuvien jakelun. He tuottavat eri lailla prosessoituja datasettejä, joista Landsat Collection-2 Level-2 datasettiä on käytetty tässä aikasarjassa. Näiden kuvien pikselien arvot ovat maanpinnan reflektanssin estimaatteja, eli kuville on tehty ilmakehäkorjaus. Kuville on tehty orto-oikaisu, eli maaston pinnanmuotojen vaikutus pikselien paikkaan on huomioitu. Kuvamosaiikit on tehty näkyvän valon B1-B4, lähi-infrapuna- B5 ja keski-infrapunakanaville B6 ja B7 (numerointi perustuu Landsat-9 Operational Land Imager-2 instrumentin kanavien numerointiin, aallonpituudet on ilmaistu mikrometreinä)
- B1: 0.433-0.453
- B2: 0.450-0.515
- B3: 0.525-0.600
- B4: 0.630-0.680
- B5: 0.845-0.885
- B6: 1.560-1.660
- B7: 2.100-2.300
Kanavien lukumäärä ja numerointi vaihtelee instrumenttien välillä; Landsat-8 OLI ja -9 OLI-2 instrumenteissa on sinisen aallonpituuden kanava B1 jota ei ole Landsat-7 ETM+ tai Landsat-5 TM-instrumenteissa. Lisäksi, kanavien aallonpituudet eri instrumenteissa vaihtelevat hieman. Lisätietoja eri instrumenttien kanavien aallonpituuksista ja numeroinnista löytyy seuraavilta WWW-sivuilta:
Yksittäiset Landsat-kuvat on ladattu Ilmatieteen laitoksen hallinnoiman Kansallisen Satelliittidatakeskuksen (NSDC) arkistoon. Varsinainen mosaiikkien tuotanto tehdään CalFin-prosessointiklusterilla. Pilvimaskaukseen käytetään kuvien mukana tulevia pilvimaskeja (menetelmä CFmask v3.3.1) ja näiden lisäksi käyttökelpoisten havaintojen valintaan käytetään myös muita testejä. Landsat-5 ja -7 kuvien tapauksessa pikseli valitaan mosaikointiin, mikäli se ei kuulu maskien cloud, cloud shadow ja fill alueeseen, pilven ja pilvivarjon todennäköisyys on alhainen, etäisyys pilven reunaan on yli 0.6 km, NDVI on vöälillä -1 ja 1, sekä kanavien reflektanssit ovat suurempia kuin 0.005. Landsat-8 ja -9 kuvien tapauksessa säännöissä on hieman muutoksia. Pikseli valitaan mosaikointiin mikäli se on maskien cloud, cloud_shadow, cirrus, dillated_cloud ja designated_fill alueeseen, etäisyys pilven reunukseen on suurempi kuin 1 km ja kanavien reflektanssit ovat suurempia kuin 0.005.
Mosaikoinnissa käytettävä kuvien aikaikkuna (ts. miltä aikaväliltä voidaan valita kuvia mosaikointiin) vaihteli etelä-pohjoissuunnassa, koska kasvukauden pituudet ovat erilaisia Suomen eri osissa. Aikaikkunat ja näiden alueet olivat:
- Etelä-Suomi (lat. 59-64 deg.): 1.5.-30.9.
- Keski-Suomi (lat. 64-68 deg.): 15.5.-15.9.
- Pohjois-Suomi (lat. 68-70 deg.): 1.6.-31.8.
Mosaikointi perustui NDVI:n 95%-persentiiliin, eli pikselin arvoksi valittiin sen päivän havainto joka oli lähinnä kyseistä presentiiliarvoa. Maksimia ei käytetty koska yritettiin välttää pilvien varjoja lopullisessa mosaiikissa.
Mosaikoinnin lopuksi suoritettiin koordinaatistomuunnos suomalaiseen TM35Fin-koordinaatistoon (EPSG 3067) 20 metrin pikselikokoon. Mosaiikin alue on
- ULE: 50000
- ULN: 7800000
- LRE: 842000
- LRN: 6600000
Seuraavat kuvaindeksimosaiikit lasketaan:
- Normalized Difference Moisture Index NDMI = ( TM_B4 – TM_B5 ) / ( TM_B4 + TM_B5 ) = ( OLI_B5 – OLI_B6 ) / ( OLI_B5 + OLI_B6 )
- Normalized Difference Tillage Index NDTI = ( TM_B5 – TM_B7 ) / ( TM_B5 + TM_B7 ) = ( OLI_B6 – OLI_B7 ) / ( OLI_B6 + OLI_B7 )
- Normalized Difference Vegetation Index NDVI = ( TM_B4 – TM_B3 ) / ( TM_B4 + TM_B3 ) = ( OLI_B5 – OLI_B4 ) / ( OLI_B5 + OLI_B4 )
- Normalized Difference Water Index NDWI = ( TM_B2 – TM_B4 ) / ( TM_B2 + TM_B4 ) = ( OLI_B3 – OLI_B5 ) / ( OLI_B3 + OLI_B5 )
- TM_Bx: Landsat-5 TM ja -7 ETM+ -kuvien kanava x
- OLI_Bx: Landsat-8 OLI ja -9 OLI-2 -kuvien kanava x
Käyttäjä voi laskea muita indeksejä lataamalla tarvittavat kanavat ja suorittamalla laskut sopivalla ohjelmistolla (esim. QGIS).
Muut julkaistavat mosaiikit ovat:
- META: Pikselin arvo on halutun havainnon päivän numero vuoden alusta laskien.
- NUMOBS: Pikselin arvo on mosaikoinnissa käytettyjen havaintojen lukumäärä. HUOM! Peräkkäiset Landsat-kuvien yhteinen alue lasketaan kahteen kertaan.
Mosaiikit ovat ladattavissa Kansallisen satelliittidatakeskuksen (https://nsdc.fmi.fi/) arkistosta:
Landsat-5 and -7 reflektanssikanavien mosaiikit: https://pta.data.lit.fmi.fi/lans/landsat57/bands/landsat57_refl_bn_year.tif
- bn: b1, b2, b3, b4, b5 or b7
- year: välillä 1984 - 2014
Landsat-5 and -7 kuvaindeksimosaiikit: https://pta.data.lit.fmi.fi/lans/landsat57/index/landsat57_ind_year.tif
- ind: ndmi, ndti, ndvi, ndwi
- year: välillä 1984 - 2014
Landsat-5 and -7 metatietomosaiikit:
Landsat-8 and -9 reflektanssikanavien mosaiikit: https://pta.data.lit.fmi.fi/lans/landsat89/bands/landsat89_refl_bn_year.tif
- bn: b1, b2, b3, b4, b5, b6 or b7
- year: välillä 2013 - 2024
Landsat-8 and -9 kuvaindeksimosaiikit: https://pta.data.lit.fmi.fi/lans/landsat89/bands/landsat89_ind_year.tif
- ind: ndmi, ndti, ndvi, ndwi
- year: välillä 2013 - 2024
Landsat-8 and -9 metatietomosaiikit:
Aikasarjan heikkoutena on että mosaiikeissa on pilvien ja kuvien puutteen aiheuttamia aukkoja, ja pilvimaskauksen virheiden vuoksi mosaiikkeihin saattaa jäädä jäänteitä pilvistä ja pilvien varjoista. Ilmakehäkorjaus näyttäisi korjaavan näkyvän valon kanavia liikaa, etenkin punaista kanavaa, jolloin kanavien reflektanssit voivat olla negatiivisia etenkin vesialueilla. Tämän vuoksi välillä pilvisten alueiden NDVI on suurempi kuin vieressä ole samanlaisen maanpeitteen pilvetön NDVI-arvo, jolloin mosaikoinnissa pikselin arvoksi tulee valittua pilvisen pikselin arvo.
Viite: Landsat Collection 2 Level-2 Science Products, https://www.usgs.gov/landsat-missions/landsat-collection-2-level-2-science-products
Alkuperäiset Landsat-satelliittikuvat: USGS, mosaikointi ja muu prosessointi SYKE. Aineisto kuuluu SYKEn avoimiin aineistoihin (CC BY 4.0).
