ArcGISPython实习指导

1、实习一 Python基础学习1.什么是python Python是一种不受局限、跨平台的开源编程语言，它功能强大且简单易学。因而得到了广泛应用和支持。 ArcGIS 9.0社区中引入了Python。此后，Python被视为可供地理处理用户选择的脚本语言并得以不断发展。每个版本都进一步增强了Python体验，从而为您提供更多的功能以及更丰富、更友好的Python体验。 ESRI已将Python完全纳入ArcGIS中，并将其视为可满足用户需求的语言。下面仅介绍Python的部分优势：· 易于学习，非常适合初学者，也特别适合专家使用· 可伸缩程度高，适于大型项目或小型的一次性程序

2、（称为脚本）· 可移植，跨平台· 可嵌入（使ArcGIS可脚本化）· 稳定成熟· 用户社区规模大 Python已延伸到ArcGIS中，成为了一种用于进行数据分析、数据转换、数据管理和地图自动化的语言，因而有助于提高工作效率。2.在ArcGIS中使用Python（1）打开ArcMap，在工具栏中点击“”按钮，位置如图1-1所示，即可打开Python编辑窗口。图1-1（2）Python程序为解释运行，输入一行代码，按“回车”键后即可运行，如图1-2。图1-23.行和缩进学习Python与其他语言最大的区别就是，Python的代码块不使用大括号（）来控制类，函

3、数以及其他逻辑判断。python最具特色的就是用缩进来写模块。（1）缩进的空白数量是可变的，但是所有代码块语句必须包含相同的缩进空白数量，这个必须严格执行。在Python窗口中输入如图1-3所示的代码块并运行，该代码块是严格缩进的。图 1-3（2）以下代码将会执行错误，如图1-4：图1-4实习二 使用ArcPy ArcPy 是一个以成功的 arcgisscripting 模块为基础并继承arcgisscripting功能进而构建而成的站点包。目的是为以实用高效的方式通过Python执行地理数据分析、数据转换、数据管理和地图自动化创建基础。 该包提供了丰富纯正的Python体验，具有代码自动完成

4、功能（输入关键字和点即可获得该关键字所支持的属性和方法的弹出列表；从中选择一个属性或方法即可将其插入），并针对每个函数、模块和类提供了参考文档。 在Python中使用ArcPy的另一个主要原因是，Python是一种通用的编程语言。Python是一种支持动态输入的解释型语言，适用于交互式操作以及为称为脚本的一次性程序快速制作原型，同时其具有编写大型应用程序的强大功能。用ArcPy 编写的ArcGIS应用程序的优势在于，可以使用由来自多个不同领域的GIS专业人员和程序员组成的众多Python小群体开发的附加模块。1.通过Python使用工具 每个地理处理工具都具有一组固定的参数，这些参数为工具提供

5、执行所需的信息。在Python中使用工具时，必须正确设置工具的参数值，以便在脚本运行时工具可以执行。下面以缓冲区分析为例练习通过Python使用工具的方法。（1）点缓冲1）在ArcMap中打开实验数据中的“point.shp”文件，如图2-1。图2-12）打开Python窗口，在其中输入图2-2所示的代码，其中workspace的路径根据数据在电脑中的实际路径设置。在调用工具（图中红色部位）的时候会在窗口右侧出现提示信息，根据提示信息，输入相关参数，运行结果如图2-3。运行后会在目录中生成一个“pointBuffer.shp”文件，这个文件的文件名和路径均在代码中设置。图2-2图2-3（2）线

6、缓冲1）在ArcMap中打开实验数据中的“line.shp”文件，选中其中的一条线，如图2-4。图2-42）打开Python窗口，在其中输入图2-5所示的代码，运行结果如图2-6。图2-5图2-63）取消选中，重新运行一次代码，运行结果如图2-7所示。从两次运行结果可以看出，地理处理工具处理对象首先应是选中的要素，如果没有选中要素，则默认为图层中所有对象。图2-7（3）多层缓冲 在输入要素周围的指定距离内创建多个缓冲区。使用缓冲距离值可随意合并和融合这些缓冲区，以便创建非重叠缓冲区。 在ArcMap中打开实验数据中的“line.shp”文件，运行如图2-8所示的程序，结果如图2-9。多个缓冲区

