I have 2 geodataframes:
- An original shapefile of water features read in as a geodataframe. (playas)
- A gdf as a result of a convex_hull operation with only polygon geometries. (polygons)
There should be intersection between these two layers, but the following code returns an empty geodataframe:
aff_playas = playas[playas.intersects(polygons)]
I am trying to think of what the problem could be with these two layers. I believe it is one of the following:
- A projection inconsistency with the created hulls_gdf
- function failure due to how detailed, varied the playas layer is.
The two layers plotted below have similar extent and are verified to intersect in Qgis, ArcGIS. The third is the output from the intersects().