{"id":1329,"date":"2022-09-13T13:09:07","date_gmt":"2022-09-13T11:09:07","guid":{"rendered":"https:\/\/wp.unil.ch\/dawn\/?page_id=1329"},"modified":"2024-11-26T16:39:38","modified_gmt":"2024-11-26T15:39:38","slug":"sig_tp7","status":"publish","type":"page","link":"https:\/\/wp.unil.ch\/dawn\/sig_tp7\/","title":{"rendered":"TP7 \u2014 Op\u00e9rations Raster"},"content":{"rendered":"\n<details class=\"wp-block-details is-layout-flow wp-block-details-is-layout-flow\"><summary>Introduction<\/summary>\n<div class=\"wp-block-group\"><div class=\"wp-block-group__inner-container is-layout-constrained wp-block-group-is-layout-constrained\">\n<p>Au cours de ce TP, tu vas te familiariser avec certains outils de g\u00e9otraitement pour les donn\u00e9es au format raster (=image). Tu travailleras avec des donn\u00e9es satellites au format raster (provenant de la <a href=\"https:\/\/www.nasa.gov\/\">NASA<\/a>), une couche th\u00e9matique raster qui identifie les classes d&#8217;occupations du sol (cr\u00e9\u00e9e par nos soins avec les donn\u00e9es satellites de la NASA) et enfin un MNT (mod\u00e8le num\u00e9rique de terrain) que tu extrapoleras (gr\u00e2ce \u00e0 la m\u00e9thode de l&#8217;IDW &#8220;Inverse Distance Weighting&#8221;) \u00e0 partir de donn\u00e9es ponctuelles disponibles sur Swisstopo. <\/p>\n\n\n\n<p>Ce TP n\u2019aurait pas \u00e9t\u00e9 possible sans les ressources list\u00e9es ci-dessous:<\/p>\n\n\n\n<ul class=\"wp-block-list\">\n<li><em>TP7 du cours \u201cG\u00e9omatique et SIG\u201d du Privat-docent Dr. Marj Tonini<\/em><\/li>\n\n\n\n<li><em>(Image de couverture) \u201cSardinia in satellite imagery<em>\u201d<\/em> by&nbsp;<a href=\"https:\/\/stock.adobe.com\/ch_de\/contributor\/210875081\/m-jurzyk?load_type=author&amp;prev_url=detail\">M. Jurzyk<\/a><\/em><\/li>\n<\/ul>\n<\/div><\/div>\n<\/details>\n\n\n\n<details class=\"wp-block-details is-layout-flow wp-block-details-is-layout-flow\"><summary>1. T\u00e9l\u00e9chargement des donn\u00e9es du TP et exploration des m\u00e9tadonn\u00e9es<\/summary>\n<div class=\"wp-block-group\"><div class=\"wp-block-group__inner-container is-layout-constrained wp-block-group-is-layout-constrained\">\n<p>1a) T\u00e9l\u00e9charge \u00e0 <a href=\"https:\/\/unils-my.sharepoint.com\/:f:\/g\/personal\/tom_beucler_unil_ch\/Eg04Oty-ty9HquHde4fcxusBrpLneBC-P8ZvDIH0VDQ6WQ?e=NIhoqQ\">cet hyperlien<\/a> un projet ArcGIS qui contient d\u00e9j\u00e0 les couches n\u00e9cessaires ainsi qu&#8217;une mise en page finale de la carte que tu pourras utiliser pour l&#8217;habillage.<\/p>\n\n\n\n<p>1b) Ouvre le projet que tu viens de t\u00e9l\u00e9charger. Il contient une g\u00e9odatabase regroupant les trois couches raster que nous allons utiliser : <\/p>\n\n\n\n<ul class=\"wp-block-list\">\n<li>MapTicino1990_8classes.tif : occupation du sol en 1990 au Tessin avec une classification en 8 classes (&#8220;For\u00eat&#8221;, &#8220;Pr\u00e9s&#8221;, &#8220;Eau&#8221;, &#8220;Neige&#8221;, &#8220;Sols nus&#8221;, &#8220;Aires urbaines&#8221;, &#8220;Nuages&#8221;, &#8220;Ombres&#8221;) et calcul\u00e9es \u00e0 partir des donn\u00e9es de Landsat4.<\/li>\n\n\n\n<li>Landsat4_1990_194028_bx.tif : les 7 bandes spectrales individuelles collect\u00e9es par le satellite Landsat 4 de la NASA (b1: Bleu ; b2:Vert ; b3: Rouge ; b4: Proche infrarouge ; b5: Infrarouge moyen -1 ; b6: Thermique ; b7: Infrarouge moyen -2). <\/li>\n\n\n\n<li>MNT_Ticino.tif : un mod\u00e8le num\u00e9rique de terrain du Canton Tessin. Nota bene : ce n&#8217;est pas ce MNT qu&#8217;il faudra utiliser pour la carte finale ! Cette couche nous servira uniquement comme base de calcul. <\/li>\n<\/ul>\n\n\n\n<p>1c) V\u00e9rifie dans les propri\u00e9t\u00e9s de la carte que le <a href=\"https:\/\/pro.arcgis.com\/fr\/pro-app\/latest\/help\/mapping\/properties\/specify-a-coordinate-system.htm\" data-type=\"link\" data-id=\"https:\/\/pro.arcgis.com\/fr\/pro-app\/latest\/help\/mapping\/properties\/specify-a-coordinate-system.htm\">syst\u00e8me de r\u00e9f\u00e9rence<\/a> est bien CH1903+ LV95. Tu peux maintenant commencer le TP.<\/p>\n\n\n\n<p>1d) Prends le temps d&#8217;explorer les m\u00e9tadonn\u00e9es des couches (on parle de m\u00e9tadonn\u00e9es, mais il s&#8217;agit, dans ArcGis Pro, des propri\u00e9t\u00e9s de la couche) et r\u00e9ponds aux premi\u00e8res questions sur le quiz Moodle. Choisis bien les m\u00e9tadonn\u00e9es de la <strong>bande 4<\/strong> pour r\u00e9pondre aux questions. <\/p>\n\n\n\n<p>Pour la suite du quiz Moodle, nous allons travailler sur les <a href=\"https:\/\/pro.arcgis.com\/fr\/pro-app\/latest\/help\/data\/imagery\/raster-bands-pro-.htm\" data-type=\"link\" data-id=\"https:\/\/pro.arcgis.com\/fr\/pro-app\/latest\/help\/data\/imagery\/raster-bands-pro-.htm\">bandes<\/a> spectrales qui apparaissent comme 7 couches s\u00e9par\u00e9es. Afin d&#8217;effectuer les op\u00e9rations d&#8217;analyse d&#8217;image, il faudra les fusionner dans un fichier multibandes unique. Pour ce faire, il existe un outil nomm\u00e9 <a href=\"https:\/\/pro.arcgis.com\/fr\/pro-app\/latest\/tool-reference\/data-management\/composite-bands.htm\" data-type=\"link\" data-id=\"https:\/\/pro.arcgis.com\/fr\/pro-app\/latest\/tool-reference\/data-management\/composite-bands.htm\">Composite Bands<\/a>.<\/p>\n\n\n\n<p>1e) Une fois le menu de l&#8217;outil ouvert, ajoute les fichiers qui correspondent \u00e0 chaque bande <strong>en ordre croissant<\/strong> (de b1 \u00e0 b7).<\/p>\n\n\n\n<p>1f) Sauvegarde la nouvelle couche dans la g\u00e9odatabase et affiche le r\u00e9sultat. <\/p>\n\n\n\n<details class=\"wp-block-details is-layout-flow wp-block-details-is-layout-flow\"><summary>Solution &#8211; Outil composite band<\/summary>\n<div class=\"wp-block-group\"><div class=\"wp-block-group__inner-container is-layout-constrained wp-block-group-is-layout-constrained\">\n<p><\/p>\n\n\n\n<figure class=\"wp-block-gallery has-nested-images columns-default is-cropped wp-block-gallery-1 is-layout-flex wp-block-gallery-is-layout-flex\">\n<figure class=\"wp-block-image size-large\"><img alt=\"\" loading=\"lazy\" decoding=\"async\" width=\"600\" height=\"400\" data-id=\"1957\" src=\"https:\/\/wp.unil.ch\/dawn\/files\/2022\/11\/Composite-Band.gif\" alt=\"\" class=\"wp-image-1957\" \/><\/figure>\n<\/figure>\n<\/div><\/div>\n<\/details>\n<\/div><\/div>\n<\/details>\n\n\n\n<details class=\"wp-block-details is-layout-flow wp-block-details-is-layout-flow\"><summary>2. Analyse des bandes spectrales<\/summary>\n<div class=\"wp-block-group\"><div class=\"wp-block-group__inner-container is-layout-constrained wp-block-group-is-layout-constrained\">\n<p>2a) Pour commencer, il faut uniquement visualiser la couche multibandes que tu viens de cr\u00e9er, en d\u00e9cochant toutes les autres couches depuis la fen\u00eatre <a href=\"https:\/\/pro.arcgis.com\/fr\/pro-app\/latest\/help\/mapping\/map-authoring\/contents-pane.htm\" data-type=\"link\" data-id=\"https:\/\/pro.arcgis.com\/fr\/pro-app\/latest\/help\/mapping\/map-authoring\/contents-pane.htm\">Contents<\/a>. Cela \u00e9vitera qu&#8217;ArcGis ne doive rafra\u00eechir toutes les couches \u00e0 chaque fois que tu bouges sur la carte.<\/p>\n\n\n\n<p>2b) Ouvre le menu de symbologie de la couche multibandes et essaye les diff\u00e9rentes combinaisons de bandes (&#8220;Band&#8221;) dans les canaux (&#8220;Channel&#8221;) correspondant aux couleurs rouge, vert et bleu. <\/p>\n\n\n\n<p>2c) Apr\u00e8s t&#8217;\u00eatre familiaris\u00e9.e avec le fonctionnement des bandes, essaye de reproduire les compositions du tableau ci-dessous, qui indique la composition recommand\u00e9e pour les images Landsat.<\/p>\n\n\n\n<figure class=\"wp-block-gallery has-nested-images columns-default is-cropped wp-block-gallery-2 is-layout-flex wp-block-gallery-is-layout-flex\">\n<figure class=\"wp-block-image size-large\"><img alt=\"\" loading=\"lazy\" decoding=\"async\" width=\"1024\" height=\"285\" data-id=\"1948\" src=\"https:\/\/wp.unil.ch\/dawn\/files\/2022\/11\/Schermata-2022-11-11-alle-16.22.28-1024x285.png\" alt=\"\" class=\"wp-image-1948\" srcset=\"https:\/\/wp.unil.ch\/dawn\/files\/2022\/11\/Schermata-2022-11-11-alle-16.22.28-1024x285.png 1024w, https:\/\/wp.unil.ch\/dawn\/files\/2022\/11\/Schermata-2022-11-11-alle-16.22.28-300x83.png 300w, https:\/\/wp.unil.ch\/dawn\/files\/2022\/11\/Schermata-2022-11-11-alle-16.22.28-768x213.png 768w, https:\/\/wp.unil.ch\/dawn\/files\/2022\/11\/Schermata-2022-11-11-alle-16.22.28-1200x334.png 1200w, https:\/\/wp.unil.ch\/dawn\/files\/2022\/11\/Schermata-2022-11-11-alle-16.22.28.png 1414w\" sizes=\"auto, (max-width: 1024px) 100vw, 1024px\" \/><\/figure>\n<\/figure>\n\n\n\n<p>Des op\u00e9rations peuvent \u00eatre effectu\u00e9es sur les diff\u00e9rentes bandes de fa\u00e7on \u00e0 mettre en \u00e9vidence certains \u00e9l\u00e9ments. Dans le menu &#8220;Analysis&#8221; en haut de la page, tu trouveras l&#8217;outil &#8220;<a href=\"https:\/\/pro.arcgis.com\/fr\/pro-app\/latest\/help\/analysis\/raster-functions\/raster-functions.htm\">Raster Functions<\/a>&#8221; qui te permet de faire diff\u00e9rentes analyse bas\u00e9es sur les images multibandes. <\/p>\n\n\n\n<p>Une des fonctions qui peuvent \u00eatre appliqu\u00e9es \u00e0 un raster multibandes est le calcul du NDVI. Cet index de v\u00e9g\u00e9tation normalis\u00e9 (&#8220;Normalized Difference Vegetation Index&#8221;) met en \u00e9vidence la couverture v\u00e9g\u00e9tale sur le territoire. Il est tr\u00e8s utilis\u00e9 en agriculture et sylviculture pour conna\u00eetre l&#8217;\u00e9tat de sant\u00e9 des plantes, comme le montre l&#8217;image ci-dessous. <\/p>\n\n\n\n<figure class=\"wp-block-gallery has-nested-images columns-default is-cropped wp-block-gallery-3 is-layout-flex wp-block-gallery-is-layout-flex\">\n<figure class=\"wp-block-image size-large\"><img alt=\"\" loading=\"lazy\" decoding=\"async\" width=\"580\" height=\"305\" data-id=\"1950\" src=\"https:\/\/wp.unil.ch\/dawn\/files\/2022\/11\/plants.jpg\" alt=\"\" class=\"wp-image-1950\" srcset=\"https:\/\/wp.unil.ch\/dawn\/files\/2022\/11\/plants.jpg 580w, https:\/\/wp.unil.ch\/dawn\/files\/2022\/11\/plants-300x158.jpg 300w\" sizes=\"auto, (max-width: 580px) 100vw, 580px\" \/><\/figure>\n<\/figure>\n\n\n\n<p>Source \u00e0 consulter pour plus d&#8217;informations : <a href=\"https:\/\/eos.com\/blog\/ndvi-faq-all-you-need-to-know-about-ndvi\/\">https:\/\/eos.com\/blog\/ndvi-faq-all-you-need-to-know-about-ndvi\/<\/a> <\/p>\n\n\n\n<p>Sur ArcGIS le calcul de cet indicateur est automatique, et effectu\u00e9 gr\u00e2ce \u00e0 la formule ci-dessous : <\/p>\n\n\n\n<figure class=\"wp-block-image size-full\"><img alt=\"\" loading=\"lazy\" decoding=\"async\" width=\"980\" height=\"328\" src=\"https:\/\/wp.unil.ch\/dawn\/files\/2022\/11\/image.png\" alt=\"\" class=\"wp-image-1978\" srcset=\"https:\/\/wp.unil.ch\/dawn\/files\/2022\/11\/image.png 980w, https:\/\/wp.unil.ch\/dawn\/files\/2022\/11\/image-300x100.png 300w, https:\/\/wp.unil.ch\/dawn\/files\/2022\/11\/image-768x257.png 768w\" sizes=\"auto, (max-width: 980px) 100vw, 980px\" \/><\/figure>\n\n\n\n<p>O\u00f9 R d\u00e9signe la r\u00e9flectance spectrale dans la bande rouge (la bande 3 dans le TP) et PIR (ou NIR, Near InfraRed en anglais) indique la r\u00e9flectance spectrale dans la bande du Proche-Infrarouge (la bande 4). <\/p>\n\n\n\n<p>2d) Applique la <a href=\"https:\/\/pro.arcgis.com\/fr\/pro-app\/latest\/help\/analysis\/raster-functions\/ndvi-function.htm\" data-type=\"link\" data-id=\"https:\/\/pro.arcgis.com\/fr\/pro-app\/latest\/help\/analysis\/raster-functions\/ndvi-function.htm\">fonction NDVI<\/a> (en noir et blanc) et dans les propri\u00e9t\u00e9s de l&#8217;op\u00e9ration, indique comme bande visible la bande 3 (correspondant au rouge) et comme bande infrarouge la bande 4 (correspondant au proche infrarouge). Assure-toi de s\u00e9lectionner l&#8217;option &#8220;Scientific Output&#8221; qui permet de maintenir les valeurs du NDVI scientifiquement correctes, comme d\u00e9crit par la formule ci-dessus (entre -1 et 1, sinon les valeurs sont comprises entre 0 et 200, plus d&#8217;infos sur le site d&#8217;<a href=\"https:\/\/pro.arcgis.com\/en\/pro-app\/latest\/help\/analysis\/raster-functions\/ndvi-function.htm#:~:text=This%20scientific%20index%20outputs%20values,from%20rock%20and%20bare%20soil.\">ArcGis Pro<\/a>). <\/p>\n\n\n\n<p><span style=\"text-decoration: underline\">Remarque<\/span> : Il est important de noter qu&#8217;apr\u00e8s avoir s\u00e9lectionn\u00e9 le rendu en noir et blanc, l&#8217;image appara\u00eet en couleurs. Ceci est normal et est fait pour faciliter la visualisation du ph\u00e9nom\u00e8ne sans toutefois compromettre les valeurs calcul\u00e9es \u00e0 partir du NDVI. Ce n&#8217;est, en revanche, pas le cas avec l&#8217;option &#8220;colorized NDVI&#8221; qui, elle, n&#8217;offre pas les valeurs scientifiques du NDVI (Elle affiche les valeurs de 0-255, plus d&#8217;infos sur le site d&#8217;<a href=\"https:\/\/pro.arcgis.com\/en\/pro-app\/latest\/help\/analysis\/raster-functions\/ndvi-colorized-function.htm\">ArcGis Pro<\/a>).<\/p>\n\n\n\n<p>2e) Explore d\u00e9sormais les autres rendus possibles dans le menu des fonctions raster, et r\u00e9ponds aux questions de la deuxi\u00e8me page du quiz Moodle.<\/p>\n\n\n\n<details class=\"wp-block-details is-layout-flow wp-block-details-is-layout-flow\"><summary>Solution NDVI<\/summary>\n<div class=\"wp-block-group\"><div class=\"wp-block-group__inner-container is-layout-constrained wp-block-group-is-layout-constrained\">\n<figure class=\"wp-block-gallery has-nested-images columns-default is-cropped wp-block-gallery-4 is-layout-flex wp-block-gallery-is-layout-flex\">\n<figure class=\"wp-block-image size-large\"><img alt=\"\" loading=\"lazy\" decoding=\"async\" width=\"600\" height=\"400\" data-id=\"1958\" src=\"https:\/\/wp.unil.ch\/dawn\/files\/2022\/11\/NDVI.gif\" alt=\"\" class=\"wp-image-1958\" \/><\/figure>\n<\/figure>\n<\/div><\/div>\n<\/details>\n\n\n\n<p>2f) Sauvegarde bien la nouvelle couche du NDVI dans ta g\u00e9odatabase ! Tu devras utiliser ce rendu \u00e0 la fin du travail.<\/p>\n<\/div><\/div>\n<\/details>\n\n\n\n<details class=\"wp-block-details is-layout-flow wp-block-details-is-layout-flow\"><summary>3. Correction du raster<\/summary>\n<div class=\"wp-block-group\"><div class=\"wp-block-group__inner-container is-layout-constrained wp-block-group-is-layout-constrained\">\n<p>Le raster &#8220;MapTicino1990_8classes.tif&#8221; est une cartographie de l&#8217;occupation du sol au Tessin en 1990. Les donn\u00e9es sont recueillies par le satellite Landsat 4 de la NASA, et ensuite analys\u00e9es de sorte \u00e0 cr\u00e9er 8 classes d&#8217;occupation du sol : &#8220;For\u00eat&#8221;, &#8220;Pr\u00e9s&#8221;, &#8220;Eau&#8221;, &#8220;Neige&#8221;, &#8220;Sols nus&#8221;, &#8220;Aires urbaines&#8221;, &#8220;Nuages&#8221;, &#8220;Ombres&#8221;. <\/p>\n\n\n\n<p>La fonction qui permet de faire cette diff\u00e9renciation en se basant uniquement sur les donn\u00e9es satellitaires (les 7 bandes spectrales vues plus haut) est par contre tr\u00e8s sensible et engendre des erreurs. <\/p>\n\n\n\n<p>L&#8217;erreur principale qu&#8217;on retrouve est due \u00e0 la similitude entre les longueurs d&#8217;ondes \u00e9mises par les sols nus et les sols anthropiques (les zones urbaines). Ainsi, on retrouve des pixels class\u00e9s comme &#8220;aires urbaines&#8221; en haute montagne \u00e0 la place d&#8217;une classification comme &#8220;sols nus&#8221;. <\/p>\n\n\n\n<p>Pour corriger ces erreurs, il existe dans ArcGIS un outil permettant d&#8217;effectuer des calculs et des modifications d&#8217;une image raster : l&#8217;outil &#8220;Raster Calculator&#8221; (dans &#8220;Spatial Analyst Tools&#8221;).<\/p>\n\n\n\n<p>3a) Cherche l&#8217;outil &#8220;Raster Calculator&#8221; dans le paquet d&#8217;outil &#8220;Spatial Analyst Tools&#8221;. <\/p>\n\n\n\n<p>\u26a0\ufe0f Fais attention \u00e0 ne pas le confondre avec un autre outil du m\u00eame nom mais contenu dans un autre paquet d&#8217;outils. <\/p>\n\n\n\n<p>3b) Ouvre l&#8217;outil en question et essaye de trouver l&#8217;expression qui te permet de modifier les valeurs erron\u00e9es dont on a parl\u00e9 avant (lis d&#8217;abord la suite avant de t&#8217;y attaquer). En pratique, il faudra cr\u00e9er une requ\u00eate permettant de modifier l&#8217;affectation des pixels class\u00e9s comme &#8220;zone urbaine&#8221; et situ\u00e9s au-dessus de 1&#8217;400 m\u00e8tres d&#8217;altitude en les classants comme &#8220;sols nus&#8221;. <\/p>\n\n\n\n<p><strong>\u26a0\ufe0f <\/strong> Utilise les variables ainsi que les op\u00e9rateurs offerts par l&#8217;outil. Tu pourrais \u00e9crire toi-m\u00eame le tout, mais il arrive que les symboles ne soient pas identiques, ce qui pourrait engendrer une erreur. En outre, si tu fais une erreur minime dans l&#8217;\u00e9criture des couches, l&#8217;outil produira un autre message d&#8217;erreur. <\/p>\n\n\n\n<p>Voici quelques exemples de requ\u00eates qui peuvent \u00eatre men\u00e9es : <\/p>\n\n\n\n<ul class=\"wp-block-list\">\n<li><em>(&#8220;MapTicino1990_8classes&#8221; == 4)<\/em> : s\u00e9lection de tous les pixels de la couche &#8220;MapTicino1990_8classes&#8221; ayant une valeur \u00e9gale \u00e0 4 (c.\u00e0.d. qui appartiennent \u00e0 la classe &#8220;neige&#8221;). <\/li>\n\n\n\n<li><em>(&#8220;MNT25_Ticino&#8221; &lt; 600)<\/em> : s\u00e9lection de tous les pixels situ\u00e9s \u00e0 une altitude inf\u00e9rieure \u00e0 600 m\u00e8tres dans le MNT &#8220;MNT25_Ticino&#8221;. <\/li>\n<\/ul>\n\n\n\n<p>Les deux requ\u00eates ci-dessus donneront comme r\u00e9sultat un nombre, puisqu&#8217;elles indiqueront combien de pixels r\u00e9pondent \u00e0 la requ\u00eate en question. Elles donnent aussi la position de ces pixels. <\/p>\n\n\n\n<p>Il existe aussi des fonctions conditionnelles, qui s&#8217;appliquent uniquement dans le cas o\u00f9 une condition pr\u00e9alablement fix\u00e9e est respect\u00e9e. <\/p>\n\n\n\n<ul class=\"wp-block-list\">\n<li><em>Con ( condition, valeur-si-vrai, valeur-si-faux)<\/em> : cette requ\u00eate appliquera la valeur &#8220;valeur-si-vrai&#8221; dans le cas o\u00f9 la condition initiale est respect\u00e9e ; et appliquera la valeur &#8220;valeur-si-faux&#8221; dans le cas o\u00f9 la condition initiale n&#8217;est pas respect\u00e9e. Il fonctionne de la m\u00eame mani\u00e8re que le &#8220;=SI()&#8221; dans Excel.<\/li>\n\n\n\n<li><em>Con ( (&#8220;MNT25_Ticino&#8221; &gt; 600) , 600 , &#8220;MNT25_Ticino&#8221;)<\/em> : cette requ\u00eate d\u00e9finit un plafond d&#8217;altitude \u00e0 600 m\u00e8tres. La condition est &#8220;si le MNT25_Ticino pr\u00e9sente une valeur <strong>sup\u00e9rieure<\/strong> \u00e0 600 m\u00e8tres d&#8217;altitude, assigne la valeur de 600m&#8221; ; &#8220;si le MNT25_Ticino pr\u00e9sente une valeur <strong>inf\u00e9rieure <\/strong> \u00e0 600 m\u00e8tres d&#8217;altitude, laisse la valeur d&#8217;origine&#8221;. <\/li>\n<\/ul>\n\n\n\n<p>Les requ\u00eates ci-dessus ont une seule condition, mais on peut tr\u00e8s bien indiquer plusieurs conditions \u00e0 remplir avec le symbole &#8220;&amp;&#8221;.<\/p>\n\n\n\n<p><em>Con ( ( &#8220;MNT25_Ticino&#8221; &gt; 600) &amp; (&#8220;MapTicino1990_8classes&#8221; == 7) , 4 , &#8220;MapTicino1990_8classes&#8221; )<\/em> : cette requ\u00eate modifie la classification des pixels 7 (&#8220;ombres&#8221;) en pixels 4 (&#8220;neige&#8221;) si l&#8217;altitude est sup\u00e9rieure \u00e0 600 m\u00e8tres. Donc : &#8220;si le MNT indique une valeur sup\u00e9rieure \u00e0 600 m, et que le pixel appartient \u00e0 la classes 7 &#8220;Ombres&#8221;, alors assigne la valeur du pixel \u00e0 4 (&#8220;neige&#8221;). Si une des condition n&#8217;est pas remplie, laisse la valeur d&#8217;origine&#8221;. <\/p>\n\n\n\n<p>\u26a0\ufe0f C&#8217;est une requ\u00eate de ce type qui te permettra d&#8217;effectuer le reclassement des pixels ayant la valeur 6 &#8220;Aires urbaines&#8221; et situ\u00e9s au-dessus de 1&#8217;400 m\u00e8tres d&#8217;altitude en pixels avec valeur 5 &#8220;Sols nus&#8221;.<\/p>\n\n\n\n<p>3c) Effectue la requ\u00eate, puis, une fois la requ\u00eate effectu\u00e9e, sauvegarde le r\u00e9sultat dans une nouvelle couche que tu nommeras avec la mention &#8220;ReClass&#8221; (ex. ReClass_MapTicino1990_8classes). <\/p>\n\n\n\n<p>Il se peut qu&#8217;il ne te sorte que 2 classes en sortie au lieu des 8: dans ce cas, relance l&#8217;outil (en changeant le nom dans &#8220;Output raster&#8221;) jusqu&#8217;\u00e0 obtenir le total des classes. Garde cette couche dans ta g\u00e9odatabase car elle te sera redemand\u00e9e plus tard pour le rendu final. <\/p>\n\n\n\n<details class=\"wp-block-details is-layout-flow wp-block-details-is-layout-flow\"><summary>Solution reclassement<\/summary>\n<div class=\"wp-block-group\"><div class=\"wp-block-group__inner-container is-layout-constrained wp-block-group-is-layout-constrained\">\n<p>Si tu ne visualises pas la vid\u00e9o ci-dessous, active le mode plein \u00e9cran si tu es sur MacOS, ou tourne ton appareil si tu es sur iOS ou iPadOS. La vid\u00e9o est toujours visible depuis Windows (sur les machines virtuelles).<\/p>\n\n\n\n<figure class=\"wp-block-video\"><video height=\"480\" style=\"aspect-ratio: 720 \/ 480;\" width=\"720\" controls src=\"https:\/\/wp.unil.ch\/dawn\/files\/2022\/11\/ReCLass.mp4\"><\/video><\/figure>\n\n\n\n<p>\u00c9tant donn\u00e9 la qualit\u00e9 de la vid\u00e9o, on t&#8217;offre un zoom sur la condition qui a \u00e9t\u00e9 utilis\u00e9e. Essaye toutefois de l&#8217;\u00e9crire toi-m\u00eame avant de regarder la solution.<\/p>\n\n\n\n<details class=\"wp-block-details is-layout-flow wp-block-details-is-layout-flow\"><summary>Solution Requ\u00eate<\/summary>\n<div class=\"wp-block-group\"><div class=\"wp-block-group__inner-container is-layout-constrained wp-block-group-is-layout-constrained\">\n<figure class=\"wp-block-gallery has-nested-images columns-default is-cropped wp-block-gallery-5 is-layout-flex wp-block-gallery-is-layout-flex\">\n<figure class=\"wp-block-image size-large\"><img alt=\"\" loading=\"lazy\" decoding=\"async\" width=\"911\" height=\"411\" data-id=\"1962\" src=\"https:\/\/wp.unil.ch\/dawn\/files\/2022\/11\/Schermata-2022-11-13-alle-14.13.27.png\" alt=\"\" class=\"wp-image-1962\" srcset=\"https:\/\/wp.unil.ch\/dawn\/files\/2022\/11\/Schermata-2022-11-13-alle-14.13.27.png 911w, https:\/\/wp.unil.ch\/dawn\/files\/2022\/11\/Schermata-2022-11-13-alle-14.13.27-300x135.png 300w, https:\/\/wp.unil.ch\/dawn\/files\/2022\/11\/Schermata-2022-11-13-alle-14.13.27-768x346.png 768w\" sizes=\"auto, (max-width: 911px) 100vw, 911px\" \/><\/figure>\n<\/figure>\n<\/div><\/div>\n<\/details>\n<\/div><\/div>\n<\/details>\n<\/div><\/div>\n<\/details>\n\n\n\n<details class=\"wp-block-details is-layout-flow wp-block-details-is-layout-flow\"><summary>4. Interpolation<\/summary>\n<div class=\"wp-block-group\"><div class=\"wp-block-group__inner-container is-layout-constrained wp-block-group-is-layout-constrained\">\n<p>Bravo, tu as presque fini&nbsp;! Dans cette derni\u00e8re partie, tu cr\u00e9eras un MNT \u00e0 partir de donn\u00e9es ponctuelles disponibles sur le site de l\u2019Office f\u00e9d\u00e9ral de topographie (Swisstopo). Il s\u2019agit du mod\u00e8le num\u00e9rique de base du terrain de la Suisse, utilis\u00e9 pour la production du MNT avec une maille de 25 m. Pour ce faire, on se basera sur une m\u00e9thode d\u00e9terministe classique, \u00e0 savoir l\u2019interpolation par moyenne mobile pond\u00e9r\u00e9e ou <a href=\"https:\/\/pro.arcgis.com\/fr\/pro-app\/2.9\/tool-reference\/spatial-analyst\/idw.htm\"><em>Inverse Distance Weighting<\/em> (IDW)<\/a>.<\/p>\n\n\n\n<p>Les ph\u00e9nom\u00e8nes spatio-continus sont d\u00e9finis en tout point de l\u2019espace g\u00e9ographique (ex. l\u2019altitude et la temp\u00e9rature) mais sont g\u00e9n\u00e9ralement \u00e9tudi\u00e9s \u00e0 travers des donn\u00e9es ponctuelles. Entre les points d\u2019\u00e9chantillonnage, les valeurs de ces ph\u00e9nom\u00e8nes ne sont pas mesur\u00e9es. L\u2019objectif des m\u00e9thodes d\u2019interpolation consiste \u00e0 pr\u00e9dire ces valeurs inconnues sur la base de l\u2019autocorr\u00e9lation spatiale : <\/p>\n\n\n\n<p>\u00ab <em>Deux objets proches ont plus de chance [d\u2019interagir] que deux objets \u00e9loign\u00e9s<\/em> \u00bb (<a href=\"https:\/\/support.esri.com\/fr-fr\/gis-dictionary\/tobler-s-first-law-of-geography\">premi\u00e8re loi de la g\u00e9ographie<\/a> de Waldo Tobler). <\/p>\n\n\n\n<p>Pour qu\u2019une mod\u00e9lisation soit satisfaisante, il est primordial qu\u2019elle soit bas\u00e9e sur une analyse exploratoire des donn\u00e9es et sur une analyse des erreurs (quelle que soit la m\u00e9thode d\u2019interpolation choisie).<\/p>\n\n\n\n<p>Le but de ce TP est d\u2019interpoler assez rapidement (au d\u00e9triment de la qualit\u00e9 du r\u00e9sultat) une surface raster \u00e0 partir de donn\u00e9es ponctuelles de l\u2019altitude.<\/p>\n\n\n\n<p>4a) Pour t\u00e9l\u00e9charger les donn\u00e9es, rends-toi sur le site de swisstopo et appuie sur \u00ab&nbsp;<a href=\"https:\/\/cms.geo.admin.ch\/ogd\/topography\/DHM25_BM_SHP.zip\">DHM25 &#8211; Mod\u00e8le de base ESRI Shapefile&nbsp;<\/a>\u00bb qui se trouve sous l\u2019onglet G\u00e9odonn\u00e9es et applications &gt; Mod\u00e8les d&#8217;altitude &gt; MNT25. <\/p>\n\n\n\n<p>4b) Ensuite, importe le fichier shape \u00ab&nbsp;dhm25_p&nbsp;\u00bb dans ton projet (sans oublier de <a href=\"https:\/\/support.microsoft.com\/fr-fr\/windows\/compresser-et-d%C3%A9compresser-des-fichiers-f6dde0a7-0fec-8294-e1d3-703ed85e7ebc\" data-type=\"URL\" data-id=\"https:\/\/support.microsoft.com\/fr-fr\/windows\/compresser-et-d%C3%A9compresser-des-fichiers-f6dde0a7-0fec-8294-e1d3-703ed85e7ebc\">d\u00e9zipper<\/a> le fichier, deux dossiers en seront extraits).<\/p>\n\n\n\n<p>Pour pr\u00e9dire l\u2019altitude en tout point de l\u2019espace g\u00e9ographique Suisse on utilisera l&#8217;outil&nbsp;<a href=\"https:\/\/pro.arcgis.com\/fr\/pro-app\/2.9\/tool-reference\/spatial-analyst\/idw.htm\" data-type=\"URL\" data-id=\"https:\/\/pro.arcgis.com\/fr\/pro-app\/2.9\/tool-reference\/spatial-analyst\/idw.htm\">IDW<\/a>. <\/p>\n\n\n\n<p>4c) Cherche cet outil dans la barre <a href=\"https:\/\/pro.arcgis.com\/fr\/pro-app\/latest\/help\/analysis\/geoprocessing\/basics\/find-geoprocessing-tools.htm\" data-type=\"URL\" data-id=\"https:\/\/pro.arcgis.com\/fr\/pro-app\/latest\/help\/analysis\/geoprocessing\/basics\/find-geoprocessing-tools.htm\">Find Tools<\/a> ou trouve-le en cliquant Tools &gt; Geostatistical Analyst Tools &gt; Interpolation &gt; IDW. <\/p>\n\n\n\n<p>4d) Gr\u00e2ce \u00e0 cet outil, estime le MNT \u00e0 partir des donn\u00e9es Swisstopo.<\/p>\n\n\n\n<p>\u26a0\ufe0f Fais bien attention \u00e0 choisir le bon champ \u00e0 interpoler, i.e., l&#8217;\u00e9l\u00e9vation Shape.Z et non le code identifiant la commune OBJECTID dans &#8220;Z value field&#8221;. Il se peut que l&#8217;outil affiche une erreur (un th\u00e8me r\u00e9current avec ce logiciel), dans ce cas : sauvegarde le projet, quitte ArcGis Pro et relance le projet.<\/p>\n\n\n\n<details class=\"wp-block-details is-layout-flow wp-block-details-is-layout-flow\"><summary>Solution<\/summary>\n<div class=\"wp-block-group\"><div class=\"wp-block-group__inner-container is-layout-constrained wp-block-group-is-layout-constrained\">\n<figure class=\"wp-block-image size-full\"><img alt=\"\" loading=\"lazy\" decoding=\"async\" width=\"1280\" height=\"686\" src=\"https:\/\/wp.unil.ch\/dawn\/files\/2022\/11\/Enregistrement-20221115_180606.gif\" alt=\"\" class=\"wp-image-2018\" \/><\/figure>\n<\/div><\/div>\n<\/details>\n\n\n\n<p>4e) Finalement, g\u00e9n\u00e8re une repr\u00e9sentation 3D de la surface du MNT gr\u00e2ce \u00e0 l\u2019outil <em>Hillshade<\/em> (<a href=\"https:\/\/pro.arcgis.com\/fr\/pro-app\/latest\/help\/analysis\/raster-functions\/hillshade-function.htm\" data-type=\"URL\" data-id=\"https:\/\/pro.arcgis.com\/fr\/pro-app\/latest\/help\/analysis\/raster-functions\/hillshade-function.htm\">Ombrage<\/a>) dans Raster Function.<\/p>\n\n\n\n<p>\u2554\u2550\u2557<br>\u2551\u2550\u256c\u2566\u2566\u2550\u2566\u2550\u2566\u2550\u2557  \u2551<br>\u2560\u2550\u2551\u2551\u2551\u256c\u2551\u2569\u2563\u2554\u255d  \u2551<br>\u255a\u2550\u2569\u2550\u2563\u2554\u2569\u2550\u2569\u255d. .  \u2550<\/p>\n<\/div><\/div>\n\n\n\n<div class=\"wp-block-group\"><div class=\"wp-block-group__inner-container is-layout-constrained wp-block-group-is-layout-constrained\"><\/div><\/div>\n<\/details>\n\n\n\n<details class=\"wp-block-details is-layout-flow wp-block-details-is-layout-flow\"><summary>5. Cr\u00e9er les trois cartes de r\u00e9sultats<\/summary>\n<p>Il ne te reste plus qu&#8217;\u00e0 rendre tes r\u00e9sultats sur Moodle. Pour ce faire, nous avons pr\u00e9par\u00e9 trois layouts : <strong>NDVI<\/strong>, <strong>occup_sol<\/strong>, et <strong>IDW<\/strong>, avec lesquels tu peux exporter les rasters respectifs au format .pdf.<\/p>\n\n\n\n<p>5a) Pour chaque layout, fais correspondre l\u2019information pr\u00e9sente sur la carte \u00e0 son titre sans oublier le fond de carte.&nbsp;<br>Conseil : Pour modifier la couche affich\u00e9e dans une <a href=\"https:\/\/pro.arcgis.com\/fr\/pro-app\/latest\/help\/layouts\/layouts-in-arcgis-pro.htm\">mises en page<\/a> (Layout), il te suffit d\u2019activer la carte et s\u00e9lectionner, d\u00e9s\u00e9lectionner ou ajouter (Map &gt; Add Data) les bonnes couches dans la <a href=\"https:\/\/pro.arcgis.com\/fr\/pro-app\/latest\/help\/mapping\/map-authoring\/contents-pane.htm\">fen\u00eatre&nbsp;contenu<\/a> (Contents).<\/p>\n\n\n\n<p>5b) Pour la carte concernant l\u2019occupation du sol, ins\u00e8re une <a href=\"https:\/\/pro.arcgis.com\/fr\/pro-app\/latest\/help\/layouts\/add-a-legend.htm\">l\u00e9gende<\/a> qui illustre les 8 classes.<\/p>\n\n\n\n<p>5c) Finalement, pour chaque carte, indique l\u2019auteur.e en utilisant ton nom, ton pr\u00e9nom, ainsi que les <strong>sources des donn\u00e9es<\/strong>.<\/p>\n<\/details>\n\n\n\n<details class=\"wp-block-details is-layout-flow wp-block-details-is-layout-flow\"><summary>6. Rendus et paquetage du projet<\/summary>\n<div class=\"wp-block-group\"><div class=\"wp-block-group__inner-container is-layout-constrained wp-block-group-is-layout-constrained\">\n<p>6a) Tu peux d&#8217;ores et d\u00e9j\u00e0 soumettre les trois rasters au format .pdf sur Moodle: <a href=\"https:\/\/moodle.unil.ch\/mod\/assign\/view.php?id=1736952\">Rendus_Cartes_TP7<\/a>. N&#8217;oublie pas le format du rendu : <strong><em>nom_prenom_cartes_TP7.pdf<\/em><\/strong> (tu peux consulter la grille d&#8217;\u00e9valuation sur Moodle).<\/p>\n\n\n\n<p>6b) Cr\u00e9e un paquetage de projet \u00e0 partager en suivant les instructions <a href=\"https:\/\/pro.arcgis.com\/fr\/pro-app\/latest\/help\/sharing\/overview\/project-package.htm\">\u00e0 ce lien<\/a>.<\/p>\n\n\n\n<p>N\u2019oublie pas d\u2019inclure un <strong>r\u00e9sum\u00e9 <\/strong>(<em>summary<\/em>) et des <strong>balises <\/strong>ou mots-cl\u00e9s (<em>tags<\/em>) contenant les m\u00e9tadonn\u00e9es identifi\u00e9es pr\u00e9c\u00e9demment. Tu peux cliquer la case \u201c<em>Share outside of organization<\/em>\u201d et le bouton \u201c<em>Analyze<\/em>\u201d pour optimiser ton paquetage pour un partage public.<\/p>\n\n\n\n<p>Il se peut que l&#8217;historique (History Items) de ton projet cr\u00e9e une erreur \u00e0 l&#8217;analyse ou pendant le paquetage, tu peux d\u00e9s\u00e9lectionner cette option si c&#8217;est le cas.<\/p>\n\n\n\n<details class=\"wp-block-details is-layout-flow wp-block-details-is-layout-flow\"><summary>Solution<\/summary>\n<div class=\"wp-block-group\"><div class=\"wp-block-group__inner-container is-layout-constrained wp-block-group-is-layout-constrained\">\n<figure class=\"wp-block-image size-full\"><img alt=\"\" loading=\"lazy\" decoding=\"async\" width=\"1915\" height=\"898\" src=\"https:\/\/wp.unil.ch\/dawn\/files\/2022\/10\/E1_3.gif\" alt=\"\" class=\"wp-image-1675\" \/><\/figure>\n<\/div><\/div>\n<\/details>\n\n\n\n<p>6c) Une fois le projet paquet\u00e9, v\u00e9rifie que le fichier <em>.ppkx<\/em> que tu viens de cr\u00e9er te permet bien de r\u00e9ouvrir le projet avec la g\u00e9odatabase compl\u00e8te (par exemple en double-cliquant dessus).<\/p>\n\n\n\n<details class=\"wp-block-details is-layout-flow wp-block-details-is-layout-flow\"><summary>Solution<\/summary>\n<div class=\"wp-block-group\"><div class=\"wp-block-group__inner-container is-layout-constrained wp-block-group-is-layout-constrained\">\n<figure class=\"wp-block-image size-full\"><img alt=\"\" loading=\"lazy\" decoding=\"async\" width=\"1553\" height=\"782\" src=\"https:\/\/wp.unil.ch\/dawn\/files\/2022\/10\/E1_4.gif\" alt=\"\" class=\"wp-image-1676\" \/><\/figure>\n<\/div><\/div>\n<\/details>\n\n\n\n<p>6d) Copie le fichier <em>.ppkx<\/em> de la machine virtuelle sur ton OneDrive et cr\u00e9e un lien de partage.<\/p>\n\n\n\n<p>6e) Tu peux maintenant te rendre \u00e0 nouveau sur Moodle pour soumettre \u00e0 ce lien : <a href=\"https:\/\/moodle.unil.ch\/mod\/quiz\/view.php?id=1736953\">Rendus_projet_TP7<\/a>. N&#8217;oublie pas le format du rendu : <strong><em>nom_prenom_TP7.ppkx<\/em><\/strong><\/p>\n\n\n\n<p>F\u00e9licitations pour avoir termin\u00e9 le TP et \u00e0 la semaine prochaine pour l&#8217;analyse de distance !<\/p>\n\n\n\n<figure class=\"wp-block-image size-full\"><img alt=\"\" loading=\"lazy\" decoding=\"async\" width=\"1238\" height=\"500\" src=\"https:\/\/wp.unil.ch\/dawn\/files\/2023\/11\/gryffindor-students-cheering-5rpou67pjsvez4uf.gif\" alt=\"\" class=\"wp-image-2746\" \/><\/figure>\n<\/div><\/div>\n\n\n\n<p>&#8220;<a href=\"https:\/\/tenor.com\/fr\/view\/gryffindor-cheering-gold-and-red-harry-potter-hogwarts-house-gif-17712184\"><em>Gryffindor cheering<\/em><\/a>&#8221; par michellesgifs<\/p>\n<\/details>\n","protected":false},"excerpt":{"rendered":"","protected":false},"author":1002254,"featured_media":1975,"parent":0,"menu_order":0,"comment_status":"closed","ping_status":"closed","template":"templates\/template-cover.php","meta":{"_seopress_robots_primary_cat":"","_seopress_titles_title":"","_seopress_titles_desc":"","_seopress_robots_index":"","footnotes":""},"class_list":["post-1329","page","type-page","status-publish","has-post-thumbnail"],"_links":{"self":[{"href":"https:\/\/wp.unil.ch\/dawn\/wp-json\/wp\/v2\/pages\/1329","targetHints":{"allow":["GET"]}}],"collection":[{"href":"https:\/\/wp.unil.ch\/dawn\/wp-json\/wp\/v2\/pages"}],"about":[{"href":"https:\/\/wp.unil.ch\/dawn\/wp-json\/wp\/v2\/types\/page"}],"author":[{"embeddable":true,"href":"https:\/\/wp.unil.ch\/dawn\/wp-json\/wp\/v2\/users\/1002254"}],"replies":[{"embeddable":true,"href":"https:\/\/wp.unil.ch\/dawn\/wp-json\/wp\/v2\/comments?post=1329"}],"version-history":[{"count":4,"href":"https:\/\/wp.unil.ch\/dawn\/wp-json\/wp\/v2\/pages\/1329\/revisions"}],"predecessor-version":[{"id":3397,"href":"https:\/\/wp.unil.ch\/dawn\/wp-json\/wp\/v2\/pages\/1329\/revisions\/3397"}],"wp:featuredmedia":[{"embeddable":true,"href":"https:\/\/wp.unil.ch\/dawn\/wp-json\/wp\/v2\/media\/1975"}],"wp:attachment":[{"href":"https:\/\/wp.unil.ch\/dawn\/wp-json\/wp\/v2\/media?parent=1329"}],"curies":[{"name":"wp","href":"https:\/\/api.w.org\/{rel}","templated":true}]}}