lat_long_distance.py

def distance(origin, destination):
    '''Haversine distance function to compute distance between 2 lat/long pairs. source: https://gist.github.com/rochacbruno/2883505'''
    lat1, lon1 = origin
    lat2, lon2 = destination
    radius = 3959  # miles

    dlat = math.radians(lat2 - lat1)
    dlon = math.radians(lon2 - lon1)
    a = math.sin(dlat / 2) * math.sin(dlat / 2) + math.cos(math.radians(lat1)) \
        * math.cos(math.radians(lat2)) * math.sin(dlon / 2) * math.sin(dlon / 2)
    c = 2 * math.atan2(math.sqrt(a), math.sqrt(1 - a))
    d = radius * c

    return d

data['lat_lon_distance_miles'] = data.apply(lambda row: distance((row['latitude1'], row['longitude1']), (row['latitude2'], row['longitude2'])), axis=1)