mapchete-eo 2025.10.1__py2.py3-none-any.whl → 2026.1.0__py2.py3-none-any.whl

mapchete_eo/settings.py

@@ -16,6 +16,7 @@ class Settings(BaseSettings):
     default_cache_location: MPathLike = MPath("s3://eox-mhub-cache/")
     default_catalog_crs: CRS = CRS.from_epsg(4326)
     blacklist: Optional[MPathLike] = None
+    lazy_load_stac_items: bool = True
 
     # read from environment
     model_config = SettingsConfigDict(env_prefix="MAPCHETE_EO_")

mapchete_eo/sort.py

@@ -5,10 +5,8 @@ This module holds all code required to sort products or slices.
 from typing import Callable, List, Optional
 
 from pydantic import BaseModel
-from pystac import Item
 
-from mapchete_eo.io.items import get_item_property
-from mapchete_eo.protocols import DateTimeProtocol
+from mapchete_eo.protocols import DateTimeProtocol, GetPropertyProtocol
 from mapchete_eo.time import timedelta, to_datetime
 from mapchete_eo.types import DateTimeLike
 
@@ -51,11 +49,11 @@ class TargetDateSort(SortMethodConfig):
 
 
 def sort_objects_by_cloud_cover(
-    objects: List[Item], reverse: bool = False
-) -> List[Item]:
+    objects: List[GetPropertyProtocol], reverse: bool = False
+) -> List[GetPropertyProtocol]:
     if len(objects) == 0:  # pragma: no cover
         return objects
-    objects.sort(key=lambda x: get_item_property(x, "eo:cloud_cover"), reverse=reverse)
+    objects.sort(key=lambda x: x.get_property("eo:cloud_cover"), reverse=reverse)
     return objects
 
 

mapchete_eo/source.py

@@ -0,0 +1,107 @@
+from functools import cached_property
+from typing import Any, Dict, List, Literal, Optional, Generator, Union, Callable
+
+from mapchete.bounds import Bounds
+from mapchete.path import MPath
+from mapchete.types import BoundsLike, CRSLike, MPathLike
+from pydantic import BaseModel, ConfigDict, model_validator
+from pystac import Item
+from shapely.geometry.base import BaseGeometry
+from shapely.errors import GEOSException
+
+from mapchete_eo.exceptions import ItemGeometryError
+from mapchete_eo.search.base import CollectionSearcher
+from mapchete_eo.search import STACSearchCollection, STACStaticCollection
+from mapchete_eo.settings import mapchete_eo_settings
+from mapchete_eo.types import TimeRange
+
+
+class Source(BaseModel):
+    """All information required to consume EO products."""
+
+    collection: str
+    catalog_crs: Optional[CRSLike] = mapchete_eo_settings.default_catalog_crs
+    query: Optional[str] = None
+    area: Optional[Union[MPathLike, dict, type[BaseGeometry]]] = None
+    bounds: Optional[BoundsLike] = None
+
+    model_config = ConfigDict(arbitrary_types_allowed=True)
+
+    @property
+    def item_modifier_funcs(self) -> List[Callable]:
+        return []
+
+    @cached_property
+    def catalog_type(self) -> Literal["search", "static"]:
+        # TODO: stupid test but probably sufficient
+        return "static" if self.collection.endswith(".json") else "search"
+
+    def _spatial_subset(
+        self,
+        bounds: Optional[BoundsLike] = None,
+        area: Optional[BaseGeometry] = None,
+    ) -> Dict[str, Any]:
+        """Combine bounds and area with bounds defined in Source if any."""
+        if self.bounds is None:
+            return {"bounds": bounds, "area": area}
+        self_bounds = Bounds.from_inp(self.bounds)
+        out = dict()
+        if bounds is not None:
+            bounds = Bounds.from_inp(bounds)
+            if bounds.intersects(self_bounds):
+                out["bounds"] = Bounds.from_inp(
+                    bounds.geometry.intersection(self_bounds.geometry)
+                )
+        if area is not None:
+            out["area"] = area.intersection(self_bounds.geometry)
+        return out
+
+    def search(
+        self,
+        time: Optional[Union[TimeRange, List[TimeRange]]] = None,
+        bounds: Optional[BoundsLike] = None,
+        area: Optional[BaseGeometry] = None,
+        base_dir: Optional[MPathLike] = None,
+    ) -> Generator[Item, None, None]:
+        for item in self.get_catalog(base_dir=base_dir).search(
+            time=time,
+            query=self.query,
+            search_kwargs=dict(query=self.query) if self.query else None,
+            **self._spatial_subset(
+                bounds=bounds,
+                area=area,
+            ),
+        ):
+            yield self.apply_item_modifier_funcs(item)
+
+    def apply_item_modifier_funcs(self, item: Item) -> Item:
+        try:
+            for modifier in self.item_modifier_funcs:
+                item = modifier(item)
+        except GEOSException as exc:
+            raise ItemGeometryError(
+                f"item {item.get_self_href()} geometry could not be resolved: {str(exc)}"
+            )
+        return item
+
+    def get_catalog(self, base_dir: Optional[MPathLike] = None) -> CollectionSearcher:
+        match self.catalog_type:
+            case "search":
+                return STACSearchCollection(self.collection)
+            case "static":
+                return STACStaticCollection(
+                    collection=MPath(self.collection).absolute_path(base_dir=base_dir)
+                )
+
+    def eo_bands(self, base_dir: Optional[MPathLike] = None) -> List[str]:
+        return self.get_catalog(base_dir=base_dir).eo_bands
+
+    @model_validator(mode="before")
+    def deprecate_max_cloud_cover(cls, values: Dict[str, Any]) -> Dict[str, Any]:
+        if "max_cloud_cover" in values:
+            raise DeprecationWarning(
+                "'max_cloud_cover' will be deprecated soon. Please use 'eo:cloud_cover<=...' in the source 'query' field.",
+            )
+        elif "area" in values:  # pragma: no cover
+            raise NotImplementedError("please use 'bounds' as spatial subset for now")
+        return values

