无法使用geopandas转换到Mercator投影

Question

我正在尝试使用Geopandas和Matplotlib绘制西班牙地图。我使用的是您在下面看到的GeoJSON文件，其中的几何体是地理坐标（EPSG 4326）

我想在墨卡托投影中绘制地图，但当我尝试使用

df.to_crs({'init': 'epsg:3395'})

我收到一个错误说：

RuntimeError: b'no arguments in initialization list'

我已按照本手册中的说明进行操作http://geopandas.org/projections.html. 其中指出，如果geopandas数据帧不与CRS的信息一起计算，则必须对其进行定义。但这不是问题所在。

这是我目前掌握的代码

import geopandas as gpd
%matplotlib inline
import matplotlib.pyplot as plt

geojson_url = 'https://raw.githubusercontent.com/codeforamerica/click_that_hood/master/public/data/spain-provinces.geojson'
df = gpd.read_file(geojson_url)

如果我运行以下命令：

df.crs

我收到以下响应（这意味着geopandas数据帧具有关于CRS的信息）

{'init': 'epsg:4326'}

然后

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

返回一个错误

我错过了什么？提前谢谢。

Answer 1

您似乎正在使用较旧版本的geopandas和pyproj。我建议升级到最新版本，这样做可以：

>>> import geopandas
>>> geopandas.__version__
'0.5.0'
>>> import pyproj
>>> pyproj.__version__
'2.2.0'
>>> geojson_url = 'https://raw.githubusercontent.com/codeforamerica/click_that_hood/master/public/data/spain-provinces.geojson'
>>> df = geopandas.read_file(geojson_url)
>>> df.crs
{'init': 'epsg:4326'}
>>> df2 = df.to_crs({'init': 'epsg:3395'})
>>> df2.crs
{'init': 'epsg:3395'}