python - Link each point in one GeoPandas dataframe to polygons in another dataframe

Question

I searched for my problem and found this question which is different from my issue.

I have two geo data frames, one contains houses locations as points (~700 points) and the other contains suburbs names and their polygon (~2973 polygons). I want to link each point to a polygon to assign each house to the correct suburb.

sample of my geo dataframe

polygon

import geopandas as gpd
from shapely.geometry import Point
from shapely.geometry.polygon import Polygon

#creating geo series
polys = gpd.GeoSeries({
    '6672': Polygon([(142.92288, -37.97886,), (141.74552, -35.07202), (141.74748, -35.06367)]),
    '6372': Polygon([(148.66850, -37.40622), (148.66883, -37.40609), (148.66920, -37.40605)]),
})

#creating geo dataframe
polysgdf = gpd.GeoDataFrame(geometry=gpd.GeoSeries(polys))
polysgdf

Which produces the following(my original geo dataframe also includes a suburb column that contains the suburb name but I couldn't add it to my sample, you can only see the suburb ID below)

        geometry
6672    POLYGON ((142.92288 -37.97886, 141.74552 -35.07202, 141.74748 -35.06367, 142.92288 -37.97886))
6372    POLYGON ((148.66850 -37.40622, 148.66883 -37.40609, 148.66920 -37.40605, 148.66850 -37.40622))

sample of the points geo dataframe

points

points=[Point(145.103,-37.792), Point(145.09720, -37.86400), 
        Point(145.02190, -37.85450)]

pointsDF = gpd.GeoDataFrame(geometry=points,
                                  index=['house1_ID', 'house2_ID', 'house3_ID'])

pointsDF

Which produces the following

            geometry
house1_ID   POINT (145.10300 -37.79200)
house2_ID   POINT (145.09720 -37.86400)
house3_ID   POINT (145.02190 -37.85450)

I would like the final output to be the pointsDF geo dataframe with each house assigned to the corresponding suburb. As a result of matching the points and the polygons.

Example:

suburbID subrubName    house_ID
6672      south apple  house1_ID
6372      water garden house2_ID

I am new to GeoPandas, I tried to explain my question in the clearest way possible. I am happy to clarify any point. Thank you.

Answer 1

I found a way to accomplish this by joining the two data frames using a spatial join

joinDF=gpd.sjoin(pointsDF, polysgdf, how='left',op="within")