使用python将经纬度几何投影到utm

Question

我有一个名为 eq 的地震数据数据框，其中列出了纬度和经度的列。使用 geopandas 我创建了一个点列，内容如下：

from geopandas import GeoSeries, GeoDataFrame
from shapely.geometry import Point
s = GeoSeries([Point(x,y) for x, y in zip(df['longitude'], df['latitude'])])
eq['geometry'] = s
eq.crs = {'init': 'epsg:4326', 'no_defs': True}
eq

现在我有一个带有经纬度坐标的几何列，但我想将投影更改为 UTM。任何人都可以帮助转型吗？

Answer 1

纬度/经度并不是真正的投影，而是一种默认的“非投影”。请参阅此页面了解更多详细信息，但这可能意味着您的数据使用WGS84或epsg:4326。

让我们构建一个数据集，在进行任何重投影之前，我们将定义crs为epsg:4326

import geopandas as gpd
import pandas as pd
from shapely.geometry import Point
df = pd.DataFrame({'id': [1, 2, 3], 'population' : [2, 3, 10], 'longitude': [-80.2, -80.11, -81.0], 'latitude': [11.1, 11.1345, 11.2]})
s = gpd.GeoSeries([Point(x,y) for x, y in zip(df['longitude'], df['latitude'])])

geo_df = gpd.GeoDataFrame(df[['id', 'population']], geometry=s)
# Define crs for our geodataframe:
geo_df.crs = {'init': 'epsg:4326'} 

我不确定你所说的“UTM投影”是什么意思。从维基百科页面我看到有 60 种不同的 UTM 投影，具体取决于世界的区域。您可以在线找到适当的epsg代码，但我只会给您一个随机epsg代码的示例。例如，这是 33N 区的区域

你如何进行重投影？您可以从投影上的 geopandas 文档轻松获取此信息。这只是一行：

geo_df = geo_df.to_crs({'init': 'epsg:3395'})

并且几何形状不再编码为纬度/经度：

    id  population  geometry
0   1   2   POINT (-8927823.161620541 1235228.11420853)
1   2   3   POINT (-8917804.407449147 1239116.84994171)
2   3   10  POINT (-9016878.754255159 1246501.097746004)