topic Re: Get line direction/orientation as a numeric field in Spatial Data Science Questions

Get line direction/orientation as a numeric field

AlanAitken — Mon, 02 May 2011 03:48:53 GMT

Should be a simple one here.

How do I get ArcGIS 10 to return the orientation (or direction) of each line in a polyline feature into a field.

what do I need to type into field calculator (or is there a better way)

Ta
Alan

Re: Get line direction/orientation as a numeric field

FrédéricPRALLY — Sat, 11 Dec 2021 23:10:15 GMT

Hi,

To calculate angle of lines fill calculate values with parser python using this python line:

180 + math.atan2((!Shape.lastpoint.Y! - !Shape.firstpoint.Y!),(!Shape.lastpoint.X! - !Shape.firstpoint.X!)) * (180 / math.pi)

Hope this help you,

Best regards,

Fred

Re: Get line direction/orientation as a numeric field

AlanAitken — Tue, 03 May 2011 00:38:43 GMT

Thanks Fred,

Trig seems obvious now, but good to get the python syntax

Works

Alan

Re: Get line direction/orientation as a numeric field

SteveFeeney — Sun, 19 Jun 2011 09:37:03 GMT

Is there any means to attribute a field of a Polyline with the line bearing within a model? I have a model I need to package and distribute to other operators who are still using 9.3. Ideally it needs as few add-ons as possible to function.

Thank you,
Steve

Re: Get line direction/orientation as a numeric field

SteveFeeney — Thu, 23 Jun 2011 13:36:24 GMT

Please disregard my last post.
Some revision and self education later and I can see it has already been answered above.

When I input the python code, I'm using the "Calculate Field" management tool with expression type as "Python 9.3". The function calculates fine but returns different values to the ones I manually measured using MA. Am I using it correctly?
Kind regards,
Steve

Re: Get line direction/orientation as a numeric field

AlexanderBeavis — Tue, 23 Aug 2011 05:11:29 GMT

So does this code provide the bearing of a polyline with reference to north?

Re: Get line direction/orientation as a numeric field

chriss_ — Wed, 07 Sep 2011 02:30:09 GMT

Hi!
I wanted to use the Code to define the deviation from north in degree. Where north would be 0 east 90 south 180 and west 270. So I changed the code to this:

180 + math.atan2((!Shape.lastpoint.X! - !Shape.firstpoint.X!),(!Shape.lastpoint.Y! - !Shape.firstpoint.Y!)) * (180 / math.pi)

So far so good, BUT for everything in western direction I get a negative value eg. -90 instead of 270.

Any mathematical genius that can enlighten me?
Thanks
Chris

Re: Get line direction/orientation as a numeric field

BradyCallahan — Sat, 11 Dec 2021 23:10:18 GMT

The radians thing didn't help me either so I used this to return a North Azimuth based reference. Hope it helps someone...

Field (type double) = NorthAzimuth( !Shape! )

Codeblock:

def NorthAzimuth(Pline):
  degBearing = math.degrees(math.atan2((Pline.lastPoint.X - Pline.firstPoint.X),(Pline.lastPoint.Y - Pline.firstPoint.Y)))
  if (degBearing < 0):
      degBearing += 360.0
  return degBearing

Re: Get line direction/orientation as a numeric field

SteveFeeney — Wed, 12 Oct 2011 18:31:29 GMT

This may be a bit basic,

There is a tool called 'Linear Directional Mean' that will create and populate a compass bearing field with the mean direction of a shapefile/featureclass or, if you use a unique field, for each feature.

As far as I can tell it is planar rather than geodesic but it's a sound, quick fix and saves coding trig.

Re: Get line direction/orientation as a numeric field

MarcoRusmini — Fri, 01 Jun 2012 12:33:02 GMT

I'm not a phyton programmer but I run into the same problem: finding the deviation from North. I solved it (or it seems) by swopping firstpoint and lastpoint in the following way:

180+math.atan2(( !Shape.firstpoint.X! - !Shape.lastpoint.X! ),( !Shape.firstpoint.Y! - !Shape.lastpoint.Y! ) ) * (180 / math.pi )

It seems to me that it returns exactly the deviation from North!!

Please check this out

Cheers

Marco Rusmini

Hi!
I wanted to use the Code to define the deviation from north in degree. Where north would be 0 east 90 south 180 and west 270. So I changed the code to this:

180 + math.atan2((!Shape.lastpoint.X! - !Shape.firstpoint.X!),(!Shape.lastpoint.Y! - !Shape.firstpoint.Y!)) * (180 / math.pi)

So far so good, BUT for everything in western direction I get a negative value eg. -90 instead of 270.

Any mathematical genius that can enlighten me?
Thanks
Chris

Re: Get line direction/orientation as a numeric field

RichardFairhurst — Wed, 09 Mar 2016 19:36:01 GMT

The formula by Brady is the only correct formula proposed above.

