python - GeoDataFrame.plot 有效但不显示在 folium GeoJson 中

Question

编辑:添加生成多边形的代码。 df 是从 csv 加载并处理成簇的基础数据:

df.head(1)

---

Requirement     PickupTime  Year    Month   DayOfWeek   Hour    Town    StreetName  StreetNr    Latitude    Longitude   GeoHash     cluster     dbcluster
0   NaN     2013-01-21 13:20:17.000     2013    1   1   13  Mikkeli     PORRASSALMENKATU    16  61.68844    27.274791   udgquw  17  29

然后通过组合属于一个簇的点来创建多边形:

polygons = []
for i in range(30):
    latlong = df[['Latitude','Longitude']].loc[df['cluster'] == i]
    polygons.append([i,MultiPoint(latlong.to_numpy()).convex_hull
                     ,centerOfRegions[i][0],centerOfRegions[i][1]])

然后将其加载到数据框中

polygondf = pd.DataFrame(polygons, columns=['cluster', 'geometry','center_lat','center_long'])
#polygondf['count'] = 1
polygondf.head(2)

---

cluster     geometry    center_lat  center_long
0   0   POLYGON ((61.6947869 27.2714365, 61.6943701 27...   61.692932   27.278602
1   1   POLYGON ((61.699156 27.2833976, 61.69188 27.28...   61.696324   27.290567

现在将此数据框对象加载到 GeoDataFrame 对象中:

import geopandas as gpd
frm = gpd.GeoDataFrame(polygondf, crs="EPSG:4326", geometry='geometry')
frm.head()

---

输出:

    cluster     geometry    center_lat  center_long     count
0   0   POLYGON ((61.69479 27.27144, 61.69437 27.27188...   61.692932   27.278602   1
1   1   POLYGON ((61.69916 27.28340, 61.69188 27.28595...   61.696324   27.290567   1
2   2   POLYGON ((62.02711 27.76096, 61.90189 27.79479...   61.843740   28.136880   1
3   3   POLYGON ((61.81762 27.06707, 61.56495 27.11654...   61.681207   27.183761   1
4   4   POLYGON ((61.69780 27.25109, 61.68140 27.25391...   61.689157   27.259389   1

frm['geometry']中的多边形是有效的，可以绘制如下图:

frm.plot(figsize=(6, 6))
plt.show()

但是当尝试在 folium map 上绘制相同的多边形时， map 显示没有绘制多边形。

m = folium.Map([61.681659, 27.271036], zoom_start=8, tiles='cartodbpositron')
for _, r in frm.iterrows():
    sim_geo = gpd.GeoSeries(r['geometry'],crs="EPSG:4326").simplify(tolerance=0.001)
    #sim_geo = sim_geo.to_crs("EPSG:4326")
    geo_j = sim_geo.to_json()
    geo_j = folium.GeoJson(data=geo_j,
                           style_function=lambda x: {'fillColor': 'orange'})
    folium.Popup(r['cluster']).add_to(geo_j)
    geo_j.add_to(m)
m

---

我能知道是什么导致了这个问题吗？Folium map 工作正常，我也能够将多边形绘制为 Choropleth。但这是我无法解决的问题。

Answer 1

frm 可能没有在 folium map 上显示的正确 crs。

即使使用 crs 参数创建 GeoDataFrame，它似乎也不会更改原始 crs。

polygondf = gpd.read_file('shapefile.shp')
frm = gpd.GeoDataFrame(polygondf, crs="EPSG:4326", geometry='geometry')
frm.crs # epsg:3006
frm = frm.to_crs(epsg=4326) # now crs is updated

编辑:

身材匀称MultiPoint有一个 (x, y[, z ]) 的构造函数，这意味着(经度，纬度)的顺序，请参见 post .以 longlat 而不是 latlong 顺序创建你的 polgyons。