bug: Collection bbox filter now works when collection bbox is a point (#41)

Collection bbox search filter was not working when the collection bbox
was defined as a single point.
This PR fixes it.

---------

Co-authored-by: Guillermo Tesoro Calvo <guillermo.tesoro@c-ssystems.de>

Author: gtesoro

pyproject.toml

@@ -24,6 +24,7 @@ dependencies = [
     "brotli-asgi",
     "starlette",
     "typing_extensions",
+    "shapely"
 ]
 
 [project.urls]
@@ -59,6 +60,7 @@ dev = [
     "stdlib-list",
     "tox",
     "tox-uv",
+    "types-shapely"
 ]
 
 [tool.mypy]

stac_fastapi/eodag/core.py

@@ -63,6 +63,7 @@
     get_sortby_to_post,
 )
 from stac_fastapi.eodag.utils import (
+    check_poly_is_point,
     dt_range_to_eodag,
     format_datetime_range,
     is_dict_str_any,
@@ -276,12 +277,15 @@ async def all_collections(
         # bbox filter
         if bbox:
             bbox_geom = get_geometry_from_various(geometry=bbox)
+
             default_extent = [[-180.0, -90.0, 180.0, 90.0]]
             collections = [
                 c
                 for c in collections
-                if get_geometry_from_various(
-                    geometry=c.get("extent", {}).get("spatial", {}).get("bbox", default_extent)[0]
+                if check_poly_is_point(
+                    get_geometry_from_various(  # type: ignore
+                        geometry=c.get("extent", {}).get("spatial", {}).get("bbox", default_extent)[0]
+                    )
                 ).intersection(bbox_geom)
             ]
 

stac_fastapi/eodag/utils.py

@@ -24,6 +24,7 @@
 from urllib.parse import unquote_plus
 
 import orjson
+from shapely.geometry import Point, Polygon
 
 if TYPE_CHECKING:
     from typing import Any, Optional, Union
@@ -119,3 +120,17 @@ def dt_range_to_eodag(
     end = end.isoformat().replace("+00:00", "Z") if end else None
 
     return start, end
+
+
+def check_poly_is_point(poly: Polygon) -> Union[Point, Polygon]:
+    """
+    Check if the polygon is a Point and returns it if so.
+
+    :param poly: A Shapely polygon.
+    :returns: Either a `Point` or the unchanged `poly`.
+    """
+
+    if poly.area == 0 and poly.area == 0 and poly.bounds[0] == poly.bounds[2] and poly.bounds[1] == poly.bounds[3]:
+        return Point(poly.exterior.coords[0])
+    else:
+        return poly