Utiliser les images Landsat (gratuites) dans QGis 2.8 pour le NDVI

Nous avons vu dans les deux articles précédents (Utiliser les images Landsat (gratuites) dans ArcGis (ArcMap) pour le NDVI ,Utiliser les images Landsat (gratuites) dans ArcGis Pro pour le NDVI) le calcul et l’affichage de l’indice NDVI avec arcMap, puis ArcGis Pro. Nous verrons dans celui-ci comment réaliser cette opération avec QGis.
Pour le téléchargement des images Landsat à partir d’Earth Explorer, consultez cet article: Utiliser les images Landsat (gratuites) dans votre SIG.

Dans QGis nous allons utiliser la calculatrice raster pour calculer le NDVI, puis nous verrons comment créer et utiliser un dégradé de couleurs spécifique à cet index.

Si vous travaillez avec une image Landsat 7, vous chargerez dans QGis les bandes 3 (rouge) et 4 (infra-rouge). Si vous travaillez avec une image Landsat 8,vous chargerez les bandes 4 (rouge ) et 5 (infra-rouge).

chargement des bandes rouges et infra-rouges dans qgis 2.8

Pour ouvrir la calculatrice raster, cliquez sur le menu Raster->Calculatrice raster

calculatrice raster de qgis 2.8

Nous allons rentrer la formule de calcul

NDVI = PIR – R / PIR + R

PIR est la bande du proche infra-rouge et R la bande du rouge

Dans notre exemple cela donne:
( « LE72020272014184ASN00_B4@1 » – « LE72020272014184ASN00_B3@1 » ) – ( « LE72020272014184ASN00_B4@1 » + « LE72020272014184ASN00_B3@1 » )

Le résultat s’affiche:

résultat du calcul du ndvi avec la calculatrice raster

Vous pouvez travailler directement sur ce résultat, mais il est toujours utile de se référer à une symbologie standard. On a l’habitude de représenter le NDVIsur une échelle complète de -1 à +1 (ici le résultat va seulement de -0.5 à +0.506), de représenter les valeurs négatives, qui ne correspondent pas à de la végétation, par un dégradé de bleu, et les valeurs positives, qui correspondent à la végétation, par un dégradé de vert au rouge passant par le jaune.

colormap pour affichage du ndvi

Seul ennui, ce dégradé n’est pas disponible dans  QGis tel quel.

Pour le créer, ouvrez le bloc-notes et tapez le texte suivant:

<?xml version= »1.0″ encoding= »UTF-8″?>
<svg xmlns= »http://www.w3.org/2000/svg » version= »1.1″ width= »300px » height= »45px » viewBox= »0 0 300 45″>
<g>
<defs>
<linearGradient id= »NDVI » gradientUnits= »objectBoundingBox » spreadMethod= »pad » x1= »0% » x2= »100% » y1= »0% » y2= »0% »>
<stop offset= »0.00% » stop-color= »rgb(0,0,255) » stop-opacity= »1.0000″/>
<stop offset= »49.90% » stop-color= »rgb(190,210,255) » stop-opacity= »1.0000″/>
<stop offset= »50.00% » stop-color= »rgb(56,168,0) » stop-opacity= »1.0000″/>
<stop offset= »66.60% » stop-color= »rgb(255,255,0) » stop-opacity= »1.0000″/>
<stop offset= »100.00% » stop-color= »rgb(255,0,0) » stop-opacity= »1.0000″/>
</linearGradient>
</defs>
<rect fill= »url(#NDVI) » x= »4″ y= »4″ width= »292″ height= »37″ stroke= »black » stroke-width= »1″/>
</g>
<metadata>
<creator name= »nasca » version= »1.0″/>
<created date= »Mon Mar 23 22:38:20 2015″/>
</metadata>
</svg>

Enregistrez ce fichier avec le nom NDVI.svg dans le répertoire \Program Files\QGIS Wien\apps\qgis\resources\cpt-city-qgis-min\cb\div

Maintenant on va appliquer ce dégradé au raster résultant du calcul du NDVI:

Double cliquez sur la légende du raster pour ouvrir la fenêtre de propriétés
Dans l’onglet Style, changez le Type de rendu en « Pseudo-couleur à bande unique »
Changez les valeurs de Min en -1 et de Max en 1.

paramétrage du style NDVI

Dans la liste déroulante de palettes, descendez jusqu’à la dernière ligne et cliquez sur Nouvelle palette de couleurs

palettes de couleurs

Dans le Type de palette de couleur, sélectionnez cpt-city et cliquez sur OK

La palette cpt-city s’ouvre

palette cpt-city

Ouvrez l’item « Diverging » et sélectionnez cb/div pour retrouver la palette ndvi que vous avez créé.
Cliquez sur celle-ci, puis sur OK. Laissez le nom par défaut et cliquez OK.
Dans Mode sélectionnez Intervalles égaux.
Dans Classe sélectionnez, par exemple, 20.
Cliquez sur le bouton Classer
Vous voyez l’échelle qui sera appliquée au raster.

classes de couleur retenues

Cliquez sur OK pour fermer la fenêtre de propriétés de la couche et voir le résultat final:

résultat final du calcul de ndvi avec qgis 2.8

15 réflexions sur « Utiliser les images Landsat (gratuites) dans QGis 2.8 pour le NDVI »

  1. Bonjour,

    Merci pour cette explication claire et précise. Il me vient une question que je n’arrive pas à percer :
    Comment pour chaque classe, arriver à connaître la surface qui lui est liée ? En effet l’image est très parlante mais avoir des données chiffrées serait encore mieux. Par exemple on pourrait connaitre la surface et ensuite la proportion de surface inerte

    bonne journée

    1. Vous avez dans l’article le mode d’emploi pour créer le fichier que vous nommerez ndvi.svg. Une fois créé, suivez le mode opératoire à la fin de l’article pour l’appliquer à votre couche.
      Pour l’endroit où mettre le fichier ndvi.svg, remplacez dans la chaîne suivante :\Program Files\QGIS Wien\apps\qgis\resources\cpt-city-qgis-min\cb\div QGIs Wien par votre version de qgis

        1. Bonjour,

          Mon probléme et que lorsque je vais pour trouver ma palette NDVI celle-ci est un dégradé de bleu et de vert uniquement… Comment pourrai-je arranger le soucis ?

          Merci !

          Bien à vous !

          1. Merci beaucoup Atilio Francois de ces conseils !

            Il ne me reste plus qu’a trouver comment enregistrer en tant qu’administrateur pour pouvoir le mettre dans le dossier div de QGIS, mais ça ne sera pas le plus dure je pense !

            Bien à vous.

  2. bonjour !
    Merci pour cette article.
    Il me vient une question! Je suis actuellement en train de faire une multitude de carte d’occupation des sols en utilisant des images landsat 7 et 8.
    Sur les images landsat 7 il y a ce problème de lignes noirs liées à la défaillance d’un des capteur et que nous pouvons voir sur votre résultat.
    J’aurai voulu savoir si vous connaissez une méthode pour pouvoir supprimer ces bandes noires ?

    merci
    Bien à vous

  3. Bonjour et merci d’avance pour ce cours intéressant. Mon seul hic est que je n’arrive pas à trouver la palette de couleur NDVI téléchargée avec le lien que vous nous avez donné sur Qgis alors que quand je me rend dans le dossier ou se trouve le NDVI.svg je le vois…

Laisser un commentaire

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