src / pysoilmap / shapeops.py

find_polygon_at_point

find_polygon_at_point(gdf, points)[source]

For each point in points look up the index of a containing polygon in a GeoDataFrame, or -1 if none is found.

Parameters:

  • gdf (GeoDataFrame) – a GeoDataFrame whose geometry is queried against

  • points (ndarray) – an (N, 2) coordinate array.

Return type:

ndarray

Returns:

an integer (N,) array