AGS Python开发-ArcPy开发基础

1、Python与ArcGIS关系
ArcGIS提供了ArcPy包，Python是通过该ArcPy包来操作ArcGIS中的功能。ArcPy可以理解为ArcGIS针对Python语言提供的开发API。
2、导入ArcPy与模块搜索路径
（1）导入ArcPy
示例1：import arcpy
示例2：import arcpy.mapping
示例3：from arcpy import env
示例4：from arcpy import env as ENV
（2）模块搜索路径
导入ArcPy时，会先搜索PYTHONPATH环境变量。如何未定义环境变量，则搜索sys.path定义的路径：Python27ArcGIS10.5Libsite-packages和Python27ArcGISx6410.5Libsite-packages。ArcGIS在安装Python时，会在这两个目录下创建Desktop10.5.pth和Server10.5.pth。这两个文件记录了ArcPy需要的路径：
*C:arcgisDesktop10.5in
C:arcgisDesktop10.5ArcPy
C:arcgisDesktop10.5ArcToolBoxScripts
•
如果出现提示找不到ArcPy模块的错误，可以手动创建该文件并添加上面的路径。
3、调用地理处理工具和地理处理服务
3.1、调用地理处理工具
所有工具可以ArcPy函数方式调用，也可以工具箱别名的模块调用。
函数方式：arcpy.GetCount_management(features)
工具箱模块方式：arcpy.management.GetCount(features)
3.2、调用自定义地理处理工具
建议在ImportToolbox()中定义自定义工具箱的别名，ArcPy使用arcpy.<工具名称>_<别名>方式调用。
arcpy.ImportToolbox("c:/mytools/geometrytools.tbx", "mytools")
arcpy.CreateRegularPolygons_mytools(params)
3.3、调用地理处理服务
地理处理服务器名称和工具箱名称使用分号隔开。
arcpy.ImportToolbox("http://flame7/arcgis/services;GP/BufferByVal","mytools")
arcpy.BufferPoints_mytools(inFeatureSet, "5 km")
3.4、调用地理处理工具的技巧
技巧1：调用ArcGIS自带的地理处理工具
ArcGIS提供了数量众多的地理处理工具，这些工具都可以通过ArcPy来调用。在调用地理处理工具时，需要清楚工具的处理流程，以及输入输出参数等信息。为了方便调用与调试，可按如下步骤进行：
1. 在桌面软件中运行地理处理工具
在工具运行窗口中可以查看工具说明，以及相关参数设置要求。同时还可以测试输入参数与输出结果。
2. 查看地理处理工具调用方式
方式1：在桌面软件中，工具拖拉到Python命令行窗口，即可查看到工具调用名称，以及相关参数设置要求。
方式2：调用help()方法查看工具调用文档。在桌面软件或者IDE环境中的Python命令行中，输入help(arcpy.toolname)查询。
方式3：直接查看桌面软件帮助中的工具参考。在参考中还包含了示例代码。
3. 在Python中调用地理处理工具
根据地理处理工具调用要求，编写Python调用代码。
技巧2：调用自定义的地理处理模型
对于自定义的地理处理模型，在Python中可以使用两种调用方式：
1. 直接调用
直接调用是将自定义处理模型当作整理看待，调用方式：
arcpy.importToolbox("toolbox_path", "mytoolbox")
arcpy.Mytool_mytoolbox("params")
1. 将模型导出为Python脚本，再复用脚本内容
ArcGIS提供了自定义地理处理模型导出为Python脚本的功能，导出的脚本可以直接在python中使用，也可以在其上修改增加自定义的功能。
这种方式可以看作是可视化编码方式，可以减少调用系统地理处理工具的复杂度。
4、ArcPy中的OMD
ArcPy由三大部分组成：ArcPy函数、ArcPy类、ArcPy模块。
（1）ArcPy函数
ArcPy函数用于执行特定的任务，方便执行地理处理工作流。ArcPy中的函数分为两种：地理处理函数和非地理处理函数。所有地理处理工具都以地理处理函数提供。
（2）ArcPy类
ArcPy类存储地理对象相关的信息，通常用于解析地理数据，或者作为地理处理工具的参数。
（3）ArcPy模块
ArcPy模块是将共性的功能以模块方式提供，方便调用。ArcPy模块包括数据访问模块（arcpy.da）、制图模块（arcpy.mapping）、空间分析扩展模块（arcpy.sa）、网络分析扩展模块（arcpy.na）。以及工具类模块：时间模块（arcpy.time）。
4.1、ArcPy中的函数
ArcPy中的函数主要包括了：描述数据的函数（Describe）、判断数据是否存在的函数（Exists）、数据迭代相关的函数（如：ListFields）、外部输入参数相关的函数（如GetParameterAsText）、许可检查相关的函数（如CheckExtension）、环境设置相关的函数（如GetSystemEnvironment）、消息错误处理相关的函数（如GetMessage）
详情参见：http://desktop.arcgis.com/zh-cn/arcmap/latest/analyze/arcpy-functions/alphabetical-list-of-arcpy-functions.htm
4.2、ArcPy中的类
分类主要类说明
要素类或表、栅格类 FeatureSet、RecordSet、Raster 代表要素图层或者表、栅格图层
要素数据内容相关类 Cursor、Row 代表要素图层或者表中记录的每行数据
要素图层字段相关类 Field、FieldInfo 代表要素图层或表的字段信息
几何类型与空间参考 Point、Multipoint、PointGeometry、Polyline、Polygon、Extent、SpatialReference
图表类 Graph、GraphTemplate 如：气泡图、直方图、饼图等
ToolValidator相关类 Filter、Parameter、Schema、Value 自定义Python工具箱所需要的参数类
异常类 ExecuteError、ExecuteWarning
其他核心类环境设置类：env；地理处理结果类：Result；索引类（空间索引和属性索引）：Index；多值参数类：ValueTable；SDE SQL操作类：ArcSDESQLExecute；辅助类：Array、RandomNumberGenerator
4.3、ArcPy中的模块
模块名称主要功能
arcpy.da 数据访问模块。提供了增、删、查、改数据的类和函数。
arcpy.mapping 地图制图模块。提供了操作mxd文档、地图渲染、地图打印输出、地图服务发布等功能。
arcpy.na 网络分析扩展模块。用于访问网络分析地理处理工具。
arcpy.sa 栅格分析模块。用于访问空间分析地理处理工具。
5、描述数据的函数
通过Describe函数可以查询所有类型数据的属性信息。如：数据几何类型、字段、索引等。不同的数据类型可查询的属性信息不同。
ftClass="D:/mapdata/test/china/chinamap.gdb/citiesbuffer"
ftDesc = arcpy.Describe(ftClass)
print(ftDesc.shapeType)
print(ftDesc.featureType)
print(ftDesc.shapeFieldName)
print(ftDesc.OIDFieldName)
print(len(ftDesc.fields))
print(ftDesc.spatialReference.name)
6、环境设置类
ArcPy中环境设置定义在env类中。
（1）环境设置读取和更改
arcpy.env.scratchWorkspace = "c:/mapdata"
arcpy.env.workspace = "c:/mapdata"
arcpy.env.XYTolerance = 2.5
workspace = arcpy.env.scratchGDB
（2）检查环境变量名称
arcpy.ListEnvironments()
（3）重设环境设置
arcpy.ResetEnvironments()
arcpy.ClearEnvironment("workspace")
7、消息处理
（1）读取消息
在地理处理工具执行过程中，会输出处理消息，消息内容包括：时间、参数、错误等。可通过如下方式获取：
方式1：直接读取
result=arcpy.GetCount_management("D:/mapdata/test/worldshp/cities.shp")
print(arcpy.GetMessageCount())
print(arcpy.GetMessages())
print(arcpy.GetMessage(0))
方式2：从result中读取
result=arcpy.GetCount_management("D:/mapdata/test/worldshp/cities.shp")
print(result.messageCount)
print(result.getMessage(0))
print(result.getMessages(2))
（2）自定义消息
如果是控制台输出，可以直接使用print语句，也可以使用ArcPy提供的消息函数：
arcpy.AddMessage("info")
arcpy.AddWarning("warning")
arcpy.AddError("error")
8、错误处理
（1）ExecuteError类
地理处理工具失败时会抛出arcpy.ExecuteError类。其他的异常可以使用sys.exc_info()获取。
import arcpy
import sys

