topic Re: How to calculate the angle of a line respect to the North in Geoprocessing Questions

How to calculate the angle of a line respect to the North

CeciliaMartinelli — Tue, 13 Sep 2022 13:57:55 GMT

Hi,

I would like to calculate the angle of the lines in the image respect to the North. Is there a tool that can do that? and eventually put the output in the Attributes table. I know that these lines aren't perfectly straight but I need the average

I'm using ArcGis Pro 3.0.0

Thanks!

Re: How to calculate the angle of a line respect to the North

JohannesLindner — Tue, 13 Sep 2022 15:06:41 GMT

~~I don't think there is a tool that does that out of the box~~ (edit: nevermind, see above). But you can easily calculate the value with Calculate Field.

angle of the lines in the image respect to the North

Do you mean the geographic angle (North=0°, East=90°)? If so, use this Arcade expression in the tool:

// get a list of line parts var paths = Geometry($feature).paths // for each part, loop through each segment and get its angle (arithmetic angle -> East=0°, North=90°) var angles = [] for(var p in paths) { var count_p = Count(paths[p]) for(var q = 1; q < count_p; q ++) { var a = Angle(paths[p][q], paths[p][q-1]) Push(angles, a) } } // get the mean arithmetic angle var arith_angle = Mean(angles) // convert to geographic and return return (450 - arith_angle) % 360

Or do you mean the arithmetic angle (rotating counter-clockwise) but with North = 0°? Then use this expression:

// get a list of line parts var paths = Geometry($feature).paths // for each part, loop through each segment and get its angle (arithmetic angle -> East=0°, North=90°) var angles = [] for(var p in paths) { var count_p = Count(paths[p]) for(var q = 1; q < count_p; q ++) { var a = Angle(paths[p][q], paths[p][q-1]) Push(angles, a) } } // get the mean arithmetic angle var arith_angle = Mean(angles) // rotate so that North = 0° var north_arith_angle = arith_angle - 90 return north_arith_angle + IIF(north_arith_angle < 0, 180, 0)

Or do you mean something else? Then please clarify.

Re: How to calculate the angle of a line respect to the North

LukeCatania — Tue, 13 Sep 2022 15:32:19 GMT

I deleted the other post as it was not quite correct, but if you extract your start and end points (https://pro.arcgis.com/en/pro-app/latest/tool-reference/data-management/feature-vertices-to-points.htm) and then use https://pro.arcgis.com/en/pro-app/latest/tool-reference/defense/coordinate-table-to-line-of-bearing.htm

You should be able to get what you need

Re: How to calculate the angle of a line respect to the North

CeciliaMartinelli — Wed, 14 Sep 2022 10:01:14 GMT

Thank you! I'm indeed calculating the geographic angle and your solution works well, but just in some cases it is odd. For example, for the lines in the image below the dashed line's angle results to be 177° and the one for the solid line is 354°, so they are both wrong... could it be due to the sense in which I drew the lines?

They should be : dashed line around 100° and solid line around 75°

Thanks

Re: How to calculate the angle of a line respect to the North

JohannesLindner — Wed, 14 Sep 2022 13:46:12 GMT

Indeed, the angle is dependent on the vertex order:

Removing that dependency is easy, but it wasn't quite as trivial as I thought. I thought I just had to replace the modulo 360 with modulo 180 (last line). That works for lines with a single segment, but polylines are still wrong:

Instead, you have to convert to geographic angle inside the loop, not at the end:

// get a list of line parts var paths = Geometry($feature).paths // for each part, loop through each segment and get its geographic angle var angles = [] for(var p in paths) { var count_p = Count(paths[p]) for(var q = 1; q < count_p; q ++) { var a_arith = Angle(paths[p][q], paths[p][q-1]) var a_geo = (450 - a_arith) % 180 // 180 instead of 360, so that vertex order doesn't matter Push(angles, a_geo) } } // return the mean return Mean(angles)