
Stack aligned rasters along a new outer dim (default time).
lazy_stack.RdAll layers must share the spatial grid (align() first otherwise);
dtypes promote to a common type. Chunks carry the stack as
(t, y, x) arrays; temporal reductions
(reduce_over(x, "median", "t")) then run chunk-locally.
Details
Element names of xs become labels on the new axis (slice dates on
t, band names on band), used by time_sel() / band_sel() for
label selection; an unnamed list leaves the axis unlabelled.