Vincent GODARD

Département de Géographie

Université de Paris 8


V.1.1.2 - Dernière mise à jour : 06/03/2015

Fiche TD du cours n°4.2. de Télédétection niveau 2 :

Rehaussement et interprétation visuelle

 

5. Exercice d'application de la fiche mémo 4.2.

Réaliser une photo-interprétation, à l'aide d'IDRISI Andes, du parc Georges Valbon à La Courneuve, puis un traitement multivarié avec SPAD.

Les étapes sont les suivantes :

 

5.1. Choix de la composition colorée

5.1.1 Réalisation de 3 compositions colorées

Réalisez les 3 compositions colorées suivantes avec une saturation à 2,5%

- COMPOSIT, où :

tab. 1 - Tableau des compositions colorées

Couleurs \ Nom fichier output

spot6-5-07xs123

spot6-5-07xs321

H87TM354

bleu

spot6-5-07xs1

spot6-5-07xs3

etm_lacour_24-8-00_L2_TM3

vert

spot6-5-07xs2

spot6-5-07xs2

etm_lacour_24-8-00_L2_TM5

rouge

spot6-5-07xs3

spot6-5-07xs1

etm_lacour_24-8-00_L2_TM4

 

5.1.2 Comment choisir la "bonne" composition colorée

- Affichez par dessus le fichier vecteur ZEA11_a.vct

Télechargement des fichiers vecteur .vct et document .vdc

A stocker dans le Working folder

- Sélectionnez la composition colorée qui vous permet de continuer la vectorisation

- Affichez le fichier EXCEL ZEA11_a.xls

Télechargement du fichier tableur .xls

A stocker deans le Working folder

- Voir avec le Curseur d'interrogation (Cursor Inquiry Mode) quelle est la relation entre la photointerprétation (fichier vecteur) et le fichier EXCEL.

 

5.2 Étape de vectorisation des ZEA

5.2.1 Consignes de vectorisation des ZEA

Poursuite de la vectorisation de chaque thème d'occupation du sol sur le fichier ortho2003_u31n_cc en complétant simultanément votre tableau (feuille EXCEL ZEA11_a.xls ) et en utilisant CartaLinx :

- avec les variables de :

- couleurs => noir, brun-rouge, marron-vert, jaune, bleu, ... ;

C'est la couleur dominante !

- texture

