检查一层的哪个多边形包含另一层的多边形

Question

我有一个包含两层的 shapefile（cpg、dbf、prj、shp、shx 文件集）。 L1 代表城市的子区域，L2 代表特定建筑物。

L1 的属性之一是子区域的名称（NAME） L2 的属性之一是建筑物的邮政编码（A_PSTCD）。

我的任务是创建一个 DataFrame，它将具有唯一对“NAME/A_PSTCD”的列表。 为此，我需要将 L1 与 L2 匹配，然后提取相关字段。

我对 shapefile 很陌生，尽管我对 python 有一些经验。 任何有关如何解决此问题的建议将不胜感激

Answer 1

因此，我设法自己找到了解决方案：

1.导入库

import geopandas as gpd

2.加载建筑物文件

buildings_polys = gpd.read_file('.../buildings.shp')

3.加载区域文件

polygons_areas = gpd.read_file('.../areas.shp')

4.找到建筑物的质心

buildings_points = buildings_polys.copy()
buildings_points['geometry'] = buildings_points['geometry'].centroid

5.查看建筑物所在的区域

i = 0
for area in areas:
    i+=1
    print('{} of {}'.format(i, len(areas)), area)
    buildings_points[area] = np.where(buildings_points['geometry'].within(dictionary[area]), 1, 0)

这将创建一个 dataframe 建筑物，其中有一堆列（区域的完整列表），其值为 1 或 0，具体取决于建筑物是否在该区域中。

您也可以执行类似的操作（而不是步骤 5 中的最后一行：

buildings_points['area_of_building'] = np.where(buildings_points['geometry'].within(dictionary[area]), area, 0)

要获得一个单独的列，您将在其中获得该区域的名称。