经纬度信息 转换 直角X-Y坐标 & 绘图

1 创建坐标转换器transformer

import pyproj
# 定义输入和输出的坐标系统
wgs84 = pyproj.CRS('EPSG:4326')  # 输入经纬度坐标系统
utm = pyproj.CRS('EPSG:32650')   # 输出坐标系统（UTM投影坐标）
# 创建一个坐标转换器
transformer = pyproj.Transformer.from_crs(wgs84, utm, always_xy=True)

 2 遍历转换

map是包含经纬度信息的pandas数据库

>>> map.head()
          LON        LAT  NAME  TYPE
0  121.234263  31.392221  嘉定北站     1
1  121.235249  31.392030   中间点     3
2  121.235416  31.391939   中间点     3
3  121.234957  31.391274   中间点     3
4  121.234909  31.391205   中间点     3

 转换并储存至position

position = []
map = np.array(map)
annotation = []
for longitude,latitude,name,type in map:
    x, y = transformer.transform(longitude, latitude)
    position+=[[x,y]]
position = np.array(position)

3 绘图

初始化