一般而言,关于点的坐标系转换非常方便,调用pyproj的Proj函数即可,如下
# 将wgs坐标系转换到utm50
utm50 = pyproj.Proj(proj='utm', zone=50, ellps='WGS84') # https://proj4.org/operations/projections/utm.html#usage 默认的ellps是GRS80
xs, ys = utm50(acc.lng.values