Geopandas 中的 Winkel 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 教程，並重新投影到墨卡托。

首先，我嘗試使用 epsg 代碼重新投影到 Winkel Tripel。 然而，這個投影似乎沒有這樣的代碼（即我似乎無法在EPSG.io上找到這個投影）。

接下來，我嘗試查找支持的 epsg 代碼列表。 這讓我意識到這個功能實際上是由 pyproj 處理的。

雖然 PROJ 似乎確實支持Winkel Tripel ，但我還沒有弄清楚在這種情況下是否以及如何使此功能起作用。

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

我似乎在v2.1.3 的源代碼中發現了wintri的提及，但是我似乎在當前穩定版本的文檔中找不到任何內容。

最后，我遇到了這種交流，這讓我懷疑在 PROJ 和擴展 GeoPandas 中不可能使用 Winkel Tripel 投影。

該線程是從 2015 年開始的，我還沒有發現任何與此相關的進展。

總之，我是否正確理解不能將 Winkel Tripel 投影與 GeoPandas 一起使用？

或者如果可能的話，我錯過了什么？

任何幫助深表感謝。

Answer 1

如果我沒記錯的話，這是正確且簡單的方法：

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