try:
result = arcpy.GetCount_management("c:/erro.shp")

except arcpy.ExecuteError:
arcpy.AddError(arcpy.GetMessages(2))
except:
e = sys.exc_info()[1]
print(e.args[0])
（2）错误跟踪
在较大复杂的脚本中，可结合sys和traceback模块来查找跟踪错误。
import arcpy
import sys
import traceback
arcpy.env.workspace = "C:/Data/myData.gdb"
try:
arcpy.CreateSpatialReference_management()
except arcpy.ExecuteError:
msgs = arcpy.GetMessages(2)
print(msgs)
except:
# Get the traceback object
#
tb = sys.exc_info()[2]
tbinfo = traceback.format_tb(tb)[0]
pymsg = "PYTHON ERRORS: Traceback info: " + tbinfo + " Error Info: " + str(sys.exc_info()[1])
msgs = "ArcPy ERRORS: " + arcpy.GetMessages(2) + " "
print(pymsg)
print(msgs)
（3）Result对象中获取错误消息
result = arcpy.Clip_analysis("roads", "urban_area", "urban_roads")
消息总数：result.messageCount
按索引号获取消息：result.getMessage(index)
按类型获取消息：result.getMessages(severity)
消息类型分：0（消息）、1（警告）、2（错误）
9、许可检查
使用arcpy.CheckExtension("3D")和arcpy.CheckOutExtension("3D")检查许可。
import arcpy
try:
if arcpy.CheckExtension("3D") == "Available":
arcpy.CheckOutExtension("3D")
else:
print("license error")
arcpy.env.workspace = "D:/GrosMorne"
arcpy.HillShade_3d("WesternBrook", "westbrook_hill", 300)
arcpy.Aspect_3d("WesternBrook", "westbrook_aspect")
except arcpy.ExecuteError:
print(arcpy.GetMessages(2))
finally:
arcpy.CheckInExtension("3D")
10、地理数据路径设置与中文处理
（1）路径设置
在访问本地数据时，常需要设置数据存放路径，设置方式有两种：
方式1：绝对路径。
直接将绝对路径传入地理处理工具函数中。
arcpy.GetCount_management(r"c: empstreams.shp")