Quelques Travaux pratique en SIG BEN ARFA Seifeddine
N.B. Toute Reproduction est Interdite !
Quelques Travaux pratique en SIG BEN ARFA Seifeddine
Mosaïquage de rasters &
fusion de shapefiles
I- Mosaïquage de rasters
Dans le cas où la zone d’étude couvre plusieurs cartes, il faut effectuer un mosaïquage de
ces cartes. Nous allons dans ce TP détailler les étapes à suivre pour effectuer ce mosaïquage
sous RasterStitch ainsi que sous ArcGIS.
1- Mosaïquage de rasters sous RasterStitch
Nous allons effectuer un mosaïquage sous RasterStitch des deux cartes topographiques
adjacentes Bizerte_NO et Bizerte_NE, cartes d’échelles 1:25000. Nous allons joindre les
points de même coordonnées sur la ligne commune entre les deux cartes topographiques.
Commencez par ouvrir RasterStitch Allez cliquez sur dans la fenêtre de
gauche Allez ouvrir à partir du répertoire de travail Allez cliquez sur
dans la fenêtre de droite Allez ouvrir .
Nous devons sélectionner un premier point sur la limite gauche de et
sélectionner le point de même coordonnées sur la limite droite de : Cliquez sur
la fenêtre d’affichage relative à pour la sélectionner Zoomez sur
la partie gauche en utilisant (atteindre un zoom de 100%) et allez par la suite vers le point
fixé Cliquez sur dans le menu et double-cliquez par la suite sur ce point Ouverture de
la fenêtre Double-cliquez sur le point zoomé Cliquez sur OK Le point 1
s’affiche sur la carte Sélectionnez par la suite la fenêtre d’affichage
relative à Zoomez sur la partie droite en utilisant (atteindre
un zoom de 100%) et allez par la suite vers le point de même coordonnées Cliquez sur
dans le menu et double-cliquez par la suite sur ce point Ouverture de la fenêtre
Double-cliquez sur le point zoomé Cliquez sur OK Le point 1 s’affiche sur la carte
.
Quelques Travaux pratique en SIG BEN ARFA Seifeddine
Prenons un deuxième point sur la limite gauche de et un point de mêmes
coordonnées sur la limite droite de : Sélectionnez la fenêtre d’affichage
relative à Cliquez sur dans le menu et double-cliquez par la suite sur le point
voulu Ouverture de la fenêtre Double-cliquez sur le point
zoomé Cliquez sur OK Le point 2 s’affiche sur la carte Sélectionnez
par la suite la fenêtre d’affichage relative à Allez vers le point de
mêmes coordonnées sur la limite gauche Cliquez sur dans le menu et double-cliquez sur ce
point Ouverture de la fenêtre Double-cliquez sur le point zoomé Cliquez
sur OK Le point 2 s’affiche sur la carte .
Reliez les deux cartes en cliquant sur :
Sauvegardez le fichier résultant issu du mosaïquage des deux cartes sous format jpg et
sous le nom de Bizerte_Nord : Allez à File choisir comme type de fichier
et mentionnez le nouveau nom. Ce nouveau fichier pourra être géoréférencé et pourra servir pour
la digitalisation.
.
Quelques Travaux pratique en SIG BEN ARFA Seifeddine
Géoréférencement d’une carte scannée,
Création de shapefiles, et Digitalisation
I- Géoréférencement d’une carte scannée (ArcMap et ArcCatalog )
Pour géoréférencer (= caler) une carte ou un plan scanné sous ArcGIS commencez
parcréer un répertoire de travail (nommé TP1-SIG par exemple) dans l’un des compartiments
de votre PC (D: par exemple). Ce répertoire doit contenir la carte ou le plan scanné à
géoréférencer (Nabeul_SE.jpg).
Lancez ArcMap : démarrer Programmes ArcGIS ArcMap. A u
démarrage, l’application ArcMap offre 3 possibilités :
- Créer une nouvelle carte (A new empty map).
- Un modèle de carte (A template) : ArcGIS propose de nombreux modèles de cartes.
L’intérêt de ces modèles est de réaliser rapidement une carte ou de permettre à l’utilisateur
d’utiliser son propre modèle de carte préalablement établi afin que toutes les cartes produites
possèdent des apparences identiques.
- Une carte existante (An existing map).
1- Choisissez dans notre exemple l’option : A new empty map
Sous ArcMap on distingue une structure hiérarchique d’organisation de
plans d’informations cartographiques qui est la suivante : Map, Frame, et Layers. Le terme
Map (Carte) est à considérer comme étant une superposition de Layers (Couches
d’informations). Par exemple, une carte d’occupation du sol peut être issue d’une
superposition de plusieurs couches : une couche ne contenant que des forêts, une couche ne
contenant que des surfaces d’eau libre, une couche ne contenant que des surfaces bâties, une
couche contenant le réseau routier et les chemins et une couche contenant les surfaces
agricoles. On peut admettre que la dernière couche est plus variable d’une année à une
autre. Les autres couches peuvent être considérées comme des thèmes relativement stables
dans le temps. Ces différentes couches sont appelées Layers sous ArcGIS. Un ensemble
de couches peuvent appartenir au même Frame (Cadre). Un ensemble de couches peut
être enregistré dans un fichier Map ayant l’extension mxd. Une couche peut être
enregistrée indépendamment des autres couches dans un fichier Layer ayant l’extension lyr.
Quelques Travaux pratique en SIG BEN ARFA Seifeddine
2- Ajoutez la carte scannée comme arrière-plan : Menu principal : File Add Data
3- Naviguez jusqu’à l’emplacement de votre image Nabeul_SE.jpg (dossier TP1-SIG),
sélectionnez-là à l’aide d’un seul clic puis cliquez sur le bouton Add. Une fenêtre Create
pyramids for Nabeul_SE.jpg s’affiche, cliquer sur Yes et cliquer par la suite sur OK dans la
fenêtre Unknown Spatial Reference qui s’affiche.
4- Vérifiez que l’outil Georeferencing est activé, sinon vous pouvez l’activer à partir du
Menu View Toolbar. La barre d’outils relative à Georeferencing est la suivante :
5- Zoomez sur la zone d’intérêt en utilisant et choisissez un premier point de contrôle
(appelé également point de calage ou point d’amer) sur le carroyage de la carte
Cliquez une seule fois sur le point de contrôle à l’aide du bouton d’ajout de points de
contrôle Cliquez enfin sur le bouton droit de la souris.
6- Introduisez les coordonnées réels en mètres de ce point de contrôle après avoir cliqué sur
Input X and Y (pour annuler le choix du point cliquez sur Cancel Point).
7- Cliquez sur pour visualiser de nouveau la carte scannée.
8- Refaites le même travail pour les 3 autres points de calage.
Notons que :
- ArcMap donne des résultats approximatifs aux coordonnées réelles à partir du 2ème
point.
- On peut introduire, sans exagérer, un nombre de points de calage supérieur à 4, mais ces
points doivent être bien répartis.
- Il faut vérifier l’erreur quadratique moyenne RMS en cliquant sur le bouton (View Link
Table). L’erreur (RMS) de chacun des points de contrôle ne doit pas dépasser 2.5m pour
toute carte d’échelle 1/25.000, et ne doit pas dépasser 5m pour une carte d’échelle 1/50.000,
etc. Les points ayant un RMS (Residual) dépassant la valeur acceptable peuvent être
supprimés en les sélectionnant dans la fenêtre Link Table et en cliquant par la suite sur
l’icône de la même fenêtre.
Quelques Travaux pratique en SIG BEN ARFA Seifeddine
9- Cliquez de nouveau sur le bouton (View Link Table) de l’outil Georeferencing pour
afficher la "Link Table". Sauvegardez le fichier de géoréférencement sous format texte en
cliquant sur Save et lui attribuez un nom (par exemple le nom Georef_Nabeul). Cliquez enfin
sur OK et enregistrer sous le répertoire TP1-SIG.
10- Allez à Georeferencing et cliquez sur Rectify Une fenêtre Save As s’ouvre, choisir
comme Resample Type : Nearest Neighbor (for discrete data), comme Format : TIFF et
enregistrez sous le nom NABEUL_SE.tif et cliquez enfin sur Save.
11- Allez à l’image NABEUL_SE.jpg affichée sous , cliquez sur par le bouton
droit de la souris et l’enlever en cliquant sur Remove. Cliquez ensuite sur le bouton Add
Data pour ajouter l’image NABEUL_SE.tif. Les unités et le système de coordonnées de
cette image sont inconnus.
Quelques Travaux pratique en SIG BEN ARFA Seifeddine
12- Allez maintenant sous ArcCatalog pour attribuer le système de coordonnées et les
unités : Lancez ArcCatalog en cliquant sous Allez dans la table d’affichage au dossier de
travail et cliquez par le bouton droit sur NABEUL_SE.tif et cliquez sur Properties
Ouverture de la fenêtre Raster Dataset Properties dans laquelle allez à Spatial Reference et
cliquez sur Edit Ouverture de la fenêtre Spatial Reference Properties dans laquelle
cliquez sur Select Ouverture de la fenêtre Browse for Coordinate System dans laquelle
cliquez sur Cliquez sur
Cliquez sur Cliquez sur dans Raster Dataset Properties.
13- Revenir à ArcMap : Affichez NABEUL_SE.tif et vérifier le système de coordonnées
et les unités : Cliquez par le bouton droit de la souris sur
Ouverture de la fenêtre dans laquelle allez à et vérifiez la
référence spatiale :
14- On peut modifier la couleur du fond de cette image du noir vers le blanc : Cliquez par le
bouton droit sur le nom de l’image NABEUL_SE.tif affichée sous , allez ensuite
cliquer sur ouverture de la fenêtre dans laquelle allez à
et cochez la fenêtre .
Quelques Travaux pratique en SIG BEN ARFA Seifeddine
15- Allez à View Toolbars et activez Effects. Vous pouvez ainsi modifier le contrase
(Contrast), la luminosité (Brightness), la transparence (Transparency) et utiliser les
fonctions Swipe Layer et Flicker Layer.
16- Terminez par enregistrer votre travail en tant que ArcMap Document sous le nom de
Nabeul_SE par exemple : Allez dans le menu principal à File Save As.
II- Création de shapefiles (ArcCatalog )
Lancez ArcCatalog : Allez à démarrer Tous les programmes
ArcGIS
ArcCatalog.
1- Créez un nouveau dossier dans lequel vous allez enregistrer vos fichiers nommé Shapefiles
par exemple (Actualisez ArcCatalog en appuyant sur la touche F5 de votre clavier) et
cliquez sur ce dossier par le bouton droit de la souris Allez à New et ensuite à Shapefile
La fenêtre Create New Shapefile apparaît.
2- Nous allons commencer dans notre cas par créer une couche d’informations pour héberger
les points cotés : Attribuez un nom (Name) à cette couche (PC_Nabeul_SE par exemple), et
sélectionnez un Feature Type (type d’entité) Point.
3- Le système de coordonnées qu’on va utiliser est le Lambert Nord Tunisie : Pour attribuer
ce système au shapefile créé cliquez sur le bouton Edit de la fenêtre Create New Shapefile
La fenêtre Spatial Reference Properties s’ouvre Cliquez sur Select Ouverture de
la fenêtre Browse for Coordinate System dans laquelle cliquez sur
Cliquez sur puis sur Cliquez
sur dans Create New Shapefile.
4- Une fois le fichier shapefile Points cotés est créé, on peut ajouter à sa table attributaire
des champs : exemple : le champ Elévat relatif à l’élévation : Sélectionnez le fichier Points
cotés et cliquez sur par le bouton droit de la souris Cliquez sur Properties La
fenêtre Shapefile Properties s’ouvre Sélectionnez Fields : L’ajout des champs selon
vos besoins tout comme on le fait sous Microsoft Access : Ajoutez un champ Elévat de
Quelques Travaux pratique en SIG BEN ARFA Seifeddine
type Short Integer : Commencez par le Data Type et sélectionnez Short Integer et allez par
la suite à la colonne Field Name pour mentionner le nom (Elévat) Cliquez sur puis sur
.
NB : Le fait d’ajouter des attributs (Fields) peut se faire également sous ArcMap . Ceci
sera évoqué dans les TPs suivants.
5- Créez dans le même dossier de travail (dossier Shapefiles) une 2ème couche qui sera
relative aux courbes de niveau : Attribuez un nom (Name) à cette couche (CN_Nabeul_SE
par exemple) et sélectionnez un Feature Type (type d’entité) Polyline. Pour attribuer le
système de coordonnées Lambert Nord Tunisie au shapefile créé cliquez sur le bouton
Edit de la fenêtre Create New Shapefile La fenêtre Spatial Reference Properties
s’ouvre Cliquez sur Import Ouverture de la fenêtre Browse for Coordinate S stem
y
dans laquelle sélectionnez dans le dossier Shapefiles le shapfile PC_Nabeul_SE préparé
précédemment et cliquez sur Cliquez enfin sur puis sur dans la fenêtre
Spatial Reference Properties.
6- Ajouter à la table attributaire du shapefile CN_Nabeul_SE un champ Elévat relatif à
l’élévation.
Quelques Travaux pratique en SIG BEN ARFA Seifeddine
III- Digitalisation (ArcMap )
Après avoir géoréférencé votre carte topographique, créé la couche d’information où vous
allez stocker vos données graphiques et attributaires (couches points cotés PC_Nabeul_SE et
couche courbes de niveau CN_Nabeul_SE), vous pouvez commencez la digitalisation sous
ArcMap .
III- 1 - D i g i t a l i sa t i o n d e l a c o u c h e p o i n t s c o t é s
1- Commencez par ouvrir ArcMap et faites appel à votre fichier
shapefile PC_Nabeul_SE : Cliquer sur Layers par le bouton droit de la souris et ensuite sur
ou directement sur .
2- Nous allons commencer par digitaliser (= numériser ou vectoriser) un premier point coté de
la carte topographique : Cliquez sur le bouton Editor Start Editing :
Les options Create New Feature et Target sont activées. Target, qui signifie cible, est la
couche PC_Nabeul_SE, couche qui contiendra les points à créer. Le bouton "Sketch Tool"
devient également actif, c’est l’outil de numérisation Cliquez sur ce bouton et
cliquez une fois sur l’emplacement du point coté Enregistrez le point numérisé en
cliquant sur sous Pour stopper la digitalisation cliquez
sous .
NB : Pour modifier la symbologie double-cliquez sur le symbole relatif à la couche
Ouverture de la fenêtre dans laquelle vous pouvez modifier
la couleur et la taille de votre symbole.
Quelques Travaux pratique en SIG BEN ARFA Seifeddine
3- Pour remplir le champ Elévat : Commencez par activez l’édition
Cliquez sur par le bouton droit de la souris
Cliquez surpour ouvrir la table attributaire relative à cette couche
Remplir le
champ relatif Elévat Finir par arrêter l’édition :
4- Faites afficher les "labels" sur la carte à partir du champ Elévat : Allez sous à
et cliquez-sur par le bouton droit de la souris Ouverture
de
la fenêtre Layer Properties Allez à Cochez Label features in this
layer : Choisir Elévat comme Label Field En cas de
besoin vous pouvez changer la police, la taille, la couleur… Cliquez sur OK.
Poursuivez la digitalisation des différents points côtés de la carte topo dans le shapefil
PC_Nabeul_SE.
III- 2 - D i g i t a l i sa t i o n d e l a c o u c h e c o u r b e s d e n i v e a u
1- Nous allons commencer par digitaliser une première polyligne de la carte topographique :
Cliquez sur le bouton Editor Start Editing :
Les options Create New Feature et Target sont activées. Target est la couche
CN_Nabeul_SE, couche qui contiendra les courbes de niveau.
2- Cliquez sur le bouton outil de numérisation "Sketch Tool" pour numériser les courbes
de niveau :
Créez des vertex en cliquant sur le bouton gauche de la souris, cliquez deux fois sur le
même bouton pour terminer.
Pour supprimer un vertex : le sélectionner et cliquer sur par le bouton droit de souris
Cliquer sur .
Pour changer la position d’un vertex : sélectionnez la polyligne pour faire apparaître les
vertex qui la forment Cliquez sur le vertex par le bouton ( )
et modifiez sa position.
Si vous voulez zoomer sur des polylignes sélectionnées cliquez par le bouton droit à
l’intérieur du cadre d’affiche et cliquez par la suite sur .
Quelques Travaux pratique en SIG BEN ARFA Seifeddine
Utilisez l’outil pour zoomer sur la couche voulue.
Si on veut joindre une nouvelle polyligne à l’extrémité d’une polyligne déjà établie :
Vérifiez que l’édition est active, sinon aller à Editor Start Editing. Allez par la suite
au menu Editor Snapping et cochez les cases relatives à Vertex, Edge, End
de la couche CN_Nabeul_SE :
Réactivez par la suite l’outil "Sketch Tool" et numérisez la deuxième polyligne.
Si on veut fusionner les deux polylignes adjacentes on doit utiliser la fonction Merge :
Sélectionnez les deux polylignes puis Ouverture de la fenêtre
Cliquez sur .
Si on veut subdiviser une polyligne : L'édition étant activée (Start Editing), sélectionner la
polyligne à subdiviser Allez par la suite au bouton Split Tool situé sur la barre
d’outilsEditor pour subdiviser la polyligne à l’endroit voulu.
Si on veut transformer des polylignes sélectionnées en un nouveau shapefile : Allez sous
Layers et cliquez par le bouton droit de la souris sur le nom de la couche Allez à Data
Export Data Attribuez un nom au nouveau shapefile.
3- Pour remplir le champ Elévat : Commencez par activez l’édition
Cliquez sur par le bouton droit de la souris
Cliquez surpour ouvrir la table attributaire relative à cette couche
Remplir le
champ relatif Elévat Finir par arrêter l’édition :
.
Poursuivez la digitalisation des différentes courbes de niveau de la carte topographique
dans le shapefile CN_Nabeul_SE.
