@motionstudies/core 0.1.0-alpha.2 → 0.1.0-alpha.21

package/domain/published-day.d.ts

@@ -0,0 +1,132 @@
+/** A published day of the national layer, as compiled by the recording host: the sea (one file of 0.1° cells per
+ * five-minute slice), the shoal (one file per slice from cell to member vehicles), and the fish (one gzipped pack
+ * per operator per day, split into parts by sample budget), with an index from operator and vehicle to part and a
+ * manifest of every file. Everything is static; opening one vehicle fetches one pack. Readers supply the transport
+ * and, where the platform lacks DecompressionStream, the inflater; nothing here touches the network itself.
+ */
+export declare const SLICE_SECONDS = 300, CELL_DEGREES = 0.1;
+export interface PublishedDayManifest {
+    readonly date: string;
+    readonly compiledAt: string;
+    readonly hours: {
+        readonly present: readonly string[];
+        readonly missing: readonly string[];
+    };
+    readonly stats: {
+        readonly samples: number;
+        readonly vehicles: number;
+        readonly operators: number;
+        readonly slices: number;
+        readonly excludedRejected: number;
+        readonly excludedOutsideDay: number;
+    };
+    readonly slices: number;
+    readonly files: number;
+    readonly bytes: number;
+    readonly list: readonly {
+        readonly path: string;
+        readonly bytes: number;
+    }[];
+}
+/** One cell of a slice: centre-less floor coordinates, samples (which sum across cells) and distinct vehicles (which do not). */
+export type SliceCell = readonly [lon: number, lat: number, samples: number, vehicles: number];
+export interface SliceFile {
+    readonly slice: string;
+    readonly cells: readonly SliceCell[];
+}
+/** Members of each cell as operator index into the day index and vehicle reference. */
+export type Member = readonly [operatorIndex: number, vehicleRef: string];
+export interface MembersFile {
+    readonly slice: string;
+    readonly cells: Readonly<Record<string, readonly Member[]>>;
+}
+export interface PackPart {
+    readonly file: string;
+    readonly bytes: number;
+    readonly vehicles: number;
+    readonly samples: number;
+    readonly refs?: readonly string[];
+}
+export interface OperatorEntry {
+    readonly ref: string;
+    readonly vehicles: number;
+    readonly samples: number;
+    readonly parts: readonly PackPart[];
+}
+export interface DayIndex {
+    readonly date: string;
+    readonly operators: readonly OperatorEntry[];
+}
+export type TrackPoint = readonly [epochSeconds: number, lon: number, lat: number];
+export type JourneyRun = readonly [epochSeconds: number, lineRef: string | null, journeyRef: string | null];
+export interface VehicleDay {
+    readonly track: readonly TrackPoint[];
+    readonly journeys: readonly JourneyRun[];
+}
+export interface Pack {
+    readonly date: string;
+    readonly operator: string;
+    readonly part: number;
+    readonly parts: number;
+    readonly vehicles: Readonly<Record<string, VehicleDay>>;
+}
+export interface DayTransport {
+    /** Text of a file under the day's base, by its relative path. */
+    readonly text: (path: string) => Promise<string>;
+    /** Bytes of a file under the day's base. */
+    readonly bytes: (path: string) => Promise<Uint8Array>;
+    /** Gunzip; defaults to the platform's DecompressionStream. */
+    readonly inflate?: (bytes: Uint8Array) => Promise<Uint8Array>;
+}
+/** The slice file name for an instant: HH-MM of the five-minute slice containing it, UTC. */
+export declare function sliceNameAt(instant: string | number): string;
+/** The ISO start of a named slice on a date. */
+export declare function sliceStart(date: string, name: string): string;
+/** The cell a position falls in, by the compiler's own rule: floor of the coordinate over the cell size, in double
+ * arithmetic, printed to one decimal. The rule is reproduced exactly rather than rounded, so a position near a
+ * cell edge lands where the compiler put it. */
+export declare function cellOf(lon: number, lat: number): readonly [number, number];
+export declare function cellKey(cellLon: number, cellLat: number): string;
+/** The pack part holding a vehicle: the only part, or the part whose reference list names it. */
+export declare function partFor(index: DayIndex, operatorRef: string, vehicleRef: string): PackPart | null;
+/** Problems with a compiled day's manifest, empty when consistent: all 24 hours of the date accounted for once, one
+ * slice and one members file per slice, and the file list agreeing with the stated file and byte totals (the manifest
+ * itself is the one unlisted file). A trimmed fixture, whose statistics describe the whole day, fails by design. */
+export declare function checkDayManifest(manifest: PublishedDayManifest): string[];
+/** Problems reconciling a slice with its members, empty when they agree: the sea and the shoal must count the same
+ * vehicles. Each cell's distinct-vehicle count equals its distinct members, no member repeats, and neither file names
+ * a cell the other lacks. Samples are not checked here; members carry identities, not sample counts. */
+export declare function reconcileSlice(slice: SliceFile, members: MembersFile): string[];
+/** Position along a track at an instant: linear between the samples around it, null outside the track or across a
+ * gap longer than the tolerance, so a vehicle that stopped reporting does not glide. */
+export declare function positionAt(track: readonly TrackPoint[], epochSeconds: number, maxGapSeconds?: number): readonly [lon: number, lat: number] | null;
+/** The line and journey a vehicle was on at an instant, from its runs. */
+export declare function journeyAt(journeys: readonly JourneyRun[], epochSeconds: number): JourneyRun | null;
+/** Minimal shapes of the platform pieces used, so the module compiles without DOM typings. */
+interface ByteSource {
+    arrayBuffer(): Promise<ArrayBuffer>;
+}
+interface ResponseLike extends ByteSource {
+    readonly ok: boolean;
+    readonly status: number;
+    text(): Promise<string>;
+}
+export type FetchLike = (url: string) => Promise<ResponseLike>;
+export declare function inflateGzip(bytes: Uint8Array): Promise<Uint8Array>;
+/** Packs are stored gzipped, but a server that declares the encoding hands the browser the decoded bytes; either form is accepted by its first two bytes. */
+export declare function inflateIfGzipped(bytes: Uint8Array, inflate: (bytes: Uint8Array) => Promise<Uint8Array>): Promise<Uint8Array>;
+export interface PublishedDay {
+    readonly manifest: PublishedDayManifest;
+    readonly index: DayIndex;
+    readonly slice: (name: string) => Promise<SliceFile>;
+    readonly members: (name: string) => Promise<MembersFile>;
+    readonly pack: (operatorRef: string, vehicleRef: string) => Promise<Pack>;
+    readonly vehicle: (operatorRef: string, vehicleRef: string) => Promise<VehicleDay | null>;
+    /** Slice names present, from the manifest's file list, in order. */
+    readonly sliceNames: readonly string[];
+}
+/** Open a published day: manifest and index up front, everything else on demand and cached per file. */
+export declare function openPublishedDay(transport: DayTransport): Promise<PublishedDay>;
+/** A transport over fetch for a day served under a base URL. */
+export declare function fetchTransport(base: string, fetchImpl?: FetchLike | undefined): DayTransport;
+export {};