7、的距离为“10，20，30”，可根据实际情况自己设定。图2-8图2-92.使用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_nam

8、es, where_clause, spatial_reference, explode_to_points, sql_clause)更新或删除行读取几何（1）读取点1）在ArcMap中打开“point.shp”文件，点的相对位置如图2-10。图2-102）在Python窗口中运行如图2-11所示的代码，“SHAPEXY”表示一个要素的质心XY坐标，点的质心为其本身。cursor为一个搜索游标，其中“for row in cursor”是用一个for循环来遍历游标中的每一行。运行结果如图2-12所示，红框中的数据为点的坐标。图2-11（2）读取折线或面1）在ArcMap中打开实验数据中的“li

9、ne.shp”文件，两个线要素相对位置如图2-13。图 2-132）在Python窗口中运行如图2-14所示的代码，结果如图2-15所示。红色方框中的数据为两个先要素的折点坐标。搜索游标中“OID”表示要素的FID，“SHAPE”表示几何要素对象。图2-14图2-15写入几何通过使用插入游标可以实现写入几何功能，（1）写入点 将“point.shp”加载到ArcMap中，在Python窗口中输入图2-16所示的代码并运行，运行结束后点击左下方的“”按钮，刷新一下地图窗口，结果如图2-17。在写入点要素时，只有单个点对象用于设置点要素几何。使用“SHAPEXY”，就可轻松地创建点。图 2-16图

10、2-17（2）写入线 将“line.shp”加载到ArcMap中，在Python窗口中输入图2-18所示的代码并运行，运行结果如图2-19。代码中“arcpy.Array”表示数组对象可包含任意数量的地理处理对象，例如点、几何或空间参考。（点坐标可自行设置）图2-18图 2-19实习三 创建ArcGIS Python Add-in 加载项是一种自定义项，它可以插入到 ArcGIS for Desktop 应用程序（即 ArcMap、ArcCatalog、ArcGlobe 和 ArcScene）中以便提供补充功能以完成自定义任务，例如工具条上的工具集合。1.创建加载项项目（1）选择加载项项目文件

11、夹：解压“addin_assistant.zip”（图3-1）文件，在解压后的文件夹处找到 addin_assistant.exe （图3-2）文件并双击。将打开浏览文件夹对话框，在其中输入用来存储加载项项目的文件夹。必须选择一个空文件夹或创建一个新文件夹（图3-3），单击确定。图 3-1图 3-2图3-3（2）输入项目设置：选择工作文件夹之后，将显示向导的第一个面板（图3-4）。可在此处输入项目设置。这些设置是项目的元数据，存储在 config.xml 文件中。使用加载项时，这些设置在部署期间十分重要。 图3-42.创建工具条 单击加载项内容选项卡，右键单击工具条，然后单击“新建工具条”，如

12、图3-5。向导的右侧面板将显示工具条的属性，如图3-6。工具条必须具有标题和ID。标题用于为工具条提供标注。ID 表示唯一的名称，用于区分不同工具条。ID不能包含空格。默认情况下，启动应用程序时初始显示工具条这一选项处于选中状态。如果不希望在启动应用程序时显示工具条，可以取消选中此选项。图3-5图3-63.创建按钮 （1）输入工具条属性后，即可创建新按钮。右键单击名为“工具条”的新工具条并选择“新建按钮”，如图3-7。图3-7（2）新建按钮后，为按钮键入相应的属性信息，如图3-8，图3-8按钮具有多个属性需要设置。以下列出了所有这些属性及其对应的说明。这些属性存储在项目的 config.xml

