ArcGISPython实习指导

实习一Python基础学习

1.什么是python

Python是一种不受局限、跨平台的开源编程语言，它功能强大且简单易学。因而得到了广泛应用和支持。

ArcGIS 9.0社区中引入了Python。此后，Python被视为可供地理处理用户选择的脚本语言并得以不断发展。每个版本都进一步增强了Python体验，从而为您提供更多的功能以及更丰富、更友好的Python体验。

ESRI已将Python完全纳入ArcGIS中，并将其视为可满足用户需求的语言。下面仅介绍Python的部分优势：

•易于学习，非常适合初学者，也特别适合专家使用

•可伸缩程度高，适于大型项目或小型的一次性程序（称为脚本）

•可移植，跨平台

•可嵌入（使ArcGIS可脚本化）

•稳定成熟

•用户社区规模大

Python已延伸到ArcGIS中，成为了一种用于进行数据分析、数据转换、数据管理和地图自动化的语言，因而有助于提高工作效率。

2.在ArcGIS中使用Python

（1）打开ArcMap，在工具栏中点击“”按钮，位置如图1-1所示，即可打开Python编辑窗口。

图1-1

（2）Python程序为解释运行，输入一行代码，按“回车”键后即可运行，如图1-2。

图1-2

3.行和缩进

学习Python与其他语言最大的区别就是，Python的代码块不使用大括号（{}）来控制类，函数以及其他逻辑判断。python最具特色的就是用缩进来写模块。（1）缩进的空白数量是可变的，但是所有代码块语句必须包含相同的缩进空白数量，这个必须严格执行。在Python窗口中输入如图1-3所示的代码块并运行，该代码块是严格缩进的。

图1-3

（2）以下代码将会执行错误，如图1-4：

图1-4

实习二使用ArcPy

ArcPy 是一个以成功的arcgisscripting 模块为基础并继承arcgisscripting功能进而构建而成的站点包。目的是为以实用高效的方式通过Python执行地理数据分析、数据转换、数据管理和地图自动化创建基础。

该包提供了丰富纯正的Python体验，具有代码自动完成功能（输入关键字和点即可获得该关键字所支持的属性和方法的弹出列表；从中选择一个属性或方法即可将其插入），并针对每个函数、模块和类提供了参考文档。

在Python中使用ArcPy的另一个主要原因是，Python是一种通用的编程语言。Python是一种支持动态输入的解释型语言，适用于交互式操作以及为称为脚本的一次性程序快速制作原型，同时其具有编写大型应用程序的强大功能。用ArcPy 编写的ArcGIS应用程序的优势在于，可以使用由来自多个不同领域的GIS 专业人员和程序员组成的众多Python小群体开发的附加模块。

1.通过Python使用工具

每个地理处理工具都具有一组固定的参数，这些参数为工具提供执行所需的信息。在Python中使用工具时，必须正确设置工具的参数值，以便在脚本运行时工具可以执行。下面以缓冲区分析为例练习通过Python使用工具的方法。（1）点缓冲

1）在ArcMap中打开实验数据中的“point.shp”文件，如图2-1。

图2-1

2）打开Python窗口，在其中输入图2-2所示的代码，其中workspace的路径根据数据在电脑中的实际路径设置。在调用工具（图中红色部位）的时候会在窗口右侧出现提示信息，根据提示信息，输入相关参数，运行结果如图2-3。运行后会在目录中生成一个“pointBuffer.shp”文件，这个文件的文件名和路径均在代码中设置。

图2-3

（2）线缓冲

1）在ArcMap中打开实验数据中的“line.shp”文件，选中其中的一条线，如图2-4。

2）打开Python窗口，在其中输入图2-5所示的代码，运行结果如图2-6。

图2-5

图2-6

3）取消选中，重新运行一次代码，运行结果如图2-7所示。从两次运行结果可以看出，地理处理工具处理对象首先应是选中的要素，如果没有选中要素，则默认为图层中所有对象。

图2-7

（3）多层缓冲

在输入要素周围的指定距离创建多个缓冲区。使用缓冲距离值可随意合并和融合这些缓冲区，以便创建非重叠缓冲区。

在ArcMap中打开实验数据中的“line.shp”文件，运行如图2-8所示的程序，结果如图2-9。多个缓冲区的距离为“[10，20，30]”，可根据实际情况自己设定。

图2-8

图2-9

2.使用Python处理游标

游标是一种数据访问对象，可用于在表中迭代一组行或者向表中插入新行。游标有三种形式：搜索、插入或更新。

arcpy.da.InsertCursor(in_table, field_names) 插入行

arcpy.da.SearchCursor(in_table, field_names, {where_clause}, {spatial_reference}, {explode_to_points}, {sql_clause}) 只读访问

arcpy.da.UpdateCursor(in_table, field_names, {where_clause}, {spatial_reference}, {explode_to_points}, {sql_clause}) 更新或删除行

读取几何

（1）读取点

1）在ArcMap中打开“point.shp”文件，点的相对位置如图2-10。

图2-10

2）在Python窗口中运行如图2-11所示的代码，“SHAPEXY”表示一个要素的质心XY坐标，点的质心为其本身。cursor为一个搜索游标，其中“for row in cursor”是用一个for循环来遍历游标中的每一行。运行结果如图2-12所示，红框中的数据为点的坐标。