topic Re: Does a GeoAccessor (SEDF) spatial join require arcpy? in ArcGIS API for Python Questions

Does a GeoAccessor (SEDF) spatial join require arcpy?

Jay_Gregory — Wed, 27 Apr 2022 17:07:31 GMT

I am using an conda environment with version 2.0.0 of the Python API installed, but not arcpy. This is running on a JupyterHub deployment.

I am trying to do a spatial join of two SEDFs (point and polygons) by running the line:

points.spatial.join(polygons[['SHAPE']], how="left", op="intersects")

Unfortunately, get no positive intersections (the length of the output is the same as the points SEDF, and the index_right field is all NaN).

If I run the exact same code on the same datasets in my ArcGIS Pro environment (arcpy is present) , using the same version of the Python API (2.0.0), I get a different result, the result I would expect (i.e. some points intersect the polygons).

There is nothing in the API documentation saying the spatial join requires arcpy, so I am wondering if something else could be wrong with the JuypterHub environment, or is the documentation incorrect.

Can anyone help verify?

Re: Does a GeoAccessor (SEDF) spatial join require arcpy?

CarlSunderman — Tue, 03 May 2022 19:56:52 GMT

from arcgis.features import GeoAccessor import geopandas as gpd import pandas as pd gdf_pts = gpd.read_file(r'C:\temp\point.shp') sdf_pts = pd.DataFrame.spatial.from_geodataframe(gdf_pts) gdf_poly = gpd.read_file(r'C:\temp\polygon.shp') sdf_poly = pd.DataFrame.spatial.from_geodataframe(gdf_poly ) sdf_pts.spatial.join(sdf_poly[['SHAPE']], how="left", op="intersects")

This seemed to work fine for me in jupyter notebooks, are your points and polys in the same projection? have you tried other featureclasses, shapefiles, etc?

Re: Does a GeoAccessor (SEDF) spatial join require arcpy?

Jay_Gregory — Tue, 03 May 2022 20:19:39 GMT

yes, everything is in the same projection - however I also do not have geopandas (and hence shapely) installed in the environment. Just the ArcGIS Python API. So possibly either shapely OR arcpy is required for intersections to work?

Re: Does a GeoAccessor (SEDF) spatial join require arcpy?

DonaldPike — Tue, 26 Sep 2023 17:24:21 GMT

Hi -
I'm currently encountering a similar problem and was wondering if you ever came up with a solution to your issue.

Thanks