0

我想在 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 string wintri,但不知道如何。

我似乎wintriv2.1.3 的源代码中找到了提及,但是我似乎在当前稳定版本的文档中找不到任何内容。

最后,我遇到了这种交流,这让我怀疑在 PROJ 以及扩展 GeoPandas 中无法使用 Winkel Tripel 投影。

该线程是从 2015 年开始的,我还没有发现任何关于此的发展。

总之,我是否正确理解不能将 Winkel Tripel 投影与 GeoPandas 一起使用?

或者,如果可能的话,我错过了什么?

任何帮助深表感谢。

4

1 回答 1

3

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

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

在此处输入图像描述

于 2020-11-17T11:34:31.120 回答