点击箭头处“蓝色字”，关注我们哦！！

问题来源

最近看到群里有小伙伴有这样的需求，大概意思是这样的。他有一个面数据，这么面数据有很多条记录(polygon),我们姑且把这个面叫A数据，还有另外一个B数据，也是有很多条记录(polygon)。现在需要A数据中每一条记录与B面求交，然后找出所有的相交面，在这些相交面中取面积最大者的一些属性填到A面的那条记录的属性表中，基本上流程就走完了。不知道我需求说明讲得清楚么。不过需要注意的是，这里需要遍历每一条记录来求交。

解决办法

好了，下面我把数据拿上来。A数据为下图的impBase.shp矢量数据。其中的DLMC为空数据，这个字段我们需要在B数据中拿过来。

而B数据，则为另外一个矢量面数据。有DLMC字段，现在需要将该字段经过上面的逻辑后取出来，用于填到上面A数据的DLMC字段。

我们来看一下实现代码。在这个例子中，我们将会使用属性表的遍历，在内存中创建要素类对象，添加字段，添加属性值，更新插入属性表。空间分析中的求交，游标遍历，令牌的使用，还有内存里数据的删除。花了一点时间调出来的。大家可以试着自己写一下。具体的逻辑可以参考我写的代码。

 1import arcpy 2imptAttrPath="D:/QQGroups/shp_inserct/impAtt.shp" 3imptBasePath="D:/QQGroups/shp_inserct/impBase.shp" 4 5arcpy.env.overwriteOutput = True 6 7 8outputSR = arcpy.Describe(imptBasePath).spatialReference 91011with arcpy.da.UpdateCursor(imptBasePath, ['FID','SHAPE@','DLMC']) as updateCursor:12    for updateRow in updateCursor:13     #createFC = "in_memory/"+str(updateRow[0])14     createFC = r"in_memory\tempUpdate"15     arcpy.CreateFeatureclass_management(16        "in_memory",17        "tempUpdate",18        "POLYGON", "", "", "",19        outputSR)2021     arcpy.AddField_management(createFC, "TDYT", "TEXT")2223     cur = arcpy.InsertCursor(createFC)24     newRow = cur.newRow()25     newRow.shape = updateRow[1]26     newRow.TDYT=updateRow[2]2728     cur.insertRow(newRow)2930     arcpy.Intersect_analysis([createFC,imptAttrPath],31                              out_feature_class="in_memory/IntersectShapeFile",32                              join_attributes="ALL", cluster_tolerance="#", output_type="INPUT")3334     area=0.035     dlmc=""36     with arcpy.da.SearchCursor("in_memory/IntersectShapeFile", ["Flag", "SHAPE@AREA", "DLMC"]) as maxCursor:37          for maxRow in maxCursor:38              tmpArea=maxRow[1]39              if area<=tmpArea:40                  area=tmpArea41                  dlmc=maxRow[2]4243     updateRow[2]=dlmc44     updateCursor.updateRow(updateRow)4546     #arcpy.Delete_management(["in_memory/"+str(updateRow[0]),"in_memory/IntersectShapeFile"])47     arcpy.Delete_management('in_memory')48print "finished"

在写代码的时候，特别要注意游标遍历方法的应用，字段在游标中是否存在。相交计算中需要注意设置要输出类型，有"LINE","INPUT"，具体可以参考ArcGIS Desktop里面的相交工具中参数。还有内存的删除，使用数组的方式一一删除内存中的数据，可能会出现错误，至少我是没有试成功。

最后，将实现得到这样的结果。

往期回顾

ArcGIS开发&读取txt坐标数据，并创建面shapefile

ArcGIS开发&创建gdb，复制shapefile与更改属性表

ArcGIS开发&多幅影像批量拼接

ArcGIS开发&数据重新分类

【ArcGIS制图】&利用字段计算器编写脚本，对数据分类

【ArcGIS制图】&ArcGIS Catalog目录误删了文件怎么办

【ArcGIS制图】&ArcGIS中导入xy经纬度坐标数据不显示字段

openlayer开发&返回Geoserver发布的wfs指定字段信息

openlayer开发&绘制点、圆、多边形

欢迎加入arcpy开发qq学习群：487352121

做一个gis公众号

长按并关注我们