leaflet-city-mapper 0.0.1 → 0.0.3

package/src/terrain/sample.js

@@ -0,0 +1,7 @@
+export function sample(tiles, px, py) {
+  const tileX = Math.floor(px / 256), tileY = Math.floor(py / 256);
+  const tile = tiles.get(`${tileX}/${tileY}`);
+  if (!tile) return null;
+  const elevation = tile[(py - tileY * 256) * 256 + px - tileX * 256];
+  return Number.isFinite(elevation) && elevation > -1000 ? elevation : null;
+}

package/src/terrain/setLayerTerrain.js

@@ -0,0 +1,5 @@
+import { setRasterTerrain } from '../raster/setRasterTerrain.js';
+
+export function setLayerTerrain(layer) {
+  if (layer._type === 'pmtiles') setRasterTerrain(layer, this);
+}

package/src/terrain/setTerrainHeight.js

@@ -0,0 +1,10 @@
+import { getElevation } from './getElevation.js';
+
+export function setTerrainHeight(layer) {
+  if (!layer.root || !layer.location) return;
+  layer._terrainBaseY ??= layer.root.position[1];
+  const terrain = layer._map?._terrain3D;
+  const elevation = terrain && getElevation(terrain, layer.location);
+  const offset = Number.isFinite(elevation) ? layer.options.terrainOffset ?? 0 : 0;
+  layer.root.position[1] = layer._terrainBaseY + (elevation ?? 0) + offset;
+}

package/src/terrain/terrainAtlas.js

@@ -0,0 +1,23 @@
+/** Float32 DEM atlas; missing samples retain the shader no-data sentinel. */
+
+export function terrainAtlas(tiles, bounds) {
+  const columns = bounds.max.x - bounds.min.x + 1, rows = bounds.max.y - bounds.min.y + 1;
+  const elevations = new Float32Array(columns * rows * 65536).fill(-1e9);
+  let low = Infinity, high = -Infinity;
+  let key, values, coordinates, x, y, start, row, column, value;
+  for ([key, values] of tiles) {
+    coordinates = key.split('/');
+    x = Number(coordinates[0]), y = Number(coordinates[1]);
+    start = (y - bounds.min.y) * 256 * columns * 256 + (x - bounds.min.x) * 256;
+    for (row = 0; row < 256; row++) {
+      for (column = 0; column < 256; column++) {
+        value = values[row * 256 + column];
+        if (!Number.isFinite(value) || value <= -1000) continue;
+        elevations[start + row * columns * 256 + column] = value;
+        low = Math.min(low, value);
+        high = Math.max(high, value);
+      }
+    }
+  }
+  return { elevations, low, high };
+}

package/src/terrain/terrainElevation.js

@@ -0,0 +1,20 @@
+import { sample } from './sample.js';
+
+/** Bilinear DEM sampling, with the original first-valid-neighbour no-data rule. */
+
+export function terrainElevation(point, tiles) {
+  const x = point.x - 0.5, y = point.y - 0.5;
+  const x0 = Math.floor(x), y0 = Math.floor(y);
+  const fx = x - x0, fy = y - y0;
+  const values = [sample(tiles, x0, y0), sample(tiles, x0 + 1, y0), sample(tiles, x0, y0 + 1), sample(tiles, x0 + 1, y0 + 1)];
+  let valid = 0, first = null;
+  let value;
+  for (value of values) {
+    if (!Number.isFinite(value)) continue;
+    if (!valid) first = value;
+    valid++;
+  }
+  if (valid < 4) return first;
+  return (values[0] * (1 - fx) + values[1] * fx) * (1 - fy) +
+    (values[2] * (1 - fx) + values[3] * fx) * fy;
+}

package/src/terrain/terrainGrid.js

