Le Calcul du NDVI (normalized difference vegetation index) a été réalisé sur des images Landsat 8 theia pour l'année 2014 sur la région Languedoc-Roussillon.

Les images Landsat 8 utilisées :

Depuis le mois de mai 2014, le pôle thématique surfaces continentales Theia (https://www.theia-land.fr/) distribue des données LANDSAT 8 de grande qualité (http://spirit.cnes.fr/resto/Landsat/). Les images de réflectance Landsat (niveau 2A) sont traitées par le pôle Theia pour l'Agence Spatiale Française (CNES). L'atelier de production développé par le CNES utilise le prototype de chaîne de Niveau 2A, MACCS, développé et conçu au CESBIO. Les images au format geoTiff contiennent les réflectances de surface, corrigées des effets atmosphériques, y compris les effets d'environnement et avec une correction supplémentaire qui prend en compte les effets des pentes. D’une résolution de 30 mètres, les réflectances de surface sont codées en entiers 16 bits signés sur 7 bandes, dans l'ordre bleu 1 (Aérosols), bleu 2, vert, rouge, proche infra-rouge (PIR) et deux moyens infra-rouge (MIR1et MIR 2) (Erreur ! Source du renvoi introuvable.). De plus, les données sont fournies projetées en Lambert 93 et découpées en tuiles de de 110x110 km et décalées de 100km les unes par rapport aux autres (Figure 1) (blog, O. Hagolle, http://www.cesbio.ups-tlse.fr ).

L’essentiel des dalles sont acquises à la date du 5 mai 2014 (Figure 2a): 05/05/2014 : D0006H0001, D0006H0002, D0006H0003, D0007H0001, D0007H0002, D0007H0003, D0007H0004, D0008H0002, D0008H0003.

Cependant, quelques extraits de dalles complémentaires ont été ajoutées afin de pallier aux zones manquantes ou nuageuses dans le secteurs extrêmes Sud-ouest et le Nord-est du Languedoc Roussillon (Figure 2b): 10/04/2014 dalles D0005H0001, D0005H0002, D0006H0002, 05/14/2014 dalle : D0008H0003. Le décalage de 36 jours entre les images principales du 5 mai et ces images complémentaires du 10 avril entraine un artéfact dans la composition du NDVI du fait du décalage phénologique de la végétation.

L’ensemble des traitements a été réalisé à l’aide du logiciel libre QGIS. Un aperçu du résultat final est donnée en figure 3. Il a consisté en la réalisation des étapes ci-dessous :

- Acquisition des données

- Mosaïquage des dalles

- Découpage de la mosaïque sur le Languedoc-Roussillon

- Calcul du NDVI (NDVI = (?PIR-?Rouge)/(?PIR+?Rouge)

- échantillonnage à 50 mètres