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).
Pour ouvrir la calculatrice raster, cliquez sur le menu Raster->Calculatrice raster
Nous allons rentrer la formule de calcul
NDVI = PIR – R / PIR + R
Où 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:
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.
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.
Dans la liste déroulante de palettes, descendez jusqu’à la dernière ligne et cliquez sur Nouvelle palette de couleurs
Dans le Type de palette de couleur, sélectionnez cpt-city et cliquez sur OK
La palette cpt-city s’ouvre
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.
Cliquez sur OK pour fermer la fenêtre de propriétés de la couche et voir le résultat final:
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
Bonjour, j’utilise qgis 2.14 mais je ne retrouve pas la couleur ndvi. Qui peut m’aider à la charger?
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
Bonjour François, merci c’était juste à côté mais le cerveau se chauffe on devient aveugle. Merci beaucoup.
Bonjour,
Je n’arrive pas à créer la palette de couleur, comment ouvre-t-on le bloc note?
Pour vous simplifier la tâche, vous pouvez télécharger le fichier ndvi.svg directement à l’adresse http://www.nasca.ovh/downloads/ndvi.svg
Bonjour, comment télécharger depuis ce liens ?
Merci d’avance :-).
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 !
Bonjour,
J’ai le même problème avec un dégradé de bleu et vert uniquement … Pourtant j’utilise votre lien que j’enregistre dans le bon dossier. Avez-vous une idée de ce qui peut poser problème ?
Merci beaucoup
A la question si j’ai une idée, pour l’instant la réponse est non. Je n’arrive pas à reproduire cette « erreur » sur ma machine. Mais je ne désespère pas de réussir (à avoir l’erreur pour vous dire comment la résoudre)
Ne cliquez pas sur le lien, faites un clic droit sur le lien et sélectionnez « enregistrer la cible du lien sous » (avec Firefox)
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.
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
Il y en a une, publiée chez Editions Universitaires Européennes : Amélioration des images Landsat SLC-off par dilatations successives, de Joseph Libar. Je n’ai pas la publication moi-même.
tuto au top, merci beaucoup
excellent travail, merci énormément !
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…
Bonjour, merci pour ce tutoriel.
j’essaye de faire une reclassification d’une carte NDVI de ma zone d’étude, mais je trouve un problème, c’est que les zones agricoles ont les mêmes valeurs NDVI que les zones forestières (> 0.5) ce qui n’est pas logique !! même la carte initiale de NDVI montrent des couleurs identiques pour ces deux types d’usages… Avez-vous une idée ?
Merci
je ne vois pas d’impossibilité logique. il y a des cultures qui peuvent avoir des NDVI de même classe que la couverture forestière. J’ai en mémoire un travail dans les landes où les valeurs des cultures de maîs en été était de même ordre que celui de la forêt des Landes.
Bref, tout dépend de la saison, du type de culture, du stress hydrique…
Et encore plus bref, à la Normande: peut-être ben que oui, peut-être ben que non
Merci pour la réponse. Mais comment peux-je obtenir une carte de NDVI qui distingue entre les cultures et les forêts ?
Sur une seule image vous n’arriverez pas. C’est en analysant des séries comprenant différentes périodes que vous pourrez voir des réponses différentes.
Bonjour,
Tout d’abord merci pour ce post très intéressant, il m’a beaucoup aidée ! Je suis actuellement en stage et j’effectue des suivis temporels de cultures avec des images Landsat 8. Je rencontre un problème : lorsque j’enregistre ma palette de couleur NDVI à l’endroit où vous le détaillez, un message s’affiche m’informant que je ne dispose pas des droits d’utilisateur nécessaires pour exécuter cette simple sauvegarde. Comment pourrais-je remédier à ce problème?
C’est un problème de droits sous Windows. Soit vous lancez QGis en mode administrateur et il n’y aura pas le message, soit vous le sauvegardez dans un autre répertoire qui ne soit pas Program Files, puis vous le copiez manuellement dans le répertoire de Qgis. Vous aurez un message qui vous préviendra que vous devez avoir des droits d’administration pour le faire, mais il le fera quand même… à condition que vous ayez les droits d’administration sur votre PC. Si ce n’est pas le cas, il faut voir avec la personne qui les a.
Bonsoir,
J’arrive à calculer NDVI par Qgis mais j’ai besoin des statistiques de chaque compte numérique (DN) en fonction des fréquences de NDVI (% NDVI)? Est ce que vous pouvez m’aider le faire s’il vous plait!!