Recent bugfix of MAJA v4, please consider upgrading to 4.5.4

A bug affecting the AOT estimates and leading to underestimated surface reflectances in some rare cases concerns all MAJA v4.x versions prior to 4.5.4. Users activating the CAMS option shouldn’t have noticed, but if you don’t use CAMS auxiliary data you may face this issue.

Therefore, we highly recommend all users to upgrade to the latest MAJA 4.5.4 available for download here: Download Maja 4.5.4

MAJA version 4 was a complete recoding of MAJA version 3’s internal scheduler (the « Main » program). It went from a monolithic C++ code, hard to maintain, to the use of OTB applications chained in a Python orchestrator.

Despite all our tests (global, algorithmic, unitary), developed on a continuous integration platform, it seems we had missed one bug (let’s hope it’s the last one). This bug appeared in the aerosol estimation method, the most complex one, and caused in some parts of rare images, the apparition of high optical thicknesses and low surface reflectances. It then could propagate to the next images, as the method is multi-temporal.

It took us some time to figure out what was happening, and Julie Brossard finally found out that it was due to the reading of a float parameter as an integer (Wdark).

Wdark is the weight of an equation in the cost function of the aerosol optical thickness (AOT) estimates. This equation increases the cost if surface reflectances are really close or below 0.  This equation is therefore useful if something goes wrong in the aerosol estimates, with a too high AOT, yielding very low surface reflectances.

In most cases this equation is not needed, but there are some cases when it is necessary: if the prescribed type of aerosol is wrong, or if there was an issue with the cloud detection. We are very conservative with clouds in MAJA and every suspicious pixel is not considered in the aerosol estimates, but still, some cases might escape to the detection, for example, thin clouds above shadows, whose reflectance stays in the expected range for a clear pixel.

In the original versions we used to use a value of 10, which is not affected by that bug. This value is still in use in our algorithmic tests. But a few years go, I figured out that Wdark was too high, and tended to bias the aerosol estimates even if the conditions were perfect for the estimate. After a lot of tests, it was set to 0.5 . And there, the truncated value has a huge effect.

Within THEIA, we process data using CAMS to decide which aerosol type has to be used. This eliminates one of the error causes which make Wdark useful. But users who do not use the CAMS data are somewhat more affected.

We thank the Sen2Agri project who raised this issue, and CS Group for the investigation and bugfix. If you are interested, you might follow the discussions here , and have an idea of the quantity of work needed to maintain a software like MAJA.

MAJA v3 MAJA v4 before bugfix
Example of MAJA v3 L2ASame product as (left) with MAJA v4 before bugfix. After bugfix, both images are identical.

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. […]

BIOMASS, the third launched satellite mission designed at CESBIO !

After SMOS in 2009, and VENµS in 2017, the CESBIO laboratory is very proud to see its third proposed mission, Biomass, reach orbit. As always, it has been a long journey from the idea, at the beginning of the century, to the selection in 2013 as the seventh Earth Explorer Mission by ESA, to the […]

Rechercher