clusterMembers = {}
for clusterID, region in tqdm(enumerate(boundaries), total=len(boundaries)):
    for point in tqdm(np.vstack([x, y]).T):
        if region.is_point_inside(point):
            if clusterID not in clusterMembers:
                clusterMembers[clusterID] = list()
            clusterMembers[clusterID].append(point)