package/domain/published-day.js

@@ -0,0 +1,162 @@
+/** A published day of the national layer, as compiled by the recording host: the sea (one file of 0.1° cells per
+ * five-minute slice), the shoal (one file per slice from cell to member vehicles), and the fish (one gzipped pack
+ * per operator per day, split into parts by sample budget), with an index from operator and vehicle to part and a
+ * manifest of every file. Everything is static; opening one vehicle fetches one pack. Readers supply the transport
+ * and, where the platform lacks DecompressionStream, the inflater; nothing here touches the network itself.
+ */
+export const SLICE_SECONDS = 300, CELL_DEGREES = 0.1;
+/** The slice file name for an instant: HH-MM of the five-minute slice containing it, UTC. */
+export function sliceNameAt(instant) {
+    const t = typeof instant === 'number' ? instant : Date.parse(instant);
+    if (!Number.isFinite(t))
+        throw new Error('Invalid instant');
+    const s = Math.floor(t / 1000 / SLICE_SECONDS) * SLICE_SECONDS, d = new Date(s * 1000);
+    return `${String(d.getUTCHours()).padStart(2, '0')}-${String(d.getUTCMinutes()).padStart(2, '0')}`;
+}
+/** The ISO start of a named slice on a date. */
+export function sliceStart(date, name) { return `${date}T${name.replace('-', ':')}:00Z`; }
+/** The cell a position falls in, by the compiler's own rule: floor of the coordinate over the cell size, in double
+ * arithmetic, printed to one decimal. The rule is reproduced exactly rather than rounded, so a position near a
+ * cell edge lands where the compiler put it. */
+export function cellOf(lon, lat) {
+    return [Number((Math.floor(lon / CELL_DEGREES) * CELL_DEGREES).toFixed(1)), Number((Math.floor(lat / CELL_DEGREES) * CELL_DEGREES).toFixed(1))];
+}
+export function cellKey(cellLon, cellLat) { return `${cellLon.toFixed(1)},${cellLat.toFixed(1)}`; }
+/** The pack part holding a vehicle: the only part, or the part whose reference list names it. */
+export function partFor(index, operatorRef, vehicleRef) {
+    const op = index.operators.find(o => o.ref === operatorRef);
+    if (!op)
+        return null;
+    if (op.parts.length === 1)
+        return op.parts[0];
+    return op.parts.find(p => p.refs?.includes(vehicleRef)) ?? null;
+}
+/** Problems with a compiled day's manifest, empty when consistent: all 24 hours of the date accounted for once, one
+ * slice and one members file per slice, and the file list agreeing with the stated file and byte totals (the manifest
+ * itself is the one unlisted file). A trimmed fixture, whose statistics describe the whole day, fails by design. */
+export function checkDayManifest(manifest) {
+    const problems = [], hours = [...manifest.hours.present, ...manifest.hours.missing], distinct = new Set(hours).size;
+    if (hours.length !== 24 || distinct !== 24)
+        problems.push(`Hours present and missing account for ${distinct} distinct of 24, listed ${hours.length} times`);
+    if (hours.some(h => !h.startsWith(`${manifest.date}T`)))
+        problems.push('An hour lies outside the manifest date');
+    const count = (prefix) => manifest.list.filter(f => f.path.startsWith(prefix)).length;
+    if (manifest.stats.slices !== manifest.slices)
+        problems.push(`Statistics count ${manifest.stats.slices} slices; the manifest ${manifest.slices}`);
+    if (count('slices/') !== manifest.slices)
+        problems.push(`${count('slices/')} slice files listed for ${manifest.slices} slices`);
+    if (count('members/') !== manifest.slices)
+        problems.push(`${count('members/')} members files listed for ${manifest.slices} slices`);
+    if (manifest.list.length + 1 !== manifest.files)
+        problems.push(`${manifest.list.length} listed files and the manifest disagree with ${manifest.files} files`);
+    const bytes = manifest.list.reduce((a, f) => a + f.bytes, 0);
+    if (bytes !== manifest.bytes)
+        problems.push(`Listed files hold ${bytes} bytes; the manifest states ${manifest.bytes}`);
+    return problems;
+}
+/** Problems reconciling a slice with its members, empty when they agree: the sea and the shoal must count the same
+ * vehicles. Each cell's distinct-vehicle count equals its distinct members, no member repeats, and neither file names
+ * a cell the other lacks. Samples are not checked here; members carry identities, not sample counts. */
+export function reconcileSlice(slice, members) {
+    if (slice.slice !== members.slice)
+        return [`Slice ${slice.slice} paired with members of ${members.slice}`];
+    const problems = [], seen = new Set();
+    for (const [lon, lat, , vehicles] of slice.cells) {
+        const key = cellKey(lon, lat), refs = members.cells[key];
+        seen.add(key);
+        if (!refs) {
+            problems.push(`Cell ${key} counts ${vehicles} vehicles and has no members`);
+            continue;
+        }
+        const unique = new Set(refs.map(member => JSON.stringify(member))).size;
+        if (unique !== refs.length)
+            problems.push(`Cell ${key} repeats a member`);
+        if (unique !== vehicles)
+            problems.push(`Cell ${key} counts ${vehicles} vehicles but lists ${unique} members`);
+    }
+    for (const key of Object.keys(members.cells))
+        if (!seen.has(key))
+            problems.push(`Members name cell ${key}, absent from the slice`);
+    return problems;
+}
+/** Position along a track at an instant: linear between the samples around it, null outside the track or across a
+ * gap longer than the tolerance, so a vehicle that stopped reporting does not glide. */
+export function positionAt(track, epochSeconds, maxGapSeconds = 180) {
+    if (!track.length || epochSeconds < track[0][0] || epochSeconds > track[track.length - 1][0])
+        return null;
+    let lo = 0, hi = track.length - 1;
+    while (hi - lo > 1) {
+        const mid = (lo + hi) >> 1;
+        if (track[mid][0] <= epochSeconds)
+            lo = mid;
+        else
+            hi = mid;
+    }
+    const a = track[lo], b = track[hi];
+    if (a[0] === epochSeconds || lo === hi)
+        return [a[1], a[2]];
+    if (b[0] === epochSeconds)
+        return [b[1], b[2]];
+    if (b[0] - a[0] > maxGapSeconds)
+        return null;
+    const f = (epochSeconds - a[0]) / (b[0] - a[0]);
+    return [a[1] + (b[1] - a[1]) * f, a[2] + (b[2] - a[2]) * f];
+}
+/** The line and journey a vehicle was on at an instant, from its runs. */
+export function journeyAt(journeys, epochSeconds) {
+    let found = null;
+    for (const run of journeys) {
+        if (run[0] <= epochSeconds)
+            found = run;
+        else
+            break;
+    }
+    return found;
+}
+const platform = globalThis;
+export async function inflateGzip(bytes) {
+    if (!platform.DecompressionStream || !platform.Blob || !platform.Response)
+        throw new Error('No DecompressionStream on this platform; supply an inflater');
+    const stream = new platform.Blob([bytes]).stream().pipeThrough(new platform.DecompressionStream('gzip'));
+    return new Uint8Array(await new platform.Response(stream).arrayBuffer());
+}
+/** Packs are stored gzipped, but a server that declares the encoding hands the browser the decoded bytes; either form is accepted by its first two bytes. */
+export async function inflateIfGzipped(bytes, inflate) {
+    return bytes.length >= 2 && bytes[0] === 0x1f && bytes[1] === 0x8b ? inflate(bytes) : bytes;
+}
+function decodeText(bytes) {
+    if (!platform.TextDecoder)
+        throw new Error('No TextDecoder on this platform');
+    return new platform.TextDecoder().decode(bytes);
+}
+/** Open a published day: manifest and index up front, everything else on demand and cached per file. */
+export async function openPublishedDay(transport) {
+    const manifest = JSON.parse(await transport.text('manifest.json'));
+    const index = JSON.parse(await transport.text('index.json'));
+    if (manifest.date !== index.date)
+        throw new Error('Manifest and index disagree on the date');
+    const cache = new Map();
+    const cached = (path, load) => { let p = cache.get(path); if (!p) {
+        p = load();
+        cache.set(path, p);
+    } return p; };
+    const inflate = transport.inflate ?? inflateGzip;
+    const sliceNames = manifest.list.filter(f => f.path.startsWith('slices/')).map(f => f.path.slice(7, -5)).sort();
+    const pack = (operatorRef, vehicleRef) => { const part = partFor(index, operatorRef, vehicleRef); if (!part)
+        return Promise.reject(new Error(`No pack for ${operatorRef} ${vehicleRef}`)); return cached(part.file, async () => JSON.parse(decodeText(await inflateIfGzipped(await transport.bytes(part.file), inflate)))); };
+    return { manifest, index, sliceNames,
+        slice: name => cached(`slices/${name}.json`, async () => JSON.parse(await transport.text(`slices/${name}.json`))),
+        members: name => cached(`members/${name}.json`, async () => JSON.parse(await transport.text(`members/${name}.json`))),
+        pack,
+        vehicle: async (operatorRef, vehicleRef) => { if (!partFor(index, operatorRef, vehicleRef))
+            return null; return (await pack(operatorRef, vehicleRef)).vehicles[vehicleRef] ?? null; } };
+}
+/** A transport over fetch for a day served under a base URL. */
+export function fetchTransport(base, fetchImpl = platform.fetch) {
+    if (!fetchImpl)
+        throw new Error('No fetch on this platform; supply one');
+    const url = (path) => `${base.replace(/\/$/, '')}/${path}`;
+    const get = async (path) => { const r = await fetchImpl(url(path)); if (!r.ok)
+        throw new Error(`${path} returned HTTP ${r.status}`); return r; };
+    return { text: async (path) => (await get(path)).text(), bytes: async (path) => new Uint8Array(await (await get(path)).arrayBuffer()) };
+}

