mean angle

meanangle.py

#!/usr/bin/python2
import sys

def mean_angle(angles):
    edges = [a - angles[i] for i, a in enumerate(angles[1:])] + [360 - angles[-1] + angles[0]]
    maxedge = max(edges)
    start = (edges.index(maxedge) + 1)  % len(edges)
    angles2 = angles[:start] + [a - 360 for a in angles[start:]]
    mean = float(sum(angles2)) / len(angles)
    if mean < 0:
        mean += 360
    return mean


if __name__ == '__main__':
    if len(sys.argv) < 2:
        print("usage: meanangle ANGLES")
        raise SystemExit(2)
    angles = sorted(map(int, sys.argv[1:]))
    mean = mean_angle(angles)
    print("mean angle: {}".format(mean))