Latest Contributions by JohannesBierer

Hello, I try to find excactly identical points in a polygon feature class with the 10.4. basic license to search for polygon errors with python 2.7. The following script can't find end points in the row[1].lastPoint in line 18? I need to exclude start and end points. What's wrong or any ideas how to do better? import arcpy infc = r"R:\XXX\LSG_Fehler.shp" outfc = r"C:\temp2\GleicheKoord.shp" CoordList = list() arcpy.env.overwriteOutput = True # Enter for loop for each feature # for row in arcpy.da.SearchCursor(infc, ["OID@", "SHAPE@"]): # Print the current multipoint's ID # print("Feature {}:".format(row[0])) partnum = 0 if row[1] == row[1].lastPoint: print row[1].lastPoint break # Step through each part of the feature else: for part in row[1]: # Print the part number # print("Part {}:".format(partnum)) # Step through each vertex in the feature for pnt in part: if pnt: # print("{}, {}".format(pnt.X, pnt.Y)) CoordList.append("{};{}".format(pnt.X, pnt.Y)) else: # If pnt is None, this represents an interior ring # print("Interior Ring:") partnum += 1 # Find identical coordinates output = set([x for x in CoordList if CoordList.count(x) > 1]) print(output) spatial_reference = arcpy.SpatialReference(25832) pt = arcpy.Point() ptGeoms = [] for p in output: coord = p.split(";") pt.X = coord[0] print(coord[0]) pt.Y = coord[1] print(coord[1]) ptGeoms.append(arcpy.PointGeometry(pt, spatial_reference)) arcpy.CopyFeatures_management(ptGeoms, outfc) arcpy.AddXY_management(outfc)‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍ ... View more

I would like to divide a shp in about 20-50 adjacent polygons. I have the following script but the while loop (line 54) becomes endless if there are less than 10 adjacent polygons to select. Any ideas for a better solution? # -*- coding: cp1252 -*- import arcpy arcpy.env.overwriteOutput = True arcpy.env.workspace = "in_memory" Alk_Sel = r"R:\Daten\sgb_Planungen\2013\Saalbachniederung\Eigentuemer\ALK_merge_E.shp" alk_sel_Layer = "alk_sel_Layer" InFile = r"C:\temp2\InFile.shp" def SicFeature(): arcpy.CopyFeatures_management(Alk_Sel, InFile) arcpy.MakeFeatureLayer_management(InFile, alk_sel_Layer) def AddField(): arcpy.AddField_management(alk_sel_Layer, "Round", "LONG") def selectMin(): sum = arcpy.GetCount_management(alk_sel_Layer) print sum i = 0 for i in range(1, int(sum.getOutput(0)), 10): i += 1 Expression = """ "Round" = {} """.format(0) arcpy.MakeFeatureLayer_management(InFile, alk_sel_Layer, Expression ) sum1 = arcpy.GetCount_management(alk_sel_Layer) print('{} has a sum of {} records'.format(alk_sel_Layer, sum1[0])) MyList = list() cursor = arcpy.da.SearchCursor(alk_sel_Layer, "FID") for row in cursor: MyList.append(int(row[0])) print 'Minimum: {0}'.format(min(MyList)) expression1 = min(MyList) arcpy.SelectLayerByAttribute_management(alk_sel_Layer, "NEW_SELECTION", """ "FID" = {} """.format(expression1)) result1 = arcpy.GetCount_management(alk_sel_Layer) while (int(result1.getOutput(0)) < 10): arcpy.SelectLayerByLocation_management (alk_sel_Layer, "INTERSECT", alk_sel_Layer) result1 = arcpy.GetCount_management(alk_sel_Layer) print(result1) arcpy.CalculateField_management(alk_sel_Layer, "Round", i) OutFile = r"C:\temp2\OutFile_{}.shp".format(i) arcpy.CopyFeatures_management(alk_sel_Layer, OutFile) SicFeature() AddField() selectMin()‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍ ... View more

Hello, I have problems to update two attribute fields with values created in the following code. The insert cursor works well but with the update cursor I get a error message. Are there any other ways to fill a table with values? Error message: line 60, in UpdateFields     cursor.updateRow([row]) StopIteration: iteration not started # -*- coding: cp1252 -*- import arcpy, os arcpy.env.overwriteOutput = True myPath = r"R:\Daten\geo_daten\Arten\Fledermaeuse\AltDaten" out_path_table = r"C:\temp4\Geodatenbestand.mdb" out_name_table = "Geodatenbestand" template_table = r"C:\temp4\Geodatenbestand.mdb\Geodatenbestand_Template" inTable = os.path.join(out_path_table, out_name_table) inTableL = "Table_Layer" MyList1 = list() MyList2 = list() # Encoding der Standardausgabe herausfinden stdout_encoding = sys.stdout.encoding or sys.getfilesystemencoding() class Geodatenbestand(object): def __init__(self): self.myPath = myPath def SearchFiles(self): SHPCounter = 0 for path, dirs, files in os.walk(myPath): for file in files: if file.lower().endswith('.shp'): print path + "\\" + file SHPCounter += 1 print SHPCounter accessTime = os.path.getatime(os.path.join(path, file)) # last access time actualTime = time.time() hoursPassed = """{}""".format(int((actualTime - accessTime)/3600)) path1 = path.decode(stdout_encoding) file1 = file.decode(stdout_encoding) print file1 cursor = arcpy.da.InsertCursor(inTable, ['PathFile']) cursor.insertRow([path1 + "\\" + file1]) MyList1.append(file1) MyList2.append(hoursPassed) def UpdateFields(self): print MyList1 with arcpy.da.UpdateCursor(inTable,["File"]) as cursor: for row in MyList1: cursor.updateRow([row]) print MyList2 with arcpy.da.UpdateCursor(inTable,["h_not_used"]) as cursor: for row in MyList2: cursor.updateRow([row]) def CreateTable(self): if arcpy.Exists(inTable): arcpy.Delete_management(inTable) arcpy.CreateTable_management(out_path_table, out_name_table, template_table) if __name__ == '__main__': MyClass = Geodatenbestand() MyClass.CreateTable() MyClass.SearchFiles() MyClass.UpdateFields() ‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍ ... View more