nav[aria-label="Primary Navigation"] { padding: 0; & ul { list-style: none; width: 100%; display: flex; flex-direction: row; justify-content: start; align-items: start; gap: 30px; padding: 0; & li { margin: 0; } & ul li { list-style: none; } } }

ArcGIS GeoAccessor Overlay Erase

TylerT

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:

Find more posts tagged with

Accepted answers

All comments

Marshal

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().

TylerT

@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

davidanderson_nifc

Do the same process for df_h

then

which yields

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

TylerT

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

TylerT

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]]]}"

TylerT

@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

davidanderson_nifc

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.

Join the Conversation

Signed in members can post, follow updates, and more. New here? Register a free account.

Quick Links

Community Help
Find useful guides, FAQs, and documents to help you navigate and make the most of Esri Community.