Add "rad" and "deg" to Units about haversine

[jdeniau]

Hi, my math is far in my memory, but if rad and deg are angle units, so how do you convert a distance to an angle ? 🤔

[jdeniau]

To rephrase what I said : the input might be in rad or deg, but the output is actually a distance,

If I look at the scikit learn implementation of haversine that you referenced in #30 their output is also in kilometers:

scikitlearn:

>>> from sklearn.metrics.pairwise import haversine_distances
>>> from math import radians
>>> bsas = [-34.83333, -58.5166646]
>>> paris = [49.0083899664, 2.53844117956]
>>> bsas_in_radians = [radians(_) for _ in bsas]
>>> paris_in_radians = [radians(_) for _ in paris]
>>> result = haversine_distances([bsas_in_radians, paris_in_radians])
>>> result * 6371000/1000  # multiply by Earth radius to get kilometers
array([[    0.        , 11099.54035582],
       [11099.54035582,     0.        ]])

haversine:

>>> bsas = (-34.83333, -58.5166646)
>>> paris = (49.0083899664, 2.53844117956)
>>> haversine(bsas, paris)
11099.555687157774

[MuellerSeb]

An angle in radians is simply the distance on the unit-sphere. So to get the result in rad, you simply don't multiply with the earth radius. The dependence between angle and distance in linear (that is why the circumference of a circle is r * 2pi (where 2pi is a 360 deg angle))

Edit:
and as you see in the sklearn example the result of the haversine function is multiplied with the earth radius. The result itself is in radians.

[jdeniau]

Fixed by #40