Skip to contents

Lazy time-sliced stack of one STAC asset on a target grid.

lazy_stac_stack.Rd

Builds a GTI (GDAL Tile Index; see stac_gti_index()) index for asset, then opens one mosaic per time slice pinned to grid (mixed source CRS is fine: the GTI driver reprojects per tile) and stacks them along t. Overlaps within a slice resolve by ascending sort_field (highest drawn on top).

Usage

lazy_stac_stack(
  sources,
  grid,
  asset,
  granularity = "day",
  sort_field = "datetime",
  nodata = NULL,
  lon = NULL,
  scale = FALSE,
  offset = NULL
)

Arguments

sources

A stac_sources() table.

grid

Target GridSpec() for every slice.

asset

Asset name to stack.

granularity

Slice granularity (see stac_time_slices()).

sort_field

Index field ordering overlaps within a slice.

nodata

Optional nodata override passed to each slice source.

lon

Longitude for granularity = "solar_day" (see stac_time_slices()).

scale, offset

Read affine, as in lazy_source(): FALSE (default) reads raw values, TRUE discovers the file's band scale/offset (probing the mosaic, then the first item), a numeric supplies it explicitly.

Value

A list: stack (LazyRaster), slices (character), index (path).

See also

collect() to materialise the stack; lazy_dataset(), the higher-level multi-band interface most users want.

Other stac helpers: stac_drop_duplicates(), stac_filter_assets(), stac_filter_cloud(), stac_filter_coverage(), stac_filter_orbit(), stac_gti_index(), stac_merge(), stac_query(), stac_rename_assets(), stac_sign_mpc(), stac_sources(), stac_time_slices()