使用geopandas +地图崩溃会话绘制热图(内核密度)

Question

最近，我在一个问题上发布了如何用地质公园绘制热图。有人建议我使用geoplot.kdeplot()，但是传递一个projection会使我的会话崩溃。

假设我的数据分别存储在df_points和df_map中：

> print(df_points)
        PointID                     geometry
0     204403876  POINT (-101.66700 21.11670)
1     204462769  POINT (-101.66700 21.11670)
2     144407530  POINT (-101.66700 21.11670)

> print(df_map)
     PolyID                                           geometry
0     01001  POLYGON ((-102.10641 22.06035, -102.10368 22.0...
1     01002  POLYGON ((-102.05189 22.29144, -102.05121 22.2...
2     01003  POLYGON ((-102.68569 22.09963, -102.69087 22.0...

这是我尝试过的，但是根据本教程，我应该将一个projection传递给热图。

# Import geoplot
import geoplot
import geoplot.crs as gcrs

# Plot points
ax = geoplot.kdeplot(df_points, shade=True, alpha=0.7)
# Plot polygons
geoplot.polyplot(df_map, ax=ax)

​

​

但是，如果我添加了projection=gcrs.AlbersEqualArea()，会话就会崩溃：

# Plot heatmap
ax = geoplot.kdeplot(df_points, cmap='Reds', shade=True,
                     projection=gcrs.AlbersEqualArea()) # This crashes my session ):
# Add polygons
geoplot.polyplot(df_map, ax=ax)

python3: geos_ts_c.cpp:3991: int GEOSCoordSeq_getSize_r(GEOSContextHandle_t，const：：geom：：协调员序列*，无符号int*)：断言‘0 != cs’失败。

如何传递projection=gcrs.AlbersEqualArea()并避免错误？

Answer 1

我是在评论前面的一个问题:我刚刚在Colab的官方地理图引用中运行了这个例子，并经历了类似的崩溃。经过一些研究后，我发现崩溃是由shaply、卸载和重新安装造成的，并按照引用显示了图表。请试试看。

pip uninstall shapely
pip install shapely --no-binary shapely

以下是帮助我的答案