- 读取几何图形
首先，让我们从阅读几何图形开始。
来自 csv
geoms = pd.read_csv('geometries.csv')想象一下文件在列几何体下包含多边形，我们现在必须将它们转换为几何类型（默认情况下它们将被读取为字符串）。通常，几何图形将以wkt格式，我们用shapely库进行转换
# if the geometries are stored as well-know text
geoms['geometry'] = geoms['geometry'].apply(lambda g: wkt.loads(g))来自形状文件
如果您要读取的文件是 shapefile 那么我们可以使用 geo-pandas 库，它将自动处理几何类型
import geopandas as gpd
gdf = gpd.read_file('geometries.shp')2 .更改坐标参考系 (CRS)
在处理空间数据时，通常需要从一种坐标系传递到另一种坐标系。地球的每个部分都有其特定的坐标参考系统或 CRS，它可以最大限度地减少在 2D 平面上处理 3D（几乎）球面坐标（例如参考地球上的点的坐标）的误差。最常见的 CRS 标有名为 epsg 的特定代码。
使用pyproj库的Proj和Transform函数（点投影）
from pyproj import Proj, transform
# the crs of the point (in this case GPS Lat/Lon)
inProj = Proj(init='epsg:4326')
# the desired crs (in this case Lambert, French crs)
outProj = Proj(init='epsg:
