Allow point comparison (#149)

Author: hoxbro

spatialpandas/geometry/_algorithms/intersection.py

@@ -511,18 +511,12 @@ def polygons_intersect_bounds(
     if y1 < y0:
         y0, y1, = y1, y0
 
-    if x0 == x1 or y0 == y1:
-        # Zero width/height rect does not intersect with anything
-        return
-
     for i in range(n):
         _perform_polygon_intersect_bounds(
             i, x0, y0, x1, y1, flat_values,
             start_offsets0, stop_offsets0, offsets1, result
         )
 
-    return result
-
 
 @ngjit
 def multipolygons_intersect_bounds(
@@ -561,10 +555,6 @@ def multipolygons_intersect_bounds(
     if y1 < y0:
         y0, y1, = y1, y0
 
-    if x0 == x1 or y0 == y1:
-        # Zero width/height rect does not intersect with anything
-        return
-
     # Populate results
     for i in range(n):
         polygon_offsets = offsets1[start_offsets0[i]:stop_offsets0[i] + 1]

spatialpandas/tests/geometry/test_cx.py

@@ -99,6 +99,18 @@ def test_polygon_cx_selection(gp_polygon, rect):
             assert all(expected == result)
 
 
+@pytest.mark.slow
+def test_polygon_cx_selection_point():
+    import geopandas as gpd
+    from shapely.geometry import shape
+
+    geometry = {"type": "Polygon", "coordinates": [[[0, 50], [20, 50], [50, 20], [50, 0]]]}
+    gp_multipolygon = gpd.GeoSeries([shape(geometry)])
+    xslice, yslice = slice(None, 0, None), slice(50, None, None)
+    expected = PolygonArray.from_geopandas(gp_multipolygon.cx[xslice, yslice])
+    result = PolygonArray.from_geopandas(gp_multipolygon).cx[xslice, yslice]
+    assert all(expected == result)
+
 @pytest.mark.slow
 @given(
     st_multipolygon_array(min_size=1, geoseries=True),
@@ -122,6 +134,18 @@ def test_multipolygon_cx_selection(gp_multipolygon, rect):
             assert all(expected == result)
 
 
+@pytest.mark.slow
+def test_multipolygon_cx_selection_point():
+    import geopandas as gpd
+    from shapely.geometry import shape
+
+    geometry = {"type": "MultiPolygon", "coordinates": [[[[0, 50], [20, 50], [50, 20], [50, 0]]]]}
+    gp_multipolygon = gpd.GeoSeries([shape(geometry)])
+    xslice, yslice = slice(None, 0, None), slice(50, None, None)
+    expected = MultiPolygonArray.from_geopandas(gp_multipolygon.cx[xslice, yslice])
+    result = MultiPolygonArray.from_geopandas(gp_multipolygon).cx[xslice, yslice]
+    assert all(expected == result)
+
 @given(st_multipoint_array(min_size=1, geoseries=True), st_bounds(orient=True))
 @hyp_settings
 def test_multipoint_cx_series_selection(gp_multipoint, rect):