I am trying to get rid of overlapping, incongruous lines representing the same road. I have done this successfully on an exported subset of my feature class (13 rows) with this process:
- Buffer FC by 15 feet
- Join FC to buffer layer via Spatial Join - "One_To_One", "Completely_Contains"
- Select features from Spatial Join with "Join_Count" > 1
- Select FC by location - has centroid in selected Spatial Join features
- Integrate selected features From there, I can use QGIS to split lines at intersections and delete duplicates
However, when I try this process on my full streets features class (~50k rows), I get different results.
As you can see, the buffer of my subset FC (top) looks identical to my full FC (bottom).
But when I select "Join_Count > 1" the road segment shown in the first set of pictures has inconsistent results (subset top, full bottom):
Why am I seeing this discrepancy and how can I fix it?
I am using ArcMap 10.2, ArcGIS Pro, and QGIS. I have tried copying the python snippet from the successful run and tweaking it to run my full FC, but I get the same results.