Optimize setbacks joins

[MRossol]

Bug Description
sjoin here is super slow:

reVX/reVX/wind_setbacks/wind_setbacks.py

Line 385 in d216964

tmp = gpd.sjoin(features, cnty,

Expected behavior
Use spatial index and bounding box to speed it up

                idx_generator = features.sindex.intersection(cnty.total_bounds)
                idx = [i for i in idx_generator]
                buildings = features.iloc[idx]
                centroids = buildings.centroid
                incentroids = [(i, c) for i, c in enumerate(centroids) if c.within(cnty["geometry"].values[0])]
                idx2 = [i[0] for i in incentroids]
                buildings = buildings.iloc[idx2]
                buildings['geometry'] = buildings.buffer(setback)
                setbacks.extend((geom, 1) for geom in buildings['geometry'])

Charge code
WEPG 10348 08.01.01

[MRossol]

fixed in #86