{"id":1062,"date":"2015-03-25T07:01:14","date_gmt":"2015-03-25T06:01:14","guid":{"rendered":"http:\/\/www.sigterritoires.fr\/?p=1062"},"modified":"2015-03-28T01:01:13","modified_gmt":"2015-03-28T00:01:13","slug":"utiliser-les-images-landsat-gratuites-dans-qgis-2-8-pour-le-ndvi","status":"publish","type":"post","link":"https:\/\/www.sigterritoires.fr\/index.php\/utiliser-les-images-landsat-gratuites-dans-qgis-2-8-pour-le-ndvi\/","title":{"rendered":"Utiliser les images Landsat (gratuites) dans QGis 2.8 pour le NDVI"},"content":{"rendered":"

Nous avons vu dans les deux articles pr\u00e9c\u00e9dents (Utiliser les images Landsat (gratuites) dans ArcGis (ArcMap) pour le NDVI<\/a> ,Utiliser les images Landsat (gratuites) dans ArcGis Pro pour le NDVI<\/a>) le calcul et l’affichage de l’indice NDVI avec arcMap, puis ArcGis Pro. Nous verrons dans celui-ci comment r\u00e9aliser cette op\u00e9ration avec QGis.
\nPour le t\u00e9l\u00e9chargement des images Landsat \u00e0 partir d’Earth Explorer, consultez cet article: Utiliser les images Landsat (gratuites) dans votre SIG<\/a>.<\/p>\n

Dans QGis nous allons utiliser la calculatrice raster pour calculer le NDVI<\/a>, puis nous verrons comment cr\u00e9er et utiliser un d\u00e9grad\u00e9 de couleurs sp\u00e9cifique \u00e0 cet index.
\n<\/p>\n

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).<\/p>\n

$\"chargement$ <\/a><\/p>\n

Pour ouvrir la calculatrice raster, cliquez sur le menu Raster->Calculatrice raster<\/p>\n

$\"calculatrice$ <\/a><\/p>\n

Nous allons rentrer la formule de calcul<\/p>\n

NDVI = PIR – R \/ PIR + R<\/strong><\/p>\n

O\u00f9 PIR<\/strong> est la bande du proche infra-rouge et R<\/strong> la bande du rouge<\/p>\n

Dans notre exemple cela donne:
\n( \u00ab\u00a0LE72020272014184ASN00_B4@1\u00a0\u00bb – \u00ab\u00a0LE72020272014184ASN00_B3@1\u00a0\u00bb ) – ( \u00ab\u00a0LE72020272014184ASN00_B4@1\u00a0\u00bb + \u00ab\u00a0LE72020272014184ASN00_B3@1\u00a0\u00bb )<\/p>\n

Le r\u00e9sultat s’affiche:<\/p>\n

$\"r\u00e9sultat$ <\/a><\/p>\n

Vous pouvez travailler directement sur ce r\u00e9sultat, mais il est toujours utile de se r\u00e9f\u00e9rer \u00e0 une symbologie standard. On a l’habitude de repr\u00e9senter le NDVIsur une \u00e9chelle compl\u00e8te de -1 \u00e0 +1 (ici le r\u00e9sultat va seulement de -0.5 \u00e0 +0.506), de repr\u00e9senter les valeurs n\u00e9gatives, qui ne correspondent pas \u00e0 de la v\u00e9g\u00e9tation, par un d\u00e9grad\u00e9 de bleu, et les valeurs positives, qui correspondent \u00e0 la v\u00e9g\u00e9tation, par un d\u00e9grad\u00e9 de vert au rouge passant par le jaune.<\/p>\n

$\"colormap$ <\/a><\/p>\n

Seul ennui, ce d\u00e9grad\u00e9 n’est pas disponible dans \u00a0QGis tel quel.<\/p>\n

Pour le cr\u00e9er, ouvrez le bloc-notes et tapez le texte suivant:
\n
\n<?xml version=\u00a0\u00bb1.0″ encoding=\u00a0\u00bbUTF-8″?>
\n<svg xmlns=\u00a0\u00bbhttp:\/\/www.w3.org\/2000\/svg\u00a0\u00bb version=\u00a0\u00bb1.1″ width=\u00a0\u00bb300px\u00a0\u00bb height=\u00a0\u00bb45px\u00a0\u00bb viewBox=\u00a0\u00bb0 0 300 45″>
\n<g>
\n<defs>
\n<linearGradient id=\u00a0\u00bbNDVI\u00a0\u00bb gradientUnits=\u00a0\u00bbobjectBoundingBox\u00a0\u00bb spreadMethod=\u00a0\u00bbpad\u00a0\u00bb x1=\u00a0\u00bb0%\u00a0\u00bb x2=\u00a0\u00bb100%\u00a0\u00bb y1=\u00a0\u00bb0%\u00a0\u00bb y2=\u00a0\u00bb0%\u00a0\u00bb>
\n<stop offset=\u00a0\u00bb0.00%\u00a0\u00bb stop-color=\u00a0\u00bbrgb(0,0,255)\u00a0\u00bb stop-opacity=\u00a0\u00bb1.0000″\/>
\n<stop offset=\u00a0\u00bb49.90%\u00a0\u00bb stop-color=\u00a0\u00bbrgb(190,210,255)\u00a0\u00bb stop-opacity=\u00a0\u00bb1.0000″\/>
\n<stop offset=\u00a0\u00bb50.00%\u00a0\u00bb stop-color=\u00a0\u00bbrgb(56,168,0)\u00a0\u00bb stop-opacity=\u00a0\u00bb1.0000″\/>
\n<stop offset=\u00a0\u00bb66.60%\u00a0\u00bb stop-color=\u00a0\u00bbrgb(255,255,0)\u00a0\u00bb stop-opacity=\u00a0\u00bb1.0000″\/>
\n<stop offset=\u00a0\u00bb100.00%\u00a0\u00bb stop-color=\u00a0\u00bbrgb(255,0,0)\u00a0\u00bb stop-opacity=\u00a0\u00bb1.0000″\/>
\n<\/linearGradient>
\n<\/defs>
\n<rect fill=\u00a0\u00bburl(#NDVI)\u00a0\u00bb x=\u00a0\u00bb4″ y=\u00a0\u00bb4″ width=\u00a0\u00bb292″ height=\u00a0\u00bb37″ stroke=\u00a0\u00bbblack\u00a0\u00bb stroke-width=\u00a0\u00bb1″\/>
\n<\/g>
\n<metadata>
\n<creator name=\u00a0\u00bbnasca\u00a0\u00bb version=\u00a0\u00bb1.0″\/>
\n<created date=\u00a0\u00bbMon Mar 23 22:38:20 2015″\/>
\n<\/metadata>
\n<\/svg>
\n<\/strong>
\nEnregistrez ce fichier avec le nom NDVI.svg<\/strong> dans le r\u00e9pertoire \\Program Files\\QGIS Wien\\apps\\qgis\\resources\\cpt-city-qgis-min\\cb\\div<\/em><\/p>\n

Maintenant on va appliquer ce d\u00e9grad\u00e9 au raster r\u00e9sultant du calcul du NDVI:<\/p>\n

Double cliquez sur la l\u00e9gende du raster pour ouvrir la fen\u00eatre de propri\u00e9t\u00e9s
\nDans l’onglet Style, changez le Type de rendu <\/em> en \u00ab\u00a0Pseudo-couleur \u00e0 bande unique\u00a0\u00bb
\nChangez les valeurs de Min en -1 et de Max en 1.<\/p>\n

$\"param\u00e9trage$ <\/a><\/p>\n

Dans la liste d\u00e9roulante de palettes, descendez jusqu’\u00e0 la derni\u00e8re ligne et cliquez sur Nouvelle palette de couleurs<\/em><\/p>\n

$\"palettes$ <\/a><\/p>\n

Dans le Type de palette de couleur<\/em>, s\u00e9lectionnez cpt-city<\/em> et cliquez sur OK<\/strong><\/p>\n

La palette cpt-city s’ouvre<\/p>\n

$\"palette$ <\/a><\/p>\n

Ouvrez l’item \u00ab\u00a0Diverging\u00a0\u00bb<\/strong> et s\u00e9lectionnez cb\/div<\/strong> pour retrouver la palette ndvi que vous avez cr\u00e9\u00e9.
\nCliquez sur celle-ci, puis sur OK. Laissez le nom par d\u00e9faut et cliquez OK.
\nDans Mode<\/strong> s\u00e9lectionnez Intervalles \u00e9gaux<\/em>.
\nDans Classe<\/strong> s\u00e9lectionnez, par exemple, 20<\/em>.
\nCliquez sur le bouton Classer<\/strong>
\nVous voyez l’\u00e9chelle qui sera appliqu\u00e9e au raster.<\/p>\n

$\"classes$ <\/a><\/p>\n

Cliquez sur OK<\/strong> pour fermer la fen\u00eatre de propri\u00e9t\u00e9s de la couche et voir le r\u00e9sultat final:<\/p>\n

$\"r\u00e9sultat$ <\/a><\/p>\n","protected":false},"excerpt":{"rendered":"

Nous avons vu dans les deux articles pr\u00e9c\u00e9dents (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…<\/p>\n","protected":false},"author":1,"featured_media":0,"comment_status":"open","ping_status":"closed","sticky":false,"template":"","format":"standard","meta":{"give_campaign_id":0,"_bbp_topic_count":0,"_bbp_reply_count":0,"_bbp_total_topic_count":0,"_bbp_total_reply_count":0,"_bbp_voice_count":0,"_bbp_anonymous_reply_count":0,"_bbp_topic_count_hidden":0,"_bbp_reply_count_hidden":0,"_bbp_forum_subforum_count":0,"sfsi_plus_gutenberg_text_before_share":"","sfsi_plus_gutenberg_show_text_before_share":"","sfsi_plus_gutenberg_icon_type":"","sfsi_plus_gutenberg_icon_alignemt":"","sfsi_plus_gutenburg_max_per_row":"","_monsterinsights_skip_tracking":false,"_monsterinsights_sitenote_active":false,"_monsterinsights_sitenote_note":"","_monsterinsights_sitenote_category":0,"_jetpack_newsletter_access":"","_jetpack_dont_email_post_to_subs":false,"_jetpack_newsletter_tier_id":0,"_jetpack_memberships_contains_paywalled_content":false,"_jetpack_memberships_contains_paid_content":false,"footnotes":"","jetpack_post_was_ever_published":false},"categories":[3,30,62],"tags":[140,139,58],"class_list":["post-1062","post","type-post","status-publish","format-standard","hentry","category-logiciel_libres","category-donnees-gratuites","category-qgis-2","tag-degrade","tag-ndvi","tag-qgis"],"aioseo_notices":[],"jetpack_featured_media_url":"","jetpack_shortlink":"https:\/\/wp.me\/p6XU0A-h8","jetpack_sharing_enabled":true,"_links":{"self":[{"href":"https:\/\/www.sigterritoires.fr\/index.php\/wp-json\/wp\/v2\/posts\/1062","targetHints":{"allow":["GET"]}}],"collection":[{"href":"https:\/\/www.sigterritoires.fr\/index.php\/wp-json\/wp\/v2\/posts"}],"about":[{"href":"https:\/\/www.sigterritoires.fr\/index.php\/wp-json\/wp\/v2\/types\/post"}],"author":[{"embeddable":true,"href":"https:\/\/www.sigterritoires.fr\/index.php\/wp-json\/wp\/v2\/users\/1"}],"replies":[{"embeddable":true,"href":"https:\/\/www.sigterritoires.fr\/index.php\/wp-json\/wp\/v2\/comments?post=1062"}],"version-history":[{"count":0,"href":"https:\/\/www.sigterritoires.fr\/index.php\/wp-json\/wp\/v2\/posts\/1062\/revisions"}],"wp:attachment":[{"href":"https:\/\/www.sigterritoires.fr\/index.php\/wp-json\/wp\/v2\/media?parent=1062"}],"wp:term":[{"taxonomy":"category","embeddable":true,"href":"https:\/\/www.sigterritoires.fr\/index.php\/wp-json\/wp\/v2\/categories?post=1062"},{"taxonomy":"post_tag","embeddable":true,"href":"https:\/\/www.sigterritoires.fr\/index.php\/wp-json\/wp\/v2\/tags?post=1062"}],"curies":[{"name":"wp","href":"https:\/\/api.w.org\/{rel}","templated":true}]}}