package/domain/station-calls.d.ts

@@ -0,0 +1,48 @@
+import type { NetworkSnapshot, NetworkTrain } from './network.ts';
+export type StationCallDirection = 'arrival' | 'departure';
+export interface StationCall {
+    readonly id: string;
+    readonly train: NetworkTrain;
+    /** Call occurrence within this train, not a global station index. */
+    readonly index: number;
+    readonly arrival: number;
+    readonly departure: number;
+    readonly origin: string;
+    readonly destination: string;
+    readonly allowsArrival: boolean;
+    readonly allowsDeparture: boolean;
+    readonly movementAvailable?: boolean;
+    readonly source?: string;
+    readonly windowStart?: number;
+    readonly windowEnd?: number;
+}
+export interface StationCallSelection {
+    readonly name: string;
+    /** When supplied, source IDs are authoritative; never fall back to a name match. */
+    readonly stopIds?: readonly string[];
+}
+export interface StationCallPolicy {
+    readonly includeTrain?: (train: NetworkTrain) => boolean;
+    readonly allows?: (train: NetworkTrain, callIndex: number, direction: StationCallDirection) => boolean;
+    readonly origin?: (train: NetworkTrain, snapshot: NetworkSnapshot) => string | undefined;
+    readonly destination?: (train: NetworkTrain, snapshot: NetworkSnapshot) => string | undefined;
+}
+/** Retain repeated visits and source-local indexes; passenger restrictions belong to the adapter. */
+export declare function stationCalls(snapshot: NetworkSnapshot, selection: StationCallSelection, policy?: StationCallPolicy): readonly StationCall[];
+export interface StationCallSource {
+    readonly id: string;
+    readonly selection: StationCallSelection;
+    readonly snapshot?: NetworkSnapshot;
+    readonly movements?: NetworkSnapshot;
+    readonly windowStart: number;
+    readonly windowEnd: number;
+    readonly policy?: StationCallPolicy;
+    readonly loading?: boolean;
+    readonly error?: boolean;
+}
+export type StationSourceStatus = 'loading' | 'error' | 'unavailable' | 'date-mismatch' | 'partial' | 'ready';
+export declare function stationSourceStatus(source: StationCallSource, serviceDate: string, requestedStart: number, requestedEnd: number): StationSourceStatus;
+/** Sources remain independently dated and retryable. Equal names/times never merge services. */
+export declare function combinedStationCalls(serviceDate: string, sources: readonly StationCallSource[]): readonly StationCall[];
+/** Clip to both requested and source windows. Never wrap to an unacquired day. */
+export declare function upcomingStationCalls(calls: readonly StationCall[], direction: StationCallDirection, time: number, windowStart: number, windowEnd: number, route?: string, maxRows?: number): readonly StationCall[];

