项目使用pyproj库进行对不同坐标系shapefile数据进行重新投影,具体实现形式是读取.prj文件,获取wkt定义的wkid,然后定影project投影函数,使用投影函数批量几何图形。根据shapely官方文档 ,使用了pyproj库定义了投影函数进行几何处理,官网例子使用了pyproj 2.X以上版本。
import pyproj
from shapely.geometry import Point
from shapely.ops import transform
wgs84_pt = Point(-72.2495, 43.886)
wgs84 = pyproj.CRS('EPSG:4326')
utm = pyproj.CRS('EPSG:32618')
project = pyproj.Transformer.from_crs(wgs84, utm, always_xy=True).transform
utm_point = transform(project, wgs84_pt)
pyproj库依赖proj库,根据pyproj官方文档指引,首先安装proj
- 下载32位或64位安装程序。
- 运行OSGeo4W安装程序。
- 选择 "高级安装 "并按下一步。
- 选择 "从互联网安装 "并按下一步。
- 选择一个安装目录。在大多数情况下,默认的建议是可以的。按 “下一步”。
- 选择 “本地软件包目录”。在大多数情况下,默认的建议是可以的。按下一步。
- 选择 "直接连接 "并按下一步。
- 选择download.osgeo.org服务器并按下一步。
- 在 "Commandline_Utilities "