Multi-box bounding_box_query works on elements but crashes on a SpatialData object (TypeError: Unsupported type list) #1238

Summary

Element-level multi-box queries return a list; the SpatialData overload passes those lists into SpatialData(points={...}) → TypeError: Unsupported type <class 'list'>.

Severity (agent's assessment): low/medium — the docstring advertises min_coordinate of shape (n_boxes, n_axes)

Where: src/spatialdata/_core/query/spatial_query.py, SpatialData overload of bounding_box_query / _dict_query_dispatcher (element results that are lists are fed to the SpatialData constructor)

Expected behaviour

A list of SpatialData objects (one per box) or a clear ValueError.

Reproduction

Save as repro.py and run uv run repro.py (the PEP 723 header pins spatialdata to the commit the bug was found on; replace the URL fragment with @main to test the current main branch).

# /// script
# requires-python = ">=3.12"
# dependencies = [
#     "spatialdata @ git+https://github.com/scverse/spatialdata.git@ccf1ea048d054b6624214bf618008a9f9ae223e0",
# ]
# ///
"""Multi-box bounding_box_query works on elements but crashes on a SpatialData object."""
import warnings
import numpy as np
import pandas as pd
from spatialdata import SpatialData, bounding_box_query
from spatialdata.models import PointsModel

warnings.simplefilter("ignore")
points = PointsModel.parse(pd.DataFrame({"x": [1.0, 5.0, 9.0], "y": [1.0, 5.0, 9.0]}, index=[100, 200, 300]))
boxes_min, boxes_max = np.array([[0, 0], [4, 4]]), np.array([[3, 3], [8, 8]])
res = bounding_box_query(points, axes=("x", "y"), min_coordinate=boxes_min, max_coordinate=boxes_max, target_coordinate_system="global")
print("element-level multi-box query ->", [r.index.compute().tolist() for r in res])
bug = False
try:
    bounding_box_query(SpatialData(points={"pts": points}), axes=("x", "y"), min_coordinate=boxes_min, max_coordinate=boxes_max, target_coordinate_system="global")
    print("SpatialData-level multi-box query -> OK")
except Exception as e:  # noqa: BLE001
    print("SpatialData-level multi-box query ->", type(e).__name__, e)
    bug = True
print("VERDICT:", "BUG REPRODUCED" if bug else "NOT REPRODUCED")

Observed output

element-level multi-box query -> [[100], [200]]
SpatialData-level multi-box query -> TypeError Unsupported type <class 'list'>
VERDICT: BUG REPRODUCED

Possible fix direction (unverified)

Either build one SpatialData per box or raise early when min_coordinate.ndim == 2.

Environment

uv run repro.py with the PEP 723 metadata in the script (fresh, isolated environment; spatialdata built from main @ ccf1ea0 (2026-08-28); Python 3.13, latest releases of the dependencies at run time: pandas 3.0, anndata 0.13, zarr 3.3, dask 2026.8, numpy 2.5, geopandas 1.1, shapely 2.1). macOS (arm64). Also reproduced in a second environment with pandas 2.3.3 / anndata 0.12.11 / numpy 2.4.4 / zarr 3.2.1.

Automatically generated; discovered by an AI agent (Claude) and not yet reviewed by a human.