[python] degrees to cardinal directions

deg_to_cardinal.py

def degrees_to_cardinal(d):
    '''
    note: this is highly approximate...
    '''
    dirs = ["N", "NNE", "NE", "ENE", "E", "ESE", "SE", "SSE",
            "S", "SSW", "SW", "WSW", "W", "WNW", "NW", "NNW"]
    ix = int((d + 11.25)/22.5)
    return dirs[ix % 16]

Used in goose-search (here). Thank you!

I would say that @umiphos is correct, however you need to round the result instead of casting it to an integer, as it will always be rounded down.

I.e. the code should be:

dirs = ['N', 'NNE', 'NE', 'ENE', 'E', 'ESE', 'SE', 'SSE', 'S', 'SSW', 'SW', 'WSW', 'W', 'WNW', 'NW', 'NNW']
ix = round(d / (360. / len(dirs)))
return dirs[ix % len(dirs)]

@RobertSudwarts and @Lauszus,
Used in my project, Thank You!

@RobertSudwarts and @Lauszus,
Thank you! Used for my weather tracking display project!

Thanks, @RobertSudwarts, and @Lauszus. Used for my project!

Haven't investigated too carefully but I needed ix to be an int ie

def calculate_bearing(d):
  dirs = ['N', 'NNE', 'NE', 'ENE', 'E', 'ESE', 'SE', 'SSE', 'S', 'SSW', 'SW', 'WSW', 'W', 'WNW', 'NW', 'NNW']
  ix = int(round(d / (360. / len(dirs))))  <--- int here
  return dirs[ix % len(dirs)]

Thanks

@mattarderne this is because the result of round in Python 2 is a float while the returned value is an int in Python 3.

See: https://docs.python.org/2/library/functions.html#round and https://docs.python.org/3/library/functions.html#round.

Python 2:

If ndigits is omitted, it defaults to zero. The result is a floating point number.

Python 3:

If ndigits is omitted or is None, it returns the nearest integer to its input.

Thank you, I am going to use this is my project

Thank you. Will use this in my project.

In case someone want's to apply for a Numpy Array of angles I've adapted as such...

def degrees_to_cardinal(d):
    '''
    note: this is highly approximate...
    '''
    dirs = np.array(["N", "NNE", "NE", "ENE", "E", "ESE", "SE", "SSE",
                     "S", "SSW", "SW", "WSW", "W", "WNW", "NW", "NNW"], dtype='U')
    ix = np.round(d / (360. / len(dirs))).astype('i')
    return dirs[ix % 16]

I did the following:

DIRECTIONS = [
    'N', 'NNE', 'NE', 'ENE', 'E', 'ESE', 'SE', 'SSE', 'S', 'SSW', 'SW', 'WSW',
    'W', 'WNW', 'NW', 'NNW', 'N'
]

x = 155 # anything 0 - 360
cardinal = DIRECTIONS[round(x / 22.5)]

Note the extra N at the end of the list.

This should probably be made a package