[GIS] Selecting features based on relation to other features within same layer using ArcGIS Desktop

arcgis-desktopselect

I want to select all the features in my Waterlines layer that share a line segment with other features in the same layer. I've done a complicated search cursor iteration in python that takes a couple hours, but I'd rather have something quicker.

I found this page that says you can do it simply with "select by location" which doesn't seem to work from the selection menu in arcmap, but kind of works in arcpy.

http://help.arcgis.com/en/arcgisdesktop/10.0/help/index.html#/Select_By_Location_Select_based_on_spatial_relationship_within_the_layer/0017000000tr000000/

When I ran this:

arcpy.SelectLayerByLocation_management("waterlines_layer", "SHARE_A_LINE_SEGMENT_WITH", "", "", "NEW_SELECTION")

…it selected things that definitely didn't share a line segment with anything else. But it didn't select everything (which would have at least made some sense).

Best Answer

The reason Select by Location is not able to do this in one fell swoop is that it considers every feature in the layer against every other feature including itself. So iteration is necessary using this tool and considering one feature at a time.

As you mentioned there are geometry on geometry checks in Data Reviewer that can do this but I'm not sure how much faster if any it would be.

Related Solutions

[GIS] Why unable to create interactive FeatureSet features for ModelBuilder tool

Open your model and right-click on Clip Features to check its Properties.

On the Parameters tab, you should see

Select Data Type set to Feature Set; and
Schema imported from set to "a feature class, a layer file (.lyr) that references a feature class, or a layer in the ArcMap table of contents" (from Setting feature set and record set schema)

I suspect that you currently have Schema imported from left blank.

[GIS] Using ArcPy to zoom to Y,X (string) coordinate pair copy/pasted from Google Maps

Below is some code that should do the basics of this.

import arcpy

yxCoordString = arcpy.GetParameterAsText(0)
yStr,xStr = yxCoordString.split(",")
xFloat = float(xStr)
yFloat = float(yStr)

mxd = arcpy.mapping.MapDocument("CURRENT")
df = arcpy.mapping.ListDataFrames(mxd)[0]
newExtent = df.extent
newExtent.XMin, newExtent.YMin = xFloat - 0.5, yFloat - 0.5
newExtent.XMax, newExtent.YMax = xFloat + 0.5, yFloat + 0.5
df.extent = newExtent

arcpy.RefreshActiveView()

You just need to Add Script it into a toolbox and then assign a parameter like below:

enter image description here

I have it zooming to a 1 x 1 degree extent by subtracting/adding 0.5 to the entered coordinates but you can set that to whatever you like.

This code takes no account of the coordinates entered being different to the coordinate system of the data frame.

Best Answer

Related Solutions

[GIS] Why unable to create interactive FeatureSet features for ModelBuilder tool

[GIS] Using ArcPy to zoom to Y,X (string) coordinate pair copy/pasted from Google Maps

Related Question