pub fn bearing(a: LonLat, b: LonLat) -> f64
Initial bearing from a to b in degrees clockwise from north.
a
b