no way to compare when less than two revisions
Differences
This shows you the differences between two versions of the page.
— | fr:int [2006/12/12 17:41] (current) – created - external edit 127.0.0.1 | ||
---|---|---|---|
Line 1: | Line 1: | ||
+ | Une fois que le projet est structuré, il faut penser à lui donner des données à afficher !!\\ | ||
+ | Pour celà, nous allons importer les cartes que nous avons sous forme de scan (format jpg ou de préférence **TIFF**). | ||
+ | ===== Création de la LOCATION temporaire ===== | ||
+ | Pour bien importer en les géorectifiant les images, il faut commencer par les importer dans une location spécifique (qui en fait va servir à mettre les images dans le format propre à **GRASS**). Cette création est créée de la façon décrite précédemment, | ||
+ | Une bonne résolution de scan est de **200dpi**, pour une carte au 1/25 000e. Une résolution supérieure ne va pas forcément être très utile et en particulier, | ||
+ | |||
+ | |||
+ | ===== Import de l' | ||
+ | Ouvrez donc **GRASS** avec la location XY (en sélectionnant dans la fenêtre de démarrage). Dans l' | ||
+ | La fenêtre qui s' | ||
+ | Vous pouvez aussi faire ça dans la fenêtre de Terminal, en tapant la commande suivante : | ||
+ | < | ||
+ | r.in.gdal input=image.tif output=imageGRASS | ||
+ | </ | ||
+ | où image.tif doit être remplacé par le chemin et le nom du fichier à traiter et imageGRASS par le nom que vous attribuez à l' | ||
+ | |||
+ | Cette manipulation prend un petit moment, variable selon la taille de votre carte et sa résolution.\\ | ||
+ | A l' | ||
+ | Il faut donc maintenant recomposer l' | ||
+ | |||
+ | ===== Recomposition de l' | ||
+ | Pour ce faire, il existe la commande r.composite... | ||
+ | Dans l' | ||
+ | < | ||
+ | Raster -> Manage map colors -> Create color image from RGB files | ||
+ | </ | ||
+ | Dans la fenêtre qui s' | ||
+ | En ligne de commande, dans le Terminal, vous saisissez : | ||
+ | < | ||
+ | r.composite red=xxx.red green=xxx.green blue=xxx.blue levels=32 output=xxx | ||
+ | </ | ||
+ | Le calcul prend un certain temps... Si vous avez plusieurs images à importer, n' | ||
+ | < | ||
+ | # | ||
+ | #Script d' | ||
+ | #dans GRASS | ||
+ | |||
+ | if test " | ||
+ | echo "Tu dois etre dans GRASS pour lancer ce programme" | ||
+ | exit 1 | ||
+ | fi | ||
+ | |||
+ | # import des cartes en XY | ||
+ | |||
+ | r.in.gdal input=/ | ||
+ | r.in.gdal input=/ | ||
+ | |||
+ | # composition couleur | ||
+ | |||
+ | r.composite red=MQ06.red green=MQ06.green blue=MQ06.blue output=MQ06R | ||
+ | r.composite red=MQ06Parc.red green=MQ06Parc.green blue=MQ06Parc.blue output=MQ06ParcR | ||
+ | |||
+ | # suppression des tempo | ||
+ | |||
+ | g.remove rast=MQ06.red, | ||
+ | g.remove group=MQ06, | ||
+ | </ | ||
+ | Ainsi, en même temps, vous importez vos images (2 ici), puis vous les recomposez en couleurs, vous supprimez les fichiers et les groupes temporaires (on verra les groupes un peu plus loin... car ils sont utiles mais ceux-ci, vous pouvez les supprimer). Pour plus de confort, j' | ||
+ | |||
+ | |||
+ | ===== Projection dans le véritable projet ===== | ||
+ | Maintenant que notre carte est tout à fait visible dans **GRASS**, on va la reprojeter dans la LOCATION que nous avons créée en premier et qui est la zone de votre étude... | ||
+ | ==== 1 : création des groupes de traitement ==== | ||
+ | Les groupes sont utilisées par les fonctions du module " | ||
+ | < | ||
+ | Image -> Develop images and group -> Create/edit imagery group | ||
+ | </ | ||
+ | Là, vous choisissez le nom du groupe et le nom du raster à inclure dans ce groupe. Pour le nom, vous êtes libre... mais restez clair tout de même !!\\ | ||
+ | En ligne de commande, ça donne : | ||
+ | < | ||
+ | i.group group=XXX input=xxx | ||
+ | </ | ||
+ | où __group__ est le nom du groupe et __input__ le nom du raster.\\ | ||
+ | Après ça, vous devez cibler le groupe de traitement.. En mode graphique : | ||
+ | < | ||
+ | Image -> Develop images and group -> Target imagery group | ||
+ | </ | ||
+ | On précise le nom du groupe à cibler, le nom de la LOCATION et du MAPSET de destination. | ||
+ | En ligne de commande : | ||
+ | < | ||
+ | i.target group=XXX location=DEST mapset=mapset | ||
+ | </ | ||
+ | où __group__ est le nom du groupe, __location__ le nom de la location de destination et __mapset__ le nom du jeu de données d' | ||
+ | |||
+ | Ca y est, vous allez pouvoir préciser les points de contrôle. | ||
+ | |||
+ | ==== 2 : Points de contrôle ==== | ||
+ | Ces points sont ceux qui servent à calculer l' | ||
+ | Pour celà, on va dans : | ||
+ | < | ||
+ | Image -> Rectify and georeference image group -> Set ground control points (GCP' | ||
+ | </ | ||
+ | A ce moment-là, 2 fenêtres s' | ||
+ | Les éléments affichés sont, à mon avis, suffisamment clairs... Si cependant vous rencontrez des difficultés, | ||
+ | Dans le cas standard (une carte scannée et la plupart des cas en fait), il faudra traiter les données avec un polynome d' | ||
+ | |||
+ | ==== 3 : La géorectification ==== | ||
+ | Le grand moment est arrivé... Il ne reste plus qu'à lancer la géorectification.\\ | ||
+ | Dans l' | ||
+ | < | ||
+ | Image -> Rectify and georeference image group -> Affine and polynomial... | ||
+ | </ | ||
+ | Vous précisez le groupe à utilisez, la carte contenue dans le groupe (c'est une option, mais ça ne coûte rien de la préciser !), l' | ||
+ | En ligne de commande, ça nous donne : | ||
+ | < | ||
+ | i.rectify group=XXX input=xxx extension=e order=1 | ||
+ | </ | ||
+ | où __group__ est le nom du groupe, __input__ le nom du raster, __extension__ la lettre qui sera rajoutée (ici, e) et __order__ l' | ||
+ | |||
+ | Cette manipulation prend "un certain temps" | ||
+ | |||
+ | Bien évidemment, | ||
+ | |||
+ | Quand votre ou vos rectifications sont terminées, vous pouvez quitter la LOCATION d' | ||
+ | < | ||
+ | File -> Exit | ||
+ | </ | ||
+ | puis, dans la fenêtre de Terminal, tapez : | ||
+ | < | ||
+ | exit | ||
+ | </ | ||
+ | C'est bon, vous pouvez maintenant retourner dans votre LOCATION géoréférencée... |