topic Re: Iterate through polygons and perform Intersect in Python Questions

Iterate through polygons and perform Intersect

KathyGillis — Thu, 25 Jun 2020 23:12:42 GMT

Hi, I'm hopping for some help as I work through the process of writing a piece of code to produce an output of polygon intersects. My shape file contains many polygons which I would like to determine if any of them intersect and if they do add that intersect results to a new shape file. The code I used is below. I'm getting and error which is shown in a pic below.

Ultimately I'd like to have a single output shape file containing all intersect results. If you know a more efficient way or know how to make what I have work I am all ears!

I have attached a sample of the polygons. Many thanks for your attention!

I'm limited to using ArcMap 10.7 and python 2.7

```

import arcpy
from arcpy import env
import os

##Set working directory and workspace
wk_dir = os.getcwd()
env.workspace = wk_dir

##Assign feature class and working fields
fc = wk_dir + '\\CBT_Composite.gdb\\example_polys'
fields = ['RecordID', 'SHAPE@']

##Assign output file path and name
out = wk_dir + '\\output.gdb\\A'
##print fc_path + fc
with arcpy.da.SearchCursor(fc, fields) as cursor:
    for row in cursor:
        record = '{0}'.format(row[0]).split('.')[0]
        polygon = row[1]
        for row2 in cursor:
            if row2 != row:
                arcpy.Intersect_arc(row2[1], polygon, out +record)

```

Re: Iterate through polygons and perform Intersect

MehdiPira1 — Sun, 12 Dec 2021 08:06:40 GMT

Hi Kathy Gillis‌,

If you change the variables and run the following code you should get the intersected shapes, all in one feature class:

arcpy.env.workspace = r"workspace path"
arcpy.analysis.Intersect("example_polys #", "example_polys_Intersect", "ALL", None, "INPUT")

‍‍‍‍‍‍‍‍‍

Re: Iterate through polygons and perform Intersect

KathyGillis — Fri, 26 Jun 2020 17:48:58 GMT

Thank you for your response Mehdi Pira!

Can you clarify what ' #' stands for in the Intersect statement?

Thank you Again - Kathy

Re: Iterate through polygons and perform Intersect

MehdiPira1 — Fri, 26 Jun 2020 18:13:23 GMT

Not a problem, Kathy Gillis‌

This sign # is for Priority Ranks. Refer to the link below for more info:

Priority ranks and geoprocessing tools—Help | Documentation

Re: Iterate through polygons and perform Intersect

KathyGillis — Fri, 26 Jun 2020 21:01:16 GMT

Thank you! @mehdi Pira It worked perfectly. Smiling really big!