Python 经纬度转高斯坐标

好的，我可以回答这个问题。你可以使用 pyproj 库来将经纬度转换为高斯坐标。下面是一个示例代码：

import pyproj

# 创建 WGS84 坐标系对象
wgs84 = pyproj.Proj('EPSG:4326')

# 创建高斯-克吕格坐标系对象
gcj02 = pyproj.Proj('EPSG:4479')

# 深圳市中心的经纬度坐标
lon, lat = 114.06667, 22.61667

# 将经纬度坐标转换为高斯坐标
x, y = pyproj.transform(wgs84, gcj02, lon, lat)

print(f'经度：{x}, 纬度：{y}')

其中，EPSG:4326 是 WGS84 坐标系的标识符，EPSG:4479 是高斯-克吕格坐标系的标识符。lon 和 lat 分别是经纬度坐标的经度和纬度值。使用 pyproj.transform() 函数将经纬度坐标转换为高斯坐标，结果保存在 x 和 y 变量中。

相关问题

使用python实现高斯-克吕格投影将经纬度转换

可以使用pyproj库来实现高斯-克吕格投影的经纬度转换，下面是一个简单的示例代码：

import pyproj

# 定义高斯-克吕格投影的参数
proj_gk = pyproj.Proj(proj='tmerc', lon_0=120.0, lat_0=0.0, k_0=1.0, x_0=500000.0, y_0=0.0, ellps='WGS84')

# 定义经纬度坐标
lng = 121.0
lat = 31.0

# 将经纬度转换为高斯-克吕格投影坐标
x, y = proj_gk(lng, lat)

# 输出结果
print("高斯-克吕格投影坐标：{}, {}".format(x, y))

在上面的代码中，使用pyproj.Proj函数定义了高斯-克吕格投影的参数，然后使用proj_gk函数将经纬度坐标转换为高斯-克吕格投影坐标。在定义高斯-克吕格投影参数时，需要指定投影类型为'tmerc'（即transverse mercator），然后指定中央经线的经度（lon_0）、标准纬度（lat_0）、比例因子（k_0）、东偏距（x_0）和北偏距（y_0）等参数。

gis经纬度转换为xy坐标

将经纬度转换为平面坐标需要进行投影变换，常用的投影方式有墨卡托投影、高斯-克吕格投影等。其中，墨卡托投影是最常用的投影方式之一。

具体步骤如下：

1. 将经纬度转换为弧度制。
2. 根据所选的投影方式计算出相应的投影系数。
3. 根据经度和纬度计算出对应的墨卡托坐标。
4. 将墨卡托坐标转换为平面坐标。

以下是一个示例代码（使用 Python 实现）：

import math

def lon_lat_to_mercator(lon, lat):
    r_major = 6378137.000
    x = r_major * math.radians(lon)
    scale = x / lon
    y = 180.0 / math.pi * math.log(math.tan(math.pi / 4.0 + lat * (math.pi / 180.0) / 2.0)) * scale
    return x, y

# 测试
print(lon_lat_to_mercator(116.4074, 39.9042))  # (12955470.331355846, 4837470.484524184)

上述代码实现的是将经度为 116.4074，纬度为 39.9042 的点转换为墨卡托坐标系中的坐标。