Annonce
Pages: 1
- Sujet précédent - Définition d'un système de coordonnées pour un ensemble de rasters. - Sujet suivant
#1 Tue 09 December 2014 09:05
- Banghs
- Participant occasionnel
- Date d'inscription: 17 May 2011
- Messages: 14
Définition d'un système de coordonnées pour un ensemble de rasters.
Bonjour,
J'ai besoin de vos explications:
Disons que j'ai un ensemble de rasters(plus de 200 coupures) dont la référence spatiale est inconnue(Unknown) et, je dois les renseigner avant de procéder à leur mosaïquage. Seulement, manuellement c'est fastidieux. Donc, j'ai pensé à l'automatisation via Python et, je suis tombé sur un code dans un document pdf de M. Hervé BOHBOT du CNRS que j'ai adapté selon mon environnement de travail:
Code:
# importer la librairie arcpy
import arcpy
#définir l'espace de travail (repertoire où se trouvent les données
arcpy.env.workspace = "P:\03Planning-FSMT\07_GIS\07-TEMPORAIRE\06-Entrepôt nouvelle données\20140602_SEPRET\20140702_Acquisition Avenant complement\Orthophoto"
#récupérer un objet liste contenant tous les fichiers TIFF de l'espace de travail
lstRasters = arcpy.ListRasters("*.tif")
#parcourir la liste
for raster in lstRasters:
#
#si le système de référence spatial du raster est inconnu...
if arcpy.Describe(raster).spatialReference.name == "Unknown":
#
#définir sa projection à l'aide du fichier PRJ ad hoc
arcpy.DefineProjection_management(raster,r'P:\03Planning-FSMT\07_GIS\04-OUTILS\03_Gabon Toolboxes\BecGIS 25608 Gabon Transverse Mercator 2010 (GTM2010).prj')
#
#afficher un message à chaque définition de système de coordonnées
print "Systeme de coordonnees defini pour",rasterLe problème est que j'ai un message d'erreur à la fin du code.
Prière donc de m'aider s'il vous plait.
Merci;
Hors ligne
#2 Tue 09 December 2014 10:08
Re: Définition d'un système de coordonnées pour un ensemble de rasters.
Bonjour,
Si vous avez un message d'erreur, il serait utile de nous le faire partager.
Hors ligne
#3 Tue 09 December 2014 10:40
- Banghs
- Participant occasionnel
- Date d'inscription: 17 May 2011
- Messages: 14
Re: Définition d'un système de coordonnées pour un ensemble de rasters.
Ok,
Voici le message d'erreur:
Code:
Parsing error SyntaxError: invalid syntax (line 1) >>>
Merci
Hors ligne
#4 Tue 09 December 2014 19:55
Re: Définition d'un système de coordonnées pour un ensemble de rasters.
Bonjour,
En ajoutant un "r" dans la définition du chemin d'accès aux fichiers images, cela devrait mieux fonctionner :
Code:
arcpy.env.workspace = r"P:\03Planning-FSMT\07_GIS\07-TEMPORAIRE\06-Entrepôt nouvelle données\20140602_SEPRET\20140702_Acquisition Avenant complement\Orthophoto"
Le problème pourrait venir du fait que le backslash est un caractère d'échappement en Python et, du coup, le chemin n'est peut-être pas interprété comme il faut...
Une solution est d'ajouter un "r" devant une chaîne est le préfix qui permet d'interpréter différemment le backslash (cf. doc Python).
C'était bien fait sur la ligne du "arcpy.DefineProjection_management"
A+
Franck
Hors ligne
#5 Wed 10 December 2014 08:08
- Banghs
- Participant occasionnel
- Date d'inscription: 17 May 2011
- Messages: 14
Re: Définition d'un système de coordonnées pour un ensemble de rasters.
Bonjour Frank,
Merci pour votre solution que je m' envais tester tout de suite et vous donner le résultat.
Bien cordialement.
Hors ligne
#6 Wed 10 December 2014 09:05
- Banghs
- Participant occasionnel
- Date d'inscription: 17 May 2011
- Messages: 14
Re: Définition d'un système de coordonnées pour un ensemble de rasters.
Merci Franck,
J'ai apporté la modification à l'endroit indiqué et, il semble tourner mais les resters référencés s'affichent en noir.
Je vais peut- être attendre jusqu'à la fin pour voir ce que ça donne.
A bientôt,
Hors ligne
Pages: 1
- Sujet précédent - Définition d'un système de coordonnées pour un ensemble de rasters. - Sujet suivant
