我正在尝试从Geopandas GeoDataFrame中的多边形制作一个Choropleth map 。我想通过GeoDataFrame列之一中的值的分位数来符号化多边形。我正在尝试找出不同的选项,并找出最适合我的需求。任何对此的建议将不胜感激。

看来Geopandas确实已经具备执行此操作的能力:http://nbviewer.ipython.org/github/geopandas/geopandas/blob/master/examples/choropleths.ipynb

tracts.plot(column='CRIME', scheme='QUANTILES', k=3, colormap='OrRd')

尽管我找不到太多文档,但这种方法行得通。我希望能够添加一个显示分位数截止值的图例,但看来Geopandas图当前仅允许基于分类数据的图例。有人为此工作吗?

另外,我希望能够调整多边形轮廓的宽度。这可能吗?

作为替代选择,我一直在使用matplotlib中的面片。这似乎涉及更多,但似乎确实提供了更多的自定义选项。如果有必要沿着这条路走以建立图例,我可以跟进另一个问题,并包括到目前为止的代码。

谢谢您的帮助。

最佳答案

以下补丁已集成在Geopandas中,因此您现在可以执行以下操作:

tracts.plot(column='CRIME', scheme='QUANTILES', k=3, colormap='OrRd', legend=True)

我对geopandas的plot_dataframe函数做了一个小补丁,以在使用PySAL方案时启用图例。您可以在这里找到它:http://nbviewer.ipython.org/gist/jorisvandenbossche/d4e6efedfa1e4e91ab65(调整仅在if scheme is not None:之后的几行中)。

这使您可以执行以下操作:
ax = plot_dataframe(tracts, column='CRIME', scheme='QUANTILES', k=3, colormap='OrRd', legend=True)

得到这样的数字:

python - Geopandas GeoDataFame的Choropleth map-LMLPHP

关于python - Geopandas GeoDataFame的Choropleth map ,我们在Stack Overflow上找到一个类似的问题:https://stackoverflow.com/questions/31755886/

10-13 06:19