libera_utils.footprint_matching.readers.base.StaticReader#
- class libera_utils.footprint_matching.readers.base.StaticReader(file_paths: Sequence[Path])#
Bases:
AncillaryReaderA reader whose source is one body of data for the staged period.
Covers the
DAILYsources (one composite or map per period), and aNEAREST_PASSsource small enough to hold whole. The subclass loads the whole source once in_load_native(); this class indexes it (a point cloud into 1 degree bins) and serves every query from memory.Methods
preload()Load and index the whole source now, rather than on the first query.
Return every stacked product variable this reader contributes.
Return every variable this reader contributes to the product definition.
samples_in(bbox, time)Return the native samples inside
bboxthat matchtime.Attributes
The staged granule paths passed at construction, in the order given.
Number of read variables (the first axis of every
valuesarray).The one granule of a reader that reads exactly one per period.
- ADDITIONAL_PRODUCT_VARIABLES: tuple[VariableSpec, ...] = ()#
- MAX_TIME_OFFSET: timedelta64 | None = None#
- PASS_WINDOW: timedelta64 | None = None#
- PRODUCT_STACKS: tuple[ProductStack, ...] = ()#
- _box_samples(native: GridSamples | _PointIndex | _BoxNative, bbox: BoundingBox) tuple[ndarray, ndarray, ndarray, ndarray | None]#
Return
(lats, lons, values, times)of the native samples inside one plain box.- Raises:
TypeError – If
nativeis none of the three native forms.
- classmethod _check_product_stacks() None#
Reject a
PRODUCT_STACKSdeclaration the product could not be built from.- Raises:
TypeError – If a stack names a spec the reader does not produce, a spec is stacked twice, a stack takes the name of a product spec, a stack’s members have standard-deviation companions for some members but not others, or a stack’s members do not all have the same
modes.
- _index_points(native: PointSamples) _PointIndex#
Validate a loaded point cloud against the reader’s contract and index it.
- abstractmethod _load_native() GridSamples | PointSamples | _BoxNative#
Open the staged granule(s) and return the whole source as native samples.
Called once per reader instance. Implementers handle the source’s own fill and valid-range conventions, returning NaN where the source has no value, and return longitudes in [-180, 180].
- Returns:
GridSamplesfor a regular latitude/longitude grid, a_BoxNativefor a source held in a projected form that geolocates the block a query selects,PointSamplesfor anything else (withtimesfor aNEAREST_PASSreader).- Return type:
GridSamples or PointSamples or _BoxNative
- property _native: GridSamples | _PointIndex | _BoxNative#
The loaded, indexed native data (loaded on first access, then held).
- Raises:
ValueError – If the loaded data carries a different number of variables than
VARIABLES.TypeError – If
_load_native()returned none of the three native forms.
- _native_for(time: datetime64) tuple[GridSamples | _PointIndex | _BoxNative, ...]#
Return the native data a query at
timedraws from; empty for no match.StaticReaderreturns its one indexed body of data regardless oftime;SteppedGridReaderthe step nearesttime, or nothing when that step is farther thanMAX_TIME_OFFSET.
- _nearest_pass_mask(times: ndarray, time: datetime64, inside: ndarray) ndarray#
Boolean mask of the samples that belong to the pass nearest
time.The nearest pass is every sample, padding included, within
PASS_WINDOWof the observation time of the sample nearesttimeamong the samplesinsidethe footprint’s own (unpadded) box, so a pass that only clips the padding, where the PSF weight is zero, cannot displace one that covers the footprint. The window is centred on that sample’s own time rather than on its offset, so two passes on either side of the footprint time are never served together. When no sample lies inside the box (a footprint smaller than one native cell), the nearest sample is taken from every served sample; on a tie, the first in array order. When its offset exceedsMAX_TIME_OFFSETno sample is kept.
- _padded(bbox: BoundingBox, resolution_km: float) BoundingBox#
Widen
bboxby half a native cell (resolution_km) on every side.A sample whose cell overlaps the footprint’s box belongs to the footprint even when its center lies just outside: without the padding a coarse grid (ERA5 at 0.25 degree) never serves the points whose cells cover the rim of the contour, and the weigher scores the rim as missing. The PSF weight of a sample outside the contour is zero, so padding changes no aggregated value, only the coverage.
- _require_valid_time(time: datetime64) None#
Raise when a footprint time a time-selecting query needs is NaT.
- Raises:
ValueError – If
timeis NaT.
- _resolution_km(natives: tuple[GridSamples | _PointIndex | _BoxNative, ...]) float#
The native cell size the served samples stand for, km.
The largest spacing a loaded point cloud measured on itself, else the reader’s declared
RESOLUTION_KM. A measured spacing wins because a swath subsampled on delivery (e.g. every eighth pixel) must be credited with the cell it really stands for, not the instrument’s nominal one.
- property file_paths: tuple[Path, ...]#
The staged granule paths passed at construction, in the order given.
- preload() None#
Load and index the whole source now, rather than on the first query.
- Raises:
ValueError – If the loaded data carries a different number of variables than
VARIABLES.TypeError – If
_load_native()returned none of the three native forms.
- classmethod product_stacks() tuple[ProductStack, ...]#
Return every stacked product variable this reader contributes.
Each stack in
PRODUCT_STACKS, followed by a<name>_Standard_Deviationstack of its members’ companions when its members are continuous. The stacked members leave the product as separate variables; every other spec inproduct_variable_specsstays a 1-D product variable.
- classmethod product_variable_specs() tuple[VariableSpec, ...]#
Return every variable this reader contributes to the product definition.
This is the full set of output variables – what appears, prefixed with
<READER_KEY>_, in the FMATCH product definition YAMLs and what the product/test cross-check is built against. It is distinct fromVARIABLES(the smaller set actually read from the source file) and is composed of, in order:each read variable in
VARIABLES;a
<name>_Standard_Deviationcompanion for every continuous (mean-aggregated) read variable, inserted right after its parent bywith_standard_deviation_companions();any reader-specific derived outputs declared in
ADDITIONAL_PRODUCT_VARIABLES(e.g. IGBP ranked scenes).
- samples_in(bbox: BoundingBox, time: datetime64) SampleSet#
Return the native samples inside
bboxthat matchtime.- Parameters:
bbox (BoundingBox) – The footprint’s geographic box, widened by half a native cell (
RESOLUTION_KM) on every side so a sample whose cell overlaps the box is served. A box that wraps the antimeridian (given in the [0, 360) form) is split and both halves served.time (np.datetime64) – The footprint’s observation time (radiometer or camera time). Ignored by a
DAILYreader.
- Returns:
The selected samples; zero samples when the box holds none, or when no sample lies within
MAX_TIME_OFFSEToftime.- Return type:
- Raises:
ValueError – If
timeis NaT and this reader selects in time; if more than one granule is staged for a reader that reads one per period; if the loaded data carries a different number of variables thanVARIABLES; if a stepped reader’s step times are empty, hold a NaT or are not strictly ascending; or if aNEAREST_PASSreader serves samples without per-sample times.TypeError – If the reader’s native data is none of
GridSamples,PointSamplesor a box-queried native.
- property single_file_path: Path#
The one granule of a reader that reads exactly one per period.
- Raises:
ValueError – If more than one granule was staged for this reader.