{mapchete_eo-2025.10.1.dist-info → mapchete_eo-2026.1.0.dist-info}/METADATA

@@ -1,6 +1,6 @@
 Metadata-Version: 2.4
 Name: mapchete-eo
-Version: 2025.10.1
+Version: 2026.1.0
 Summary: mapchete EO data reader
 Project-URL: Homepage, https://gitlab.eox.at/maps/mapchete_eo
 Author-email: Joachim Ungar <joachim.ungar@eox.at>, Petr Sevcik <petr.sevcik@eox.at>
@@ -21,6 +21,7 @@ Requires-Dist: mapchete[complete]>=2025.10.0
 Requires-Dist: opencv-python-headless
 Requires-Dist: pillow
 Requires-Dist: pydantic
+Requires-Dist: pygeofilter
 Requires-Dist: pystac-client>=0.7.5
 Requires-Dist: pystac[urllib3]>=1.12.2
 Requires-Dist: retry
@@ -32,9 +33,9 @@ Provides-Extra: docs
 Requires-Dist: sphinx; extra == 'docs'
 Requires-Dist: sphinx-rtd-theme; extra == 'docs'
 Provides-Extra: test
+Requires-Dist: pytest; extra == 'test'
 Requires-Dist: pytest-coverage; extra == 'test'
-Requires-Dist: pytest-lazy-fixture; extra == 'test'
-Requires-Dist: pytest<8; extra == 'test'
+Requires-Dist: pytest-lazy-fixtures; extra == 'test'
 Description-Content-Type: text/x-rst
 
 .. image:: logo/mapchete_eo.svg

mapchete_eo-2026.1.0.dist-info/RECORD

