使用pyproj库实现地理坐标向投影坐标转换,一开始用的是:
from pyproj import Proj
p1 = pyproj.Proj(init="epsg:4326")
p2 = pyproj.Proj(init="epsg:32649")
x2, y2= pyproj.transform(p1,p2,x1,y1)
发生报错:pyproj FutureWarning: '+init=<authority>:<code>' syntax is deprecated
查了一下官方文档,发现不推荐使用此功能。请参阅:从pyproj 1升级到pyproj 2
参考:https://pyproj4.github.io/pyproj/stable/gotchas.html#upgrade-transformer
系统安装的pyproj版本是:3.0.1
修改后:
from pyproj import Transformer
transformer = Transformer.from_crs("epsg:4326", "epsg:32649")
x2, y2=transformer.transform(x1,y1)