Correction des cirrus sur les images de Sentinel-2 et Landsat 8

=>Grâce à notre collaboration avec le DLR, qui a mis au point la méthode [1], nous avons mis en place dans MAJA une correction des cirrus fins pour LANDSAT 8 et Sentinel-2. Cette correction utilise la bande cirrus située à 1.38 µm, qui permet d’estimer la réflectance des cirrus qui est ensuite soustraite des autres bandes avec un facteur de proportionnalité calculé empiriquement. 

La chaîne MAJA est faite pour traiter automatiquement de grandes séries de données, il n’est donc pas possible de contrôler visuellement ses sorties pour décider d’appliquer ou non, la correction des cirrus. Or, même si cette correction donne souvent de bons résultats, qui paraissent même parfois magiques, comme dans l’image ci-contre, elle n’est pas toujours fiable, pour plusieurs raisons :

  • certains cirrus sont trop épais pour être corrigés,
  • d’autres sont trop inhomogènes (car la bande cirrus n’observe pas sous le même angle que les autres bandes),
  • lorsqu’on corrige le cirrus, l’ombre du nuage apparaît souvent en dessous
  • le facteur empirique peut être imprécis

 Nous avons donc décidé de n’appliquer la correction qu’aux pixels classés comme cirrus, que les traitements automatiques n’utilisent pas.  Il nous a donc fallu lisser un peu le résultat au niveau des contours pour éviter que des discontinuités n’apparaissent [2]. Si vous souhaitez utiliser ces pixels correspondant à des cirrus après correction, nous recommandons de ne le faire qu’après avoir visualisé l’image. Voici donc quelques exemples de résultats. Les animations ci-dessous alternent la sortie actuelle de MAJA, sans correction des cirrus (RUN2), avec celle que nous aurons dans quelques jours (RUN3). Vous pouvez remarquer que seuls les pixels situés à l’intérieur des contours verts changent.

Quelques autres exemples issus de Sentinel-2. N’hésitez pas à cliquer pour voir les images à meilleure résolution. Sur l’image au centre, l’algo a tenté de corriger un nuage plutôt épais, et la correction laisse apparaître une zone sombre, qui n’est pas une réussite sur le plan esthétique. Mais pas de problème pour l’utilisation des données, la zone est bien marquée comme nuageuse.

 Quelques autres exemples issus de Landsat-8 :

References

[1] Richter, Rudolf, Xingjuan Wang, Martin Bachmann, and Daniel Schläpfer. « Correction of cirrus effects in Sentinel-2 type of imagery. » International journal of remote sensing 32, no. 10 (2011): 2931-2941.[2] Olivier Hagolle, Mireille Huc, Camille Descardins, Stefan Auer, Rudolf Richter, MAJA Algorithm Theoretical Basis Document, https://doi.org/10.5281/zenodo.1209633

Plus d'actualités

Dates de fauche en France

En 2024, des collègues du Cesbio ont publié un article* sur la cartographie de la date de fauche en France en 2022 à partir des données satellitaires Sentinel-2.   Leur magnifique figure 10 a attiré mon attention. La version à haute résolution fournie avec l’article (2861 × 2911 pixels) correspond à une image dont les pixels font […]

Évolution du jour de déneigement dans les Alpes françaises et les Pyrénées

Les socio-écosystèmes des Alpes et des Pyrénées dépendent étroitement des fluctuations annuelles du manteau neigeux. En particulier, le moment de l’année où la neige disparait détermine le début de la saison de croissance de la végétation de montagne et donc la période des estives. Le changement climatique est en train de bouleverser ce rythme saisonnier. […]

Sentinel-2 reveals the surface deformation after the 2025 Myanmar earthquake

Sentinel-2 captured several clear-sky images of Myanmar before and after the 28 March 2025 earthquake. The animation below shows a 5-day apart sequence of images captured by Sentinel-2B and Sentinel-2C (10 m resolution) near the epicenter located close to Mandalay. The surface slip due to the earthquake follows the Sagaing Fault, a major fault in […]

Rechercher