13、 文件中属性描述标题（必填）定义按钮标题。下图显示了 ArcGIS 加载项管理器中使用的标题，有助于标识不同类型的可用自定义内容（加载项类型在括号中标识）：类（必填）在桌面应用程序中单击按钮后执行的 Python 类。Python 类是写入按钮业务逻辑的位置。构造类时应使用 Python 单词首字母大写命名约定，例如应使用“ZoomToSelectedFeatures”而不是“zoomtoselectedfeatures”。ID（必填）用于标识按钮的唯一名称。您可以为给定项目创建多个按钮，此 ID 便用来区分不同的按钮。理想情况下，应使用更有意义的值替换默认 ID。ID 不应包含任何空格。可以

14、使用下划线连接单词。不应使用 Python 关键字。有关保留字信息，请查阅 Python 文档。工具提示（可选）将鼠标指针悬停在桌面应用程序中的按钮上方时出现的简短描述。消息（可选）按钮功能的详细描述。将鼠标指针悬停在按钮上方时，此消息会显示在工具提示下方。图像（可选）此图像应为 16 * 16 像素的图像，用于符号化按钮。图像格式应为任意一种常用图片格式（即.bmp、.jpg等）。图像将复制到加载项项目中创建的图像文件夹中。4.编辑 Python 脚本 此时已为定义加载项按钮所需的属性添加了值。此工作流的下一步是编辑 Python 脚本并更新 Python 类，以包含“缩放至所选要素”这一功

15、能。如果未选择任何要素，则会缩放至所有图层的全图。要将此功能添加到自定义按钮，执行以下步骤：（1） 在arcgis软件目录下找到如图3-9所示的Python编辑器，打开通过向导创建的工作文件夹下的“install”文件夹中Python脚本。 图3-9（2） 在onClick(self)函数中键入相关代码，如图3-10.图3-10 （3） 保存脚本，如果出现图3-11所示的警告信息，则需要在脚本头上加上“# coding gbk”如图3-12.图 3-11图3-125.测试按钮（1）创建加载项文件： 转到工作文件夹，双击该文件夹以运行 makeaddin.py （图3-13）脚本。此脚本将加载项

16、所需的全部文件和文件夹都复制到工作文件夹中的加载项压缩文件中。文件名包括工作文件夹名称和.esriaddin 扩展名（图3-14）。图3-13图3-14（2）安装加载项：双击新的加载项文件以启动 ESRI ArcGIS 加载项安装工具。调用此工具后，它分析加载项文件内的 config.xml，并显示创建加载项时输入的项目设置以及加载项是否包含受信数字签名（如图3-15）。图3-15（3）使用加载项：此时可使用为 ArcGIS for Desktop 应用程序设计的自定义功能。启动桌面应用程序 (ArcMap)。自定义工具条（）可能已经处于可见状态并可用于测试。如果不可见，可转到自定义菜单，然后

17、单击加载项管理器。加载项管理器对话框将列出针对当前应用程序安装的加载项。应显示作为项目设置而输入的加载项信息，如名称、描述和图像。如果加载项列于管理器中，可单击自定义按钮打开自定义 对话框。要将工具条添加到应用程序，可单击工具条选项卡，然后选择您所创建的工具条。（4）将示例数据加载到ArcMap中，开始编辑后任意选中一个要素（如图3-16），单击自定义工具条中的“”按钮，地图会自动缩放到选中要素的位置（如图3-17）。图3-16图3-17实习四 宗地四至提取1. 四至提取步骤 界址线表示地块的边界线，界址点表示界址线的转折点或交点。界址点、线是确定地块位置的重要资料。界址点、线属性信息要采集的

18、重要内容。属性信息中，规定的界标类型有九种，分别是钢钉、水泥桩、石灰桩、喷涂、瓷标志、无标识和其他界标；界址线类别有九种，分别是田埂、沟渠、道路、行树、围墙、墙壁、栅栏、两点连线；界址线位置分为内中外三种，分别对应界址线地物的自有、共有和他有；毗邻地物权利人代表界址线两侧对应也是需四的承包地块方名称，中间用逗号隔开。四至提取就是提取宗地四邻（东南西北）的地块权属，或者相邻道路。 1）提取界址点的步骤：循环读取地块要素类的地块，取出每个地块的所有节点，若点要素类中没有坐标相同的点，则将节点存入点要素类中。界址点位置由农户种地时自行约定，实际上表现为田埂的拐角，无明显地物标志，故界标类型默认设为“

