text: Title: ENH: add dwithin predicate and associated distance to sindex query / spatial join
  Body: PyGEOS and Shapely >= 2.0 support spatial index queries with `dwithin` predicate, but in order to use it we need to add support for the distance parameter to `sindex.query()` and `sjoin()`.
  
  
  