Geospatial Modelling Environment

listintersectingfeatures (List Intersecting Features)

For each feature in an input layer this tool creates a list of other features (in the same layer or a different layer) that intersect it


This tool creates a delimited text file that summarizes intersections among features. For each feature (point, line or polygon) in a dataset the tool identifies which features from the same or another dataset intersect it, and writes the unique ID numbers to the output file.

By default the output file contains one line per feature of the input data source (not including features with no intersections at all), followed by the count of intersecting features, and then a list of the ID numbers. The alternative ‘linear’ format consists of one pair of intersecting features per line (only the two ID numbers are shown).

The delimiting character can be specified. See the ‘delimiter’ command for instructions on how to change the delimiting character on all tabular text output (this is a separate command that you would issue before this command, not a parameter for this command).

Although the default spatial relationship is intersect, you can actually specify the type of relationship to use. The options (and their abbreviations in parentheses) are: INTERSECTS (I), CROSSES (CR), OVERLAPS (O), TOUCHES (T), WITHIN (W), CONTAINS (CO), ENVELOPEINTERSECTS (EI), INDEXINTERSECTS (II). You may also specify a custom relationship using: RELATExxxxxxxxx, where the x’s define the custom relationship and can assume values of F, T, or * (e.g. RELATEFFFTTT***). Refer to the ESRI Help documentation for descriptions of these relationships.

Intersections can be computationally expensive, therefore it is recommended that you consider taking steps to increase processing speed, especially for large datasets. Ensuring that your feature dataset has a spatial index can improve search speed (see the Add Spatial Index tool in ArcToolbox). You might also consider eliminating any features that are known to have no overlap with the features in another layer. For instance, the ArcMap Select By Location tool allows you to easily select overlapping features between two layers, which you can then export into two new (temporary) layers that are the input to this tool. Reducing the number of features in the layers will thus also improve search speeds.


listintersectingfeatures(in, uidfield, in2, uidfield2, out, [relationship], [linear], [where]);

inthe input feature data source
uidfieldthe unique ID field of the first vector data source
in2the second vector data source
uidfield2the unique ID field of the second vector data source
outthe full path to the output delimited text file
[relationship]the spatial relationship to detect features (default=INTERSECTS; options: INTERSECTS, CROSSES, OVERLAPS, TOUCHES, WITHIN, CONTAINS, ENVELOPEINTERSECTS, INDEXINTERSECTS)
[linear](TRUE/FALSE) if true, one pair of ID numbers is written per line (default=FALSE)
[where]the filter/selection statement that will be applied to the first (‘in’, not ‘in2’) feature class to identify a subset of features to process


listintersectingfeatures(in=”C:dataplots.shp”, uidfield=”PLOTID”, in2=”C:datatransects.shp”, uidfield2=”TID”, out=”C:datatransinplots.csv”, linear=TRUE);

listintersectingfeatures(in=”C:datazones.shp”, uidfield=”ZNID”, in2=”C:dataroads.shp”, uidfield2=”RDID”, out=”C:dataisectfeat.csv”);

listintersectingfeatures(in=”C:dataparcels.shp”, uidfield=”PID”, in2=”C:dataschoolzones.shp”, uidfield2=”SCHID”, out=”C:dataisectfeat.csv”, relationship=”OVERLAPS”, linear=TRUE);


Please consider making a purchase to support the continued development of these tools  Read more...

Tips on how to use this interface efficiently 


Open Source GIS

Copyright © 2001-2014 Hawthorne L. Beyer, Ph.D., Spatial Ecology LLC    Connect on LinkedIn