topic Re: Creating polygons from closed polylines with arcpy in Python Questions

Creating polygons from closed polylines with arcpy

RafaeldePedro — Tue, 28 May 2019 11:57:10 GMT

Does anyone have any advice on how to invoke the "construct polygon" (I believe this lives in the Advance editing toolbar in Arcmap) on a python script?

I am currently trying to write a piece of code that will take in closed polylines which represent tree canopies, and convert them to polygons. I had tried using the tool Feature to Polygon in the Management toolbox, but it does not react well to closed polylines when they overlap each other as tree canopies do in real life. I am trying to get one polygon for each closed polyline, and the "construct polygon" tool seems to do the job, but I don't know how to call it on python.

If anyone here can suggest a way (with "construct polygon" or otherwise) to achieve this it would be very much appreciated.

Thanks.

Re: Creating polygons from closed polylines with arcpy

JoshuaBixby — Sat, 11 Dec 2021 10:06:34 GMT

I assume when you state, "they overlap each other," a given polyline does not overlap itself.

There is no direct way to run the toolbar command you want from Python. If you were so inclined, you could use comtypes and call the underlying ArcObjects, but I think the simpler path would be to use a Search Cursor to retrieve the polyline objects from the feature class and then convert them to polygon objects and insert them into a new feature class.

Once you have the polyline object, converting to a polygon is straightforward using ArcPy.

>>> import arcpy
>>> pl = arcpy.FromWKT('LINESTRING(0 0, 10 0, 10 10, 0 10, 0 0)', arcpy.SpatialReference(3857))
>>> pl
<Polyline object at 0x19159987dd8[0x19159c5a210]>
>>>
>>> pg = arcpy.Polygon(pl.getPart(), pl.spatialReference)
>>> pg
<Polygon object at 0x19162933908[0x19159ef0aa8]>
>>>‍‍‍‍‍‍‍‍‍

Re: Creating polygons from closed polylines with arcpy

LukeWebb — Sat, 11 Dec 2021 10:06:37 GMT

Joshuas answer seems the best way.

Another take, is to use the tool you were doing, but only run it on 1 feature at a time, something like:

src_shp = r'a_shapfile'
result_fc = 'a_different_shapefile'


for row in arcpy.SearchCursor(src_shp):
    row_oid = row.getValue('OBJECTID')
    row_where = """ObjectID = %s""" % row_oid
    active_lyr = arcpy.MakeFeatureLayer(src_shp, where_clause=row_where)
    if arcpy.Exists(row_result):
        arcpy.Delete(row_result)
    row_result = Feature to Polygon (active_lyr,  temp_layer)
    append_management(row_result, result_fc)
‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍

I imagine Joshuas approach above will perform much faster though!

Re: Creating polygons from closed polylines with arcpy

RafaeldePedro — Wed, 29 May 2019 10:30:24 GMT

Thank you both. Yes when I said 'overlap' I should have said 'cross', but since the polylines represent an actual area I used that word. I will experiment and see how far I get. Thanks again!

Re: Creating polygons from closed polylines with arcpy

ThomasBlankinship — Wed, 19 Jun 2019 15:09:41 GMT

Did this work? I'm looking for a way to do the same thing.

Re: Creating polygons from closed polylines with arcpy

RafaeldePedro — Wed, 26 Jun 2019 10:25:26 GMT

Hi Thomas,

The responses above set me on the right path to my solution. However my approach was in the end different. Here is a sample of the code that I used to "read" the original polyline's geometry, make it an array for each feature (in my case they were tree canopies) and use it to create a polygon feature in the destination feature class.

Sorry for the formatting of my code I simply copied and pasted and amended certain variable names. I hope it helps though

for row in arcpy.da.SearchCursor(infc, ["OID@", "SHAPE@"]):
vertexArray = arcpy.Array()

      for part in row[1]:
            for pnt in part:
                  vertexArray.add(pnt)

      with arcpy.da.InsertCursor(outfc, ("SHAPE@",)) as cursor:
            feature = arcpy.Polygon(vertexArray, spatialRef)
            cursor.insertRow((feature,))

Re: Creating polygons from closed polylines with arcpy

MatthewSackrison — Thu, 21 Nov 2024 00:18:30 GMT

Thank you! With the "Preserve attributes" functionality having been removed from Feature To Polygon (Data Management)—ArcGIS Pro | Documentation this is very helpful in creating a workaround.