Haversine Distance Calculator

haversine.py

"""Calculate the difference between two GPS coordinates using the Haversine formula"""
from math import radians, cos, sin, sqrt, atan2


def haversine(coord1, coord2):
    # Radius of the Earth in kilometers
    R = 6378.137

    # Converting coordinates from degrees to radians
    lat1, lon1 = radians(coord1[0]), radians(coord1[1])
    lat2, lon2 = radians(coord2[0]), radians(coord2[1])

    # Differences in coordinates
    dlat = lat2 - lat1
    dlon = lon2 - lon1

    # Haversine formula
    a = sin(dlat / 2)**2 + cos(lat1) * cos(lat2) * sin(dlon / 2)**2
    c = 2 * atan2(sqrt(a), sqrt(1 - a))
    distance = R * c

    return distance


brian = 49.2602538, -123.2488872
me = 49.2643602, -123.2522186

print("Distance:", haversine(brian, me), "km")