@@ -0,0 +1,89 @@
+mapchete_eo/__init__.py,sha256=Q2R77nEGEKjxOBz3qMm2JGtkZn-nVFmReBAVH_YeZ4M,25
+mapchete_eo/base.py,sha256=_6wbR9IZBBBdRBbZY4FEQ0uu603C-tMPYK4MGQ5xSYY,21953
+mapchete_eo/blacklist.txt,sha256=6KhBY0jNjXgRpKRvKJoOTswvNUoP56IrIcNeCYnd2qo,17471
+mapchete_eo/eostac.py,sha256=M-EL8y8pPw8H9hb85FPxSl-1FL_kIYOzxsTIhkrbBZY,570
+mapchete_eo/exceptions.py,sha256=ul7_9o6_SfJXQPsDOQqRBhh09xv2t4AwHl6YJ7yWN5s,2150
+mapchete_eo/product.py,sha256=0GJDk4pPvdCLxuJkZCMQnM7LIjvOBRUhIzchcQCqMjI,9644
+mapchete_eo/protocols.py,sha256=KUbteFTOCSYEiSUJZj2MeckmrJHfJFrXCvZE70tByP8,1617
+mapchete_eo/settings.py,sha256=lbwYoqiQXBRexxL7Y524Ry0dk4vBKpjN1u09UOP1_CM,734
+mapchete_eo/sort.py,sha256=DfYAGEM6vzHJOFhaKxcKr-GwH5xM5ar4_mlp0uiVUK0,1787
+mapchete_eo/source.py,sha256=Tfo5IIYdRWpS06K-M7mzesl4WIUEt4onGmc1uKJ0F14,4158
+mapchete_eo/time.py,sha256=mCIrn6C-B7CDoFIIM4u8pG15nVV8KWQixF2fuhdDTdE,1799
+mapchete_eo/types.py,sha256=yIHKZHlGCSerJzDBS2AH7yYLun4rY3rZZ73SCKiN26U,1814
+mapchete_eo/array/__init__.py,sha256=47DEQpj8HBSa-_TImW-5JCeuQeRkm5NMpJWZG3hSuFU,0
+mapchete_eo/array/buffer.py,sha256=GeeM7RVQ-Z_SKlkaeztWUrxEsa_yck-nXwXRIM3w0_I,483
+mapchete_eo/array/color.py,sha256=ArJ3-0qjJJGAlRQodLios6JwD0uHLSRRb9BwfB31v-U,810
+mapchete_eo/array/convert.py,sha256=oimklMfyS2XxQLd73zQ2cr9aE-39z1UIt6ry_2ZQI10,5390
+mapchete_eo/cli/__init__.py,sha256=SRRilPKUtwsUCi1GdMNPxjIrCM-XyYM2hK1T9vvKUI8,946
+mapchete_eo/cli/bounds.py,sha256=QcOXOXMcnd8ecwiCFVWukTLUFIkpj9SbB3ot-GA4heA,551
+mapchete_eo/cli/options_arguments.py,sha256=6VqiwswEKpr5bS0w1ugRSf2z7rMKA-wiEFZdcmri94c,6916
+mapchete_eo/cli/s2_brdf.py,sha256=siUGau8HzXNyQZkDUx1zAeaLcCfonJG5kWHSCvwJOhA,2729
+mapchete_eo/cli/s2_cat_results.py,sha256=tR5Qm6kiCAcDO5V9AkYhoUbHyEx7OItSf6IqtIjhRLg,4648
+mapchete_eo/cli/s2_find_broken_products.py,sha256=4nhtMLkys75enPV5lETw4iOn3CxxNnc-sMyhk_kaKx0,2946
+mapchete_eo/cli/s2_jp2_static_catalog.py,sha256=_j_rfFnYHtE8k79TS9KFuX_IxOx_e-r7u6beoLYyxPY,5911
+mapchete_eo/cli/s2_mask.py,sha256=1K5BiyTKyClGBdh6-Vh_ORTcgJOr6R5JfXtpEpDEQvk,2243
+mapchete_eo/cli/s2_mgrs.py,sha256=EJxbrg4I8IttWIcLUEmmqkjsTUbVLufgrkDVrRbbkYE,1276
+mapchete_eo/cli/s2_rgb.py,sha256=nzK5m4s0m2VOXgtit0uUOOhNJ14wquikWmbEQON5WJw,3747
+mapchete_eo/cli/s2_verify.py,sha256=atfJA5luHZjDSqiM2udaufmnQ2B4uW2Tur_OgtYLa2g,4080
+mapchete_eo/cli/static_catalog.py,sha256=O7oCJB6mxWwnEMowPBo6UMtWOu52mSNUBB9tpSw91a4,3005
+mapchete_eo/image_operations/__init__.py,sha256=c36Vs53mT5UQR9Iwd2Ln_vnT4K_Ni8g-K7J9Fl6UsJo,432
+mapchete_eo/image_operations/blend_functions.py,sha256=YaKLOEapMbsUR-_mkBatXwPcmim8GKpIAVua_GOdaHs,24131
+mapchete_eo/image_operations/color_correction.py,sha256=vMXxcdyH3qOJDSKcXGU3FgXWIcgQxMy-QVe2fGQTbdk,4447
+mapchete_eo/image_operations/compositing.py,sha256=K0bGCZqJuS39IBsjrRiJz5zowItI6RsiA-EBUVlbboo,8449
+mapchete_eo/image_operations/dtype_scale.py,sha256=_CCdTBkMLDkZdjV_PsVEkqLWT6jWuTJPeMzCcRwBOns,1243
+mapchete_eo/image_operations/fillnodata.py,sha256=Xdd-kacfo7VKqX7HRP0hqnnvV5OwNqnHLRTW64QhG24,4856
+mapchete_eo/image_operations/filters.py,sha256=sa_Igv0zMIANHLEALVi8mxohZMoSrdQc6XVjXDYwzy8,7201
+mapchete_eo/image_operations/linear_normalization.py,sha256=-eQX3WLCUYWUv-um3s1u33rgjAwCz84Ht8gcPsMWtJE,2468
+mapchete_eo/image_operations/sigmoidal.py,sha256=IKU8F89HhQJWGUVmIrnkuhqzn_ztlGrTf8xXZaVQWzU,3575
+mapchete_eo/io/__init__.py,sha256=1-1g4cESZZREvyumEUABZhDwgVuSxtxdqbNlLuVKlow,950
+mapchete_eo/io/assets.py,sha256=OVXLxrwCCUa6frP-xpm2CL30MR8eMABnwVD6EpN1cqU,17215
+mapchete_eo/io/items.py,sha256=5H2xWAyLptDuuHh0rjY1MXReJomV0XBZgMr_rpdpYXU,6244
+mapchete_eo/io/levelled_cubes.py,sha256=ZF7BLn9MHnJCCDjAoR9D7MNbHdusjJ9EACdM7rKNlyM,9018
+mapchete_eo/io/path.py,sha256=fCjvowd9rGJbGio5Z-jDLxa49M2BVurYM9hJy53Im48,4831
+mapchete_eo/io/products.py,sha256=2gELjHoEe4QH-Yjk82tA2nZUnBgfPbaGvk5GMHUWs4w,14642
+mapchete_eo/io/profiles.py,sha256=l9YiXbKYs674xz6NLy5EAq3fBEvHTVSf_gopXD-CuSY,935
+mapchete_eo/platforms/sentinel2/__init__.py,sha256=dviB4f5k45t_iOScTtnPc1kaWcohWKRq-8RaQzdggmU,401
+mapchete_eo/platforms/sentinel2/_mapper_registry.py,sha256=3kuyW2g_xJoFUzZrAmhKn9ASHl4gBKuOlHsSZnpt21w,2746
+mapchete_eo/platforms/sentinel2/bandpass_adjustment.py,sha256=DA0cQtjr8UH7r_kizAD-EmBsZvEtVVThf0hD0SCY7b4,3202
+mapchete_eo/platforms/sentinel2/config.py,sha256=ifmzxDEuOFfoGAvbGyWyUWwF3Dw0NxEUEA3bWNByNR0,8482
+mapchete_eo/platforms/sentinel2/driver.py,sha256=hNc9jmLQMmnA2ciVBPSTmVW6H8clL6sydkxcy6BjO40,1166
+mapchete_eo/platforms/sentinel2/masks.py,sha256=6ig8sQhXkm1u6Bwbe7n8ewW8gTdbVJp-iaDCk780Qxg,11220
+mapchete_eo/platforms/sentinel2/preprocessing_tasks.py,sha256=DJhbvYhNDGZ2Er5polK2LXWcg-MASZMFparNOw56yDk,1645
+mapchete_eo/platforms/sentinel2/processing_baseline.py,sha256=uTE7zgzgi_W3NENLIAnJbWjyUD5_znpGa3cV-VoaEws,5179
+mapchete_eo/platforms/sentinel2/product.py,sha256=UUbngHI2qPri5GkMPdmaw-ZS0y6aIQUtgZ5oZXUcIJc,28095
+mapchete_eo/platforms/sentinel2/source.py,sha256=qSziCNTpInGy2LWIMiC3M4EP9G9_VtbCRNXB1XAO1P4,4269
+mapchete_eo/platforms/sentinel2/types.py,sha256=EV7TleKZnCNuDlEMJ0UtJjDsvOUGbvh_pZToeOULQ18,2122
+mapchete_eo/platforms/sentinel2/brdf/__init__.py,sha256=W7zAJZ5F5Xla7j4ua2opRANUbcD3jZbhTwhW5_cpOUI,242
+mapchete_eo/platforms/sentinel2/brdf/config.py,sha256=v3WCEu4r-tEPQgWBEkYNAPLgCQobvYtQ2YiDeAdImMk,1133
+mapchete_eo/platforms/sentinel2/brdf/correction.py,sha256=vi-BwAA6Cbayzk6WOuOESLhRkGLhiUAbVYfJUSfYhXA,8651
+mapchete_eo/platforms/sentinel2/brdf/hls.py,sha256=pzrRTNfJhuqRLidLQvKJpSnaz816x7jUesXt_aRgFjc,8818
+mapchete_eo/platforms/sentinel2/brdf/models.py,sha256=BBwGX1OLMRLncK1jD2reidaF82IKR4sjgNGHdTgBXaw,1446
+mapchete_eo/platforms/sentinel2/brdf/protocols.py,sha256=fxLDFeaL44eIPC2UwcbEQ-AXPREWgBv6Iytc74AQ5eA,752
+mapchete_eo/platforms/sentinel2/brdf/ross_thick.py,sha256=OvuuQxqWN8qrKpz8kDhmCd8R1uMLa8h5LISJSJ_NvYk,4685
+mapchete_eo/platforms/sentinel2/brdf/sun_angle_arrays.py,sha256=RxhQ8vBU0r5vsQ3FnyfDdYS0EGDGKCwAXND0Gn1m5mQ,2148
+mapchete_eo/platforms/sentinel2/metadata_parser/__init__.py,sha256=Vio87daAQ1eBQfDfA97c9byVnAUEPCidjTBg6jnZhXk,167
+mapchete_eo/platforms/sentinel2/metadata_parser/base.py,sha256=AvPfsQEAyXOrwOQ_MQgWLUFFJL0VhTzzWHpuSe14IS0,1520
+mapchete_eo/platforms/sentinel2/metadata_parser/default_path_mapper.py,sha256=Gt9jUZwO3u_p6iH3hkOGNS5Tad5yn-Y1F7uADtalBSo,5344
+mapchete_eo/platforms/sentinel2/metadata_parser/models.py,sha256=FOdEeIKk_Sa_O3iGRqoR-5YLQuwjPG6pNbdFRgQVYm8,2325
+mapchete_eo/platforms/sentinel2/metadata_parser/s2metadata.py,sha256=i6O9t0Bbb5BCiktKmXgDBKWdzdtVoGJ-Nt-xib-w3qw,23786
+mapchete_eo/platforms/sentinel2/preconfigured_sources/__init__.py,sha256=W-PFNuC0zXfF07uTUH_-r2Macv-_5qlz_YsP04Jxh8Y,1821
+mapchete_eo/platforms/sentinel2/preconfigured_sources/guessers.py,sha256=lkftO9NCJqwV663K3YqUFkNQxQWrCuwp_enbfTkjdbQ,3746
+mapchete_eo/platforms/sentinel2/preconfigured_sources/item_mappers.py,sha256=b_lvs6-rKQjb-8vBjOLV0Zi4a2ZXZxMqinyeRljA5lY,5957
+mapchete_eo/platforms/sentinel2/preconfigured_sources/metadata_xml_mappers.py,sha256=8UbF5Riox1FZIYGuDMBiGgxNfAQkOmu07ZMc2Eqiel0,8647
+mapchete_eo/processes/__init__.py,sha256=47DEQpj8HBSa-_TImW-5JCeuQeRkm5NMpJWZG3hSuFU,0
+mapchete_eo/processes/config.py,sha256=WorMzHJzJaZPh-aEyEEwcXN08ALxn6zAhVMrjkDUWDU,1745
+mapchete_eo/processes/dtype_scale.py,sha256=7hJlIe2DsE6Kmk1c3PLwLr8rnn_4_9S6Cz7bmRuzP9M,4176
+mapchete_eo/processes/eo_to_xarray.py,sha256=Gcp4qju2C9S8KeUnVY5f3nAsrdckPhGRguzsgRtWBFs,514
+mapchete_eo/processes/merge_rasters.py,sha256=jSXUI8pEWXbQWrM9iMtrH4ACH2VkCaBUH9CQg4z6zLA,8486
+mapchete_eo/search/__init__.py,sha256=bHvaqoHrtoNyO8_3kyJfvhh1yCz0hYvU2fzXhjOXNhU,528
+mapchete_eo/search/base.py,sha256=B8vlwgvDYRqxHre5aYWfPZ0tmUic-mb1VBnzl0llUSQ,9757
+mapchete_eo/search/config.py,sha256=L3DBEuKNGh_WXi8_Ha-qj6HW9ksRioYeFyPYGGN09qA,3696
+mapchete_eo/search/s2_mgrs.py,sha256=x12m-tJZNc41lFHNR2h7eSyEZ-ceBO4Gdz785btcVN8,10859
+mapchete_eo/search/stac_search.py,sha256=gWMSmDaXZTfH_GV1YiQ8ThbnJYsxHB2j_y_SgMX41Bg,10079
+mapchete_eo/search/stac_static.py,sha256=e53nKu8QRmanGvPmrvIDavHXm8v_4aT5V-wvU8wnkiw,6979
+mapchete_eo/search/utm_search.py,sha256=ZbNJCf461H_709Bpqzuwiu4_iwXC4M-wkemxXKmGenk,9489
+mapchete_eo-2026.1.0.dist-info/METADATA,sha256=0qagoc1eCiRIjxn78e_jJIieQkjV018I_0ylh-1leQQ,3262
+mapchete_eo-2026.1.0.dist-info/WHEEL,sha256=aha0VrrYvgDJ3Xxl3db_g_MDIW-ZexDdrc_m-Hk8YY4,105
+mapchete_eo-2026.1.0.dist-info/entry_points.txt,sha256=edU29CA3j8hwBafEd6jO2Zcou1sGkyjNro-y4_LnAnw,320
+mapchete_eo-2026.1.0.dist-info/licenses/LICENSE,sha256=gmZNf2EjyzaNbUyhfsBJlS0yQeG7GEmd-29miNXtEwA,1089
+mapchete_eo-2026.1.0.dist-info/RECORD,,