=> lisse (pas de nuances de ton réellement perceptibles : effet miroir comme sur les plans d'eau, ...) ;

=> ponctuée (petites nuances de ton : variation de teintes dans le couvert végétal, différence de hauteurs d'eau, ...) ;

=> points grossiers (ronds, "choux-fleurs" ...).

- structure

=> homogène régulier : l'organisation des texture est régulière, pas de "trous" dans les tiretés ou les pontués ; il n'y a pas d'inclusions d'objets "étrangers" dans la ZEA.

=> homogène irrégulier : l'organisation des texture est irrégulière, il y a des "trous" dans les tiretés ou les pontués ; il n'y a pas ou très peu d'inclusions d'objets "étrangers" dans la ZEA.

ou irrégulière

=> hétérogène : l'échelle fait que l'on inclut des objets "étrangers" dans la ZEA.

C'est une suggestion. On peut faire plus détaillé !

A quelle échelle réaliser la vectorisation de la carte d'occupation du sol ?

La couverture du MOS est réalisée au 1/5 000 * avec des objets ne pouvant pas être plus petits que 625 m2 **. Cependant, comme on le voit sur la figure 1, la couche d'occupation du sol 2003 n'est pas très détaillée à l'intérieur du parc. Il doit être possible de faire mieux !

* HENIN (M.), AWADA (F.), FOULARD (S.) - 2001 - L’Occupation du Sol en Ile-de-France aujourd’hui, 1er bilan du MOS 1999. Notes rapide. Mode d'occupation du sol., num. 282, juillet, 1-6.

** BIASI (L. de) et NASCIMENT (I.) - 2000 - Éléments méthodologiques pour le suivi de la consommation de l’espace agricole et naturel Notes rapide sur l'occupation du sol., num. 2, février, 1-4.

fig. 1 - Comparatif de la BD Ortho de l'IGN et du MOS de l'IAU

Sources : BD Ortho 2003, IGN et MOS 2003, IAU

 On se contentera d'objets faisant au moins un hectare (environ 100 m x 100 m).

 

5.2.2 Paramétrer CartaLinx

- Lancez CartaLinx.

- Indiquez le chemin d'accès aux données :

[FILE]/[SET DEFAULT DATA PATH]. Indiquez où se trouve votre fichier image "La Courneuve".

- Convertir les images Idrisi en BMP géoréférencées :

[FILE]/ [IMAGE CONVERSIONS]. Sélectionnez ortho2003_u31n_cc

Donner un nom à l'image en .BMP qui servira de fond pour vectoriser.

On peut lui donner le même nom que l'image d'Idrisi (les formats étant différents, il n'y aura pas d'interférence), mais pour la différencier facilement ajouter : ortho2003_u31n_cc_clx.

- Affichez le fond à l'écran :

[FILE]/[NEW COVERAGE]/[COVERAGE BASED UPON BITMAP].

Entrez le nom de l'image en .BMP => ortho2003_u31n_cc_clx.bmp

L'image s'affiche !

Si le fond n'apparaissait pas, sélectionnez :

[FILE]/[PREFERENCES], onglet [BACKDROP], puis [BROWSE)] !

=> le fond voulu doit apparaître dans le répertoire de travail, le sélectionner.

- Paramétrez CartaLinx pour commencer la vectorisation :

[FILE]/[PREFERENCES],

- onglet [DISPLAY], cochez => tous les "Show..."

pour "Show map grid", portez les incréments à 100.

- onglet [TOLERANCES], cochez => tous les "Show..."

portez les "... tolerance" à 10.

 

5.2.3 La vectorisation sous CartaLinx

- Poursuite de la vectorisation de chaque thème d'occupation du sol sur le fichier ortho2003_u31n_cc en complétant simultanément votre tableau (feuille EXCEL ZEA11_a.xls ) et en utilisant CartaLinx.

Rappel des icônes sous CartaLinx => voir ici !

 

5.3 Création d'une image résultat de votre photo-interprétation

Une fois votre couche vectorielle Phot_int1 complète, il faut la rasteriser.

- Si vous n'avez pas fini, vous pouvez télécharger une vectorisation réalisée par votre camarade Vincent Choffrut. C'est un fichier zippé qui contient les 4 fichiers CartaLinx et le fichier Excel de cette vectorisation :

Il va falloir les dézipper et les stocker dans votre Working folder

- Quelque soit l'origine de votre vectorisation, il faut l'exporter en un fichier compatible avec Idrisi Andes.

[FILE]/[EXPORT],

- choisir Idrisi 32

Tous les fichiers exportés sont également à stocker dans le Working folder

- Chaînez les étapes à l'aide des fonctions suivantes :

- INITIAL : pour créer une carte vide qui va accueillir votre photo-interprétation

- Laissez cocher => Copy spatial parameters

- Output image => ZEA1

- Image to copy parameters => spot6-5-07xs123 (une des images qui sert à la photo-interprétation)

- Output documentation => mettez le titre de votre choix (cela va bientôt accueillir votre photo-interprétation)

Le reste par défaut !

- POLYRAS : pour rasteriser un fichier vecteur polygone

- Vector polygone file => Phot_int1

- Image file to be updated => ZEA1

Constatez que le thème zéro est celui que vous avez laissé par défaut en dernier, mais :

- il déborde éventuellement du parc da La Courneuve ;

- vous avez rempli un enregistrement de votre feuille d'enqête EXCEL !

Le cas échéant, il faut faire un masque qui englobe le parc et votre vectorisation.

- RECLASS : pour attribuer à la dernière ZEA (celle non vectorisée) un numéro différent de zéro

- Input file => ZEA1

- Output file => ZEA2

- Assign a new value of => n + 1 (par rapport à la dernière ZEA vectorisée)

- To all values from => 0

- To just less than => 1

Le reste par défaut !

- FILTER : pour harmoniser

- Filtre modal => combler les trous !

- POLYVEC : pour vecteuriser le fichier raster UEA2 reclassé et pouvoir le lier à la future base de données

- Output option

- Cocher => Idrisi polygon file

- Input file => ZEA2

- Output Idrisi polygone file => Pi2_vc_vg

Le fichier vecteur Phot_int2 est superposable au fichier raster ZEA2 (Composer / Add layer)

 

5.4 Regroupement des ZEA en UEA sous SPAD

- Sous EXCEL transformation des "x" en 1, le reste en 0

- Onglet Feuil1 (renommé en photo-int)

- remplacer la ZEA "0" par la valeur "n+1" utilisée précédemment dans RECLASS

- Onglet Feuil2 (renommé en codage)

- recopier les lignes de titre de la feuille photo-int (col. A et ligne 1)

- remplacer à l'aide de la formule suivante les "x" de la feuille photo-int par des 1 et le reste par des 0

Se mettre en B2 et taper =>

=SI('photo-int'!B2="x";1;0)

- recopier sur l'ensemble du tableau

- enregistrer le classeur

- puis enregistrer en format texte tabulé => codage.txt

- fermer sans réenregistrer

- Sous Windows

- Créez un répertoire temporaire "ZEA" dans Votre_Nom \ LaCourneuve

- Sous SPAD7 regroupement des ZEA en UEA

a) Import de la Base