package/domain/station-calls.js

@@ -0,0 +1,73 @@
+/** Retain repeated visits and source-local indexes; passenger restrictions belong to the adapter. */
+export function stationCalls(snapshot, selection, policy = {}) {
+    const stops = new Set(snapshot.stops.flatMap((stop, index) => (selection.stopIds ? Boolean(stop[4] && selection.stopIds.includes(stop[4])) : stop[2] === selection.name) ? [index] : []));
+    const calls = [], ids = new Map();
+    for (const train of snapshot.trains) {
+        if (policy.includeTrain && !policy.includeTrain(train))
+            continue;
+        const prior = ids.get(train.id);
+        if (prior) {
+            if (prior === train || JSON.stringify(prior) === JSON.stringify(train))
+                continue;
+            throw new Error('Ambiguous duplicate journey identity in station source');
+        }
+        ids.set(train.id, train);
+        train.stops.forEach(([stop, arrival, departure], index) => {
+            if (!stops.has(stop) || !Number.isFinite(arrival) || !Number.isFinite(departure) || departure < arrival)
+                return;
+            const allowsArrival = index > 0 && (policy.allows?.(train, index, 'arrival') ?? true);
+            const allowsDeparture = index < train.stops.length - 1 && (policy.allows?.(train, index, 'departure') ?? true);
+            if (!allowsArrival && !allowsDeparture)
+                return;
+            calls.push({ id: `${train.id}:call:${index}`, train, index, arrival, departure,
+                origin: policy.origin?.(train, snapshot) ?? snapshot.stops[train.stops[0][0]]?.[2] ?? '',
+                destination: policy.destination?.(train, snapshot) ?? (train.headsign || snapshot.stops[train.stops.at(-1)[0]]?.[2] || ''),
+                allowsArrival, allowsDeparture, });
+        });
+    }
+    return calls;
+}
+export function stationSourceStatus(source, serviceDate, requestedStart, requestedEnd) {
+    if (source.error)
+        return 'error';
+    if (source.loading)
+        return 'loading';
+    if (!source.snapshot || ![source.windowStart, source.windowEnd, requestedStart, requestedEnd].every(Number.isFinite) ||
+        source.windowEnd <= source.windowStart || requestedEnd < requestedStart)
+        return 'unavailable';
+    if (source.snapshot.metadata.serviceDate !== serviceDate)
+        return 'date-mismatch';
+    return source.windowStart > requestedStart || source.windowEnd < requestedEnd ? 'partial' : 'ready';
+}
+/** Sources remain independently dated and retryable. Equal names/times never merge services. */
+export function combinedStationCalls(serviceDate, sources) {
+    const result = [], ids = new Set();
+    for (const source of sources) {
+        if (!source.id || ids.has(source.id))
+            throw new Error('Station sources require distinct identities');
+        ids.add(source.id);
+        if (!source.snapshot || source.error || source.loading || source.snapshot.metadata.serviceDate !== serviceDate ||
+            !Number.isFinite(source.windowStart) || !Number.isFinite(source.windowEnd) || source.windowEnd <= source.windowStart)
+            continue;
+        const movementIds = source.movements?.metadata.serviceDate === serviceDate
+            ? new Set(source.movements.trains.map(train => train.id)) : undefined;
+        for (const call of stationCalls(source.snapshot, source.selection, source.policy)) {
+            result.push({ ...call, id: `${encodeURIComponent(source.id)}:${call.id}`, source: source.id,
+                movementAvailable: source.movements ? (movementIds?.has(call.train.id) ?? false) : undefined,
+                windowStart: source.windowStart, windowEnd: source.windowEnd,
+            });
+        }
+    }
+    return result;
+}
+/** Clip to both requested and source windows. Never wrap to an unacquired day. */
+export function upcomingStationCalls(calls, direction, time, windowStart, windowEnd, route = '', maxRows = 4) {
+    if (![time, windowStart, windowEnd].every(Number.isFinite) || time < windowStart || time >= windowEnd || !Number.isSafeInteger(maxRows) || maxRows < 1)
+        return [];
+    const end = Math.min(windowEnd, time + 3600);
+    return calls.filter(call => (direction === 'arrival' ? call.allowsArrival : call.allowsDeparture)
+        && call[direction] >= Math.max(windowStart, call.windowStart ?? windowStart)
+        && call[direction] < Math.min(windowEnd, call.windowEnd ?? windowEnd)
+        && (!route || call.train.route === route) && call[direction] >= time && call[direction] < end)
+        .sort((a, b) => a[direction] - b[direction] || a.id.localeCompare(b.id)).slice(0, maxRows);
+}

