topic Re: Select points within polygons i can't for the life of me get it to work in ArcGIS API for Python Questions

Select points within polygons i can't for the life of me get it to work

StuartMoore — Wed, 10 Apr 2024 09:04:44 GMT

help, i have two datasets a FGDB with a polygon feature class containing over 350k polygons and a csv file with an xy point of a little over 275k records.

using Arcgis api for python (on a computer that doesn't have any ESRI software like Pro or Desktop) i've successfully read the polygons into a spatial dataframe and read the csv file also into a spatial dataframe

but when i try and identify which points intersect the polygons i get this error:

i have tried a few other ways but non seem to work, even the documentation doesn't seem to help

any ideas?

Re: Select points within polygons i can't for the life of me get it to work

DanPatterson — Wed, 10 Apr 2024 10:50:57 GMT

Most likely that the points and polygons are in different coordinate systems and the extents don't overlap (last error message) . If this is the case, the csv needs to be converted with a know coordinate system to something that the coordinate systems are know and/or the same

Re: Select points within polygons i can't for the life of me get it to work

StuartMoore — Wed, 10 Apr 2024 12:01:09 GMT

Hi Dan, thanks for that, i checked and both have the same coordinate system with then wkid of 27700.

from manually doing the intersection in Pro on another computer around 260k of the points fall within polygons

i did also try this method i found on another post and this looks more promising but flags some geometry issues which i need to work out how to fix...

Re: Select points within polygons i can't for the life of me get it to work

Clubdebambos — Sun, 14 Apr 2024 06:55:57 GMT

Hi @StuartMoore

You could use spatial join.

# Load point and polygon spatial dataframes points_gdf = gpd.read_file('path_to_points.shp') polygons_gdf = gpd.read_file('path_to_polygons.shp') # Perform spatial join intersection = gpd.sjoin(points_gdf, polygons_gdf, how="inner", op='intersects') # intersection dataframe will contain points that intersect with polygons print(intersection)

Re: Select points within polygons i can't for the life of me get it to work

StuartMoore — Mon, 22 Apr 2024 06:48:47 GMT

thanks i tried that but thats where i was getting the geometry issues