The basic formula to calculate degrees from two points is only correct in Brady's formula:

math.degrees(math.atan2((!Shape.lastPoint.X! - !Shape.firstPoint.X!),(!Shape.lastPoint.Y! - !Shape.firstPoint.Y!)))

math.atan2((!Shape.lastpoint.X! - !Shape.firstpoint.X!),(!Shape.lastpoint.Y! - !Shape.firstpoint.Y!)) * (180 / math.pi)

This formula starts with 0 degrees pointing north and produces positive angles in the positive X axis (east half) and negative angles in the negative x axis (west half) of the quadrants: To make all angles positive so that it is clearer that they are all increasing in a clockwise direction, the formula needs to only add 360 to all of the negative angles in the West half of the quadrants. Brady's codeblock does this and using the codeblock this way is the only efficient way to do this calculation.

Although the expression below runs much slower than Brady's codeblock, it is possible to do this calculation correctly without using a codeblock by using a single in-line if python expression:

360 + math.degrees(math.atan2((!Shape.lastpoint.X! - !Shape.firstpoint.X!),(!Shape.lastpoint.Y! - !Shape.firstpoint.Y!))) if math.degrees(math.atan2((!Shape.lastpoint.X! - !Shape.firstpoint.X!),(!Shape.lastpoint.Y! - !Shape.firstpoint.Y!))) < 0 else math.degrees(math.atan2((!Shape.lastpoint.X! - !Shape.firstpoint.X!),(!Shape.lastpoint.Y! - !Shape.firstpoint.Y!)))

This code is slower than Brady's code since it has to extract the end points from the line and do the degree calculation multiple times, while his calculation only does those 2 steps once.

By adding 180 degrees to the basic formula, Criss has made 0 and 360 degrees represent a line that is pointed South and messed up the quadrants of the angles. Marco's formula defaults to a South (180 degree) line orientation, but by swapping the points he has reversed the orientation of the line relative to its real orientation, which I believe restores the real line orientation to show north as 0 or 360 degrees and fixes the messed up quadrants of the angles. I think Marco's formula may actually work except when it is used with a closed line (both ends are the same point) which will be reported as 180 degrees and not as 0 degrees. However, closed lines are a special case, and none of the formulas make it possible to distinguish that situation from a line that is actually oriented due north (for Brady's formula) or due south (for Marco's formula).

The Linear Directional Mean tool does produce angles that are compass oriented, and takes into consideration more than just the two end points of the line. So it provides a fuller analysis of the line direction for complex polylines and populates additional fields that make it possible to identify closed polylines and polylines where the angle between the end points does not really correspond to the overall directional trend of the polyline. The trig formula is really only accurate for two point lines that fall within the bounds that a Projected Coordinate System has been optimized to cover. It does not really work very well for polylines with many vertices that have severe changes in directions between the two end points or lines that extend beyond the appropriate bounds of a Projected Coordinate System.

Re: Get line direction/orientation as a numeric field

TedKowal — Sat, 11 Dec 2021 23:10:21 GMT

I use Radians for determining Azimuth's and convert it back to 360 degrees.... here is a VB snippet of code I use from an interactive addIn I created. Me.lng/Me.Lat is the user's selected start point of the line destination is the lat/lng of the endpoint. Some Python guru can easily convert this to Python...

Public Function GetAzimuth(destination As LatLng) As Double
  Dim longitudinalDifference = destination.Lng - Me.Lng
  Dim latitudinalDifference = destination.Lat - Me.Lat
  Dim azimuth = (Math.PI * 0.5) - Math.Atan(latitudinalDifference / longitudinalDifference)
  If longitudinalDifference > 0 Then
   Return azimuth
  ElseIf longitudinalDifference < 0 Then
   Return azimuth + Math.PI
  ElseIf latitudinalDifference < 0 Then
   Return Math.PI
  End If
  Return 0.0
End Function


Public Function GetDegreesAzimuth(destination As LatLng) As Double
  Return RadiansToDegreesConversionFactor * GetAzimuth(destination)
End Function

Re: Get line direction/orientation as a numeric field

DanPatterson_Retired — Thu, 10 Mar 2016 15:56:20 GMT

this thread shows the basics of a numpy/python implementation... getting the start and endpoints can be done with FeatureclassToNumPyArray, then distances and angles from a source to a destination(s) is quick

https://community.esri.com/thread/66222#comment-570366

Re: Get line direction/orientation as a numeric field

DiegoNegron — Mon, 06 Jun 2016 20:02:27 GMT

I do want to note that if you have a logitudinalDifference that equals 0 you will generate a zero division error with this code.

Re: Get line direction/orientation as a numeric field

Anonymous User — Thu, 05 Sep 2019 13:16:31 GMT

Thank you! it worked !

Re: Get line direction/orientation as a numeric field

MarkusGeier — Sun, 26 Jan 2020 14:11:58 GMT

You, sir, are my personal hero of the day. Thank you so much for your help.