{mapchete_eo-2025.10.1.dist-info → mapchete_eo-2026.1.0.dist-info}/WHEEL

@@ -1,5 +1,5 @@
 Wheel-Version: 1.0
-Generator: hatchling 1.27.0
+Generator: hatchling 1.28.0
 Root-Is-Purelib: true
 Tag: py2-none-any
 Tag: py3-none-any

{mapchete_eo-2025.10.1.dist-info → mapchete_eo-2026.1.0.dist-info}/entry_points.txt

@@ -2,7 +2,7 @@
 eo = mapchete_eo.cli:eo
 
 [mapchete.formats.drivers]
-eostac_dev = mapchete_eo.eostac
+eostac = mapchete_eo.eostac
 sentinel2 = mapchete_eo.platforms.sentinel2
 
 [mapchete.processes]

{mapchete_eo-2025.10.1.dist-info → mapchete_eo-2026.1.0.dist-info}/licenses/LICENSE

@@ -1,6 +1,6 @@
 The MIT License (MIT)
 
-Copyright (c) 2022 - 2025 EOX IT Services
+Copyright (c) 2022 - 2026 EOX IT Services
 
 Permission is hereby granted, free of charge, to any person obtaining a copy
 of this software and associated documentation files (the "Software"), to deal

mapchete_eo/archives/base.py