package/domain/station-departures.d.ts

@@ -0,0 +1,15 @@
+import type { NetworkSnapshot, ServiceCategory } from './network.ts';
+export interface StationDeparture {
+    readonly id: string;
+    readonly trainId: string;
+    /** Time from the supplied snapshot: realtime-adjusted snapshots contain adjusted times. */
+    readonly time: number;
+    readonly service: string;
+    readonly category: ServiceCategory;
+    readonly destination: string;
+    readonly platform?: string;
+    readonly via: readonly string[];
+    readonly status: 'scheduled' | 'adjusted' | 'cancelled';
+}
+/** Resolve the station against this snapshot, retaining repeat calls but excluding terminators. */
+export declare function stationDepartures(snapshot: NetworkSnapshot, stationName: string): readonly StationDeparture[];

package/domain/station-departures.js

@@ -0,0 +1,25 @@
+/** Resolve the station against this snapshot, retaining repeat calls but excluding terminators. */
+export function stationDepartures(snapshot, stationName) {
+    const stopIndexes = new Set();
+    snapshot.stops.forEach((stop, index) => { if (stop[2] === stationName)
+        stopIndexes.add(index); });
+    if (!stopIndexes.size)
+        return [];
+    const departures = [];
+    for (const train of snapshot.trains) {
+        train.stops.forEach(([stopIndex, , departure], callIndex) => {
+            if (!stopIndexes.has(stopIndex) || !Number.isFinite(departure) || callIndex === train.stops.length - 1)
+                return;
+            const onward = train.stops.slice(callIndex + 1).map(([index]) => snapshot.stops[index]?.[2]).filter((name) => Boolean(name));
+            if (!onward.length)
+                return;
+            const destination = train.headsign.trim() || onward.at(-1);
+            departures.push({ id: `${train.id}:${callIndex}`, trainId: train.id, time: departure,
+                service: train.route || train.shortName, category: train.category, destination,
+                platform: snapshot.stops[stopIndex]?.[3]?.trim() || undefined,
+                via: [...new Set(onward.slice(0, -1).filter((name) => name !== stationName && name !== destination))].slice(0, 3),
+                status: train.realtime?.status ?? 'scheduled' });
+        });
+    }
+    return departures.sort((a, b) => a.time - b.time || a.id.localeCompare(b.id));
+}

