V.1.1.2 - Dernière mise à jour : 06/03/2015
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 :
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
- 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.
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).
- 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.
- 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 !
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 :
à télécharger ici
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)
- 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
- 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 !
NB : les mots suivis de "*" font partie du vocabulaire géographique, donc leur définition doit être connue. Faites-vous un glossaire.