@@ -1,65 +0,0 @@
-from abc import ABC
-import logging
-from typing import Any, Callable, Dict, Generator, List, Optional, Union
-
-from mapchete.io.vector import IndexedFeatures
-from mapchete.types import Bounds
-from pystac import Item
-from shapely.errors import GEOSException
-from shapely.geometry.base import BaseGeometry
-
-from mapchete_eo.exceptions import ItemGeometryError
-from mapchete_eo.search.base import CatalogSearcher
-from mapchete_eo.types import TimeRange
-
-logger = logging.getLogger(__name__)
-
-
-class Archive(ABC):
-    """
-    An archive combines a Catalog and a Storage.
-    """
-
-    time: Union[TimeRange, List[TimeRange]]
-    area: BaseGeometry
-    catalog: CatalogSearcher
-    search_kwargs: Dict[str, Any]
-    _items: Optional[IndexedFeatures] = None
-    item_modifier_funcs: Optional[List[Callable[[Item], Item]]] = None
-
-    def __init__(
-        self,
-        time: Union[TimeRange, List[TimeRange]],
-        bounds: Optional[Bounds] = None,
-        area: Optional[BaseGeometry] = None,
-        search_kwargs: Optional[Dict[str, Any]] = None,
-        catalog: Optional[CatalogSearcher] = None,
-    ):
-        if bounds is None and area is None:
-            raise ValueError("either bounds or area have to be provided")
-        elif area is None:
-            area = Bounds.from_inp(bounds).geometry
-        self.time = time
-        self.area = area
-        self.search_kwargs = search_kwargs or {}
-        if catalog:
-            self.catalog = catalog
-
-    def get_catalog_config(self):
-        return self.catalog.config_cls(**self.search_kwargs)
-
-    def apply_item_modifier_funcs(self, item: Item) -> Item:
-        try:
-            for modifier in self.item_modifier_funcs or []:
-                item = modifier(item)
-        except GEOSException as exc:
-            raise ItemGeometryError(
-                f"item {item.get_self_href()} geometry could not be resolved: {str(exc)}"
-            )
-        return item
-
-    def items(self) -> Generator[Item, None, None]:
-        for item in self.catalog.search(
-            time=self.time, area=self.area, search_kwargs=self.search_kwargs
-        ):
-            yield self.apply_item_modifier_funcs(item)

