Hi everyone,

I have to identify overlapping polygons, with one of the datasets containing thousands of polygons. I am using the sf package and its st_intersects function, as:

dataframe1 %>%

st_filter(y = dataframe2, .predicate = st_intersects)

which takes about 6sec to compare each polygon of the 1st dataframe, and so, days for my current dataframes.

The only way I have encountered so far to make it possible is to first remove some manually and then split the dataframe to run the intersecting.

Would anyone have an advice on how to make it faster?

thanks a lot!