### 使用Python实现经纬度坐标转换 #### 地址与经纬度互转 为了在Python中完成地址和经纬度之间的相互转换，可以利用`geopy`库来简化这一过程。该库提供了多种地理编码服务接口的选择，比如Nominatim，它能够依据输入的地址返回对应的经纬度信息，反之亦然。 ```python from geopy.geocoders import Nominatim def address_to_coords(address): geolocator = Nominatim(user_agent="geoapiExercises") location = geolocator.geocode(address) return (location.latitude, location.longitude) def coords_to_address(latitude, longitude): geolocator = Nominatim(user_agent="specify_your_app_name_here") location = geolocator.reverse(f"{latitude}, {longitude}") return location.address ``` 以上代码展示了两个简单的函数定义，分别用于执行从地址到经纬度以及从经纬度到地址的转换操作[^1]。 #### 经纬度至平面坐标(X,Y)转换 对于某些应用场景而言，可能还需要进一步将球面坐标下的经纬度映射到二维平面上作为X-Y形式存在。这里提供了一个基于简单投影模型的方法： ```python import math EARTH_RADIUS = 6378137.0 # 半径单位为米 def lnglat_to_xy(lng, lat): x = EARTH_RADIUS * math.radians(lng) y = EARTH_RADIUS * math.log(math.tan((math.pi / 4) + (math.radians(lat)/2))) return x, y ``` 需要注意的是这种方法适用于较小范围内的近似计算；如果涉及到更大尺度的空间变换，则应该考虑采用更精确的地图投影方法[^3]。 #### WGS84坐标系下经纬度间的距离计算 当需要衡量两点间沿地球表面的实际直线距离时，可以通过Haversine公式来进行估算: ```python def haversine(lon1, lat1, lon2, lat2): """ Calculate the great circle distance between two points on the earth. Parameters: lon1(float), lat1(float): Longitude and Latitude of point A. lon2(float), lat2(float): Longitude and Latitude of point B. Returns: float: Distance in kilometers. """ # 将十进制度数转化为弧度制 lon1, lat1, lon2, lat2 = map(math.radians, [lon1, lat1, lon2, lat2]) d_lon = lon2 - lon1 d_lat = lat2 - lat1 a = math.sin(d_lat/2)**2 + math.cos(lat1) * math.cos(lat2) * math.sin(d_lon/2)**2 c = 2 * math.asin(math.sqrt(a)) r = 6371 # 地球平均半径，单位为公里 return c * r ``` 这段程序实现了给定两组经纬度参数后求解它们之间最短路径长度的功能，结果以千米计[^2]。