I have a geodataframe with points, each with a scalar value. I also have a geodataframe with shapes. I want to add the value field of all points that intersect shapes to their respective shape. So if point 1 has value 10 and point 2 has value 3 and both intersect shape 1, shape 1 then has a field called value with the value 13.
I looked at using sjoin in geopandas, but I can only see it matching one point to one shape, not if there are multiple points in a shape. I suppose I could do a spatial inner join with both and have the separate entries for each point matched with a shape and then flatten it, but I don't know how to combine rows like that either.