Combined exploitation of VENμS, Sentinel-2 and Landsat-8: the spectral bands

=>

The combined use of VENμS, Sentinel-2 and Landsat-8 data can increase the likelihood of obtaining cloud-free images or may allow detailed tracking of rapidly evolving phenomena.

In order to facilitate this combination, the table below summarizes the correspondences between the spectral bands of the instruments. VENμS does not have a spectral band in the middle infrared.

The figure below shows the spectral bands of VENμS and Sentinel-2 in the 400 to 1000 nm range. The SWIR bands of Sentinel-2 are not included.The table below shows the usual band combinations


The figure below makes it possible to assess the degree of similarity of the spectral responses of these usual bands.

The detailed spectral responses of each instrument are available via the following web pages:

VENµS

http://www.cesbio.ups-tlse.fr/multitemp/?page_id=14229

SENTINEL-2

https://earth.esa.int/web/sentinel/user-guides/sentinel-2-msi/document-library/-/asset_publisher/Wk0TKajiISaR/content/sentinel-2a-spectral-responses

LANDSAT

https://landsat.usgs.gov/spectral-characteristics-viewer

https://landsat.usgs.gov/landsat/spectral_viewer/bands/Ball_BA_RSR.xlsx

 

 

 

Exploitation combinée de VENµS, Sentinel-2 and Landsat-8 : les bandes spectrales

=>

L’utilisation combinée des données de VENµS, Sentinel-2 et Landsat-8 peut permettre d’augmenter la probabilité d’obtenir des images sans nuage ou de suivre de manière détaillée des phénomènes à évolution rapide.

Afin de faciliter cette combinaison, le tableau ci-dessous présente de manière résumée les correspondances entre les bandes spectrales des instruments. VENµS ne comporte pas de bande spectrale dans le moyen infrarouge.

La figure ci-dessous présente les bandes spectrales de VENµS et Sentinel-2 dans le domaine 400 à 1000 nm. Les bandes SWIR de Sentinel-2 ne sont incluses.Le tableau ci-dessous présente les combinaisons de bandes usuelles

La figure ci-après permet d'apprécier le degré de similarité des réponses spectrales de ces bandes usuelles.

Les réponses spectrales détaillées de chaque instrument sont disponibles via les pages web suivantes :

VENµS

http://www.cesbio.ups-tlse.fr/multitemp/?page_id=14229

SENTINEL-2

https://earth.esa.int/web/sentinel/user-guides/sentinel-2-msi/document-library/-/asset_publisher/Wk0TKajiISaR/content/sentinel-2a-spectral-responses

LANDSAT

https://landsat.usgs.gov/spectral-characteristics-viewer

https://landsat.usgs.gov/landsat/spectral_viewer/bands/Ball_BA_RSR.xlsx

 

 

 

MAJA 3.3 is available, with a LOT of improvements

What's new ?

Pfew ! It has been quite long, but MAJA 3.3 is available, and it improves a LOT of things !

  • Some bugs have been fixed, like the one which caused detection of cloud or cloud shadows on the edges of the images
  • It seems we have finally solved the bugs that plagued the CAMS option since we released MAJA V3.0. Since V3.0, this option uses the Copernicus Atmosphere aerosol forecasts to set the aerosol type before retrieving the aerosol optical thickness (AOT) from Sentinel-2 data
  • We now also use CAMS AOT as a default value, when it is not possible to estimate AOT using the images, for instance above a snow covered landscape or for small gaps in a large cloud cover. Before that, we used 0.1 everywhere as a default value. The default value is used in the cost function with a very low weight, it has no impact when conditions for AOT estimates are good, but a large impact in bad conditions.
  • The cirrus correction module was over correcting the impact of thick cirrus clouds, providing images with dark clouds. We have limited the correction in order to get more realistic values
  • We have improved the cloud detection, with a better compromise between false positives and false negatives. We also handle better the variation against altitude of the cirrus cloud detection with band 10 (1.38 µm). MAJA 3.3 is the version with which we obtained the results of our recently published article. This paper shows that MAJA has slightly better performances than FMask 4.0, and much better performances than Sen2Cor.

Moreover, when we validated the results, we figured out that one of the parameters in our settings had a wrong value (10 instead of 1). It is easy to make such errors, because there are about 150 parameters in MAJA, and it's easy to make an error. We have set up a version management of MAJA settings since 2017, but the erroneous value was already there before that. And this value has a big impact ! The standard deviation of errors in AOT estimates is reduced by 30 to 40% !!

