to_polygons(labels) turns a label containing a hole into an invalid MultiPolygon with the hole counted as area
Maintainer antworten meist innerhalb von 1 Tag
Dieses Issue hat noch niemand übernommen.
Bewertung
- Schwierigkeit
- 3/5
- Geschätzter Aufwand
- 1-2 Tage
- Anfängerfreundlichkeit
- 64/100
Rechercherichtung
Beginne mit src/spatialdata/_core/operations/vectorize.py::_region_props_to_polygon und führe die bereitgestellte repro.py mit uv aus. Vergleiche die resultierende Geometrie mit dem erwarteten gültigen Polygon, das einen inneren Ring und die Fläche 24.0 enthält, und füge einen Regressionstest für den Fall eines beschrifteten Rings hinzu, wenn sich der Testort im Repository ermitteln lässt.
Vom Indexierungsmodell aus dem Issue-Text verfasst.
Beschreibung
[!NOTE]
This whole message is AI-generated. The issue was automatically discovered and reported by an AI agent (Claude) during an autonomous bug hunt on thespatialdatacode base. It has not been verified or triaged by a human yet; theneeds: triagelabel is set so that a maintainer can confirm it. The reproduction script below was executed by the agent in an isolated environment (see Environment) and its output is pasted verbatim.
Summary
A 5×5 label with a 1-pixel hole (true area 24) becomes MultiPolygon, is_valid=False, area 25.0. Inner contours should become interiors of the enclosing polygon.
Severity (agent's assessment): medium — invalid geometries break downstream sjoin/overlay/clip; 24 of 3,556 polygons are invalid on the CosMx sandbox dataset
Where: src/spatialdata/_core/operations/vectorize.py::_region_props_to_polygon (every contour from skimage.measure.find_contours becomes an independent Polygon wrapped in a MultiPolygon)
Expected behaviour
A valid Polygon with one interior ring and area 24.
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",
# ]
# ///
"""to_polygons(labels) turns a label with a hole into an invalid MultiPolygon whose area includes the hole."""
import warnings
import numpy as np
from spatialdata import to_polygons
from spatialdata.models import Labels2DModel
warnings.simplefilter("ignore")
arr = np.zeros((7, 7), dtype=np.uint16)
arr[1:6, 1:6] = 1
arr[3, 3] = 0 # 5x5 square with a one-pixel hole -> true area 24
polys = to_polygons(Labels2DModel.parse(arr))
g = polys.geometry.iloc[0]
print(f"geometry type={g.geom_type}, is_valid={g.is_valid}, area={g.area} (expected: valid Polygon with one interior, area 24.0)")
bug = (not g.is_valid) or abs(g.area - 24.0) > 1e-9
print("VERDICT:", "BUG REPRODUCED" if bug else "NOT REPRODUCED")
Observed output
geometry type=MultiPolygon, is_valid=False, area=25.0 (expected: valid Polygon with one interior, area 24.0)
VERDICT: BUG REPRODUCED
Possible fix direction (unverified)
Build polygons with shapely.polygonize/build_area from the contour rings, or classify contours by containment (Polygon(outer, holes=[...])); add make_valid() as a safety net and a test with a labelled ring.
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.
- Vorherrschende Sprache
- Python
- Sterne
- 404
- Forks
- 106
- Ø Merge
- 3 T. 14 Std.
- Gemergte PRs (30 T.)
- 3
Entwicklungsumgebung
- Enthält ein Dockerfile oder eine Docker-Compose-Datei
- Keine Pull-Request-Vorlage
- Beitragsleitfaden lesen
Erste Schritte
- Lesen Sie das ganze Issue und danach den Beitragsleitfaden des Projekts.
- Schreiben Sie ins Issue, dass Sie es übernehmen — das erspart doppelte Arbeit.
- Forken Sie das Repository und arbeiten Sie in einem Branch.
- Öffnen Sie einen Pull Request, der die Issue-Nummer nennt.
Mehr aus scverse/spatialdata
-
Schwierigkeit 2/5 1-3 Stunden Anfängerfreundlichkeit 67/100
scverse/spatialdata#1293 · 1 Kommentar ·
Maintainer antworten meist innerhalb von 1 Tag
-
SpatialData 0.8.0 permits Zarr versions that do not provide ArrayNotFoundErrorEvtl. vergeben @Darlokt hat das vor 38 Tagen übernommen. Offen
Schwierigkeit 2/5 1-3 Stunden Anfängerfreundlichkeit 74/100
scverse/spatialdata#1256 ·
Maintainer antworten meist innerhalb von 1 Tag
-
bug 🚨 element: labels 🏷️ method: aggregation 🔢 needs: triage priority: medium
Schwierigkeit 2/5 1-3 Stunden Anfängerfreundlichkeit 68/100
scverse/spatialdata#1249 ·
Maintainer antworten meist innerhalb von 1 Tag
-
bug 🚨 element: images 🌌 element: labels 🏷️ needs: triage
Schwierigkeit 2/5 1-3 Stunden Anfängerfreundlichkeit 78/100
scverse/spatialdata#1239 ·
Maintainer antworten meist innerhalb von 1 Tag
-
bug 🚨 element: shapes ▲ models needs: triage priority: medium
Schwierigkeit 2/5 1-3 Stunden Anfängerfreundlichkeit 78/100
scverse/spatialdata#1234 ·
Maintainer antworten meist innerhalb von 1 Tag
Alle Issues in scverse/spatialdata
Ähnliche Issues
-
Schwierigkeit 2/5 1-3 Stunden Anfängerfreundlichkeit 68/100
-
Schwierigkeit 2/5 1-3 Stunden Anfängerfreundlichkeit 70/100
Maintainer antworten meist innerhalb von 2 Tagen
-
Schwierigkeit 2/5 1-3 Stunden Anfängerfreundlichkeit 62/100
Maintainer antworten meist innerhalb von 2 Tagen
-
bug
Schwierigkeit 2/5 1-3 Stunden Anfängerfreundlichkeit 82/100
developmentseed/titiler-cmr#278 ·
Maintainer antworten meist innerhalb von 1 Tag
-
Schwierigkeit 2/5 1-3 Stunden Anfängerfreundlichkeit 82/100
Maintainer antworten meist innerhalb von 2 Tagen