@@ -0,0 +1,24 @@
+/** Unit tile grid; the vertex shader supplies the elevation. */
+
+export function terrainGrid(divisions) {
+  const positions = [], normals = [], uvs = [], indices = [];
+  let row, column, a, b, c, d;
+  for (row = 0; row <= divisions; row++) {
+    for (column = 0; column <= divisions; column++) {
+      positions.push(column / divisions, 0, row / divisions);
+      normals.push(0, 1, 0);
+      uvs.push(column / divisions, row / divisions);
+      if (row < divisions && column < divisions) {
+        a = row * (divisions + 1) + column;
+        b = a + 1, c = a + divisions + 1, d = c + 1;
+        indices.push(a, c, b, b, c, d);
+      }
+    }
+  }
+  return {
+    position: { size: 3, data: new Float32Array(positions) },
+    normal: { size: 3, data: new Float32Array(normals) },
+    uv: { size: 2, data: new Float32Array(uvs) },
+    index: { data: new Uint16Array(indices) },
+  };
+}

package/src/traffic/addTraffic.js

@@ -0,0 +1,7 @@
+import { placeTraffic } from './placeTraffic.js';
+import { startTraffic } from './startTraffic.js';
+
+export function addTraffic(state) {
+  state.layer._beforeRaycast = placeTraffic.bind(null, state);
+  if (state.ready) startTraffic(state);
+}

package/src/traffic/along.js

@@ -0,0 +1,14 @@
+/** Point and unit direction at distance `s` along `points` ([x, z, y?]) with cumulative `lengths`. */
+
+export function along(points, lengths, s) {
+  let i = 1, end = lengths.length - 1;
+  let middle;
+  while (i < end) {
+    middle = i + Math.floor((end - i) / 2);
+    if (lengths[middle] < s) i = middle + 1;
+    else end = middle;
+  }
+  const a = points[i - 1], b = points[i], d = lengths[i] - lengths[i - 1] || 1;
+  const t = (s - lengths[i - 1]) / d;
+  return { x: a[0] + (b[0] - a[0]) * t, y: (a[2] ?? 0) + ((b[2] ?? 0) - (a[2] ?? 0)) * t, z: a[1] + (b[1] - a[1]) * t, dx: (b[0] - a[0]) / d, dz: (b[1] - a[1]) / d };
+}

package/src/traffic/clockMinutes.js

@@ -0,0 +1,6 @@
+
+
+export function clockMinutes(clockTime) {
+  const parts = clockTime.split(':');
+  return Number(parts[0]) * 60 + Number(parts[1]);
+}

package/src/traffic/compareDepartures.js

@@ -0,0 +1,5 @@
+
+
+export function compareDepartures(a, b) {
+  return a - b;
+}

package/src/traffic/createTrafficMesh.js

@@ -0,0 +1,7 @@
+import * as createMeshApi from '../rendering/createMesh.js';
+
+export function createTrafficMesh(state, data, featureFaces, material, usage = 0) {
+  const result = createMeshApi.createMesh(state.layer, data, material, usage);
+  result.extras.featureFaces = featureFaces;
+  return result;
+}

package/src/traffic/destroyTrafficGpu.js

@@ -0,0 +1,5 @@
+export function destroyTrafficGpu(state) {
+  let buffer;
+  for (buffer of [state.gpu?.params, state.gpu?.template, state.gpu?.points, state.gpu?.vehicles]) buffer?.destroy();
+  state.gpu = null;
+}

package/src/traffic/disposeTrafficObject.js

@@ -0,0 +1,13 @@
+import * as disposeMesh from '../rendering/disposeMesh.js';
+import * as setHighlight from '../picking/setHighlight.js';
+import * as showFeature from '../picking/showFeature.js';
+
+export function disposeTrafficObject(state, object) {
+  if (!object) return;
+  if (state.layer.highlightedFeature?.mesh === object) {
+    setHighlight.setHighlight(state.layer, null);
+    state.layer.hoveredFeature = null;
+    showFeature.showFeature(state.layer, null);
+  };
+  disposeMesh.disposeMesh(state.layer, object);
+}

package/src/traffic/drawTrafficRoute.js

@@ -0,0 +1,6 @@
+import { routeGeometry } from './routeGeometry.js';
+import { createTrafficMesh } from './createTrafficMesh.js';
+
+export function drawTrafficRoute(state) {
+  state.route = createTrafficMesh(state, routeGeometry(state.lanes), [], { unlit: true });
+}

package/src/traffic/failTraffic.js

