迭代2个数据帧以查找多边形中的点的最快方法

Question

我有两个数据帧，一个包含点坐标，另一个包含多边形边界（区域）。 迭代点列表并计算每个多边形中有多少点的最快方法是什么

以下是我目前的执行情况：

    for index, hotspot in hotspots.iterrows():
        for index, zone in zones.iterrows():
            polygon = shapely.wkt.loads(zone['the_geom'])
            point = shapely.wkt.loads(hotspot['the_geom'])
            if polygon.contains(point) == True:
                zonesNum[zone["OBJECTID"]] += 1

Answer 1

我最终使用了“应用”，这样运行速度快了13倍：

zones.apply(lambda row: row['polygon'].contains(point), axis=1)