by
by
by
XanderBakker
Apart from the indentation there are some funcional errors in your script. You only process the multipart polygons, while the final polygon is created from a arr_pol, which will be empty when the feature is not a multipart. You should process all features the same way regardless of it being multi or single part (all the features in your sample data have only 1 part). If you don't do that and you write the original polygons when it is single part, you will get this:
And yes rounding the coordinates on 2 decimales is visible in you data. I didn't your code, but the da cursor (ArcGIS 10.1 and higher) which is far more efficient.
Code used (to produce the "problem" above):
def main():
import arcpy
fc = r'D:\Xander\GeoNet\Donuts\shp\donutTesting.shp'
fc_out = r'D:\Xander\GeoNet\Donuts\shp\donutTesting_out.shp'
sr = arcpy.Describe(fc).spatialReference
lst_feats = []
with arcpy.da.UpdateCursor(fc, ('SHAPE@', 'OID@')) as curs:
for row in curs:
polygon = row[0]
oid = row[1]
print oid
if polygon.isMultipart:
print "multipart, # of parts:", polygon.partCount
lst_pol = []
for i in range(polygon.partCount):
part = polygon.getPart(i)
print" - part", i
lst_part = []
for pnt in part:
pnt2 = roundCoords(pnt, 2)
print " - pnt", pnt, pnt2
lst_part.append(pnt2)
lst_pol.append(lst_part)
polygon2 = arcpy.Polygon(arcpy.Array(lst_pol), sr)
lst_feats.append(polygon2)
else:
print "singlepart, # of parts:", polygon.partCount
lst_feats.append(polygon)
arcpy.CopyFeatures_management(lst_feats, fc_out)
def roundCoords(pnt, decimals):
if pnt:
return arcpy.Point(round(pnt.X, decimals), round(pnt.Y, decimals))
else:
return pnt
if __name__ == '__main__':
main()
Output messages:
0
singlepart, # of parts: 1
1
multipart, # of parts: 1
- part 0
- pnt 48371,4859999996 42834,2699999996 NaN NaN 48371,49 42834,27 NaN NaN
- pnt 48372,1289999997 42966,5089999996 NaN NaN 48372,13 42966,51 NaN NaN
- pnt 48371,9040000001 42977,0130000003 NaN NaN 48371,9 42977,01 NaN NaN
- pnt 48400,7130000005 42992,7489999998 NaN NaN 48400,71 42992,75 NaN NaN
- pnt 48620,4069999997 42952,8029999994 NaN NaN 48620,41 42952,8 NaN NaN
- pnt 48627,1799999997 42891,8430000003 NaN NaN 48627,18 42891,84 NaN NaN
- pnt 48481,5530000003 42841,0429999996 NaN NaN 48481,55 42841,04 NaN NaN
- pnt 48371,4859999996 42834,2699999996 NaN NaN 48371,49 42834,27 NaN NaN
- pnt None None
- pnt 48435,8629999999 42969,7170000002 NaN NaN 48435,86 42969,72 NaN NaN
- pnt 48466,034 42961,6860000007 NaN NaN 48466,03 42961,69 NaN NaN
- pnt 48468,5360000003 42963,7890000008 NaN NaN 48468,54 42963,79 NaN NaN
- pnt 48470,3729999997 42965,7080000006 NaN NaN 48470,37 42965,71 NaN NaN
- pnt 48435,8629999999 42969,7170000002 NaN NaN 48435,86 42969,72 NaN NaN
- pnt None None
- pnt 48391,9639999997 42944,5050000008 NaN NaN 48391,96 42944,51 NaN NaN
- pnt 48408,0039999997 42942,8200000003 NaN NaN 48408 42942,82 NaN NaN
- pnt 48416,9890000001 42950,9759999998 NaN NaN 48416,99 42950,98 NaN NaN
- pnt 48391,4349999996 42957,2050000001 NaN NaN 48391,43 42957,21 NaN NaN
- pnt 48391,9639999997 42944,5050000008 NaN NaN 48391,96 42944,51 NaN NaN
2
singlepart, # of parts: 1Doing this (the right way) in ArcGIS 10.0, you could use this code:
def main():
import arcpy
fc = r'D:\Xander\GeoNet\Donuts\shp\donutTesting.shp'
fc_out = r'D:\Xander\GeoNet\Donuts\shp\donutTesting_out2.shp'
sr = arcpy.Describe(fc).spatialReference
fld_shp = arcpy.Describe(fc).ShapeFieldName
fld_oid = arcpy.Describe(fc).OIDFieldName
curs = arcpy.UpdateCursor(fc)
for row in curs:
polygon = row.getValue(fld_shp)
oid = row.getValue(fld_oid)
print "oid", oid
lst_pol = []
for i in range(polygon.partCount):
part = polygon.getPart(i)
print" - part", i
lst_part = []
for pnt in part:
pnt2 = roundCoords(pnt, 2)
print " - pnt", pnt, pnt2
lst_part.append(pnt2)
lst_pol.append(lst_part)
polygon2 = arcpy.Polygon(arcpy.Array(lst_pol), sr)
row.setValue(fld_shp, polygon2)
curs.updateRow(row)
del curs, row
def roundCoords(pnt, decimals):
if pnt:
return arcpy.Point(round(pnt.X, decimals), round(pnt.Y, decimals))
else:
return pnt
if __name__ == '__main__':
main()Zoom in to coordinates 48391.435 42957.205 Meters to see the effect.
by
by
