=>Voici une description rapide, succincte et simpliste des produits de niveau 3A distribuées par THEIA. La description officielle du format des données est ici xxx

Niveau 3A : Synthèses mensuelles de réflectances de surface

Le fichier zip distribué par THEIA contient plusieurs répertoires et plusieurs fichiers.

Les fichiers du répertoire principal sont les suivants :
  • les réflectances de surface sont fournies sous la forme d’un fichier GeoTiff par bande, pour les bandes à 10m (B2,B3,B4,B8) et les bandes à 20m (B5, B6, B7, B8A, B11 et B12) de Sentinel-2. Les bandes à 60m (B1, B9 et B10) ne sont pas fournies.
  • Les réflectances de surface intègrent une correction des variations d’éclairement dûes au relief
  • Les données sont codées sur 16 bits et il faut diviser par 10000 pour obtenir des réflectances de surface.
    • <REFLECTANCE_QUANTIFICATION_VALUE>10000</REFLECTANCE_QUANTIFICATION_VALUE>
    • La valeur No_Data (en dehors de l’image) vaut -10000

Le répertoire MASK contient les informations suivantes :
  • Un masque de flags, à résolution 10m (FLG_R1.tif) et 20m (FLG_R2.tif), codé sur 8 bits, avec 5 valeurs utiles
    • 0 = Pas de données

    • 1 = Nuage

    • 2 = Neige

    • 3 = Eau

    • 4 = Terre émergée

  • Une image, codée sur 16 bits, fournissant le total des poids utilisés dans la moyenne pondérée à résolution 10m (WGT_R1.tif) et 20m (WGT_R2.tif)
    • Cette information sera assez difficile à utiliser, mais en gros, plus la valeur observée est élevée, plus on aura confiance dans la valeur du pixel. Cette information nous a été très utile lors de la phase de validation du produit.
  • Une image de la moyenne pondérée des dates utilisées, sur 16 bits, codée en jours depuis le 1er janvier 2018, à résolution 10m (DTS_R1.tif) et 20m (DTS_R2.tif)
    Composite Date (Weighted average) Weight

=> This page describes the data format of Sentinel-2 L3A generated by THEIA MUSCATE processing center at CNES. ESA Sen2Cor format is different.The official description of the data format est ici xxx

Level3A : Monthly Syntheses of cloud free surface reflectances

The ZIP file distributed by Theia contains several files and directories.

The files are as follows :
  • the surface reflectances are provided as a GeoTiff per spectral band, for the 10m bands (B2, B3, B4, B8) and the 20m bands (B5, B6, B7, B8A, B11 et B12) of Sentinel-2. The bands at 60m resolution are not provided.
  • The surface reflectances have been corrected from illumination variation due to topography
  • The data are coded on 16bits signed integers and you have to divide by 10000 to obtain reflectances.
    • <REFLECTANCE_QUANTIFICATION_VALUE>10000</REFLECTANCE_QUANTIFICATION_VALUE>
    • The No_Data value (outside the image) is -10000

 

The MASKS directory provides :
  • A classification mask, at 10m resolution (FLG_R1.tif) and 20m (FLG_R2.tif), coded as 8 bits integers with 5 useful values
    • 0 = No Data

    • 1 = Cloud

    • 2 = Snow

    • 3 = Water

    • 4 = Land

  • An image, Coded as 16 bit integers, that provided the sum of the weights for each pixel at 10m resolution (WGT_R1.tif) and 20m resolution (WGT_R2.tif)
    • This information is very useful for our processing as the weighted average is constructed iteratively, updating the L3A product each time a new L2A is available. The higher the weight, the more confidence we have in the surface reflectance.
  • An image of the weighted average of the dates used for the weighted average, coded as Int16, expressed in days since January 1st, 2018, at 10m resolution (DTS_R1.tif) and 20m  resolution 20m (DTS_R2.tif)

8 thoughts on “THEIA’s L3A product format

  1. Bravo, c’est impressionnant.
    Cependant, je ne vois pas bien ce que représente les couches DTS:

    « Une image de la moyenne pondérée des dates utilisées, sur 16 bits, codée en jours depuis le 1er janvier 2018, à résolution 10m (DTS_R1.tif) et 20m (DTS_R2.tif) »

    Les moyennes pondérées des dates en chaque pixels ne devraient-elles pas correspondre à la date de synthèse voulue ? Ou est-ce les moyennes « non-pondérées » qui sont représentées dans ces couches ?

    1. Bonjour,
      WASP calcule ne moyenne pondérée des réflectances sur 45 jours centrés sur le 15 du mois (ou le premier pour le sahel). Mais supposons que la première moitié de la période ait été complètement nuageuse, les données utilisées correspondront plutôt à la fin de la période. Pour donner une idée de la date à laquelle chaque pixel correspond, nous calculons la moyenne pondérée des dates sans nuages, avec les même poids que ceux utilisés pour calculer la moyenne pondérée des réflectances. Le 15 du mois est donc la date que nous utilisons comme objectif, et la date dans le plan DTS, la date moyenne réelle à laquelle correspond chaque pixel.
      Olivier

  2. Hello Olivier,
    First of all, thankyou for providing such a beautiful compilation of datasets.
    However, i am facing a little difficulties while extracting the dataset.

    I have downloaded, Sentinel-2 level 3A product for these tiles (T32ULU, T32ULV and T32UMV) for August 2017, August 2019 and January 2020.
    The extraction of these zip files fails and return data error. I have tried using 7zip and winrar.
    Do i have to use specific software for extraction or am i missing something?

    It would be really helpful if you shed some light on it.

    Thank you
    Ankit

    1. Hi,It should be possible to unzip them with any unzip software. I guess something must have happened during the download. I have forwarded the messages to my colleagues, did you get an answer ?

      1. Thank you for your response. I found that it was an issue while downloading these files.
        I solved the problem by downloading these tiles again.

  3. Hi, I downloaded Sentinel-2 level 3A product for the tiles 32/U/PV/2022/. The quality seems great, however there are so many cells that have the value 0, practically all water surfaces but also other parts of the image. Only few of the « holes » are correlating with the clouds or water from the masks. Most of them I cant explain. I handle this zeros as nodata or does this make any sense for classification? Maybe you could give me a hint, why there is so much missing data.

Laisser un commentaire

Votre adresse e-mail ne sera pas publiée. Les champs obligatoires sont indiqués avec *

Ce site utilise Akismet pour réduire les indésirables. En savoir plus sur comment les données de vos commentaires sont utilisées.