topic Re: ArcGIS GeoAccessor Overlay Erase in ArcGIS API for Python Questions

ArcGIS GeoAccessor Overlay Erase

TylerT — Sun, 13 Mar 2022 22:55:58 GMT

I have two SEDFs. df_w has 3 polygons. df_h has 2 polygons both within one of the df_w polygons. Ultimately I need to cut holes in df_w with the df_h polygon. The subject method is throwing an error. I’m missing something about the functionality of the GeoAccessor overlay although it seems straightforward. Anyway, any insight into the issue would be appreciated. Thank you. Tyler

Polygons Visualized:

Method and Error:

Re: ArcGIS GeoAccessor Overlay Erase

DavidAnderson_1701 — Mon, 14 Mar 2022 14:14:44 GMT

This worked for me

`temp_df.spatial.overlay(temp_df_2,'erase')`

It appears the overlay function is not set up with key word arguments.

Re: ArcGIS GeoAccessor Overlay Erase

TylerT — Mon, 14 Mar 2022 15:56:15 GMT

@DavidAnderson_1701

hmm. Removing the 'op=' did not do the trick. Same result. Would you mind posting a screenshot of your dataframes with the index in view. I'm curious to how many features you have and if they share indices. I am sensing overlay is trying to match index to index (from df1 to df2), but not sure.

Alternatively, try it with dataframes of different sizes, and different indices. The reason I'm wondering is because when I run op='union' I get this error.

Thank you.

Tyler

Re: ArcGIS GeoAccessor Overlay Erase

TylerT — Mon, 14 Mar 2022 16:07:39 GMT

Here are the polygons for testing. Thx. Tyler

df_w ,OID,SHAPE_PG 241,242,"{'spatialReference': {'wkid': None}, 'rings': [[[15291.5, 6192.937499999999], [15332.5, 6192.9375], [15332.5, 6341.9375], [15291.5, 6341.9375], [15291.5, 6192.937499999999]]]}" 242,243,"{'spatialReference': {'wkid': None}, 'rings': [[[15037.5, 6301.321927246139], [15037.499999999998, 6192.937499999999], [15285.5, 6192.937499999999], [15285.5, 6347.9375], [15332.5, 6347.9375], [15332.499999999998, 6591.937499999997], [14864.5, 6591.937499999999], [14864.5, 6192.937499999999], [14914.781651165493, 6192.937499999999], [14914.781651165496, 6301.321927246138], [15037.5, 6301.321927246139]]]}" 243,244,"{'spatialReference': {'wkid': None}, 'rings': [[[14994.781651165495, 6293.321927246138], [14922.781651165495, 6293.321927246138], [14922.781651165495, 6200.9375], [14951.5, 6200.937499999999], [14951.5, 6186.937499999999], [14987.5, 6186.937499999999], [14987.5, 6200.9375], [15029.5, 6200.9375], [15029.5, 6293.321927246138], [14994.781651165495, 6293.321927246138]]]}"

df_h ,OID,SHAPE_PG 245,246,"{'spatialReference': {'wkid': None}, 'rings': [[[14903.25, 6381.9375], [14969.25, 6381.9375], [14969.25, 6453.9375], [14903.25, 6453.9375], [14903.25, 6381.9375]]]}" 246,247,"{'spatialReference': {'wkid': None}, 'rings': [[[15119.264718869708, 6351.838419508824], [15168.764718869708, 6351.838419508824], [15168.764718869708, 6405.838419508824], [15119.264718869708, 6405.838419508824], [15119.264718869708, 6351.838419508824]]]}"

Re: ArcGIS GeoAccessor Overlay Erase

TylerT — Wed, 16 Mar 2022 17:35:42 GMT

@DavidAnderson_1701, circling back to see if you have screen shots for your jupiter lab notebook dataframe outputs. Thx. Tyler

Re: ArcGIS GeoAccessor Overlay Erase

DavidAnderson_1701 — Wed, 16 Mar 2022 21:28:27 GMT

Do the same process for df_h

then

which yields

Short story is it worked for me. Perhaps try recreating the spatial data frames.

Re: ArcGIS GeoAccessor Overlay Erase

TylerT — Wed, 16 Mar 2022 21:42:08 GMT

@DavidAnderson_1701 Interesting. Should the two holes be white? What you have shown looks like the holes over the top of the second larger feature. At least you didn't get errors like I got. Let me study this and work on my SEDFs.

Much appreciated.

Tyler

Re: ArcGIS GeoAccessor Overlay Erase

DavidAnderson_1701 — Wed, 16 Mar 2022 22:00:16 GMT

Probably something to do with how I'm displaying the data. I wrestled with how to get the dataframe to show the geometry.

The area values show that the erase did happen.

Not a very satisfying answer I know. Saying it worked for me is an old programmer excuse/cop out. You might double check that the dataframe geometries are indeed polygons using the spatial.area property.

Start a clean Jupyter project, try the code I show to see if that works. Isolate the problem code to just this little bit.

Re: ArcGIS GeoAccessor Overlay Erase

TylerT — Wed, 16 Mar 2022 22:06:01 GMT

To display try:

df.spatial

Example.

I'll keep wrestling with it, but I'm still getting a positional index error.

Thanks for helping.

Tyler

Re: ArcGIS GeoAccessor Overlay Erase

DavidAnderson_1701 — Wed, 16 Mar 2022 22:37:14 GMT

Could you post the full traceback stack text. It might be informative to see exactly which module is throwing the error. Is it a Python function? Pandas? Arcgis API?

Re: ArcGIS GeoAccessor Overlay Erase

TylerT — Thu, 17 Mar 2022 13:59:26 GMT

@DavidAnderson_1701, I finally made SOME progress. I switched machines and lo and behold, it worked. No error. Original machine has ArcGIS Pro 2.8.3 and the new machine has 2.9.1. I'm wondering if there's a python module bug. Do arcgis/arcpy modules get updated when ArcGIS Pro gets updated? Thanks for all your help troubleshooting. I'll try updating ArcGIS Pro and see what happens.

Side note: Jupyter displays holes oddly. They look more like feature over feature, so I checked in ArcGIS Pro and the output feature looks just fine.

Jupyter output feature with holes:

Same SEDF exported to gdb fc and displayed in ArcGIS Pro:

Cheers,

Tyler

Re: ArcGIS GeoAccessor Overlay Erase

Marshal — Tue, 07 Nov 2023 15:05:10 GMT

Just wanted to leave the solution that worked for me to help future visitors.

Using .reset_index() method on the spatial dataframes before executing spatial.overlay() eliminated this error.

This seems to occur when taking a slice of a dataframe and using that slice in spatial.overlay().