我創建了一個帶有 Pointsgeometry列的 GeoDataFrame:1. 創建dfdf = pd.DataFrame([[51.502687, -3.538329, 2242, 1, 47], [52.699185, -0.050122, 870, 2, 35], [51.574387, 0.397882, 651, 3, 47], [51.43874, 0.395791, 625, 4, 35], [51.23965, 0.561919, 614, 5, 36]], columns = ["lat","long","num_of_trucks","performance","num_of_routes"] )df 2. 創建gdf自dfgdf = gpd.GeoDataFrame(df, geometry=gpd.points_from_xy(df["lat"],df["long"]), crs={"init": "epsg:4326"})gdf3. 重新投影到使用米的 CRS 并計算buffer_radius每個點周圍10KMgdf.to_crs(epsg=3395,inplace=True)#gdf.to_crs(epsg=3857,inplace=True)#gdf.to_crs(epsg=27700,inplace=True)gdf["buffer_radius"] = gdf.geometry.buffer(10000)4. 將幾何列更改為新buffer_radius列gdf = gdf.set_geometry('buffer_radius')gdf.geometry.nameOut: 'buffer_radius'上述過程似乎產生了我想要的 GeoDataFrame和一個名為包含多邊形的gdf新“幾何”列。buffer_radius然后我想繪制這些新創建的多邊形,所以首先我轉換了列的 CRS buffer_radius:5. 重新投影到 CRS,讓我可以生成繪圖:gdf.to_crs(epsg=4326,inplace=True)6. 生成最終情節:然后我嘗試繪制多邊形,但它返回一個空圖:gv.Polygons(gdf)我的繪圖返回空這一事實讓我想知道這是否可能是投影問題?有誰知道我在這里做錯了什么?任何想法或指示將不勝感激。謝謝
生成緩沖區半徑多邊形 - 可能的投影問題
ibeautiful
2023-10-06 11:02:50