The W_dark parameter controls the weight of the dark pixel method in the AOT estimation. This method is just supposed to be used as a safeguard in case the multi-temporal or multi-spectral methods provide wrong results. It should therefore have a low weight, but with a weight of 10, it was in fact the method which had the highest weight in our estimates. As this method provides a maximum value of the AOT, it tended to reduce the provide too low AOTs. This improvement is therefore a great piece of news, but it comes with some shame not to have found this error before.

 

AOT validation against AERONET for 10 sites with the wrong W_dark value AOT validation against AERONET for 10 sites with the correct W_dark value

The blue dots correspond to validation obtained in good conditions, while red dots correspond to less reliable validation points.

These are the results of comparison of version 3.3 with CAMS activated, changing only the W_dark parameter, but if we compare with the results of version 3.1, the improvement is even more impressive :

AOT validation against AERONET for 10 sites, version 3.1/td> AOT validation against AERONET for 10 sites version 3.3, with the correct W_dark value

How to access MAJA ?

Here is how to access MAJA 3.3 :

  • MAJA 3.3 is distributed as a free software for non commercial purposes from CNES free software site (select the 3.3 TM version in the download tab). If you need it for commercial purposes, you just have to ask me for a different licence, but it still will be free of charge.
  • The best way to use MAJA is to run it with Start_MAJA, which is a simple python code that runs MAJA for a whole time series for a given Sentinel-2 tile. The Strart MAJA readme also explains how to get the good settings, with the good W_dark value, how to prepare the DEM or how to get CAMS data.
  • PEPS on-demand processing facility will be updated soon, but it is still working with MAJA 3.2 so far
  • THEIA is also running MAJA 3.2. We will update first the wdark parameter, and then MAJA 3.3 and then start production with CAMS, hopefully before summer. If everything goes well, we will then start a reprocessing of all our data set. So, stay tuned on this information channel.

Les séries temporelles de niveau 3A de Sentinel-2

=>

Le beau temps est revenu avant la fin de l'hiver cette année, et les synthèses de février et mars 2019 sont donc très belles, avec bien peu d'artefacts, même si les forts angles solaires ne simplifient pas le travail. Comme chaque mois, Peter Kettig du CNES a produit les synthèses de niveau 3A à partir des données Sentinel-2 du mois précédent. D'ailleurs, c'est avec un grand plaisir que je peux annoncer que, grâce à son brillant travail, Peter a été recruté au CNES sur un poste permanent !

 
Les données à pleine résolution, avec leurs masques de qualité, peuvent être téléchargées depuis le serveur de distribution Theia au CNES.
 
Si vous n'avez pas peur d'y passer trop de temps, alors que de nombreuses urgences vous attendent, vous pouvez jeter un œil aux mosaïques de ces produits disponibles sur la France depuis Juillet. Chaque mosaïque est accessible à partir des liens ci-dessous :

 
Une chouette interface de visualisation (merci à Michel Lepage !), est aussi disponible ci-dessous, pour comparer la synthèse d'octobre avec celle des mois précédents.

voir en plein écran

Continue reading

Sentinel-2 Level-3A time series

=>

This year, good weather came back earlier, and the syntheses of February and March are quite nice ! As every month, Peter Kettig from CNES processed the Sentinel-2 L3A composites of France from the Month before. And I am very pleased to announce that thanks to his brilliant work, Peter has now been hired by CNES on a permanent position.

 
The full resolution data, and the corresponding data quality masks, can be downloaded from Theia's distribution server at CNES.
 
If you are not afraid to spend too much time while you have urgent things to do, you may have a look to the mosaic of Sentinel-2 monthly syntheses for each month since July over France. Each monthly synthesis is accessible using the following links :

Or you may also use the nice viewer below (merci Michel Lepage !) to compare with the previous months.

See it full screen
Continue reading

Spot the odd one out

=>

There is an odd image in this time series of L2A products of the 31TCJ Sentinel-2 tile (Toulouse region). Can you guess which one ?

 

Yes it is the last one, acquired on the 26th of February. But what's odd with it ?

  • The black South East corner ? No, just a different orbit
  • The date ? Close enough...
  • Two images separated by one day (25 th, 26 th of January) ? You're boiling !

Continue reading

Trouvez l'intrus

=>

Parmi cette série temporelle d'images de niveau 2A de la tuile 31TCJ (Toulouse), il y a une intruse. Saurez vous déviner laquelle ?

 

Oui, c'est la dernière,, acquise le 26 Février. Mais qu'a t'elle de spécial ?

  • Le coin noir au sud-est ? Oui, mais c'est juste le résultat d'une orbite diffférente. Je ne vais quand même pas écrire tout un article là dessus.
  • La date ? Oui, vous chauffez.
  • Deux images séparées par un seul jour (les 25 et 26 février ?) ? Ah, là, vous brulez !

Continue reading