fix performance degration of sjoin when using pygeos

geopandas/tools/sjoin.py

@@ -146,24 +146,35 @@ def sjoin(
             r_idx = np.concatenate([[i] * len(v) for i, v in idxmatch.iteritems()])
 
     if len(r_idx) > 0 and len(l_idx) > 0:
-        # Vectorize predicate operations
-        def find_intersects(a1, a2):
-            return a1.intersects(a2)
+        if compat.USE_PYGEOS:
+            import pygeos
+
+            predicate_d = {
+                "intersects": pygeos.intersects,
+                "contains": pygeos.contains,
+                "within": pygeos.contains,
+            }
+            check_predicates = predicate_d[op]
+        else:
+            # Vectorize predicate operations
+            def find_intersects(a1, a2):
+                return a1.intersects(a2)
 
-        def find_contains(a1, a2):
-            return a1.contains(a2)
+            def find_contains(a1, a2):
+                return a1.contains(a2)
 
-        predicate_d = {
-            "intersects": find_intersects,
-            "contains": find_contains,
-            "within": find_contains,
-        }
+            predicate_d = {
+                "intersects": find_intersects,
+                "contains": find_contains,
+                "within": find_contains,
+            }
 
-        check_predicates = np.vectorize(predicate_d[op])
+            check_predicates = np.vectorize(predicate_d[op])
 
         if compat.USE_PYGEOS:
             res = check_predicates(
-                left_df.geometry[l_idx], right_df[right_df.geometry.name][r_idx]
+                left_df.geometry[l_idx].array.data,
+                right_df[right_df.geometry.name][r_idx].array.data,
             )
         else:
             res = check_predicates(