mapchete_eo/geometry.py

@@ -1,271 +0,0 @@
-import logging
-import math
-from functools import partial
-from typing import Callable, Iterable, Tuple
-
-from fiona.crs import CRS
-from fiona.transform import transform as fiona_transform
-from mapchete.geometry import reproject_geometry
-from mapchete.types import Bounds, CRSLike
-from shapely.geometry import (
-    GeometryCollection,
-    LinearRing,
-    LineString,
-    MultiLineString,
-    MultiPoint,
-    MultiPolygon,
-    Point,
-    Polygon,
-    box,
-    shape,
-)
-from shapely.geometry.base import BaseGeometry
-from shapely.ops import unary_union
-
-CoordArrays = Tuple[Iterable[float], Iterable[float]]
-
-
-logger = logging.getLogger(__name__)
-
-
-def transform_to_latlon(
-    geometry: BaseGeometry, src_crs: CRSLike, width_threshold: float = 180.0
-) -> BaseGeometry:
-    """Transforms a geometry to lat/lon coordinates.
-
-    If resulting geometry crosses the Antimeridian it will be fixed by moving coordinates
-    from the Western Hemisphere to outside of the lat/lon bounds on the East, making sure
-    the correct geometry shape is preserved.
-
-    As a next step, repair_antimeridian_geometry() can be applied, which then splits up
-    this geometry into a multipart geometry where all of its subgeometries are within the
-    lat/lon bounds again.
-    """
-    latlon_crs = CRS.from_epsg(4326)
-
-    def transform_shift_coords(coords: CoordArrays) -> CoordArrays:
-        out_x_coords, out_y_coords = fiona_transform(src_crs, latlon_crs, *coords)
-        if max(out_x_coords) - min(out_x_coords) > width_threshold:
-            # we probably have an antimeridian crossing here!
-            out_x_coords, out_y_coords = coords_longitudinal_shift(
-                coords_transform(coords, src_crs, latlon_crs), only_negative_coords=True
-            )
-        return (out_x_coords, out_y_coords)
-
-    return custom_transform(geometry, transform_shift_coords)
-
-
-def repair_antimeridian_geometry(
-    geometry: BaseGeometry, width_threshold: float = 180.0
-) -> BaseGeometry:
-    """
-    Repair geometry and apply fix if it crosses the Antimeridian.
-
-    A geometry crosses the Antimeridian if it is at least partly outside of the
-    lat/lon bounding box or if its width exceeds a certain threshold. This can happen
-    after reprojection if the geometry coordinates are transformed separately and land
-    left and right of the Antimeridian, thus resulting in a polygon spanning almost the
-    whole lat/lon bounding box width.
-    """
-    # repair geometry if it is broken
-    geometry = geometry.buffer(0)
-    latlon_bbox = box(-180, -90, 180, 90)
-
-    # only attempt to fix if geometry is too wide or reaches over the lat/lon bounds
-    if (
-        Bounds.from_inp(geometry).width >= width_threshold
-        or not geometry.difference(latlon_bbox).is_empty
-    ):
-        # (1) shift only coordinates on the western hemisphere by 360°, thus "fixing"
-        # the footprint, but letting it cross the antimeridian
-        shifted_geometry = longitudinal_shift(geometry, only_negative_coords=True)
-
-        # (2) split up geometry in one outside of latlon bounds and one inside
-        inside = shifted_geometry.intersection(latlon_bbox)
-        outside = shifted_geometry.difference(latlon_bbox)
-
-        # (3) shift back only the polygon outside of latlon bounds by -360, thus moving
-        # it back to the western hemisphere
-        outside_shifted = longitudinal_shift(
-            outside, offset=-360, only_negative_coords=False
-        )
-
-        # (4) create a MultiPolygon out from these two polygons
-        geometry = unary_union([inside, outside_shifted])
-
-    return geometry
-
-
-def buffer_antimeridian_safe(
-    footprint: BaseGeometry, buffer_m: float = 0
-) -> BaseGeometry:
-    """Buffer geometry by meters and make it Antimeridian-safe.
-
-    Safe means that if it crosses the Antimeridian and is a MultiPolygon,
-    the buffer will only be applied to the edges facing away from the Antimeridian
-    thus leaving the polygon intact if shifted back.
-    """
-    if footprint.is_empty:
-        return footprint
-
-    # repair geometry if it is broken
-    footprint = footprint.buffer(0)
-
-    if not buffer_m:
-        return footprint
-
-    if isinstance(footprint, MultiPolygon):
-        # we have a shifted footprint here!
-        # (1) unshift one part
-        subpolygons = []
-        for polygon in footprint.geoms:
-            lon = polygon.centroid.x
-            if lon < 0:
-                polygon = longitudinal_shift(polygon)
-            subpolygons.append(polygon)
-        # (2) merge to single polygon
-        merged = unary_union(subpolygons)
-
-        # (3) apply buffer
-        if isinstance(merged, MultiPolygon):
-            buffered = unary_union(
-                [
-                    buffer_antimeridian_safe(polygon, buffer_m=buffer_m)
-                    for polygon in merged.geoms
-                ]
-            )
-        else:
-            buffered = buffer_antimeridian_safe(merged, buffer_m=buffer_m)
-
-        # (4) fix again
-        return repair_antimeridian_geometry(buffered)
-
-    # UTM zone CRS
-    utm_crs = latlon_to_utm_crs(footprint.centroid.y, footprint.centroid.x)
-    latlon_crs = CRS.from_string("EPSG:4326")
-
-    return transform_to_latlon(
-        reproject_geometry(
-            footprint, src_crs=latlon_crs, dst_crs=utm_crs, clip_to_crs_bounds=False
-        ).buffer(buffer_m),
-        src_crs=utm_crs,
-    )
-
-
-def longitudinal_shift(
-    geometry: BaseGeometry, offset: float = 360.0, only_negative_coords: bool = False
-) -> BaseGeometry:
-    """Return geometry with either all or Western hemisphere coordinates shifted by some offset."""
-    return custom_transform(
-        geometry,
-        partial(
-            coords_longitudinal_shift,
-            by=offset,
-            only_negative_coords=only_negative_coords,
-        ),
-    )
-
-
-def latlon_to_utm_crs(lat: float, lon: float) -> CRS:
-    min_zone = 1
-    max_zone = 60
-    utm_zone = (
-        f"{max([min([(math.floor((lon + 180) / 6) + 1), max_zone]), min_zone]):02}"
-    )
-    hemisphere_code = "7" if lat <= 0 else "6"
-    return CRS.from_string(f"EPSG:32{hemisphere_code}{utm_zone}")
-
-
-def bounds_to_geom(bounds: Bounds) -> BaseGeometry:
-    # TODO: move into core package
-    if bounds.left < -180:
-        part1 = Bounds(-180, bounds.bottom, bounds.right, bounds.top)
-        part2 = Bounds(bounds.left + 360, bounds.bottom, 180, bounds.top)
-        return unary_union([shape(part1), shape(part2)])
-    elif bounds.right > 180:
-        part1 = Bounds(-180, bounds.bottom, bounds.right - 360, bounds.top)
-        part2 = Bounds(bounds.left, bounds.bottom, 180, bounds.top)
-        return unary_union([shape(part1), shape(part2)])
-    else:
-        return shape(bounds)
-
-
-def custom_transform(geometry: BaseGeometry, func: Callable) -> BaseGeometry:
-    # todo: shapely.transform.transform maybe can make this code more simple
-    # https://shapely.readthedocs.io/en/stable/reference/shapely.transform.html#shapely.transform
-    def _point(point: Point) -> Point:
-        return Point(zip(*func(point.xy)))
-
-    def _multipoint(multipoint: MultiPoint) -> MultiPoint:
-        return MultiPoint([_point(point) for point in multipoint])
-
-    def _linestring(linestring: LineString) -> LineString:
-        return LineString(zip(*func(linestring.xy)))
-
-    def _multilinestring(multilinestring: MultiLineString) -> MultiLineString:
-        return MultiLineString(
-            [_linestring(linestring) for linestring in multilinestring.geoms]
-        )
-
-    def _linearring(linearring: LinearRing) -> LinearRing:
-        return LinearRing(((x, y) for x, y in zip(*func(linearring.xy))))
-
-    def _polygon(polygon: Polygon) -> Polygon:
-        return Polygon(
-            _linearring(polygon.exterior),
-            holes=list(map(_linearring, polygon.interiors)),
-        )
-
-    def _multipolygon(multipolygon: MultiPolygon) -> MultiPolygon:
-        return MultiPolygon([_polygon(polygon) for polygon in multipolygon.geoms])
-
-    def _geometrycollection(
-        geometrycollection: GeometryCollection,
-    ) -> GeometryCollection:
-        return GeometryCollection(
-            [_any_geometry(subgeometry) for subgeometry in geometrycollection.geoms]
-        )
-
-    def _any_geometry(geometry: BaseGeometry) -> BaseGeometry:
-        transform_funcs = {
-            Point: _point,
-            MultiPoint: _multipoint,
-            LineString: _linestring,
-            MultiLineString: _multilinestring,
-            Polygon: _polygon,
-            MultiPolygon: _multipolygon,
-            GeometryCollection: _geometrycollection,
-        }
-        try:
-            return transform_funcs[type(geometry)](geometry)
-        except KeyError:
-            raise TypeError(f"unknown geometry {geometry} of type {type(geometry)}")
-
-    if geometry.is_empty:
-        return geometry
-
-    # make valid by buffering
-    return _any_geometry(geometry).buffer(0)
-
-
-def coords_transform(
-    coords: CoordArrays, src_crs: CRSLike, dst_crs: CRSLike
-) -> CoordArrays:
-    return fiona_transform(src_crs, dst_crs, *coords)
-
-
-def coords_longitudinal_shift(
-    coords: CoordArrays,
-    by: float = 360,
-    only_negative_coords: bool = False,
-) -> CoordArrays:
-    x_coords, y_coords = coords
-    x_coords = (
-        (
-            x_coord + by
-            if (only_negative_coords and x_coord < 0) or not only_negative_coords
-            else x_coord
-        )
-        for x_coord in x_coords
-    )
-    return x_coords, y_coords

