Python Scripting for ArcGIS - Paul Zandbergen

118 Pages·2011·2.05 MB·English

by  Paul Zandbergen

Python Scripting for ArcGIS
Paul Zandbergen
Deppartment of Geoggrapphyy
University of New Mexico
OOuuttlliinnee ooff TTooppiiccss
• Introdduction
– Examples, Python and ArcGIS, Python versions
• FFuunnddaammeennttaallss ooff ggeeoopprroocceessssiinngg iinn AArrccGGIISS
• Python language fundamentals
– Where to run Python code
– Data types: numbers, strings, lists
– Functions and modules
– Controllingg workflow
• ArcPy: Geoprocessing using Python
– Using tools, functions, classes
– Describbing ddata, llisting ddata, workking withh llists
• Creating custom tools
– Scrippt tools,, tool pparameters
• Resources
WWoorrkksshhoopp MMaatteerriiaallss PPoosstteedd
ppoosstteedd uunnttiill OOccttoobbeerr 2244
hhttp:////www.paullzanddbbergen.com//workkshhops
FFoorrtthhccoommiinngg BBooookk
• Python Scripting for ArcGIS
• EEssrrii PPrreessss
• Sometime in 2012
• Updated for ArcGIS 10.1
• Sample exercises posted (for 10.0)
Introduction
PPrriioorr KKnnoowwlleeddggee aanndd EExxppeerriieennccee
• Using ArcGIS 9.3 or 10.0?
– Workshopp is for 10.0
• Prior Python experience?
– I’m not assumingg anyy
• OOthher programmiing experiience??
– I’m not assuming any
EExxaammppllee 11
• Script to copy all shapefiles in a folder into a geodatabase
import arcpy
from arcpy import env
env.overwriteOutput = True
env.workkspace = ""c://workkshhop//ex0011""
fclist = arcpy.ListFeatureClasses()
for fc in fclist:
fcdesc = arcppyy.Describe((fc))
arcpy.CopyFeatures_management(fc, "c:/workshop/ex01/study.mdb/"
+ fcdesc.basename)
EExxaammppllee 22
• Script tool to generate a k‐nearest neighbor table
• Runs an existing ArcGIS tool multiple times, writes the result
import arcpy
from arcpy import env
env.overwriitteouttputt = TTrue
infc = arcpy.GetParameterAsText(0)
output = arcpy.GetParameterAsText(1)
kk == aarrccppyy.GGeettPPaarraammeetteerr((22))
n = 1
f = open(output, "w")
wwhhiillee nn <<= kk::
result = arcpy.CalculateDistanceBand_stats(infc, n)
f.write(str(n) + " " + str(result[1])+ "\n")
nn = nn ++ 11
f.close()
EExxaammppllee 33
• Script tool to run Huff model
• Sophisticated analysis not available in ArcGIS
EExxaammppllee 33