Spacial Join closest line with a certain angle. Is it possible?

Question

Hi,

Is it possible to spacial join a line feature set to another line feature set by proximity (closest) AND angle (ignoring all lines that have an angle/direction difference higher than X)?

From what I read about the Spacial Join Tool it seems it can only do the first, but wanted to check with those with more ArcGIS experience.

If it's indeed not possible, I 'm having doubts what the best solution would be:

1- Join each line with all that are within a certain distance (keeping the angle field of all) and then calculate which one has the most similar angle - but since I want to know the closest line that's in a certain angle interval in relation to the target line, I think this could give me the wrong results (especially as the distance radius increases)

2- Make and merge several Spacial Joins (closest) filtering the target feature ( angle=n ) and join feature ( angle=n+-x ) and just chose a x and go from n=0 to n=359? But I don't how to code this.

Cheers

JúlioLuta · Answer

Edit: It worked.

At first I thought it didn't but after I reloaded the file the AngleMin column had the correct values.

Thanks!

ClintonDow1 · Answer

Are you just missing the 'rows.updateRow(row)' in your for loop in the last snippet?

JúlioLuta · Answer

Hi, I wanted to share my progress in case it's useful to someone in the future. I also have a question because I can't understand why one step is not working. I calculated the angle of the two feature classes using the following code I found online and which uses the !shape! and the field calculator: # Pre-Loci Script Code import math def GetGeographicalDegrees(shape): radian = math.atan2(shape.lastpoint.y - shape.firstpoint.y, shape.lastpoint.x - shape.firstpoint.x) radian = radian - (math.pi /2 ) # turn minus 90° if (radian > 0): degrees = 360 - ( radian * 360) / ( 2 * math.pi ) else: degrees = 360 - ((2* math.pi + radian ) * 360) / ( 2 * math.pi ) return degrees # angle = GetGeographicalDegrees( !SHAPE! )‍‍‍‍‍‍‍‍‍‍‍‍‍‍ Now both line features have an angle field and after the spatial join (join one to many – within a distance of XXm) I created a new field that has the difference between angles (AngleDif): abs(!AngleArcs! - !AngleSegments!)‍ Next using the python console I used code I found as an answer to a simular problem, which needed to find for every field A with the same value the minimum value of field B: import arcpy fc = "C:\Dropbox\NewThesis\databases\Final.gdb\ArcosSolSombra_VGA2M_Segments" data = [row for row in arcpy.da.SearchCursor(fc, ["TARGET_FID", "AngleDif"])] import collections ddict = collections.defaultdict(list) for x,y in data: ddict[x].append(y) minvals = dict((key, min(ddict[key])) for key in ddict.keys()) with arcpy.da.UpdateCursor(fc, ["TARGET_FID", "AngleMin"]) as rows: for row in rows: row[1] = minvals[row[0]] ‍‍‍‍‍‍‍‍‍‍‍ The problem I had here is after getting the minvals list the updatecursor is not writing anything in the AngleMin field. To finish I'll create and calculate a field to identify when AngleMin is equal to AngleDif. I know this solution isn't very elegant but as I mentioned before I didn't practice enough after I learnt a bit of Python and now I barely remember anything.