- Projet / Nouveau Projet (pour paramétrer les chemins d'accès)

- indiquer Répertoire => "ZEA" dans Votre_Nom \ LaCourneuve \ ZEA

- indiquer Nom du projet => "ZEA"

Validez

- Faire glisser depuis la fenêtre Méthodes vers la fenêtre Diagramme :

- Imports / Exports des données / Imports / Feuille Excel

- Double cliquer sur la "Feuille Excel" pour la paramétrer

- Indiquer dans l'onglet :

- Import

le nom du fichier => ZEA11.xls

- Métadonnées

par un double clic dans Rôle pour "zea" => Identifiant

les autres restent à => Continue

La validation exécute l'import.

 

b) Traitement multivarié (Amado)

- Faire glisser depuis la fenêtre Méthodes vers la fenêtre Diagramme :

- Amado - Graphique de Bertin / Sélection des données d'une matrice pour Amado

- Lier les deux icônes en "cherchant" les triangles bleus puis en "tirant" un lien fléché !

- Double cliquer sur la "Amado" pour la paramétrer

- Sélectionner toutes les variables

- Valider

 

- Double-clic sur le graphique dans la filière

- Permutation des lignes et des colonnes jusqu'à obtenir des regroupements

- Traitement / Insérer des séparateurs (ou Ctrl+I)

Pour isoler visuellement les groupes

- Les placer entre des groupes homogènes

- Une fois le classement achevé, noter le groupe d'appartenance (UEA) de chaque individu (ZEA)

Car Amado n'exporte pas les séparateurs ni les numéro des groupes !

- Export vers IDRISI

- Fichier + Enregistrer sous ...

=> Classement.txt

 

5.5 Cartographie des UEA

- Sous EXCEL saisie des UEA d'appartenance pour chaque ZEA

- Ouverture du fichier issu d'Amado => Classement.txt

Attention au Type de fichier => Fichiers texte

- Onglet Classement

- réécrire ZEA en A1

- Insertion / Feuille

- Onglet Feuil2 (renommé en UEA)

- recopier la colonne A1 de la feuille classement

- écrire UEA en B1

- saisir en colonne B, en face de chaque ZEA son numéro d'UEA

- sélectionner les colonnes A et B

- Données / Trier

- Trier par ZEA croissantes

- ligne de titres => oui

- enregistrer le classeur => UEA.xls

Attention au Type de fichier => Classeur Microsoft Excel (d'avant 2007)

 

- Sous IDRISI regroupement des ZEA en UEA

a) Création de la base par import du .xls

- Lancez le Database Workshop (import du fichier dBASE)

- File / New

- lui donner un nom par enregistrer sous => UEA

- File / Import / Table / from External File

- Fichiers de type => Excel (*.xls)

- Select sheets (choisir la feuille)

Source List => UEA11

Destination List => UEA11

- Import the first row of "UEA11" as field names? => YES

 

b) Établissement du lien entre la Base et le fichier vecteur

- Database Workshop (Établissement du lien)

- Query / Establish Display Link

- Vector link file (*.vlx) => par défaut Idrisi 15 met l'en-tête de la première colonne

- Vector file => Pi2_vc_vg (réalisé à l'étape POLYVEC)

- Link field name => ZEA

En cas de problèmes, vérifiez que tous les fichiers sont bien dans le même répertoire !

- Affichez une carte vecteur pour voir si le lien fonctionne

- Sélectionnez une cellule

- Cliquez sur l'icône d'affichage (la carte)

 

c) Création d'une image résultante de la photo-interprétation avec Idrisi 15 (Andes) et ultérieur

- File / Export / Field / to Raster Image

- Output raster image => zea11_uea

- Baclground value => 0

- Reference Parameters

Copy from eisting file => spot6-5-07xscc321 (par exemple)

 

Admirez le résultat !

 

Retour au début de la fiche TD

 

Fiche Mémo suivante

 

NB : les mots suivis de "*" font partie du vocabulaire géographique, donc leur définition doit être connue. Faites-vous un glossaire.