Quickstart#
Five lines of Python — every NEXRAD KLOT (Chicago, IL) scan ever published to the public archive, in one xarray.DataTree. No downloads, no decoders, no waiting.
import xarray as xr
import icechunk
storage = icechunk.s3_storage(
bucket="nexrad-arco",
prefix="KLOT",
region="us-east-1",
anonymous=True,
)
session = icechunk.Repository.open(storage).readonly_session("main")
dt = xr.open_datatree(session.store, engine="rustytree", chunks=None)
print(sorted(dt.children)) # → ['VCP-12', 'VCP-212', 'VCP-34', ...]
That’s it. dt is a hierarchical, time-indexed view of the entire KLOT archive on the AWS Open Data Registry, grouped by Volume Coverage Pattern and sweep. Every variable is lazy — only the Zarr v3 chunks you slice into are fetched from S3.
What’s next#
Notebook 1: NEXRAD KLOT Demo — the full walkthrough, including a 2×2 polarimetric visualization (Z, ZDR, RhoHV, PhiDP).
Tutorials — paper reproduction and large-scale rainfall accumulation.
Installation — if
import icechunkfailed, start here.
Note
engine="rustytree" is a Rust-backed xarray.DataTree backend (rustytree-xarray on PyPI) recommended for radar-datatree archives. Drop-in replacement for engine="zarr", ~10× faster on icechunk repos served from object storage.