首頁猿問生成緩沖區(qū)半徑多邊形 - 可能的投影問題

生成緩沖區(qū)半徑多邊形 - 可能的投影問題

Python

ibeautiful 2023-10-06 11:02:50

我創(chuàng)建了一個(gè)帶有 Pointsgeometry列的 GeoDataFrame：1. 創(chuàng)建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. 創(chuàng)建gdf自dfgdf = gpd.GeoDataFrame(df, geometry=gpd.points_from_xy(df["lat"],df["long"]), crs={"init": "epsg:4326"})gdf3. 重新投影到使用米的 CRS 并計(jì)算buffer_radius每個(gè)點(diǎn)周圍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列g(shù)df = gdf.set_geometry('buffer_radius')gdf.geometry.nameOut: 'buffer_radius'上述過程似乎產(chǎn)生了我想要的 GeoDataFrame和一個(gè)名為包含多邊形的gdf新“幾何”列。buffer_radius然后我想繪制這些新創(chuàng)建的多邊形，所以首先我轉(zhuǎn)換了列的 CRS buffer_radius：5. 重新投影到 CRS，讓我可以生成繪圖：gdf.to_crs(epsg=4326,inplace=True)6. 生成最終情節(jié)：然后我嘗試?yán)L制多邊形，但它返回一個(gè)空?qǐng)D：gv.Polygons(gdf)我的繪圖返回空這一事實(shí)讓我想知道這是否可能是投影問題？有誰知道我在這里做錯(cuò)了什么？任何想法或指示將不勝感激。謝謝

查看完整描述

1 回答

倚天杖

TA貢獻(xiàn)1828條經(jīng)驗(yàn) 獲得超3個(gè)贊

這顯然是一個(gè)問題geoviews，它不使用活動(dòng)幾何列，而是使用名為的默認(rèn)列geometry。如果您只是使用緩沖區(qū) ( gdf["geometry"] = gdf.geometry.buffer(10000)) 覆蓋幾何列，或者確保要繪制的數(shù)據(jù)幀在幾何列中具有您想要的幾何圖形，那么它應(yīng)該可以工作。