19、其他界标”。2）界址线提取的步骤：循环读取地块，将面状的地块转换为多段线，再将多段线转换为单个的线段，若线要素类中没有相同的线段，则将线段作为界址线存入界址线要素类中。查找与界址线相交的地块并取出地块的“CBFMC”属性，将属性赋值给界址线的毗邻地物权利人属性。界址线类别和界址线位置的判断标准是，如果与界址线相交的地块有两个则界址线类别为“田埂”，界址线位置为“中”；如果与界址线相交的地块只有一个则界址线类别为“道路”，界址线位置为“内”。3）四至提取中心法：图4-1，中间的点为地块的中心或质心，四至区域通过与x轴成45°和135°的两条直线划分，周围地块的中心落入哪个四至

20、区域则选择其作为地块对应的四至地块。图 4-1具体步骤为： 循环读取地块图层，获取地块的中心点 对当前地块作缓冲区，根据叠置关系找到当前地块周围的地块，并求得周围地块的中心； 连接当前地块的中心与周围地块的中心，根据斜率判断地块的四至； 对地块的四至信息进行更新。 4）四至提取界址线法：地块对应的四至地块与该地块东、南、西、北四个方向最远四个界址点相接。而且其中的每个界址点对应该地块的两条界址线，斜率绝对值较小的界址线为该地块与其南至和北至地块共有的界址线；斜率绝对值较大的界址线则为该地块与其东至和西至地块共有的界址线。界址线中的毗邻地物权利人属性中有两个名字，一个是该地块的承包方名称，另一个

21、则为其毗邻承包地块权利人姓名；如果毗邻地物权利人中只有一个名字，该地块在这条界址线处没有毗邻地块，四至信息默认设置为“道路”。具体步骤是： 循环读取地块图层，取出一个地块的所有节点，将其存入一个链表中； 通过Y坐标值的比较找到地块南北方向上两个最远界址点；通过X坐标值比较地块东西方向两个最远界址点，将这些点以东南西北的顺序进行存储； 通过九交模型查找与每个方向最远界址点相接的两条界址线，然后通过斜率判断与地块四至有关的四条界址线； 根据这四条界址线寻找本地块毗邻地块得到四至信息并写入地块四至属性之中。2.代码实现1）设置默认地理数据库：在Python中如果没有设置工作环境（env.worksp

22、ace），默认的环境为“默认地理数据库”。在目录中右键单击“test.mdb”，选择“设为默认地理数据库”，如图4-2。图 4-22）根据第一小节中的界址点、线与四至提取步骤在Python窗口中进行代码编写和调试；3）代码成功运行后，根据实习三中的步骤，自定义一个工具条和四个按钮如图4-3，将Python窗口中的代码复制到与按钮对应的函数中。图4-34）提取完成后打开“jzd”图层属性表，如图4-4，图4-45)提取完成后打开“jzx”图层属性表，如图4-5图4-55)提取完成后打开“dk”图层属性表，如图4-6,图4-6附：参考代码# coding:gbkimport arcpyimport

