围绕带标签的坐标创建形状优美的多边形

Question

我有一个城市建筑物的框架。这些建筑有一个纬度和经度列，并分组到单独的管理位置。

我已经成功地创建了一个geopandas数据框，每个建筑都有一个“几何”列，其中包含形状优美的点，并在地图上用不同的颜色为每个管理位置绘制这些点。但是，是否有办法按管理位置对geodataframe进行分组，并将“几何体”列作为一个形状良好的多边形，包含每个管理位置内每个建筑对应的点？

我目前有：

geo_poly['geometry'] = geo_poly['geometry'].apply(lambda x: x.coords[0])
geo_poly = geo_poly.groupby('management location')['geometry'].apply(lambda x: Polygon(x.tolist())).reset_index()
geo_poly = gpd.GeoDataFrame(geo_poly, geometry = 'geometry')

但当我绘制它时：

geo_poly.plot(ax = ax)

它显然通过遍历地理平台数据框中的每个点来绘制多边形。

我本以为geopandas会有（或者在某处存在）一种更好的方法来围绕已定义簇边缘的点绘制多边形，但我很难找到一个。

如有任何帮助，我们将不胜感激，谢谢

Answer 1

您可以使用GeoDataFrame的溶解函数来融合组中的所有点，然后使用convex_hull属性来提取围绕所有分组/融合/溶解点的多边形。

下面是一个可重复的小示例：

# Importing libraries used
import numpy as np
import geopandas as gpd
import shapely

# Setting the random seed
np.random.seed(6)

# Determining the number of points in the example
npoints = 15

# Generating the GeoDataFrame with random Points
mydf = gpd.GeoDataFrame({'id':range(npoints),
                         'geometry':(pd.Series(list(np.random.rand(npoints,2)))
                                     .apply(shapely.geometry.Point)),
                         'group':np.random.randint(1,3,npoints)})

# Plotting the GeoDataFrame according to the "group" column
mydf.plot(column='group', cmap='Set1')

# Fusing all of the Points based on the "group" column and 
# generating a Polygon that envelops all of the Points in
# each group
new_df = (mydf.dissolve(by='group')
          .convex_hull
          .reset_index()
          .rename({0:'geometry'},axis=1))

# Plotting the original Point and the newly-generated Polygons
ax = new_df.plot(column='group',cmap='tab20c', alpha=.7)
mydf.plot(ax = ax, column='group', cmap='Set1')

以下是原始点的分组图：

这是包围点组的多边形图：

所以，对于你的具体例子，你可以这样做：

from shapely.geometry import Point

# Creating a new column with Point geometries
geo_poly['geometry_pt'] = geo_poly['geometry'].apply(lambda x: Point(x.coords[0]))

# Transforming into a GeoDataFrame
geo_poly = gpd.GeoDataFrame(geo_poly, geometry = 'geometry_pt')

# Dissolving and extracting the new Polygons
new_geo_poly = (geo_poly.dissolve(by='management location')
                .convex_hull
                .reset_index()
                .rename({0:'geometry'},axis=1))

在创建新的“geometry_pt”列时，可能需要调整该行。但一旦你修好了，剩下的就可以用了