nav[aria-label="Primary Navigation"] { padding: 0; & ul { list-style: none; width: 100%; display: flex; flex-direction: row; justify-content: start; align-items: start; gap: 30px; padding: 0; & li { margin: 0; } & ul li { list-style: none; } } }

Intersect with a specific feature and one unknown from gdb

DirkBrunken

Hello,

i want to find matching data from geodatabase. Thought at first, that arcpy.da.Walk is best way to get result, but maybe executing the tool Intersect is better.

Is it possible to set a feature and say, intersect with everything you find from gdb?

I am just a python beginner and want to shortly get results 😉

Find more posts tagged with

Accepted answers

JakeSkinner

Here is some code that you can work with. What the code does is iterate through a feature class and searches all features in every feature class within the geodatabase to see if any features intersect (see the disjoint method here). If features do, it adds the feature class name to a list. I would recommend creating a single polygon from the input feature class's extent. The code will then just have to iterate through one input feature rather than many.

import arcpy from arcpy import env env.workspace = r"C:\temp\python\test.gdb"  list = []  dataset = "Polygon" spatial_ref = arcpy.Describe(dataset).spatialReference  with arcpy.da.SearchCursor(dataset, ["SHAPE@"], "") as cursor:     for row in cursor:         for fc in arcpy.ListFeatureClasses("*"):             with arcpy.da.SearchCursor(fc, ["SHAPE@"], "", spatial_ref) as cursor2:                 for row2 in cursor2:                     if not row[0].disjoint(row2[0]):                         list.append(fc)  del cursor, cursor2  #remove duplicates from list list = dict.fromkeys(list) list = list.keys()  print list

All comments

JakeSkinner

import arcpy from arcpy import env env.workspace = r"C:\temp\python\test.gdb"  list = []  dataset = "Polygon" spatial_ref = arcpy.Describe(dataset).spatialReference  with arcpy.da.SearchCursor(dataset, ["SHAPE@"], "") as cursor:     for row in cursor:         for fc in arcpy.ListFeatureClasses("*"):             with arcpy.da.SearchCursor(fc, ["SHAPE@"], "", spatial_ref) as cursor2:                 for row2 in cursor2:                     if not row[0].disjoint(row2[0]):                         list.append(fc)  del cursor, cursor2  #remove duplicates from list list = dict.fromkeys(list) list = list.keys()  print list

DirkBrunken

and more specific asked..

how could i find data from a geodatabase that intersects with a single feature (e.g. a pointfeature).
Because, when i search for extent, its not that it is directly shown what really intersects or is ranked by percentage of intersection (what also would be nice to see and get)

DirkBrunken

Hi Jake,

thanks for your answer. That option is indeed really helpful. It's maybe indeed more effort to build a script that automatically finds overlapping features/feature_classes. But in the other hand, i still want to know how that task could be solved 😄

Greets, Dirk

JakeSkinner

Hi Dirk,

Does this have to be accomplished using python? It may be easiest to use Search in ArcMap. You can add data to ArcMap that you would like to find which data intersects with it and specify an extent option:

[ATTACH=CONFIG]32229[/ATTACH]

screen.png

Join the Conversation

Signed in members can post, follow updates, and more. New here? Register a free account.

Quick Links

Community Help
Find useful guides, FAQs, and documents to help you navigate and make the most of Esri Community.