用SelectLayerByLocation_management对图层进行选择,该函数需要传入选择图层,被选择图层,以及选择的方式,是否相邻之类的。
关于该函数的使用,大家可以查阅相关的资料,这里就不一一说明。
首先在编写这个程序时,需要对一个区域所有的面进行遍历,然后再与另外区域的所有面进行选择。因此,这中间涉及到面遍历后创建新元素,然后将新元素放到选择函数中,之后还需要删除该图层。这中间比较浪费时间,特别是在创建元素和删除元素的时候,以及进行选择操作时。经过比较后找到相应的id然后更新原来的传入的数据。具体的实现看一下源代码。
# coding:utf-8
import arcpy
import os
import sys
reload(sys)
sys.setdefaultencoding('utf-8')
#交换id
def switchID(selectLayer,selectID,selectedLayer,selectedID):
#处理选择图层
with arcpy.da.UpdateCursor(selectLayer,["ID","JBX"]) as cursor:
for row in cursor:
tmpSelectID=row[0].encode("utf-8")
if tmpSelectID==selectID:
row[1]=selectedID
cursor.updateRow(row)
break
del cursor
#处理被选择图层
with arcpy.da.UpdateCursor(selectedLayer,["ID","JBX"]) as cursor:
for row in cursor:
tmpSelectID=row[0].encode("utf-8")
if tmpSelectID==selectedID:
row[1]=selectID
cursor.updateRow(row)
break
del cursor
def delShpFile(root):
extends=[".cpg",".dbf",".prj",".sbn",".sbx",".shp",".shp.xml",".shx"]
files=[]
for extend in extends:
tmpFile=root+extend
if os.path.isfile(tmpFile):
files.append(tmpFile)
for file in files:
os.remove(file)
selectLayer="C:\Users\Desktop\STProject\statistics\测试数据\xx1遥感统计数据库.gdb\STM"
selectedLayer="C:\Users\Desktop\STProject\statistics\测试数据\xx2遥感统计数据库.gdb\STM"
arcpy.MakeFeatureLayer_management(selectLayer,"STM")
arcpy.MakeFeatureLayer_management(selectedLayer,"STM2")
selection=arcpy.SelectLayerByLocation_management("STM",
"BOUNDARY_TOUCHES",
selectedLayer,"","NEW_SELECTION")
selection2=arcpy.SelectLayerByLocation_management("STM2",
"BOUNDARY_TOUCHES",
selectLayer,"","NEW_SELECTION")
cnt=arcpy.GetCount_management(selection);
rootPath="C:\Users\qrb_PC\Desktop\mestShp"
# with arcpy.da.SearchCursor(selectLayer, ["SHAPE@",'ID']) as cursor:
# for row in cursor:
# geometry = row[0]
# arcpy.MakeFeatureLayer_management(geometry,"TmpSTM")
# itemSelect=arcpy.SelectLayerByLocation_management("TmpSTM",selectedLayer,"","NEW_SELECTION")
# tmpCNT = arcpy.GetCount_management(itemSelect);
# del cursor
selectResultPath=rootPath+"\Metd.shp"
arcpy.CopyFeatures_management(selection, selectResultPath)
selectedResultPath=rootPath+"\Metd2.shp"
finalResPath=rootPath+"\Final.shp"
cnt2=arcpy.GetCount_management(selection2)
arcpy.CopyFeatures_management(selection2, selectedResultPath)
# 遍历选择面
with arcpy.da.SearchCursor(selectResultPath, ["SHAPE@",'ID']) as selectCursor:
for selectRow in selectCursor:
out_name = "Select"+selectRow[1] + '.shp'
tmpSelectOutName="Select"+selectRow[1];
selectID=selectRow[1].encode("utf-8")
arcpy.FeatureClassToFeatureClass_conversion(selectRow[0], rootPath, out_name)
# 遍历被选择面
with arcpy.da.SearchCursor(selectedResultPath, ["SHAPE@",'ID']) as selectedCursor:
for selectedRow in selectedCursor:
out_name2 = "SelectED" + selectedRow[1] + '.shp'
tmpSelectedOutName = "SelectED" + selectedRow[1]
creatName="RSelectED" + selectedRow[1]
creatName.encode("utf-8")
arcpy.FeatureClassToFeatureClass_conversion(selectedRow[0], rootPath, out_name2)
tmpResultPath=rootPath +"\\"+ out_name
tmpResultPath2 = rootPath + "\\" + out_name2
#创建临时图层
arcpy.MakeFeatureLayer_management(tmpResultPath, creatName)
resSelection = arcpy.SelectLayerByLocation_management(creatName,"BOUNDARY_TOUCHES",tmpResultPath2)
fCNT = int(arcpy.GetCount_management(resSelection).getOutput(0))
selectedID=selectedRow[1].encode("utf-8")
mypath=rootPath+"\\"+tmpSelectedOutName
if fCNT==1:
print "已经找到"
switchID(selectLayer,selectID,selectedLayer,selectedID)
delShpFile(rootPath+"\\"+tmpSelectedOutName)
else:
delShpFile(rootPath+"\\"+tmpSelectedOutName)
#删除临时图层
arcpy.Delete_management(creatName)
del selectedCursor
delShpFile(rootPath + "\\" + tmpSelectOutName)
del selectCursor