R查找两个美国邮政编码列之间的距离

import requests
import math

# 腾讯云地理编码API密钥
api_key = 'YOUR_API_KEY'

def get_coordinates(zip_code):
    url = f'https://apis.map.qq.com/ws/geocoder/v1/?address={zip_code}&key={api_key}'
    response = requests.get(url)
    data = response.json()
    if data['status'] == 0:
        return data['result']['location']
    else:
        raise Exception(f'Geocode failed for {zip_code}')

def haversine(coord1, coord2):
    R = 6371  # 地球半径，单位为公里
    lat1, lon1 = coord1['lat'], coord1['lng']
    lat2, lon2 = coord2['lat'], coord2['lng']
    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

# 示例邮政编码
zip_code1 = '90001'
zip_code2 = '90002'

# 获取坐标
coord1 = get_coordinates(zip_code1)
coord2 = get_coordinates(zip_code2)

# 计算距离
distance = haversine(coord1, coord2)
print(f'Distance between {zip_code1} and {zip_code2} is {distance:.2f} km')