生成缓冲区半径多边形 - 可能的投影问题
Producing Buffer Radius Polygons - Possible Projection Issue
我创建了一个带有点 geometry 列的 GeoDataFrame：
1。创建df
df = 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。从 df
创建 gdf
gdf = gpd.GeoDataFrame(df, geometry=gpd.points_from_xy(df["lat"],df["long"]),
crs={"init": "epsg:4326"})
gdf
3。重新投影到使用米的 CRS 并计算每个点周围 10KM buffer_radius
gdf.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.name
Out: 'buffer_radius'
上述过程似乎生成了我想要的 GeoDataFrame gdf，其中包含名为 buffer_radius 的新“几何”列，其中包含多边形。
然后我想绘制这些新创建的多边形，所以首先我转换了 buffer_radius 列的 CRS：
5。重新投影到允许我制作情节的 CRS：
gdf.to_crs(epsg=4326,inplace=True)
6.制作最终情节：
然后我尝试绘制多边形，但它返回了一个空图：
gv.Polygons(gdf)
事实上我的地块 returns 是空的，这让我想知道这是否可能是投影问题？
有人知道我在这里做错了什么吗？任何想法或指示将不胜感激。
谢谢
这显然是 geoviews 的一个问题，它没有使用活动几何列，而是使用名为 geometry 的默认列。如果你只是用你的缓冲区覆盖几何列（gdf["geometry"] = gdf.geometry.buffer(10000)）或者确保你想要绘制的数据框在几何列中有你想要的几何，它应该可以工作。
我创建了一个带有点 geometry 列的 GeoDataFrame：
1。创建df
df = 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。从 df
创建 gdf
gdf = gpd.GeoDataFrame(df, geometry=gpd.points_from_xy(df["lat"],df["long"]),
crs={"init": "epsg:4326"})
gdf
3。重新投影到使用米的 CRS 并计算每个点周围 10KM buffer_radius
gdf.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.name
Out: 'buffer_radius'
上述过程似乎生成了我想要的 GeoDataFrame gdf，其中包含名为 buffer_radius 的新“几何”列，其中包含多边形。
然后我想绘制这些新创建的多边形，所以首先我转换了 buffer_radius 列的 CRS：
5。重新投影到允许我制作情节的 CRS：
gdf.to_crs(epsg=4326,inplace=True)
6.制作最终情节：
然后我尝试绘制多边形，但它返回了一个空图：
gv.Polygons(gdf)
事实上我的地块 returns 是空的，这让我想知道这是否可能是投影问题？
有人知道我在这里做错了什么吗？任何想法或指示将不胜感激。
谢谢
这显然是 geoviews 的一个问题，它没有使用活动几何列，而是使用名为 geometry 的默认列。如果你只是用你的缓冲区覆盖几何列（gdf["geometry"] = gdf.geometry.buffer(10000)）或者确保你想要绘制的数据框在几何列中有你想要的几何，它应该可以工作。