Fun fact: the EXPLAIN VERBOSE of reading Overture buildings is > 65536 characters wide (because this will list all the file paths that make up the dataset and there are hundreds of those...)
import sedona.db
sd = sedona.db.connect()
buildings_url = (
"s3://overturemaps-us-west-2/release/2025-11-19.0/theme=buildings/type=building/"
)
target_wkt = (
"POLYGON ((-73.21 44.03, -73.21 43.98, -73.11 43.97, -73.12 44.03, -73.21 44.03))"
)
sd.read_parquet(
buildings_url,
options={"aws.skip_signature": True, "aws.region": "us-west-2"},
).to_view("buildings")
sd.sql(f"""
SELECT count(*) FROM buildings
WHERE ST_Intersects(geometry, ST_SetSRID(ST_GeomFromText('{target_wkt}'), 4326))
""").explain("extended").show()
#> thread '<unnamed>' panicked at rust/sedona/src/show.rs:425:46:
#> PanicException: attempt to add with overflow