23、 pythonaddinsclass Bound1(object): """Implementation for bound1_addin.button (Button)""" def _init_(self): self.enabled = True self.checked = False def onClick(self): curdk = arcpy.da.SearchCursor("dk","OID", "SHAPE","DKBM","

24、CBFMC","") for row in curdk: mdk = row1 cbfmc = row3 for part in row1: arrpoint = for pnt in part: arrpoint.append(pnt) boundpoints = #存储四个方向上的最远点 eastPoint = arrpoint0 southPoint = arrpoint0 westPoint = arrpoint0 northPoint = arrpoint0 #根据坐标判求得四个方向最远点 for mpnt in arrpoint: if northPo

25、int.Y <= mpnt.Y: northPoint = mpnt if southPoint.Y >= mpnt.Y: southPoint = mpnt if eastPoint.X <= mpnt.X: eastPoint = mpnt if westPoint.X >= mpnt.X: westPoint = mpnt boundpoints.append(eastPoint) boundpoints.append(southPoint) boundpoints.append(westPoint) boundpoints.append(northPoint)

26、containedjzxs = #存储与四个点相交的8条界址线 cjzxqlrs = #存储与四个点相交的8条界址线的毗邻地物权利人 for boundPoint in boundpoints: njzx = arcpy.da.SearchCursor("jzx","OID", "SHAPE","PLDWQLR","") for crowjzx in njzx: jzxpoint = cjzxqlr = crowjzx2 jzx = crowjzx1 for partjzx in crowjzx

27、1: for jzxpnt in partjzx: jzxpoint.append(jzxpnt) ptstart = jzxpoint0 ptend = jzxpoint1 isStartTouch = ptstart.touches(mdk) isEndTouch = ptend.touches(mdk) isContain = (isStartTouch and isEndTouch) isTouch = jzx.touches(boundPoint) isBoundjzx = (isContain and isTouch) if isBoundjzx : containedjzxs.a

28、ppend(jzx) cjzxqlrs.append(cjzxqlr) del njzx dksz = default = "道路" #根据界址线的斜率确定哪个是与四至相关的界址线 #地块东至 if abs(containedjzxs0.firstPoint.Y - containedjzxs0.lastPoint.Y)/(containedjzxs0.firstPoint.X - containedjzxs0.lastPoint.X) >= abs(containedjzxs1.firstPoint.Y - containedjzxs1.lastPoint.Y)/(

29、containedjzxs1.firstPoint.Y - containedjzxs1.lastPoint.Y): mboundStrings = cjzxqlrs0.split(",") if len(mboundStrings) = 1: dksz.append(default) else: for bString in mboundStrings: if not(bString = cbfmc): dksz.append(bString) else: mboundStrings = cjzxqlrs1.split(",") if len(mbou

30、ndStrings) = 1: dksz.append(default) else: for bString in mboundStrings: if not(bString = cbfmc): dksz.append(bString) #地块南至 if abs(containedjzxs2.firstPoint.Y - containedjzxs2.lastPoint.Y)/(containedjzxs2.firstPoint.X - containedjzxs2.lastPoint.X) >= abs(containedjzxs3.firstPoint.Y - containedjz

31、xs3.lastPoint.Y)/(containedjzxs3.firstPoint.Y - containedjzxs3.lastPoint.Y): mboundStrings = cjzxqlrs3.split(",") if len(mboundStrings) = 1: dksz.append(default) else: for bString in mboundStrings: if not(bString = cbfmc): dksz.append(bString) else: mboundStrings = cjzxqlrs2.split(",&

32、quot;) if len(mboundStrings) = 1: dksz.append(default) else: for bString in mboundStrings: if not(bString = cbfmc): dksz.append(bString) #地块西至 if abs(containedjzxs4.firstPoint.Y - containedjzxs4.lastPoint.Y)/(containedjzxs4.firstPoint.X - containedjzxs4.lastPoint.X) >= abs(containedjzxs5.firstPoi

33、nt.Y - containedjzxs5.lastPoint.Y)/(containedjzxs5.firstPoint.Y - containedjzxs5.lastPoint.Y): mboundStrings = cjzxqlrs4.split(",") if len(mboundStrings) = 1: dksz.append(default) else: for bString in mboundStrings: if not(bString = cbfmc): dksz.append(bString) else: mboundStrings = cjzxql

