如何在python中使用经纬度数据计算最小距离

import math
from geopy.distance import geodesic

lat1 = 40.7128  # 纬度1
lon1 = -74.0060  # 经度1
lat2 = 34.0522  # 纬度2
lon2 = -118.2437  # 经度2

def haversine(lat1, lon1, lat2, lon2):
    R = 6371  # 地球平均半径（单位：千米）

    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))
    distance = R * c

    return distance

min_distance = haversine(lat1, lon1, lat2, lon2)

start = (lat1, lon1)
end = (lat2, lon2)
min_distance = geodesic(start, end).kilometers