{"id":1877,"date":"2015-06-09T08:30:06","date_gmt":"2015-06-09T06:30:06","guid":{"rendered":"http:\/\/www.sigterritoires.fr\/?p=1877"},"modified":"2015-06-09T08:30:06","modified_gmt":"2015-06-09T06:30:06","slug":"arcgis-et-python-3-votre-premier-script","status":"publish","type":"post","link":"https:\/\/www.sigterritoires.fr\/index.php\/arcgis-et-python-3-votre-premier-script\/","title":{"rendered":"ArcGis et Python (3) : votre premier script"},"content":{"rendered":"

Si vous avez lu les deux articles pr\u00e9c\u00e9dents\u00a0 ArcGis et Python :Avant de faire vos premiers pas<\/a> et\u00a0 ArcGis et Python : pour bien d\u00e9marrer<\/a>, vous \u00eates pr\u00eat pour vous lancer.<\/p>\n

Il y a deux choses que vous devez apprendre : le langage Python et comment utiliser la biblioth\u00e8que de g\u00e9otraitement Python d’ArcGis. En ce qui concerne le langage Python, des tutoriaux ne manquent pas. Mais n’oubliez pas que vous en aurez besoin que des bases du langage. Vous allez essentiellement encha\u00eener des outils ArcGis d\u00e9j\u00e0 tous faits.<\/p>\n

Nous allons voir ici un premier exemple de script et son int\u00e9gration dans votre toolbox personnelle.Ne faisons pas dans l’originalit\u00e9. On va encore utiliser l’outil Clip<\/strong> <\/em>(d\u00e9coupe) d’une couche. Disons que nous voulons adapter l’outil Clip<\/strong> <\/em>pour que l’outil d\u00e9coupe automatiquement toutes les couches situ\u00e9es das un espace de travail (r\u00e9pertoire, geodatabase), en stockant les r\u00e9sultats dans un autre espace de travail.<\/p>\n

On pourrait faire un mod\u00e8le avec Model Builder en utilisant un it\u00e9rateur en entr\u00e9e de donn\u00e9es. Et d’ailleurs c’est surement ce que vous feriez actuellement si vous ne faites que d\u00e9marrer avec Python. Voyons donc comment faire la m\u00eame chose en Python.<\/p>\n

Nous allons utiliser pour l’exemple PyScripter,<\/em> mais pour cet exemple, vous n’avez m\u00eame pas besoin d’IDE. Un \u00e9diteur de texte peut tr\u00e8s bien faire l’affaire. Par contre, sachez que dans Python, l’indentation des lignes fait partie de la syntaxe et que taper une ligne de code en d\u00e9but de ligne ou apr\u00e8s une ou deux tabulations ne veut pas dire la m\u00eame chose.<\/p>\n

Un IDE vous positionnera automatiquement l\u00e0 o\u00f9 il faut.<\/p>\n

Lancez votre IDE ou ouvrez votre \u00e9diteur de texte.<\/p>\n

$\"ouverture$ <\/a><\/p>\n

Vous allez taper votre ligne incontournable<\/p>\n

import arcpy<\/strong><\/em><\/p>\n

sans laquelle Python ne conna\u00eet pas du tout les traitements ArcGis.<\/p>\n

$\"d\u00e9tail$ <\/a><\/p>\n

Vous allez en rajouter une deuxi\u00e8me, qui appelle le module os<\/span><\/strong><\/em> qui offre un acc\u00e8s rapide aux outils fondamentaux du syst\u00e8me d’exploitation. Certaines m\u00e9thodes de gestion des noms de fichiers du module os<\/span> <\/strong><\/em>seront utilis\u00e9es dans ce script.<\/p>\n

import os<\/strong><\/em><\/p>\n

$\"d\u00e9tail$ <\/a><\/p>\n

Le script aura deux parties:<\/p>\n

\n
la premi\u00e8re o\u00f9 vous allez d\u00e9finir les donn\u00e9es en entr\u00e9e et en sortie<\/li>\n
le deuxi\u00e8me o\u00f9 vous allez ex\u00e9cuter les traitements<\/li>\n<\/ul>\n
D\u00e9finition des donn\u00e9es en entr\u00e9e et en sortie<\/h2>\n
Cette partie a un double int\u00e9r\u00eat. Bien s\u00fbr pour l’ex\u00e9cution du traitement il faut que l’outil sache sur quelle donn\u00e9e travailler. Mais aussi parce que pour pouvoir ex\u00e9cuter le script en tant qu’outil de la toolbox il va falloir cr\u00e9er la fen\u00eatre de param\u00e9trage qui sera affich\u00e9 lors du lancement de l’outil. Il faudra alors reprendre, une par une, les d\u00e9finitions de notre script pour cr\u00e9er un champ d’entr\u00e9e dans la fen\u00eatre. Vous avez donc tout int\u00e9r\u00eat \u00e0 commenter vos lignes de code pour bien vous retrouver au bout du processus.<\/p>\n
Nous allons utiliser 4 variables en entr\u00e9e\/sortie<\/p>\n
\n
l’espace de travail dans lequel se trouvent les couches \u00e0 d\u00e9couper<\/li>\n
la couche qui servira \u00e0 d\u00e9couper les autres<\/li>\n
l’espace de travail o\u00f9 mettre les couches d\u00e9coup\u00e9es r\u00e9sultantes<\/li>\n
la tol\u00e9rance \u00e0 utiliser par l’outil D\u00e9coupe<\/li>\n<\/ul>\n
Entrez les lignes de code suivantes:<\/p>\n
# Espace de travail en entr\u00e9e : Type WORKSPACE<\/strong><\/em>
\n arcpy.env.workspace = arcpy.GetParameterAsText(0)<\/strong><\/em><\/p>\n
#Classe d’entit\u00e9s servant \u00e0 d\u00e9couper Type FEATURECLASS<\/strong><\/em>
\n clipFeatures = arcpy.GetParameterAsText(1)<\/strong><\/em><\/p>\n
# Espace de travail en sortie : Type WORKSPACE<\/strong><\/em>
\n outWorkspace = arcpy.GetParameterAsText(2)<\/strong><\/em><\/p>\n
# Tol\u00e9rance XY de l’outil D\u00e9coupe : Type DOUBLE<\/strong><\/em>
\n clusterTolerance = arcpy.GetParameterAsText(3)<\/strong><\/em><\/p>\n
$\"d\u00e9finition$ <\/a><\/p>\n
La fonction GetParameterAsText<\/strong> <\/em>r\u00e9cup\u00e8re les param\u00e8tres en entr\u00e9e du script. Si vous voulez l’ex\u00e9cuter en ligne de commande (dans une fen\u00eatre de commande) vous utiliserez un commande de type<\/p>\n
python.exe \u00ab\u00a0patram\u00e8tre1\u00a0\u00bb \u00ab\u00a0param\u00e8tre2\u00a0\u00bb \u00ab\u00a0param\u00e8tre3\u00a0\u00bb Nombre<\/strong><\/em><\/p>\n
Si, comme nous le verrons plus loin, vous utilisez un outil de la toolbox, vous devrez faire correspondre l’ordre des champs de la fen\u00eatre de dialogue avec l’ordre d\u00e9fini dans les variables du script: le premier champ de la fen\u00eatre de dialogue devra \u00eatre l’espace de travail en entr\u00e9e, le deuxi\u00e8me la couche de d\u00e9coupe, etc.<\/p>\n
Dans nos quatre d\u00e9finitions, la premi\u00e8re est une d\u00e9finition d’une variable ArcGis, c’est l’espace de travail courant, tandis que les trois autres sont des variables propres \u00e0 ce script.<\/p>\n
Le corps du script : le traitement<\/h2>\n
Avant de commencer \u00e0 rentrer les lignes de code du traitement en lui-m\u00eame, nous allons rentrer un instruction bizarre. Nous allons taper simplement<\/p>\n
try:<\/strong><\/em><\/p>\n
A la fin de notre traitement nous allons rentrer l’instruction compl\u00e9mentaire de celle-ci : except<\/strong><\/em><\/p>\n
Elles ne sont pas obligatoires. Si vous ne les entrez pas, le script fonctionnera sans probl\u00e8me, sauf… quand il y aura un probl\u00e8me. Si une erreur survient lors du traitement et que vous n’avez pas ces instructions, vous aurez une belle fen\u00eatre \u00ab\u00a0ArcGis a rencontr\u00e9 une erreur et doit \u00eatre ferm\u00e9…\u00a0\u00bb.
\nPar contre, si vous avez mi try\/except, si une erreur survient dans le corps du traitement, au lieu de planter le syst\u00e8me, le script ex\u00e9cutera les instructions qu’il trouvera dans le paragraphe except (en g\u00e9n\u00e9ral l’affichage d’un message d’erreur et d’une explication de celui-ci), tout en vous rendant la main de mani\u00e8re propre.<\/p>\n
La premi\u00e8re instruction que nous allons rentrer correspond \u00e0 la cr\u00e9ation de la liste des couches pr\u00e9sentes dans l’espace de travail en entr\u00e9e et qu’il faudra d\u00e9couper.<\/p>\n
fcs = arcpy.ListFeatureClasses()<\/strong><\/em><\/p>\n
$\"gestion$ <\/a><\/p>\n
Notez que les lignes ne commencent pas au d\u00e9but de la ligne. L’indentation indique \u00e0 Python que cette commande est \u00e0 l’int\u00e9rieur de la commande try<\/strong><\/em>.
\nLa variable fcs<\/strong> <\/em>contiendra la liste des couches \u00e0 d\u00e9couper. L’instruction ListFeatureClasses<\/strong><\/em> g\u00e9n\u00e8re la liste des classes d’entit\u00e9s de l’espace de travail courant. Rappelez vous, nous l’avons d\u00e9fini dans la premi\u00e8re instruction du script, le param\u00e8tre 0.<\/p>\n
Maintenant que nous avons la liste des couches \u00e0 traiter, nous allons mettre une instruction de type boucle pour que le traitement se r\u00e9p\u00e8te sur chacune des couches de la liste<\/p>\n
for fc in fcs:<\/strong><\/em><\/p>\n
Dans cette boucle nous allons d\u00e9couper une couche en entr\u00e9e et cr\u00e9er une couche en sortie. Avant de la d\u00e9couper, nous devons g\u00e9rer certains choses dont nous nous sommes jusqu’\u00e0 pr\u00e9sent d\u00e9charg\u00e9s sur ArcGis, soit sur les commandes soit sur les outils de la toolbox. Par exemple, si notre espace de travail en entr\u00e9e est une g\u00e9odatabase et que nous d\u00e9finissons notre espace de sortie comme un r\u00e9pertoire, les couches en sortie seront des fichiers shapefile<\/a>. Pas de probl\u00e8me, sauf que les noms des shapefile sont beaucoup plus contraignants que ceux de la g\u00e9odatabase. Nous risquons d’avoir des noms trop longs ou avec des points ou des tirets. Dans un script, c’est \u00e0 nous de g\u00e9rer ce type de probl\u00e8me. Rassurez vous. Ce n’est pas compliqu\u00e9, mais il faut le faire.<\/p>\n
Il suffit d\u2019appeler la fonction ValidateTableName<\/strong> <\/em>qui transforme les noms en entr\u00e9e dans des noms corrects en sortie, tout en s’assurant qu’il n’y a pas de noms r\u00e9p\u00e9t\u00e9es en sortie.<\/p>\n
featureClassName = arcpy.ValidateTableName(fc, outWorkspace)<\/strong><\/em>
\n outFeatureClass = os.path.join(outWorkspace, featureClassName)<\/strong><\/em><\/p>\n
$\"boucle$ <\/a><\/p>\n
Maintenant nous allons d\u00e9couper la couche en entr\u00e9e avec la couche de d\u00e9coupe. Mais si la couche de d\u00e9coupe se trouve dans l’espace de travail en entr\u00e9e, nous ne voulons pas la d\u00e9couper! Nous mettons un test pour faire la d\u00e9coupe \u00e0 condition que le nom de la couche \u00e0 d\u00e9couper soit diff\u00e9rent de celui qui d\u00e9coupe:<\/p>\n
if fc != os.path.basename(clipFeatures):<\/strong><\/em><\/p>\n
Cette fonction r\u00e9cup\u00e8re juste le nom du fichier, sans le chemin, car dans notre liste en entr\u00e9e nous n’avons que les noms sans chemin.
\nPuis la fonction de d\u00e9coupe:<\/p>\n
arcpy.Clip_analysis(fc, clipFeatures, outFeatureClass, clusterTolerance)<\/strong><\/em><\/p>\n
$\"corps$ <\/a><\/p>\n
Voil\u00e0! C’est fini pour le traitement.<\/p>\n
Nous devons maintenant renseigner quoi faire en cas d’erreur:<\/p>\n
except Exception as err:<\/strong><\/em>
\n arcpy.AddError(err)<\/strong><\/em>
\n print err<\/strong><\/em><\/p>\n
Comme le script peut \u00eatre ex\u00e9cut\u00e9 dans un outil toolbox ou tout seul, nous avons deux affichages de l’erreur. AddError<\/strong> <\/em>va ajouter le message d’erreur dans la fen\u00eatre de r\u00e9sultat de l’outil toolbos. Print<\/strong> <\/em>va afficher l’erreur dans la fen\u00eatre de commande si le script est ex\u00e9cut\u00e9 tout seul.<\/p>\n
$\"gestion$ <\/a><\/p>\n
Notez l’indentation de except<\/strong><\/em>. Il se place au m\u00eame niveau que try<\/strong><\/em>.<\/p>\n
Enregistrer votre fichier. Il est pr\u00eat pour passer dans la Toolbox.<\/p>\n
Ajout du script dans la toolbox<\/h2>\n
Dans la Toolbox, chargez une bo\u00eete \u00e0 outils personnelle (ici Mes Scripts<\/strong><\/em>).
\nCliquez droit sur la bo\u00eete \u00e0 outils et s\u00e9lectionnez Ajouter->Script<\/strong><\/em><\/p>\n
$\"ajouter$ <\/a>La premi\u00e8re fen\u00eatre permet de d\u00e9finir le Nom de l’outil (sans espaces ni caract\u00e8res sp\u00e9ciaux) et l’\u00e9tiquette qui appara\u00eetra dans l’affichage de la bo\u00eete \u00e0 outils (avec espaces et caract\u00e8res sp\u00e9ciaux si vous voulez).<\/p>\n
$\"nom$ <\/a>La deuxi\u00e8me fen\u00eatre permet de faire le lien avec notre script.<\/p>\n
$\"lien$ <\/a>La derni\u00e8re fen\u00eatre c’est celle qui demande le plus de travail. Vous allez d\u00e9finir lrs c=hamps d’entr\u00e9e de param\u00e8tres de la fen\u00eatre de l’outil.<\/p>\n
Pour chaque param\u00e8tre du script (nous en avons quatre) il faudra une ligne. Dans la colonne Nom complet <\/strong><\/em>vous tapez le texte que vous voulez voir \u00e0 c\u00f4t\u00e9 du champ d’entr\u00e9e dans la fen\u00eatre de dialogue.<\/p>\n
Dans la colonne Type de donn\u00e9es<\/strong><\/em> vous devez indiquer le type de donn\u00e9e \u00e0 entrer dans le champ de la fen\u00eatre de dialogue.<\/p>\n
$\"d\u00e9finition$ <\/a>Pour le premier param\u00e8tre du script (param\u00e8tre 0) nous indiquons \u00ab\u00a0Espace de travail en entr\u00e9e<\/strong><\/em>\u00a0\u00bb et le type de donn\u00e9es \u00ab\u00a0Espace de travail<\/strong><\/em>\u00ab\u00a0.<\/p>\n
Voici la capture d’\u00e9cran avec l’ensemble des d\u00e9finitions pour notre script.<\/p>\n
$\"entr\u00e9es$ <\/a>L’ordre des lignes est tr\u00e8s importante. La premi\u00e8re ligne sera le premier param\u00e8tre d\u00e9fini dans le script, param\u00e8tre 0. La deuxi\u00e8me ligne correspondra ai=u deuxi\u00e8me param\u00e8tre, param\u00e8tre 1. Et ainsi de suite.<\/p>\n
Une fois termin\u00e9, fermez la fen\u00eatre d’ajout de script. Vous pouvez maintenant l’ex\u00e9cuter comme tous les autres outils, en double cliquant dessus dans la Toolbox.<\/p>\n
$\"fen\u00eatre$ <\/a>Vous retrouvez les quatre param\u00e8tres \u00e0 rentrer.<\/p>\n
Non seulement vous avez un nouvel outil qui peut s’ex\u00e9cuter \u00e0 partir de la toolbox mais vous pourrez l’inclure dans un mod\u00e8le de Model Builder.<\/p>\n","protected":false},"excerpt":{"rendered":"
Si vous avez lu les deux articles pr\u00e9c\u00e9dents\u00a0 ArcGis et Python :Avant de faire vos premiers pas et\u00a0 ArcGis et Python : pour bien d\u00e9marrer, vous \u00eates pr\u00eat pour vous lancer. Il y a deux choses…<\/p>\n","protected":false},"author":1,"featured_media":0,"comment_status":"open","ping_status":"closed","sticky":false,"template":"","format":"standard","meta":{"give_campaign_id":0,"_bbp_topic_count":0,"_bbp_reply_count":0,"_bbp_total_topic_count":0,"_bbp_total_reply_count":0,"_bbp_voice_count":0,"_bbp_anonymous_reply_count":0,"_bbp_topic_count_hidden":0,"_bbp_reply_count_hidden":0,"_bbp_forum_subforum_count":0,"sfsi_plus_gutenberg_text_before_share":"","sfsi_plus_gutenberg_show_text_before_share":"","sfsi_plus_gutenberg_icon_type":"","sfsi_plus_gutenberg_icon_alignemt":"","sfsi_plus_gutenburg_max_per_row":"","_monsterinsights_skip_tracking":false,"_monsterinsights_sitenote_active":false,"_monsterinsights_sitenote_note":"","_monsterinsights_sitenote_category":0,"jetpack_post_was_ever_published":false,"_jetpack_newsletter_access":"","_jetpack_dont_email_post_to_subs":false,"_jetpack_newsletter_tier_id":0,"_jetpack_memberships_contains_paywalled_content":false,"_jetpack_memberships_contains_paid_content":false,"footnotes":""},"categories":[145,222],"tags":[15,131,223,116],"class_list":["post-1877","post","type-post","status-publish","format-standard","hentry","category-arcgis","category-python","tag-arcgis","tag-python","tag-script","tag-toolbox"],"aioseo_notices":[],"jetpack_featured_media_url":"","jetpack_shortlink":"https:\/\/wp.me\/p6XU0A-uh","jetpack_sharing_enabled":true,"_links":{"self":[{"href":"https:\/\/www.sigterritoires.fr\/index.php\/wp-json\/wp\/v2\/posts\/1877","targetHints":{"allow":["GET"]}}],"collection":[{"href":"https:\/\/www.sigterritoires.fr\/index.php\/wp-json\/wp\/v2\/posts"}],"about":[{"href":"https:\/\/www.sigterritoires.fr\/index.php\/wp-json\/wp\/v2\/types\/post"}],"author":[{"embeddable":true,"href":"https:\/\/www.sigterritoires.fr\/index.php\/wp-json\/wp\/v2\/users\/1"}],"replies":[{"embeddable":true,"href":"https:\/\/www.sigterritoires.fr\/index.php\/wp-json\/wp\/v2\/comments?post=1877"}],"version-history":[{"count":0,"href":"https:\/\/www.sigterritoires.fr\/index.php\/wp-json\/wp\/v2\/posts\/1877\/revisions"}],"wp:attachment":[{"href":"https:\/\/www.sigterritoires.fr\/index.php\/wp-json\/wp\/v2\/media?parent=1877"}],"wp:term":[{"taxonomy":"category","embeddable":true,"href":"https:\/\/www.sigterritoires.fr\/index.php\/wp-json\/wp\/v2\/categories?post=1877"},{"taxonomy":"post_tag","embeddable":true,"href":"https:\/\/www.sigterritoires.fr\/index.php\/wp-json\/wp\/v2\/tags?post=1877"}],"curies":[{"name":"wp","href":"https:\/\/api.w.org\/{rel}","templated":true}]}}