@@ -0,0 +1,6 @@
+export function failTraffic(state, error) {
+  state.error = error;
+  state.status.textContent = `Error: ${error.message}`;
+  state.status.setAttribute('role', 'alert');
+  console.error('Traffic layer', error);
+}

package/src/traffic/failedTrafficLoad.js

@@ -0,0 +1,5 @@
+import { failTraffic } from './failTraffic.js';
+
+export function failedTrafficLoad(state, signal, error) {
+  if (!signal.aborted && error.name !== 'AbortError') failTraffic(state, error);
+}

package/src/traffic/fetchTrafficLanes.js

@@ -0,0 +1,10 @@
+export async function fetchTrafficLanes(state, time, signal) {
+  const href = new URL(state.url, document.baseURI).href;
+  if (/\.m?js(?:[?#]|$)/i.test(state.url)) return (await import(href)).default(time, state.layer, signal);
+  if (!state.file) {
+    const response = await fetch(href, { signal });
+    if (!response.ok) throw new Error(`HTTP ${response.status}`);
+    state.file = await response.json();
+  }
+  return state.file;
+}

package/src/traffic/groundTraffic.js

@@ -0,0 +1,15 @@
+import { getElevation } from '../terrain/getElevation.js';
+
+export function groundTraffic(state) {
+  const { layer, lanes } = state;
+  const terrain = layer._map._terrain3D?._terrainReady && layer._map._terrain3D, base = layer.root.position[1];
+  let lane, i, elevation;
+  for (lane of lanes) {
+    for (i = 0; i < lane.points.length; i++) {
+      elevation = terrain && (lane.elevation ?? getElevation(terrain, lane.latlngs[i]));
+      lane.points[i][2] = Number.isFinite(elevation) && Number.isFinite(base) ? elevation - base : 0;
+    }
+  }
+  state.lanesDirty = true;
+  state.nextCheck = -Infinity;
+}

package/src/traffic/initTraffic.js

@@ -0,0 +1,54 @@
+/**
+ * Vehicles moving along GeoJSON lanes; a compute pass places them on the GPU every frame.
+ * Source options:
+ * - `url`: GeoJSON (Multi)LineStrings, or a JS module whose default export `(time, layer, signal)` returns them (time-dependent data);
+ * - `vehicle`: 'car' (default), 'boat', 'bus', 'tram', 'metro', 'bicycle' or 'pedestrian'; `side`: offset right of the lane centreline in metres (vehicle default);
+ * - `speed`: m/s of `flow` lanes; `outline`: draw the lanes.
+ * Lane properties: `flow` (vehicles/h) or `timetable` (see scheduledTrips.js), optional fixed `elevation`
+ * in metres, `name`/`description` for the default tooltip. The status line (count and vehicle type) sits under the layer's entry in the layers control.
+ * L.Control.Time on the same map drives the time; without one, vehicles run live.
+ */
+import { build } from '../geometry/build.js';
+import { init3D } from '../scene/init3D.js';
+import { tooltip } from './tooltip.js';
+import { speed } from '../constants/speed.js';
+import { vehicles } from '../constants/vehicles.js';
+import { trafficShader } from '../constants/trafficShader.js';
+import { addTraffic } from './addTraffic.js';
+import { placeTraffic } from './placeTraffic.js';
+import { removeTraffic } from './removeTraffic.js';
+import { readyTraffic } from './readyTraffic.js';
+import { tickTraffic } from './tickTraffic.js';
+import { updateTraffic } from './updateTraffic.js';
+import { updateTrafficTerrain } from './updateTrafficTerrain.js';
+
+export function initTraffic(layer, { url, vehicle = 'car', side, outline, speed: velocity = speed, ...options } = {}) {
+  const { model, lights, side: defaultSide } = vehicles[vehicle];
+  side ??= defaultSide;
+  init3D(layer, {
+    tooltip,
+    nightLight: Boolean(lights),
+    ...options,
+  });
+  const status = L.DomUtil.create('small', 'traffic-status');
+  status.style.cssText = 'display:block;max-width:200px;margin:0 0 4px 20px;color:#555;';
+  const template = build(model(0)), stride = template.position.length, vertices = stride / 3;
+  const params = new ArrayBuffer(16), paramFloats = new Float32Array(params), paramInts = new Uint32Array(params);
+  const state = {
+    layer, url, model, lights, side, outline, velocity, vehicle, status, template, stride, vertices,
+    params, paramFloats, paramInts,
+    ready: false, frame: undefined, request: undefined, lanes: [], tripsKey: undefined,
+    mesh: undefined, route: undefined, litAt: -Infinity,
+    trips: [], rebasedAt: 0, nextCheck: -Infinity, elapsed: 0, gpu: null,
+    lanesDirty: true, stale: false, loaded: undefined, file: undefined,
+  };
+  // Lane points are [x, z, ground, cumulative length]. Mirrors along() and place().
+  state.shader = trafficShader;
+  state.tick = tickTraffic.bind(null, state);
+  state.update = updateTraffic.bind(null, state);
+  state.terrain = updateTrafficTerrain.bind(null, state);
+  layer._beforeRaycast = placeTraffic.bind(null, state);
+  layer.on('ready', readyTraffic.bind(null, state));
+  layer.on('add', addTraffic.bind(null, state));
+  layer.on('remove', removeTraffic.bind(null, state));
+}

package/src/traffic/loadTraffic.js

@@ -0,0 +1,15 @@
+import { trafficNow } from './trafficNow.js';
+import { fetchTrafficLanes } from './fetchTrafficLanes.js';
+import { loadedTrafficLanes } from './loadedTrafficLanes.js';
+import { failedTrafficLoad } from './failedTrafficLoad.js';
+
+export function loadTraffic(state) {
+  const { time } = trafficNow(state);
+  if (!Number.isFinite(time?.getTime())) return;
+  state.request?.abort();
+  const { signal } = state.request = new AbortController();
+  fetchTrafficLanes(state, time, signal).then(
+    loadedTrafficLanes.bind(null, state, signal),
+    failedTrafficLoad.bind(null, state, signal),
+  );
+}

package/src/traffic/loadedTrafficLanes.js

@@ -0,0 +1,16 @@
+import { toLanes } from './toLanes.js';
+import { prepareTraffic } from './prepareTraffic.js';
+
+export function loadedTrafficLanes(state, signal, data) {
+  if (signal.aborted || !state.layer._map) return;
+  state.error = null;
+  if (data === state.loaded) {
+    state.nextCheck = -Infinity;
+    if (!state.frame) state.tick();
+    return;
+  }
+  state.loaded = data;
+  state.lanes = toLanes(data);
+  prepareTraffic(state);
+  if (!state.frame) state.tick();
+}

package/src/traffic/locateTraffic.js

@@ -0,0 +1,9 @@
+import { along } from './along.js';
+
+export function locateTraffic(state, trip) {
+  const { points, lengths, length } = trip.lane;
+  let distance = trip.distance + state.elapsed * trip.speed;
+  distance -= length * Math.floor(distance / length);
+  const { x, y, z, dx, dz } = along(points, lengths, distance);
+  return { x: x - dz * state.side, y, z: z + dx * state.side, dx, dz };
+}

package/src/traffic/moveTraffic.js

@@ -0,0 +1,15 @@
+export function moveTraffic(state, seconds) {
+  const { device } = state.gpu;
+  state.paramFloats[0] = state.elapsed = seconds;
+  state.paramFloats[1] = state.side;
+  state.paramInts[2] = state.vertices;
+  state.paramInts[3] = state.trips.length * state.vertices;
+  device.queue.writeBuffer(state.gpu.params, 0, state.params);
+  const encoder = device.createCommandEncoder(), pass = encoder.beginComputePass();
+  pass.setPipeline(state.gpu.pipeline);
+  pass.setBindGroup(0, state.gpu.bindGroup);
+  pass.dispatchWorkgroups(state.trips.length);
+  pass.end();
+  device.queue.submit([encoder.finish()]);
+  state.stale = true;
+}

package/src/traffic/placeTraffic.js

@@ -0,0 +1,21 @@
+import { locateTraffic } from './locateTraffic.js';
+
+export function placeTraffic(state) {
+  if (!state.mesh || !state.stale) return;
+  state.stale = false;
+  const { position, normal } = state.mesh.geometry.attributes;
+  let i, trip, cx, y, cz, dx, dz, k, o, lx, lz, nx, nz;
+  for (i = 0; i < state.trips.length; i++) {
+    trip = state.trips[i];
+    ({ x: cx, y, z: cz, dx, dz } = locateTraffic(state, trip));
+    for (k = 0, o = i * state.stride; k < state.stride; k += 3, o += 3) {
+      lx = state.template.position[k], lz = state.template.position[k + 2], nx = state.template.normal[k], nz = state.template.normal[k + 2];
+      position.data[o] = cx - lx * dz - lz * dx;
+      position.data[o + 1] = y + state.template.position[k + 1];
+      position.data[o + 2] = cz + lx * dx - lz * dz;
+      normal.data[o] = -nx * dz - nz * dx;
+      normal.data[o + 1] = state.template.normal[k + 1];
+      normal.data[o + 2] = nx * dx - nz * dz;
+    }
+  }
+}

package/src/traffic/prepareTraffic.js

@@ -0,0 +1,25 @@
+import { disposeTrafficObject } from './disposeTrafficObject.js';
+import { groundTraffic } from './groundTraffic.js';
+
+export function prepareTraffic(state) {
+  const { layer, lanes } = state;
+  const crs = layer._map.options.crs, origin = crs.project(layer.location);
+  const metres = Math.cos(layer.location.lat * Math.PI / 180);
+  let lane, lat, lng, point, i;
+  for (lane of lanes) {
+    lane.points = [];
+    for ([lat, lng] of lane.latlngs) {
+      point = crs.project(new L.LatLng(lat, lng));
+      lane.points.push([(point.x - origin.x) * metres, -(point.y - origin.y) * metres, 0]);
+    }
+    lane.lengths = [0];
+    for (i = 1; i < lane.points.length; i++) {
+      lane.lengths.push(lane.lengths[i - 1] + Math.hypot(lane.points[i][0] - lane.points[i - 1][0], lane.points[i][1] - lane.points[i - 1][1]));
+    }
+    lane.length = lane.lengths.at(-1);
+  }
+  groundTraffic(state);
+  disposeTrafficObject(state, state.route);
+  state.route = null;
+  state.tripsKey = null;
+}

package/src/traffic/readyTraffic.js

@@ -0,0 +1,6 @@
+import { startTraffic } from './startTraffic.js';
+
+export function readyTraffic(state) {
+  state.ready = true;
+  if (state.layer._map) startTraffic(state);
+}

package/src/traffic/rebaseTraffic.js

@@ -0,0 +1,58 @@
+import { destroyTrafficGpu } from './destroyTrafficGpu.js';
+import { trafficStorage } from './trafficStorage.js';
+
+export function rebaseTraffic(state) {
+  const { device } = state.layer.renderer;
+  if (state.gpu?.device !== device) {
+    destroyTrafficGpu(state);
+    state.gpu = {
+      device,
+      pipeline: device.createComputePipeline({ layout: 'auto', compute: { module: device.createShaderModule({ code: state.shader }), entryPoint: 'main' } }),
+      params: device.createBuffer({ size: 16, usage: GPUBufferUsage.UNIFORM | GPUBufferUsage.COPY_DST }),
+      template: trafficStorage(device, new Float32Array([...state.template.position, ...state.template.normal])),
+    };
+    state.lanesDirty = true;
+  }
+  let lane, i, point, trip, binding;
+  if (state.lanesDirty) {
+    let count = 0;
+    for (lane of state.lanes) count += lane.points.length;
+    const data = new Float32Array(count * 4);
+    let offset = 0;
+    for (lane of state.lanes) {
+      lane.start = offset / 4;
+      for (i = 0; i < lane.points.length; i++) {
+        point = lane.points[i];
+        data[offset++] = point[0];
+        data[offset++] = point[1];
+        data[offset++] = point[2];
+        data[offset++] = lane.lengths[i];
+      }
+    }
+    state.gpu.points?.destroy();
+    state.gpu.points = trafficStorage(device, data);
+    state.lanesDirty = false;
+  }
+  const data = new Float32Array(state.trips.length * 4), index = new Uint32Array(data.buffer);
+  for (i = 0; i < state.trips.length; i++) {
+    trip = state.trips[i];
+    index[i * 4] = trip.lane.start;
+    index[i * 4 + 1] = trip.lane.points.length;
+    data[i * 4 + 2] = trip.distance;
+    data[i * 4 + 3] = trip.speed;
+    state.mesh.extras.featureFaces[i].properties = trip.properties;
+  }
+  if (state.gpu.vehicles?.size >= data.byteLength) device.queue.writeBuffer(state.gpu.vehicles, 0, data);
+  else {
+    state.gpu.vehicles?.destroy();
+    state.gpu.vehicles = trafficStorage(device, data);
+  }
+  const { position, normal } = state.mesh.geometry.attributes;
+  const buffers = [state.gpu.params, state.gpu.template, state.gpu.points, state.gpu.vehicles, position.buffer, normal.buffer];
+  const entries = [];
+  for (binding = 0; binding < buffers.length; binding++) entries.push({ binding, resource: { buffer: buffers[binding] } });
+  state.gpu.bindGroup = device.createBindGroup({
+    layout: state.gpu.pipeline.getBindGroupLayout(0),
+    entries,
+  });
+}

package/src/traffic/rebuildTrafficMesh.js

@@ -0,0 +1,28 @@
+import { build } from '../geometry/build.js';
+import { createTrafficMesh } from './createTrafficMesh.js';
+import { disposeTrafficObject } from './disposeTrafficObject.js';
+import { trafficBounds } from './trafficBounds.js';
+
+export function rebuildTrafficMesh(state) {
+  const { trips, model, vertices } = state;
+  disposeTrafficObject(state, state.mesh);
+  state.mesh = null;
+  if (!trips.length) return;
+  const color = [], featureFaces = [];
+  let i, trip, component;
+  for (i = 0; i < trips.length; i++) {
+    trip = trips[i];
+    for (component of build(model(i, trip.properties)).color) color.push(component);
+    featureFaces.push({
+      id: trip.id, properties: trip.properties, faceStart: i * vertices / 3, faceEnd: (i + 1) * vertices / 3,
+      vertexStart: i * vertices, vertexEnd: (i + 1) * vertices
+    });
+  }
+  const count = color.length / 3;
+  state.mesh = createTrafficMesh(state, {
+    position: new Array(color.length).fill(0), normal: new Array(color.length).fill(0), color,
+  }, featureFaces, {}, GPUBufferUsage.STORAGE | GPUBufferUsage.COPY_SRC);
+  state.mesh.geometry.drawRange.count = count;
+  state.mesh.extras.dynamicLights = true;
+  trafficBounds(state);
+}

package/src/traffic/removeTraffic.js

@@ -0,0 +1,17 @@
+import { disposeTrafficObject } from './disposeTrafficObject.js';
+import { destroyTrafficGpu } from './destroyTrafficGpu.js';
+
+export function removeTraffic(state) {
+  state.request?.abort();
+  cancelAnimationFrame(state.frame);
+  state.frame = undefined;
+  state.layer._map.off('simulationtimechange layerremove', state.update).off('terrainchange', state.terrain);
+  disposeTrafficObject(state, state.mesh);
+  disposeTrafficObject(state, state.route);
+  state.mesh = state.route = null;
+  delete state.layer._beforeRaycast;
+  destroyTrafficGpu(state);
+  state.lanes = [];
+  state.trips = [];
+  state.tripsKey = state.loaded = null;
+}

package/src/traffic/routeColor.js

@@ -0,0 +1,6 @@
+export function routeColor(properties) {
+  const color = properties?.route_color;
+  if (!/^[\da-f]{6}$/i.test(color)) return [0.1, 0.55, 0.65];
+  const value = Number.parseInt(color, 16);
+  return [(value >> 16 & 255) / 255, (value >> 8 & 255) / 255, (value & 255) / 255];
+}

package/src/traffic/routeGeometry.js

@@ -0,0 +1,31 @@
+/** Dashed lane outline, in layer-local metres. */
+
+export function routeGeometry(lanes) {
+  const position = [], normal = [], color = [];
+  let points, lengths, properties, raw, rgb, tint, i, a, b, dx, dz, length, x, z, distance, start, end;
+  let p, q, point;
+  for ({ points, lengths, properties } of lanes) {
+    raw = properties?.route_color;
+    rgb = /^[\da-f]{6}$/i.test(raw) ? Number.parseInt(raw, 16) : null;
+    tint = rgb === null ? [0.08, 0.7, 0.85]
+      : [(rgb >> 16 & 255) / 255, (rgb >> 8 & 255) / 255, (rgb & 255) / 255];
+    for (i = 1; i < points.length; i++) {
+      a = points[i - 1], b = points[i];
+      dx = b[0] - a[0], dz = b[1] - a[1], length = lengths[i] - lengths[i - 1];
+      if (!length) continue;
+      x = -dz / length * 2.5, z = dx / length * 2.5;
+      for (distance = Math.floor(lengths[i - 1] / 30) * 30; distance < lengths[i]; distance += 30) {
+        start = Math.max(0, (distance - lengths[i - 1]) / length);
+        end = Math.min(1, (distance + 18 - lengths[i - 1]) / length);
+        if (end <= start) continue;
+        p = [a[0] + dx * start, a[1] + dz * start, a[2] + (b[2] - a[2]) * start + 0.4];
+        q = [a[0] + dx * end, a[1] + dz * end, a[2] + (b[2] - a[2]) * end + 0.4];
+        for (point of [[p[0] - x, p[2], p[1] - z], [p[0] + x, p[2], p[1] + z], [q[0] + x, q[2], q[1] + z],
+        [p[0] - x, p[2], p[1] - z], [q[0] + x, q[2], q[1] + z], [q[0] - x, q[2], q[1] - z]]) {
+          position.push(...point); normal.push(0, 1, 0); color.push(...tint);
+        }
+      }
+    }
+  }
+  return { position, normal, color };
+}

package/src/traffic/scheduledTrips.js

@@ -0,0 +1,40 @@
+import { clockMinutes } from './clockMinutes.js';
+import { compareDepartures } from './compareDepartures.js';
+
+/**
+ * Trips of a frequency-based timetable (like GTFS frequencies.txt) underway at `date`, in the timetable's local time.
+ * @param {{timeZone: string, duration: number, services: Array<{from: string, to: string, every?: number, months?: number[], weekdays?: number[]}>}} timetable
+ * `duration` and `every` are minutes; `from`/`to` are 'HH:MM' first/last departures; months are 1-12, weekdays 0 (Sunday)-6.
+ * @returns {Array<{id: string, departure: string, progress: number}>} `progress` is 0-1 along the lane.
+ */
+
+export function scheduledTrips(date, { timeZone, duration, services }) {
+  const time = new Date(date);
+  if (!Number.isFinite(time.getTime())) throw new Error('Invalid date');
+  const clock = new Intl.DateTimeFormat('en-GB', {
+    timeZone, year: 'numeric', month: '2-digit', day: '2-digit',
+    hour: '2-digit', minute: '2-digit', second: '2-digit', hourCycle: 'h23',
+  });
+  const parts = {};
+  let part, from, to, every, months, weekdays, departure;
+  for (part of clock.formatToParts(time)) {
+    if (part.type !== 'literal') parts[part.type] = Number(part.value);
+  }
+  const minute = parts.hour * 60 + parts.minute + (parts.second + time.getUTCMilliseconds() / 1000) / 60;
+  const weekday = new Date(Date.UTC(parts.year, parts.month - 1, parts.day)).getUTCDay();
+  const departures = [];
+  for ({ from, to, every = Infinity, months, weekdays } of services) {
+    if (months?.includes(parts.month) === false || weekdays?.includes(weekday) === false) continue;
+    for (departure = clockMinutes(from); departure <= clockMinutes(to); departure += every) departures.push(departure);
+  }
+  departures.sort(compareDepartures);
+  const trips = [];
+  for (departure of departures) {
+    if (minute >= departure && minute < departure + duration) trips.push({
+      id: `${parts.year}-${parts.month}-${parts.day}-${departure}`,
+      departure: `${String(Math.floor(departure / 60)).padStart(2, '0')}:${String(departure % 60).padStart(2, '0')}`,
+      progress: (minute - departure) / duration,
+    });
+  }
+  return trips;
+}

package/src/traffic/startTraffic.js

@@ -0,0 +1,8 @@
+import { underEntry } from '../scene/underEntry.js';
+
+export function startTraffic(state) {
+  state.layer._map.off('simulationtimechange layerremove', state.update).on('simulationtimechange layerremove', state.update);
+  state.layer._map.off('terrainchange', state.terrain).on('terrainchange', state.terrain);
+  underEntry(state.layer, state.status);
+  state.update();
+}

package/src/traffic/tickTraffic.js

@@ -0,0 +1,57 @@
+import { update as updateScene } from '../scene/update.js';
+import { tripsAt } from './tripsAt.js';
+import { trafficNow } from './trafficNow.js';
+import { failTraffic } from './failTraffic.js';
+import { drawTrafficRoute } from './drawTrafficRoute.js';
+import { rebuildTrafficMesh } from './rebuildTrafficMesh.js';
+import { rebaseTraffic } from './rebaseTraffic.js';
+import { moveTraffic } from './moveTraffic.js';
+import { locateTraffic } from './locateTraffic.js';
+
+export function tickTraffic(state) {
+  state.frame = undefined;
+  if (!state.ready || !state.layer._map) return;
+  if (!state.layer.renderer) { state.frame = requestAnimationFrame(state.tick); return; }
+  if (state.mesh && !(state.mesh.geometry.attributes.position.buffer.usage & GPUBufferUsage.STORAGE)) {
+    state.tripsKey = null;
+    state.nextCheck = -Infinity;
+  }
+  const { live, time } = trafficNow(state), ms = time?.getTime(), valid = Number.isFinite(ms);
+  let changed = false;
+  let trip, x, y, z, dx, dz, lx, ly, lz, light;
+  if (!live || !valid || !(ms >= state.rebasedAt && ms < state.nextCheck)) {
+    let text;
+    try {
+      state.trips = valid ? tripsAt(time, state.lanes, state.velocity) : [];
+      text = !valid ? 'Select a valid simulation date' : !state.loaded ? 'Loading…' :
+        `${state.trips.length} ${state.vehicle === 'boat' ? 'boats' : state.vehicle}`;
+    } catch (error) { return failTraffic(state, error); }
+    state.rebasedAt = ms;
+    state.nextCheck = ms + 1000;
+    const ids = [];
+    for (trip of state.trips) ids.push(trip.id);
+    const key = ids.join();
+    changed = key !== state.tripsKey || (state.outline && !state.route && state.lanes.length);
+    if (state.outline && !state.route && state.lanes.length) drawTrafficRoute(state);
+    if (key !== state.tripsKey) { rebuildTrafficMesh(state); state.tripsKey = key; }
+    if (state.mesh) rebaseTraffic(state);
+    if (!state.error) {
+      if (state.status.textContent !== text) state.status.textContent = text;
+      state.status.removeAttribute('role');
+    }
+  }
+  if (state.mesh) moveTraffic(state, (ms - state.rebasedAt) / 1000);
+  if (state.mesh && state.lights && (!live || !(performance.now() - state.litAt < 100))) {
+    state.litAt = performance.now();
+    const lights = [];
+    for (trip of state.trips) {
+      ({ x, y, z, dx, dz } = locateTraffic(state, trip));
+      for ({ position: [lx, ly, lz], ...light } of state.lights) {
+        lights.push({ ...light, position: [x - lx * dz - lz * dx, y + ly, z + lx * dx - lz * dz] });
+      }
+    }
+    state.mesh.extras.lights = lights;
+  }
+  if (state.mesh || changed) updateScene(state.layer);
+  if (live && valid) state.frame = requestAnimationFrame(state.tick);
+}