Remove all wrapping of stackstac and odc-stac

Author: jsignell

README.md

@@ -96,23 +96,6 @@ xr.open_dataset(asset, patch_url=planetary_computer.sign)
 
 Note that this zarr asset uses the xarray-assets extension to store `open_kwargs` and `storage_options` which xpystac can then pass along to `xr.open_dataset`.
 
-### Open a single item
-
-A single item containing many COGs:
-
-```python
-import pystac
-import xarray as xr
-
-
-item = pystac.Item.from_file(
-    "https://earth-search.aws.element84.com/v1/collections/landsat-c2-l2/items/LC09_L2SR_081108_20250311_02_T2"
-)
-
-xr.open_dataset(item)
-```
-This takes advantage of a stacking library (either
-[odc-stac](https://github.com/opendatacube/odc-stac) or [stackstac](https://github.com/gjoseph92/stackstac) - configurable via the `stacking_library` option)
 
 ### Open many items
 
@@ -152,9 +135,8 @@ xr.open_dataset(item_collection)
 ## How it works
 
 When you call ``xarray.open_dataset(object, engine="stac")`` this library maps that `open` call to the correct library.
-Depending on the ``type`` of ``object`` that might be a stacking library (either
-[odc-stac](https://github.com/opendatacube/odc-stac) or [stackstac](https://github.com/gjoseph92/stackstac))
-or back to ``xarray.open_dataset`` itself but with the engine and other options pulled from the pystac object.
+Depending on the ``type`` of ``object`` that might be
+back to ``xarray.open_dataset`` itself but with the engine and other options pulled from the pystac object.
 
 ## Prior Art
 

environment.yaml

@@ -28,20 +28,21 @@ Repository = "https://github.com/stac-utils/xpystac"
 dev = [
     "adlfs>=2024.12.0",
     "aiohttp>=3.12.15",
+    "dask>=2025.9.1",
     "fsspec>=2025.3.0",
     "icechunk>=1.1.5",
+    "ipdb>=0.13.13",
     "ipython>=9.5.0",
     "kerchunk>=0.2.7",
-    "odc-stac>=0.3.11",
     "planetary-computer>=1.0.0",
     "pre-commit>=4.0.1",
     "pystac<1.12",
     "pystac-client>=0.8.6",
     "pytest>=8.3.3",
     "pytest-cov>=6.0.0",
+    "pytest-xdist>=3.8.0",
     "requests>=2.32.3",
     "rioxarray>=0.18.2",
-    "stackstac>=0.5.1",
     "urllib3>=2.3.0",
     "zarr>2",
 ]

pyproject.toml

@@ -1,7 +1,5 @@
-import dask.array
 import pystac_client
 import pytest
-import xarray as xr
 
 from tests.utils import STAC_URLS, requires_icechunk, requires_planetary_computer
 from xpystac.core import to_xarray
@@ -12,30 +10,6 @@ def test_to_xarray_with_cog_asset(simple_cog):
     assert ds
 
 
-def test_to_xarray_with_pystac_client_search(simple_search):
-    ds = to_xarray(simple_search)
-    assert ds
-
-
-def test_to_xarray_returns_dask_backed_object(simple_search):
-    ds = to_xarray(simple_search)
-    assert isinstance(ds.blue.data, dask.array.Array)
-    assert ds.blue.data.npartitions > 1
-
-
-def test_to_xarray_with_pystac_client_search_passes_kwargs_through(simple_search):
-    ds = to_xarray(simple_search, bands=["red", "green", "blue"], chunks={})
-    assert list(ds.data_vars) == ["red", "green", "blue"]
-    assert ds.blue.data.npartitions == 1
-
-
-@pytest.mark.parametrize("stacking_library", ["odc.stac", "stackstac"])
-def test_to_xarray_with_different_stacking_library(simple_search, stacking_library):
-    ds = to_xarray(simple_search, stacking_library=stacking_library)
-    assert isinstance(ds, xr.Dataset)
-    assert "band" not in ds.dims
-
-
 @requires_planetary_computer
 def test_to_xarray_with_pystac_client_search_with_patch_url():
     import planetary_computer as pc
@@ -47,24 +21,13 @@ def test_to_xarray_with_pystac_client_search_with_patch_url():
         collections=["sentinel-2-l2a"],
         datetime="2020-05-01",
     )
-
-    ds = to_xarray(search, assets=["B4", "B3", "B2"], stacking_library="odc.stac")
+    item = next(search.items())
+    asset = item.assets["B04"]
 
     with pytest.raises(RasterioIOError, match="HTTP response code: 409"):
-        ds.B01.max().compute()
+        to_xarray(asset)
 
-    ds = to_xarray(
-        search,
-        assets=["B4", "B3", "B2"],
-        stacking_library="odc.stac",
-        patch_url=pc.sign,
-    )
-    assert ds.B01.max().compute() == 11080
-
-
-def test_to_xarray_with_drop_variables_raises(simple_search):
-    with pytest.raises(KeyError, match="not implemented for pystac items"):
-        to_xarray(simple_search, drop_variables=["blue"])
+    to_xarray(asset, patch_url=pc.sign)
 
 
 def test_to_xarray_with_bad_type():
@@ -145,20 +108,6 @@ def test_to_xarray_zarr_with_zarr_extension():
     to_xarray(zarr_asset)
 
 
-@pytest.mark.skip(reason="not yet supported with kerchunk >=0.2.8")
-def test_to_xarray_with_item_collection_with_kerchunk_attrs_in_data_cube(
-    data_cube_kerchunk,
-):
-    ds = to_xarray(data_cube_kerchunk)
-    assert ds
-
-
-@pytest.mark.skip(reason="not yet supported with kerchunk >=0.2.8")
-def test_to_xarray_with_list_with_kerchunk_attrs_in_data_cube(data_cube_kerchunk):
-    ds = to_xarray([i for i in data_cube_kerchunk])
-    assert ds
-
-
 @pytest.mark.skip(reason="not yet supported with kerchunk >=0.2.8")
 def test_to_xarray_with_item_with_kerchunk_attrs_in_data_cube(data_cube_kerchunk):
     ds = to_xarray([i for i in data_cube_kerchunk][-1])

tests/test_core.py

@@ -1,17 +1,6 @@
-import pytest
 import xarray
 
 
 def test_xarray_open_dataset_can_guess_for_pystac_objects(simple_cog):
     ds = xarray.open_dataset(simple_cog)
     assert ds
-
-
-def test_xarray_open_dataset_can_guess_for_pystac_client_searchs(simple_search):
-    ds = xarray.open_dataset(simple_search, assets=["blue", "green", "red"])
-    assert ds
-
-
-def test_xarray_open_dataset_with_drop_variables_raises(simple_search):
-    with pytest.raises(KeyError, match="not implemented for pystac items"):
-        xarray.open_dataset(simple_search, engine="stac", drop_variables=["B0"])