Integration with `pystac_client`

[robbibt]

Hi all, this looks really fantastic! I'd love to try this out, but I'm unfortunately being blocked by rustac issues:

     14 # Search a STAC API and cache results to a local stac-geoparquet file.
---> 15 await rustac.search_to(
     16     "items.parquet",
     17     "https://explorer.dea.ga.gov.au/stac",
     18     collections=["ga_s2ls_intertidal_cyear_3"],

RustacError: builder error

Is it possible to run lazycogs directly on the list of items returned by pystac_client instead? For example, the output of something like this:

client = pystac_client.Client.open("https://explorer.dea.ga.gov.au/stac")

search = client.search(
    collections=["ga_ls8c_ard_3"],
    bbox=[150.15, -35.75, 150.25, -35.65], 
    datetime="2023-01-01/2023-01-31",
)

# Retrieve all matched items as a pystac.ItemCollection
items = search.item_collection()

[hrodmn]

@robbibt thanks for reaching out with your issue! lazycogs depends pretty heavily on rustac's ability to search through stac-geoparquet files. Can you share your runtime environment? I am curious about why rustac is not working for you.

[robbibt]

Hey @hrodmn, thanks for the quick reply! I'm currently installing lazycogs via uv, based on this pyproject.toml:

[project]
name = "lazycogs-dev"
version = "0.1.0"
description = "Add your description here"
readme = "README.md"
requires-python = ">=3.12"
dependencies = [
    "ipykernel>=7.2.0",
    "lazycogs>=0.1.1",
]

I've uploaded a uv.lock here - I've been testing in Python 3.13 so far: https://gist.github.com/robbibt/48a84fd723356f158520013bcf1e3e53

uv pip list:

Package                 Version
----------------------- -----------
affine                  2.4.0
arro3-core              0.8.0
asttokens               3.0.1
async-geotiff           0.4.0
async-tiff              0.7.1
certifi                 2026.4.22
comm                    0.2.3
cql2                    0.5.5
debugpy                 1.8.20
decorator               5.2.1
defusedxml              0.7.1
executing               2.2.1
ipykernel               7.2.0
ipython                 9.13.0
ipython-pygments-lexers 1.1.1
jedi                    0.19.2
jupyter-client          8.8.0
jupyter-core            5.9.1
lazycogs                0.1.2
matplotlib-inline       0.2.1
nest-asyncio            1.6.0
numpy                   2.4.4
obspec                  0.1.0
obstore                 0.9.4
packaging               26.2
pandas                  3.0.2
parso                   0.8.6
pexpect                 4.9.0
platformdirs            4.9.6
prompt-toolkit          3.0.52
psutil                  7.2.2
ptyprocess              0.7.0
pure-eval               0.2.3
pygments                2.20.0
pyproj                  3.7.2
python-dateutil         2.9.0.post0
pyzmq                   27.1.0
rustac                  0.9.10
six                     1.17.0
stack-data              0.6.3
tornado                 6.5.5
traitlets               5.14.3
wcwidth                 0.6.0
xarray                  2026.4.0

The only change I made to the Example in the readme was changing alwaysxy=True (which returned an error) to always_xy=True.

[hrodmn]

What operating system are you using?

[robbibt]

I'm running it via Jupyter Hub in a Linux system:

jovyan@jupyter-robbi-2ebishoptaylor-40ga-2egov-2eau:~$ cat /etc/lsb-release
DISTRIB_ID=Ubuntu
DISTRIB_RELEASE=22.04
DISTRIB_CODENAME=jammy
DISTRIB_DESCRIPTION="Ubuntu 22.04.3 LTS"

[hrodmn]

@gadomski have you seen RustacError: builder error before?

[gadomski]

IIRC it often has to do with building the SQL query. @robbibt can you share your whole script so I can try to reproduce?

If you turn up the logging (e.g. w/ debug) you might get some meaningful output, too.