Geopandas 中的 Winkel Triplel 投影

Question

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

考虑以下示例。

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，但我还没有弄清楚在这种情况下是否以及如何让这个功能发挥作用。

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

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

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

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

总之，我是否正确理解不能将 Winkel Triplel 投影与 GeoPandas 一起使用？

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

非常感谢任何帮助。

Answer 1

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

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