### 实现空间坐标转换的Python程序 以下是基于提供的参考资料和专业知识编写的用于实现多种空间坐标系之间相互转换的Python代码。此代码涵盖了WGS84、GCJ02（火星坐标系）、BD09（百度坐标系）以及Web墨卡托之间的转换。 #### 转换逻辑概述 为了完成不同坐标系间的转换，需了解每种坐标系统的定义及其差异。例如，GCJ02 是在中国范围内应用的一种加密算法，而 BD09 则是在 GCJ02 的基础上再次偏移的结果。这些转换可以通过一系列数学公式来实现[^1]。 #### Python代码示例 以下是一个完整的 `lnglattransfer` 类，支持上述四种坐标系统之间的互转： ```python import math class LngLatTransfer: def __init__(self): self.x_pi = 3.14159265358979324 * 3000.0 / 180.0 @staticmethod def wgs84_to_gcj02(lng, lat): """ WGS84 -> GCJ02 (火星坐标系) :param lng: 经度 :param lat: 纬度 :return: 转换后的经度和纬度 """ if not LngLatTransfer.out_of_china(lng, lat): delta_lng, delta_lat = LngLatTransfer.delta(lng, lat) gcj_lng = lng + delta_lng gcj_lat = lat + delta_lat return round(gcj_lng, 6), round(gcj_lat, 6) else: return round(lng, 6), round(lat, 6) @staticmethod def gcj02_to_wgs84(lng, lat): """ GCJ02 -> WGS84 :param lng: 经度 :param lat: 纬度 :return: 转换后的经度和纬度 """ d_lng, d_lat = LngLatTransfer.delta(lng, lat) wgs_lng = lng - d_lng wgs_lat = lat - d_lat return round(wgs_lng, 6), round(wgs_lat, 6) @staticmethod def gcj02_to_bd09(lng, lat): """ GCJ02 -> BD09 (百度坐标系) :param lng: 经度 :param lat: 纬度 :return: 转换后的经度和纬度 """ z = math.sqrt(lng * lng + lat * lat) + 0.00002 * math.sin(lat * LngLatTransfer.x_pi) theta = math.atan2(lat, lng) + 0.000003 * math.cos(lng * LngLatTransfer.x_pi) bd_lng = z * math.cos(theta) + 0.0065 bd_lat = z * math.sin(theta) + 0.006 return round(bd_lng, 6), round(bd_lat, 6) @staticmethod def bd09_to_gcj02(lng, lat): """ BD09 -> GCJ02 :param lng: 经度 :param lat: 纬度 :return: 转换后的经度和纬度 """ x = lng - 0.0065 y = lat - 0.006 z = math.sqrt(x * x + y * y) - 0.00002 * math.sin(y * LngLatTransfer.x_pi) theta = math.atan2(y, x) - 0.000003 * math.cos(x * LngLatTransfer.x_pi) gcj_lng = z * math.cos(theta) gcj_lat = z * math.sin(theta) return round(gcj_lng, 6), round(gcj_lat, 6) @staticmethod def out_of_china(lng, lat): """ 判断给定经纬度是否位于中国境内 :param lng: 经度 :param lat: 纬度 :return: 是否超出范围 """ return not (lng > 73 and lng < 136 and lat > 3 and lat < 54) @staticmethod def delta(lng, lat): """ 计算偏移量 :param lng: 经度 :param lat: 纬度 :return: 偏移量 """ ret = [] if LngLatTransfer.out_of_china(lng, lat): return 0, 0 d_lng = (lng - 105.0) * 0.0066 + (lng - 105.0) * (lat - 35.0) * 0.00002 d_lat = (lat - 35.0) * 0.0066 + (lng - 105.0) * (lat - 35.0) * 0.00002 ret.append(d_lng) ret.append(d_lat) return ret[0], ret[1] # 使用实例 if __name__ == "__main__": converter = LngLatTransfer() # 示例：WGS84 to GCJ02 wgs_lng, wgs_lat = 116.397128, 39.907494 gcj_lng, gcj_lat = converter.wgs84_to_gcj02(wgs_lng, wgs_lat) print(f"WGS84({wgs_lng}, {wgs_lat}) -> GCJ02({gcj_lng}, {gcj_lat})") # 示例：GCJ02 to BD09 bd_lng, bd_lat = converter.gcj02_to_bd09(gcj_lng, gcj_lat) print(f"GCJ02({gcj_lng}, {gcj_lat}) -> BD09({bd_lng}, {bd_lat})") ``` 以上代码实现了主流地理坐标系之间的转换功能，并提供了详细的注释以便于理解其工作原理。 #### 关于Excel坐标的处理 如果需要从Excel文件中读取并转换坐标数据，则可结合pandas库进行操作。具体而言，先加载Excel文件中的DMS格式数据，再将其转化为小数度（DD），最后调用前述函数执行坐标转换[^2]。 ---