Average bearing: one wrong way and two correct ones

circular-mean.py

import numpy as np
import scipy.stats

# compass bearings in degrees:
data = [2, 356]

# average bearing:
print(np.mean(data)) # == 179.0, aka "wrong"
print(scipy.stats.circmean(data, high=360, low=0)) # == 359.0, aka "correct"

print(np.rad2deg(np.angle(np.mean(np.cos(np.deg2rad(data)) + np.sin(np.deg2rad(data)) * 1j))%(2*math.pi)))
# == 359.0, aka also correct, but achieving the result in an overly complex way