Skip to contents

A GTI index is a vector layer of raster footprints read by GDAL's GTI (GDAL Tile Index) raster driver, available from GDAL 3.9, which mosaics the indexed rasters on the fly. Footprints are stored in crs (transformed from the table's EPSG:4326 bboxes), so the index layer SRS matches the grid the GTI dataset will be pinned to. Each footprint is padded by 1% of its width and height on every side: GDAL's GTI driver (3.12.3 onwards) silently drops a reprojected tile whose warped extent falls outside its footprint, and both STAC bboxes and exactly transformed bounds can fall slightly short of that extent. The footprint only selects tiles, so the padding costs at most a tile open for windows just beyond a tile's edge. Most users reach this through lazy_dataset(), which builds the index internally.

Usage

stac_gti_index(
  sources,
  asset,
  path = tempfile(fileext = ".gti.fgb"),
  crs = "EPSG:4326"
)

Arguments

sources

A stac_sources() table with a slice column (see stac_time_slices()).

asset

Which asset's rows to index.

path

Index path; defaults to a tempfile.

crs

CRS for the index footprints (use the target grid's CRS).

Value

The index path, invisibly.