package/domain/supported-field.d.ts

@@ -0,0 +1,78 @@
+/** A continuous estimate between scattered instruments. Values are reconstructed from weighted
+ * neighbouring stations; `support` reports how near the closest contributing evidence is, so a
+ * rendering can fade where its evidence weakens. The estimate is a model choice, not an observation
+ * at the sampled place, and missing station values stay missing rather than becoming zero.
+ */
+export interface FieldStation {
+    readonly latitude: number;
+    readonly longitude: number;
+}
+/** Kilometre scales of the neighbour kernel and the support fade. */
+export interface FieldSupportOptions {
+    /** Gaussian scale of neighbour weights. */
+    readonly scaleKm: number;
+    /** Weights taper smoothly from `reachStartKm` to zero at `reachKm`; stations beyond `reachKm` do not contribute. */
+    readonly reachStartKm: number;
+    readonly reachKm: number;
+    /** Support is 1 within `supportNearKm` of the nearest contributing station and fades to 0 at `supportFarKm`. */
+    readonly supportNearKm: number;
+    readonly supportFarKm: number;
+}
+export declare const DEFAULT_FIELD_SUPPORT: FieldSupportOptions;
+export interface FieldNeighbour {
+    readonly index: number;
+    readonly distanceKm: number;
+    readonly weight: number;
+}
+/** One station's contribution at one instant. A null value excludes the station; an incomplete vector excludes it from the vector estimate only. */
+export interface FieldSample {
+    readonly value: number | null;
+    readonly vector?: readonly [east: number, north: number] | null;
+}
+export interface FieldEstimate {
+    readonly value: number | null;
+    readonly vector: readonly [east: number, north: number] | null;
+    /** Fade from the nearest station that supplied a value. */
+    readonly support: number;
+    /** Fade from the nearest station that supplied a complete vector; 0 when none did. */
+    readonly vectorSupport: number;
+}
+export interface FieldGridBounds {
+    readonly west: number;
+    readonly east: number;
+    readonly south: number;
+    readonly north: number;
+    readonly step: number;
+}
+export interface FieldGridCell {
+    readonly latitude: number;
+    readonly longitude: number;
+    readonly neighbours: readonly FieldNeighbour[];
+}
+export interface FieldGrid {
+    readonly width: number;
+    readonly height: number;
+    readonly bounds: FieldGridBounds;
+    readonly cells: readonly FieldGridCell[];
+}
+/** A bilinear reading of a sampled grid. The vector is null when any contributing corner lacks one. */
+export interface FieldGridSample {
+    readonly value: number;
+    readonly vector: readonly [east: number, north: number] | null;
+    readonly support: number;
+}
+/** Great-circle distance on a spherical Earth. */
+export declare function distanceKm(a: FieldStation, b: FieldStation): number;
+/** Stations within reach of a place, weighted by distance. `exclude` withholds one station, for example to test the estimate against it. */
+export declare function fieldNeighbours(stations: readonly FieldStation[], latitude: number, longitude: number, options?: FieldSupportOptions, exclude?: number): FieldNeighbour[];
+/** Weighted estimate from the supplied station samples. Vector components blend as east/north components, never as bearings. */
+export declare function estimateField(neighbours: readonly FieldNeighbour[], sample: (index: number) => FieldSample | null | undefined, options?: FieldSupportOptions): FieldEstimate;
+/** Precompute neighbours for a regular grid. Rows run south to north; the grid includes both bounds. */
+export declare function createFieldGrid(stations: readonly FieldStation[], bounds: FieldGridBounds, options?: FieldSupportOptions): FieldGrid;
+/** Four channels per cell and level: value, east, north, support; levels are the outer dimension.
+ * Unavailable channels are NaN. With `movement`, a cell without a complete vector is entirely
+ * unavailable and its support also fades from the nearest station with a vector.
+ */
+export declare function sampleFieldGrid(grid: FieldGrid, sample: (index: number, level: number) => FieldSample | null | undefined, levels: number, options?: FieldSupportOptions, movement?: boolean): Float32Array;
+/** Bilinear reading of `sampleFieldGrid` output. Null outside the grid or where any contributing corner lacks a value. */
+export declare function sampleFieldGridAt(data: Float32Array, grid: FieldGrid, latitude: number, longitude: number, level: number): FieldGridSample | null;

