topic Re: Python script for convert adjacent polygons to multipart polygons without dissolving in Python Questions

Python script for convert adjacent polygons to multipart polygons without dissolving

AmrHassan — Tue, 24 Feb 2015 10:16:25 GMT

I have a polygon feature class contains groups of adjacent polygons and the question is how to create a multipart polygon for each group of the adjacent polygons without dissolving it using python scripts.

thanks

Re: Python script for convert adjacent polygons to multipart polygons without dissolving

DanPatterson_Retired — Tue, 24 Feb 2015 14:09:48 GMT

Is this a task that you are going to repeat multiple times?

Is it part of a workflow for something much larger?

Is it safe to assume that this is just a sample of a huge dataset making the manual method of doing this inappropriate?

Do you have anything within your attribute table that indicates adjacency?

Just need to know how far along you are in your script development. If you have some code you are working with, posting would be useful since people won't have to reinvent the wheel.

Re: Python script for convert adjacent polygons to multipart polygons without dissolving

XanderBakker — Tue, 24 Feb 2015 14:27:10 GMT

Might be best to dissolve the polygons not creating multiparts (yep the opposite of what you want) to a tmp featureclass and then loop through these features, select the original polygons, and add the polygons to parts of the new polygon to be created.

Re: Python script for convert adjacent polygons to multipart polygons without dissolving

XanderBakker — Tue, 24 Feb 2015 14:33:38 GMT

You may also want to consider moving this thread to the Python place:

(see: What are Places? )

Re: Python script for convert adjacent polygons to multipart polygons without dissolving

AmrHassan — Tue, 24 Feb 2015 21:15:37 GMT

This task only for time and i do it rarely,

Not a in a workflow,

The manual method of doing this inappropriate as it dissolve them and i want each group of the adjacent polygons to be one feature one record as a multipart polygon but not dissolved,

There are common attributes between the adjacent polygons.

I didn't used a script so far, I am asking about it

Re: Python script for convert adjacent polygons to multipart polygons without dissolving

XanderBakker — Tue, 24 Feb 2015 21:25:45 GMT

... and the common attributes need to be summed in some way or can you just to the attributes of the first part and assign it to the output multipart feature?

Re: Python script for convert adjacent polygons to multipart polygons without dissolving

AmrHassan — Tue, 24 Feb 2015 21:31:28 GMT

All attributes should be assigned to the output multipart feature.

Re: Python script for convert adjacent polygons to multipart polygons without dissolving

XanderBakker — Tue, 24 Feb 2015 21:35:34 GMT

... the question remains: you are combining multi features (multi records) into a single record. Is it valid that a single (random) record is used to assign to the output multipart feature?

Re: Python script for convert adjacent polygons to multipart polygons without dissolving

XanderBakker — Tue, 24 Feb 2015 21:36:14 GMT

maybe it would be easier to attach the features you used in the picture. This will enable us to see what problems you may run into ...

Re: Python script for convert adjacent polygons to multipart polygons without dissolving

XanderBakker — Sun, 12 Dec 2021 00:24:22 GMT

I hacked some code to create the multipart as described before, but I noticed a strange effect. I started a new thread to confirm this undesired behavior: Re: Create a multipart polygon with arcpy and it appears that what you are looking for is not possible as Darren Wiens‌ indicated in the other thread.

...any parts that share an edge will be merged into a single part (source: ArcGIS Help 10.1)

Below the python code that I had so far, but this does not give you the result you want, since it is not possible.

import arcpy, os
arcpy.env.overwriteOutput = True

fc_in = r"D:\Xander\GeoNet\HouseMultiPart\data.gdb\houses"
fc_dis = r"D:\Xander\GeoNet\HouseMultiPart\data.gdb\housesDissolved3"
fc_out = r"D:\Xander\GeoNet\HouseMultiPart\data.gdb\houses_MPs13"

lst_flds = arcpy.ListFields(fc_in)

ws, name = os.path.split(fc_out)
arcpy.CreateFeatureclass_management(ws, name, "POLYGON", fc_in, spatial_reference=fc_in)

##for fld in lst_flds:
##    if len(arcpy.ListFields(fc_out, wild_card=fld.name)) == 0:
##        arcpy.AddField_management(fc_out, )
arcpy.MakeFeatureLayer_management(fc_in, "lyr_in")

# loop through dissolved features
flds = ("SHAPE@")
with arcpy.da.InsertCursor(fc_out, flds) as curs_out:
    with arcpy.da.SearchCursor(fc_dis, flds) as curs_dis:
        for row_dis in curs_dis:
            pol_dis = row_dis[0]

            # select features that overlap with dissolved feature
            arcpy.SelectLayerByLocation_management("lyr_in", "INTERSECT", pol_dis)
            with arcpy.da.SearchCursor("lyr_in", ("SHAPE@", "OID@")) as curs_in:
                lstParts = []
                for row_in in curs_in:
                    pol_in = row_in[0]
                    oid = row_in[1]
                    for part in pol_in:
                        lstPnt = []
                        for pnt in part:
                            lstPnt.append([pnt.X, pnt.Y])
                        # construct multipart
                        lstParts.append(lstPnt)

            # write the multipart to output fc
            pol_out = arcpy.Polygon(arcpy.Array(lstParts))
            curs_out.insertRow((lstParts[0], ))

del curs_dis, curs_in, curs_out, row_dis, row_in