34、rs5.split(",") if len(mboundStrings) = 1: dksz.append(default) else: for bString in mboundStrings: if not(bString = cbfmc): dksz.append(bString) #地块北至 if abs(containedjzxs6.firstPoint.Y - containedjzxs6.lastPoint.Y)/(containedjzxs6.firstPoint.X - containedjzxs6.lastPoint.X) >= abs(conta

35、inedjzxs7.firstPoint.Y - containedjzxs7.lastPoint.Y)/(containedjzxs7.firstPoint.Y - containedjzxs7.lastPoint.Y): mboundStrings = cjzxqlrs7.split(",") if len(mboundStrings) = 1: dksz.append(default) else: for bString in mboundStrings: if not(bString = cbfmc): dksz.append(bString) else: mbou

36、ndStrings = cjzxqlrs6.split(",") if len(mboundStrings) = 1: dksz.append(default) else: for bString in mboundStrings: if not(bString = cbfmc): dksz.append(bString) #更新四至 curupdk = arcpy.da.UpdateCursor("dk","DKDZ","DKXZ","DKNZ","DKBZ","

37、SHAPE","") for rowup in curupdk: if rowup4.equals(mdk): rowup0 = dksz0 rowup1 = dksz2 rowup2 = dksz1 rowup3 = dksz3 curupdk.updateRow(rowup) del curupdk del curdk passclass Bound2(object): """Implementation for bound2_addin.button (Button)""" def _init_(s

38、elf): self.enabled = True self.checked = False def onClick(self): #中心法 curdk = arcpy.da.SearchCursor("dk","OID", "SHAPE","CBFMC","") for row in curdk: overdk = #存储当前地块周围地块 overcbfmc = #存储当前地块周围地块承包方名称 dkdz = dknz = dkxz = dkbz = mdk = row1 cbfmc = ro

39、w2 bmdk = mdk.buffer(0.5) #对当前地块作缓冲区分析 #获取当前地块周围地块 alldk = arcpy.da.SearchCursor("dk","OID", "SHAPE","CBFMC","") for arow in alldk: if not mdk.equals(arow1): isOverlap = bmdk.overlaps(arow1) if isOverlap : overdk.append(arow1) overcbfmc.append(arow2)

40、 del alldk i = 0 #获取四个方向的地块 while i < len(overdk) : k = abs(mdk.centroid.Y - overdki.centroid.Y)/(mdk.centroid.X - overdki.centroid.X) if k <= 1 : if mdk.centroid.X > overdki.centroid.X : dkxz.append(overcbfmci) else: dkdz.append(overcbfmci) else : if mdk.centroid.Y > overdki.centroid.Y

41、: dknz.append(overcbfmci) else: dkbz.append(overcbfmci) i = i + 1 #如果没有则默认为道路 if len(dknz) = 0 : dknz.append("道路") if len(dkdz) = 0 : dkdz.append("道路") if len(dkxz) = 0 : dkxz.append("道路") if len(dkbz) = 0 : dkbz.append("道路") #更新四至 upCur = arcpy.da.UpdateCurso

42、r("dk","DKDZ","DKXZ","DKNZ","DKBZ","SHAPE","") for rowup in upCur: if rowup4.equals(mdk): rowup0 = dkdz0 rowup1 = dkxz0 rowup2 = dknz0 rowup3 = dkbz0 upCur.updateRow(rowup) del upCur del curdk passclass JZD(object): ""&q

43、uot;Implementation for JZD_addin.button (Button)""" def _init_(self): self.enabled = True self.checked = False def onClick(self): #提取界址点 jzd = #循环读取每个地块获取所有的界址点 curdk = arcpy.da.SearchCursor("dk","OID", "SHAPE","CBFMC","") for row in c

44、urdk: cbfmc = row2 for part in row1: arrpoint = for pnt in part: arrpoint.append(pnt) for cpoint in arrpoint: isSaved = False for cjzd in jzd : if cjzd.equals(cpoint): isSaved = True if not isSaved : jzd.append(cpoint) del curdk #插入界址点 curInsert = arcpy.da.InsertCursor("jzd","SHAPE","YSDM","JZDLX","JBL