package/domain/supported-field.js

@@ -0,0 +1,95 @@
+export const DEFAULT_FIELD_SUPPORT = { scaleKm: 90, reachStartKm: 180, reachKm: 240, supportNearKm: 120, supportFarKm: 240 };
+const rad = Math.PI / 180;
+const smooth = (t) => { t = Math.max(0, Math.min(1, t)); return t * t * (3 - 2 * t); };
+const fade = (distance, near, far) => 1 - smooth((distance - near) / Math.max(far - near, 1e-9));
+/** Great-circle distance on a spherical Earth. */
+export function distanceKm(a, b) {
+    const p = (b.latitude - a.latitude) * rad, q = (b.longitude - a.longitude) * rad;
+    const h = Math.sin(p / 2) ** 2 + Math.cos(a.latitude * rad) * Math.cos(b.latitude * rad) * Math.sin(q / 2) ** 2;
+    return 12742.0176 * Math.asin(Math.min(1, Math.sqrt(h)));
+}
+/** Stations within reach of a place, weighted by distance. `exclude` withholds one station, for example to test the estimate against it. */
+export function fieldNeighbours(stations, latitude, longitude, options = DEFAULT_FIELD_SUPPORT, exclude = -1) {
+    const place = { latitude, longitude }, neighbours = [];
+    stations.forEach((station, index) => {
+        const distance = distanceKm(station, place);
+        if (index === exclude || distance >= options.reachKm)
+            return;
+        neighbours.push({ index, distanceKm: distance, weight: Math.exp(-.5 * (distance / options.scaleKm) ** 2) * fade(distance, options.reachStartKm, options.reachKm) });
+    });
+    return neighbours;
+}
+/** Weighted estimate from the supplied station samples. Vector components blend as east/north components, never as bearings. */
+export function estimateField(neighbours, sample, options = DEFAULT_FIELD_SUPPORT) {
+    let weight = 0, value = 0, vectorWeight = 0, east = 0, north = 0, nearest = Infinity, nearestVector = Infinity;
+    for (const n of neighbours) {
+        const s = sample(n.index);
+        if (!s || s.value === null || !Number.isFinite(s.value))
+            continue;
+        weight += n.weight;
+        value += n.weight * s.value;
+        nearest = Math.min(nearest, n.distanceKm);
+        const v = s.vector;
+        if (v && Number.isFinite(v[0]) && Number.isFinite(v[1])) {
+            vectorWeight += n.weight;
+            east += n.weight * v[0];
+            north += n.weight * v[1];
+            nearestVector = Math.min(nearestVector, n.distanceKm);
+        }
+    }
+    if (weight < 1e-12)
+        return { value: null, vector: null, support: 0, vectorSupport: 0 };
+    return { value: value / weight, vector: vectorWeight ? [east / vectorWeight, north / vectorWeight] : null,
+        support: fade(nearest, options.supportNearKm, options.supportFarKm), vectorSupport: vectorWeight ? fade(nearestVector, options.supportNearKm, options.supportFarKm) : 0 };
+}
+/** Precompute neighbours for a regular grid. Rows run south to north; the grid includes both bounds. */
+export function createFieldGrid(stations, bounds, options = DEFAULT_FIELD_SUPPORT) {
+    const width = Math.round((bounds.east - bounds.west) / bounds.step) + 1, height = Math.round((bounds.north - bounds.south) / bounds.step) + 1, cells = [];
+    for (let y = 0; y < height; y++)
+        for (let x = 0; x < width; x++) {
+            const latitude = bounds.south + y * bounds.step, longitude = bounds.west + x * bounds.step;
+            cells.push({ latitude, longitude, neighbours: fieldNeighbours(stations, latitude, longitude, options) });
+        }
+    return { width, height, bounds, cells };
+}
+/** Four channels per cell and level: value, east, north, support; levels are the outer dimension.
+ * Unavailable channels are NaN. With `movement`, a cell without a complete vector is entirely
+ * unavailable and its support also fades from the nearest station with a vector.
+ */
+export function sampleFieldGrid(grid, sample, levels, options = DEFAULT_FIELD_SUPPORT, movement = false) {
+    const data = new Float32Array(grid.cells.length * levels * 4);
+    data.fill(NaN);
+    for (let level = 0; level < levels; level++)
+        for (let i = 0; i < grid.cells.length; i++) {
+            const e = estimateField(grid.cells[i].neighbours, (index) => sample(index, level), options), o = (level * grid.cells.length + i) * 4;
+            if (e.value === null || (movement && !e.vector))
+                continue;
+            data[o] = e.value;
+            data[o + 3] = movement ? Math.min(e.support, e.vectorSupport) : e.support;
+            if (e.vector) {
+                data[o + 1] = e.vector[0];
+                data[o + 2] = e.vector[1];
+            }
+        }
+    return data;
+}
+/** Bilinear reading of `sampleFieldGrid` output. Null outside the grid or where any contributing corner lacks a value. */
+export function sampleFieldGridAt(data, grid, latitude, longitude, level) {
+    const { width, height, bounds } = grid, x = (longitude - bounds.west) / bounds.step, y = (latitude - bounds.south) / bounds.step;
+    if (!(x >= 0 && x <= width - 1 && y >= 0 && y <= height - 1) || level < 0 || (level + 1) * width * height * 4 > data.length)
+        return null;
+    const x0 = Math.floor(x), y0 = Math.floor(y), fx = x - x0, fy = y - y0, out = [0, 0, 0, 0];
+    let vector = true;
+    for (const [xx, yy, w] of [[x0, y0, (1 - fx) * (1 - fy)], [Math.min(x0 + 1, width - 1), y0, fx * (1 - fy)], [x0, Math.min(y0 + 1, height - 1), (1 - fx) * fy], [Math.min(x0 + 1, width - 1), Math.min(y0 + 1, height - 1), fx * fy]]) {
+        if (w < 1e-10)
+            continue;
+        const o = (level * width * height + yy * width + xx) * 4;
+        if (!Number.isFinite(data[o]) || !Number.isFinite(data[o + 3]))
+            return null;
+        if (!Number.isFinite(data[o + 1]) || !Number.isFinite(data[o + 2]))
+            vector = false;
+        for (let c = 0; c < 4; c++)
+            out[c] += data[o + c] * w;
+    }
+    return { value: out[0], vector: vector ? [out[1], out[2]] : null, support: out[3] };
+}

package/domain/vehicle-calls.d.ts

@@ -0,0 +1,11 @@
+import type { NetworkSnapshot, NetworkTrain } from './network.ts';
+export interface VehicleCall {
+    readonly id: string;
+    readonly name: string;
+    readonly time: number;
+    readonly platform?: string;
+    readonly atStop: boolean;
+    readonly isDestination: boolean;
+}
+/** Remaining calls from this snapshot's train. Times may already include realtime adjustments. */
+export declare function vehicleCalls(snapshot: NetworkSnapshot, train: NetworkTrain, time: number): readonly VehicleCall[];