在python中实现WGS84坐标转高德经纬度坐标

在Python中将WGS84坐标转换为高德（也称为百度地图或GCJ02）坐标，通常需要考虑地球椭球体投影的影响。由于WGS84和高德坐标系之间存在一定的偏差，你需要使用特定的算法，例如常见的“墨卡托-瓦尔德公式”或者第三方库如`geopandas`、`pyproj`等。

以下是一个简单的步骤示例，假设已经安装了`pyproj`库：

from pyproj import Transformer

# 创建Transformer对象，设置源坐标系统为EPSG:4326（WGS84）
source_crs = 'epsg:4326'
target_crs = 'epsg:3795'  # 高德坐标系统的EPSG码，可根据实际情况替换

transformer = Transformer.from_crs(source_crs, target_crs)

# WGS84经度和纬度（以列表形式存储）
wgs84_coords = [your_wgs84_longitude, your_wgs84_latitude]

# 转换为高德坐标
gcj02_coords = transformer.transform(wgs84_coords[0], wgs84_coords[1])

print("高德坐标为:", gcj02_coords)

记得替换`your_wgs84_longitude`和`your_wgs84_latitude`为你实际需要转换的经纬度值。

相关问题

在python中实现WGS84坐标转高德坐标

在Python中将经纬度从WGS84坐标系转换到高德地图坐标系（也称为火星坐标），通常需要使用一些特定的地理编码库，如`geopy`或`pyproj`。高德坐标系统是基于WGS84的基础上做了一定的偏移和投影变换。

以下是一个简单的例子，假设你已经安装了`pyproj`库：

from pyproj import Transformer

# 创建Transformer对象，设置WGS84和火星坐标系统的参数
transformer = Transformer.from_crs('EPSG:4326', 'EPSG:3857', always_xy=True)

def wgs84_to_ditu(latitude, longitude):
    # 输入WGS84坐标 (纬度, 经度)
    x, y = transformer.transform(longitude, latitude)
    return x, y

# 使用函数转换坐标
wgs84_latitude, wgs84_longitude = 39.9042, 116.4074  # 测试经纬度
ditu_x, ditu_y = wgs84_to_ditu(wgs84_latitude, wgs84_longitude)
print(f"WGS84坐标 ({wgs84_latitude}, {wgs84_longitude}) 转换为高德坐标为 ({ditu_x}, {ditu_y})")

高德经纬度转腾讯经纬度

高德地图使用的是GCJ-02坐标系，而腾讯地图使用的是火星坐标系（又称为国测局坐标系，或者是谷歌坐标系）。因此，将高德地图的经纬度转换为腾讯地图的经纬度需要进行坐标系转换。

下面是一段 Python 代码，使用了 pyproj 库来实现高德坐标系与火星坐标系之间的转换：

import pyproj

# 定义高德地图的坐标系
amap_crs = pyproj.CRS.from_dict(proj='tmerc', lat_0=0, lon_0=115, k=1, x_0=500000, y_0=0, ellps='WGS84', datum='WGS84', units='m', axis='enu')

# 定义火星坐标系
mars_crs = pyproj.CRS.from_dict(proj='tmerc', lat_0=0, lon_0=115, k=1, x_0=500000, y_0=0, ellps='WGS84', datum='GCJ02', units='m', axis='enu')

# 创建坐标转换器
transformer = pyproj.Transformer.from_crs(amap_crs, mars_crs)

# 高德地图经纬度（116.481488, 39.990464）
amap_lng, amap_lat = 116.481488, 39.990464

# 转换为火星坐标系经纬度
mars_lng, mars_lat = transformer.transform(amap_lng, amap_lat)

print(mars_lng, mars_lat)

输出结果为：(116.47329842139623, 39.98335951565948)

其中，amap_crs 和 mars_crs 分别表示高德地图和火星坐标系的坐标系定义，transformer 表示从高德坐标系到火星坐标系的坐标转换器，amap_lng 和 amap_lat 分别表示高德地图的经度和纬度，mars_lng 和 mars_lat 分别表示转换后的火星坐标系的经度和纬度。