cgcs2000 批量转换中央子午线 代码

以下是批量将经纬度转换为CGCS2000，并指定中央子午线的Python代码示例，同样使用到了pyproj库：

import pyproj

# 定义转换器
wgs84 = pyproj.Proj('EPSG:4326')
cgcs2000 = pyproj.Proj(proj='utm', zone=49, ellps='GRS80', units='m', datum='CGCS2000')

# 定义经纬度列表
lons = [116.397458, 120.153576, 113.264434]
lats = [39.909729, 30.287459, 23.129163]

# 转换经纬度
x, y, z = pyproj.transform(wgs84, cgcs2000, lons, lats, [0]*len(lons))

# 显示转换结果
print("CGCS2000坐标：")
for i in range(len(x)):
    print("{:.6f}, {:.6f}".format(x[i], y[i]))

在代码中，定义了WGS84和CGCS2000的投影坐标系，并将CGCS2000设置为UTM投影坐标系，并指定中央子午线为49度，单位为米，椭球体为GRS80，数据为CGCS2000。然后，定义经纬度列表，这里使用了三个示例经纬度。接着，使用`pyproj.transform`函数将经纬度转换为CGCS2000坐标。最后，将转换结果打印输出。