For my process I make a spatial join between a polygon grid and points :
join = gpd.sjoin(rotate_grid, gdf, op='contains')
It's look like this.
the result is a grid where each square of my grid is multiplied by as many points that are contained.
My goal is to have only one square and mean of a particular column, NdviM. First I just tried to dissolve square without mean calcul with
dissol = join.dissolve(by='geometry')
But it returns me
Somebody knows how to do this ?