mapchete_eo/known_catalogs.py

@@ -1,42 +0,0 @@
-"""
-Catalogs define access to a search interface which provide products
-as pystac Items.
-"""
-
-from typing import List
-
-from mapchete_eo.search import STACSearchCatalog, UTMSearchCatalog
-
-
-class EarthSearchV1S2L2A(STACSearchCatalog):
-    """Earth-Search catalog for Sentinel-2 Level 2A COGs."""
-
-    endpoint: str = "https://earth-search.aws.element84.com/v1/"
-
-
-class CDSESearch(STACSearchCatalog):
-    """Copernicus Data Space Ecosystem (CDSE) STAC API."""
-
-    endpoint: str = "https://stac.dataspace.copernicus.eu/v1"
-
-
-class PlanetaryComputerSearch(STACSearchCatalog):
-    """Planetary Computer Search."""
-
-    endpoint: str = "https://planetarycomputer.microsoft.com/api/stac/v1/"
-
-
-class AWSSearchCatalogS2L2A(UTMSearchCatalog):
-    """
-    Not a search endpoint, just hanging STAC collection with items separately.
-    Need custom parser/browser to find scenes based on date and UTM MGRS Granule
-
-    https://sentinel-s2-l2a-stac.s3.amazonaws.com/sentinel-s2-l2a.json
-    """
-
-    id: str = "sentinel-s2-l2a"
-    endpoint: str = "s3://sentinel-s2-l2a-stac/"
-    day_subdir_schema: str = "{year}/{month:02d}/{day:02d}"
-    stac_json_endswith: str = "T{tile_id}.json"
-    description: str = "Sentinel-2 L2A JPEG2000 archive on AWS."
-    stac_extensions: List[str] = []