在 Geopandas 中购买 Tripel 投影

Question

我想在 GeoPandas 0.8.1 中使用 Winkel Tripel 投影

考虑以下示例。

import geopandas as gpd

world = gpd.read_file(gpd.datasets.get_path('naturalearth_lowres'))
world = world[(world.name != "Antarctica") & (world.name != "Fr. S. Antarctic Lands")]
world['gdp_per_cap'] = world.gdp_md_est / world.pop_est
world = world.to_crs(epsg=3395)
world.plot(column = 'gdp_per_cap')

此示例基于 GeoPandas tutorial on managing projections，并重新投影到墨卡托。

首先，我尝试使用 epsg 代码重新投影到 Winkel Tripel。但是，这个投影好像没有这样的代码（即我在EPSG.io上好像找不到这个投影）。

接下来，我尝试查找受支持的 epsg 代码列表。这让我意识到这个功能实际上是由 pyproj 处理的。

虽然 PROJ 似乎确实支持 Winkel Tripel，但我还没有弄清楚在这种情况下是否以及如何让此功能发挥作用。

我怀疑我可以调用项目字符串 wintri，但不知道如何调用。

我似乎在 the source code of v2.1.3 中找到了 wintri 的提及，但是我似乎无法在当前稳定版本的文档中找到任何内容。

最后，我遇到了 this exchange，这让我怀疑无法在 PROJ 中使用 Winkel Tripel 投影，并扩展为 GeoPandas。

那个帖子是 2015 年的，我还没有发现这方面的任何进展。

总之，我是否正确理解了 GeoPandas 无法使用 Winkel Tripel 投影？

或者如果可能的话，我错过了什么？

非常感谢任何帮助。

Answer 1

如果我没记错的话，这是正确且简单的方法：

world = world.to_crs('+proj=wintri')
world.plot()

在 Geopandas 中购买 Tripel 投影

Winkel Tripel projection in Geopandas

python

geopandas

pyproj