经纬度批量转cgcs2000 代码
时间: 2023-09-05 15:09:03 浏览: 218
经纬度转换代码
以下是经纬度批量转CGCS2000的Python代码示例,需要使用到pyproj库:
```python
import pyproj
# 定义转换器
wgs84 = pyproj.Proj('EPSG:4326')
cgcs2000 = pyproj.Proj('EPSG:4490')
# 定义经纬度列表
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的投影坐标系。然后,定义经纬度列表,这里使用了三个示例经纬度。接着,使用`pyproj.transform`函数将经纬度转换为CGCS2000坐标。最后,将转换结果打印输出。
阅读全文