借助python构建自定义工具之推断并删除重复要素.docVIP

FOR U I ONLY(仅供参考) 借助python构建自定义工具之判断并删除重复要素 tj051181@ /文 问题： 对输入的数据，按照图形或者指定字段进行判断重复，并由用户来决定是否删除 粗略分析，工具界面大致需要输入项： 1.Input Feature Class：可能是shapefile，也可能是GDB要素集。为了不影响源数据，统一对源数据拷贝一份，对拷贝后数据操作，作为本工具的输出结果。 2.Duplicate Field(s)：可能需要根据图形，也可能根据某个（些）字段。不管是那种，字段都是输入文件或者要素集含的字段，所以这里设想能够关联Input Feature Class获取其字段，让用户勾选，类似Delete Field 3.Delete Duplicate Features：用户定义是否删除，所以这项定义为Boolean类型，界面上为勾选（check） 界面预览： 编写脚本： 1.python编写环境ArcGIS Desktop Help里介绍的是 pythonwin，但是目前自带的只有IDLE，可以根据当前python版本down一个适合的pythonwin，也可以就用IDLE（Python GUI） 2.参考  HYPERLINK /arcgisdesktop/9.3/index.cfm?TopicName=Creating_a_new_Python_script_module ArcGIS Desktop Help\Geoprocessing\Automating your work with scripts\Getting started with writing geoprocessing scripts\Creating a new Python script module 此部分介绍如何创建一个python模块 3.具体函数 在这个脚本里面，主要用到了这几个函数 3.1文本处理函数 us，作用是去掉s 前后空格，并将其转为大写输出 #trim and ucase a string def us(s): return string.upper(strip(str(s))) 3.2判断要素集fc中是否存在某字段field，返回True或者False #check if field exists in FeatureClass fc def blnExistField(fc,field): flag = False fds = gp.ListFields(fc) fds.Reset() fd = fds.Next() while fd: if us() == us(field): flag = True fd = None else: fd = fds.Next() del fds,fd return flag 3.3 对要素集遍历，计算每一个要素的四至值（Geometry.Extent MinX,MinY,MaxX,MaxY） #calculate MinX;MinY;MaxX;MaxY def cal(fc): #add MinX;MinY;MaxX;MaxY fds = MinX;MinY;MaxX;MaxY for i in range(0,len(fds.split(;))): fd = fds.split(;)[i] if blnExistField(fc,fd) == False: gp.addfield_management(fc,fd,Double) #calculate MinX;MinY;MaxX;MaxY rows = gp.UpdateCursor(fc) rows.Reset() row = rows.Next() while row: Geo = row.GetValue(shapeFieldName) for i in range(0,len(fds.split(;))): fd = fds.split(;)[i] row.SetValue(fd,Geo.extent.split()[i]) rows.UpdateRow(row) row = rows.Next() del rows, row, fds, fd return MinX;MinY;MaxX;M