Fixed fallback to shapely spatial select
@@ -109,7 +109,6 @@ def spatial_select_columnar(xvals, yvals, geometry):
poly = Polygon([np.concatenate([geometry, geometry[:1]]).flatten()])
geom_mask = points.intersects(poly)
except Exception:
- pass
try:
from shapely.geometry import Point, Polygon
points = (Point(x, y) for x, y in zip(masked_xvals, masked_yvals))
