La France sans nuages et sans coutures, vue par Sentinel-2 et Theia
=>
UPDATE: résolution améliorée à 20m !Voilà, après quelques jours de bidouille, j’ai réussi à produire la mosaïque de la France à partir:
- des synthèses mensuelles de produits Sentinel-2 de Niveau 3A fournis par Theia en juillet 2018; Ces produits sont réalisés à partir de données de niveau 2A, corrigé des effets atmosphériques, et avec une bonne détection des nuages, grâce à la chaîne MAJA. Ces données seront produites tous les mois
- du script mis au point par Simon Gascoin
- des conseils avisés de Simon et de Michel Le Page
- et de Gdal, (Merci Gdal !)
La mosaïque vous est présentée à 20m de résolution, elle occupe déjà 8GO. Il aurait été possible de la produire à 10m, mais il aurait fallu 32 GO et quelques dizaines d’heures de plus. Voir en plein écran
La sécheresse qui a atteint le Nord-Est de la France est déjà visible, avec une couleur beaucoup plus brune que le reste du pays. Notre produit de niveau 3A cherche à minimiser les artefacts, mais sur une base purement physique, sans homogénéisation cosmétique sur les lignes de raccord, car il n’y a pas de lignes de raccord . Vous trouverez très peu d’artefacts, mais vous en trouverez :
- des bords de tuiles horizontaux ou verticaux, qui sont dus au fait que les dates centrales des produits de niveau 3A ne sont pas toutes les mêmes (du 7 au 27 juillet). A partir du mois d’août, nous utiliserons la date centrale du mois dans sa production.
- quelques zones montagneuses sont restées nuageuses pendant toutes la période. Elles sont marquées nuageuses dans les masques qui accompagnent le produit
- des bords d’orbite, peut-être, mais je n’en ai pas encore vu
- des problèmes sur l’eau et sur la neige, car ces surfaces varient avec le temps
- quelques petits nuages par ci par là, oubliés par MAJA
Si les artefacts que vous trouvez ne rentrent pas dans ces catégories, ils m’intéressent !La génération de la mosaïque fut un peu plus compliquée que pour la Belgique, parce que la France, ce grand pays ;), s’étale sur 3 fuseaux UTM (31, 31, 32). Il faut donc rassembler les tuiles sous une même projection. On s’en sort quand même en quelques lignes, comme vous pouvez le voir ici .