Skip to main content
15 events
when toggle format what by license comment
Mar 6, 2023 at 16:31 answer added Hope timeline score: 1
Nov 12, 2019 at 2:02 answer added koshmaster timeline score: 1
Apr 26, 2018 at 15:54 comment added Stephen This is the one. @Jon I posted my workaround with this below
Apr 26, 2018 at 15:53 vote accept Stephen
Apr 26, 2018 at 15:53 answer added Stephen timeline score: 3
Apr 25, 2018 at 14:12 comment added Jon Ah, I guess you're using shapely's intersects method. You could try int_gdf = gpd.sjoin(hulls_gdf, water_gdf, op='intersects').
Apr 25, 2018 at 14:06 comment added Jon What is the variable "polygons" in your script? I ask because you refer to the hulls as hulls_gdf but you write "polygons" in your code. How are you checking the CRS and what are the epsgs? If you type "hulls_gdf['geometry'][0]" does it plot the expected geometry of the first polygon in the gdf? (I'm assuming you're using an IDE) @user30184 If you intersect two geodataframes, you get a geodataframe as an output. Not clear if OP is intersecting geometries or geodataframes.
Apr 25, 2018 at 13:59 comment added Stephen Good point, the following code also returns no results: aff_playas = playas[playas.intersects(polygons)]
Apr 25, 2018 at 13:54 history edited Stephen CC BY-SA 3.0
added 36 characters in body
Apr 24, 2018 at 6:14 comment added user30184 Doesn't "Intersects" return a boolean value? What do you expect to get as a sum of Trues and Falses?
Apr 23, 2018 at 21:35 history edited Stephen CC BY-SA 3.0
deleted 378 characters in body
Apr 23, 2018 at 19:17 comment added Stephen I added some more information about my code in the question, and I checked the crs for both layers and the seem to be consistent.
Apr 23, 2018 at 19:15 history edited Stephen CC BY-SA 3.0
added 471 characters in body
Apr 23, 2018 at 18:42 comment added user30184 Good idea to check the projections of the both layers. Add the result to your question.
Apr 23, 2018 at 18:29 history asked Stephen CC BY-SA 3.0