I have a set of candidate points and a set of reference points in lon/lat coordinate space. I have created a bounding box around each of my candidate points, and I would like to check the reference points to see what reference points fall within each box.
What is the advantage of using a geopandas spatial join over something like slicing the geopandas frame with .cx[lon1:lon2, lat1:lat2] ? IS there a difference?
.cx
is non-indexed spatial operation, which can be relatively expensive. – martinfleis