@macrostrat/column-views 3.11.0 → 3.12.0

package/dist/age-model/types.cjs

@@ -0,0 +1,229 @@
+"use strict";
+Object.defineProperty(exports, Symbol.toStringTag, { value: "Module" });
+const columnComponents = require("@macrostrat/column-components");
+const surfaceStatuses = [
+  "absolute",
+  "relative",
+  "spike",
+  "imposed",
+  "modeled",
+  "",
+  "derived"
+];
+const TIE_POINT_STATUSES = [
+  "absolute",
+  "relative",
+  "spike",
+  "imposed"
+];
+const surfaceStatusLabels = {
+  absolute: "Absolute",
+  relative: "Relative",
+  spike: "Spike",
+  imposed: "Imposed",
+  modeled: "Modeled",
+  "": "Unspecified",
+  derived: "Derived"
+};
+const surfaceStatusDescriptions = {
+  absolute: "Independently dated; the age is authoritative and the interval position is derived from it",
+  relative: "Calibrated to a position within a geologic interval; the age follows the timescale",
+  spike: "Tied to a formally defined boundary point",
+  imposed: "Interval-relative representation imposed on a boundary with no chronostratigraphic tie of its own (for instance, thickness interpolation)",
+  modeled: "Interpolated by the age model between calibrated surfaces",
+  "": "No status recorded",
+  derived: "Computed from unit tops and bottoms; not stored in the column's age model"
+};
+function surfaceToken(value) {
+  if (value == null || value === "") return "unspecified";
+  return value.replace(/[^a-z0-9]+/gi, "-").toLowerCase();
+}
+function nullifyUnitID(id) {
+  if (id == null || id === 0) return null;
+  return id;
+}
+const INTERVAL_BOUND_TOLERANCE = 1e-3;
+function inferTiePointStatuses(surfaces, tolerance = INTERVAL_BOUND_TOLERANCE) {
+  const extents = sectionAgeExtents(surfaces);
+  return surfaces.map((surface) => {
+    if (surface.status !== "modeled") return surface;
+    if (surface.calibration == null) return surface;
+    const { proportion } = surface;
+    const onIntervalBound = proportion != null && (Math.abs(proportion) <= tolerance || Math.abs(proportion - 1) <= tolerance);
+    if (!onIntervalBound && !boundsPackage(surface, extents)) {
+      return surface;
+    }
+    return { ...surface, status: "relative", statusInferred: true };
+  });
+}
+function boundsPackage(surface, extents) {
+  const extent = extents.get(surface.section_id ?? null);
+  if (extent == null) return false;
+  const [youngest, oldest] = extent;
+  if (surface.unitsAbove.length === 0 && surface.age === youngest) return true;
+  return surface.unitsBelow.length === 0 && surface.age === oldest;
+}
+function sectionAgeExtents(surfaces) {
+  const extents = /* @__PURE__ */ new Map();
+  for (const surface of surfaces) {
+    const key = surface.section_id ?? null;
+    const extent = extents.get(key);
+    if (extent == null) {
+      extents.set(key, [surface.age, surface.age]);
+      continue;
+    }
+    extent[0] = Math.min(extent[0], surface.age);
+    extent[1] = Math.max(extent[1], surface.age);
+  }
+  return extents;
+}
+function boundaryToSurface(boundary) {
+  const above = nullifyUnitID(boundary.unit_above);
+  const below = nullifyUnitID(boundary.unit_below);
+  let calibration = null;
+  if (boundary.interval_id != null) {
+    calibration = {
+      id: boundary.interval_id,
+      name: boundary.interval_name,
+      b_age: boundary.age_bottom,
+      t_age: boundary.age_top
+    };
+  }
+  return {
+    id: boundary.boundary_id,
+    age: boundary.model_age,
+    position: boundary.boundary_position,
+    status: boundary.boundary_status ?? "",
+    type: boundary.boundary_type ?? "",
+    calibration,
+    proportion: boundary.rel_position,
+    unitsAbove: above == null ? [] : [above],
+    unitsBelow: below == null ? [] : [below],
+    section_id: boundary.section_id,
+    ref_id: boundary.ref_id,
+    boundary
+  };
+}
+function surfacesFromBoundaries(boundaries, options = {}) {
+  const { inferTiePoints = true } = options;
+  if (boundaries == null) return [];
+  let surfaces = boundaries.map(boundaryToSurface);
+  if (inferTiePoints) {
+    surfaces = inferTiePointStatuses(surfaces);
+  }
+  return surfaces.sort((a, b) => a.age - b.age);
+}
+function surfacesFromUnits(units, options = {}) {
+  const { tolerance = 1e-3, axisType = columnComponents.ColumnAxisType.AGE } = options;
+  if (units == null) return [];
+  const byPosition = isPositionAxis(axisType);
+  const surfaces = [];
+  function keyOf(surface) {
+    if (byPosition) return surface.position ?? null;
+    return surface.age;
+  }
+  function addBound(unit, top) {
+    const age = top ? unit.t_age : unit.b_age;
+    const position = asNumber(top ? unit.t_pos : unit.b_pos);
+    const key = byPosition ? position : age;
+    if (key == null) return;
+    let surface = surfaces.find((d) => {
+      const k = keyOf(d);
+      return k != null && Math.abs(k - key) < tolerance;
+    });
+    if (surface == null) {
+      surface = {
+        id: "",
+        age,
+        position,
+        status: "derived",
+        type: "",
+        calibration: null,
+        proportion: null,
+        unitsAbove: [],
+        unitsBelow: [],
+        section_id: unit.section_id
+      };
+      surfaces.push(surface);
+    }
+    if (top) {
+      surface.unitsBelow.push(unit.unit_id);
+    } else {
+      surface.unitsAbove.push(unit.unit_id);
+    }
+  }
+  for (const unit of units) {
+    addBound(unit, true);
+    addBound(unit, false);
+  }
+  surfaces.sort((a, b) => compareAlongAxis(a, b, axisType));
+  return surfaces.map((s, i) => ({ ...s, id: `derived-${i}` }));
+}
+function surfacePosition(surface, axisType) {
+  if (isPositionAxis(axisType)) {
+    return surface.position ?? null;
+  }
+  return surface.age ?? null;
+}
+function isPositionAxis(axisType) {
+  return axisType === columnComponents.ColumnAxisType.HEIGHT || axisType === columnComponents.ColumnAxisType.DEPTH;
+}
+function compareAlongAxis(a, b, axisType) {
+  const pa = surfacePosition(a, axisType);
+  const pb = surfacePosition(b, axisType);
+  if (pa == null || pb == null) return 0;
+  if (axisType === columnComponents.ColumnAxisType.HEIGHT) {
+    return pb - pa;
+  }
+  return pa - pb;
+}
+function formatProportion(proportion) {
+  if (proportion == null) return null;
+  return `${Math.round(proportion * 100)}%`;
+}
+function surfaceLabel(surface) {
+  const { calibration } = surface;
+  if (calibration != null) {
+    const prop = formatProportion(surface.proportion);
+    if (prop != null) {
+      return `${calibration.name} · ${prop}`;
+    }
+    return calibration.name;
+  }
+  return formatAge(surface.age);
+}
+function formatAge(age) {
+  if (age == null) return "";
+  let value = age;
+  let unit = "Ma";
+  if (Math.abs(age) < 1) {
+    value = age * 1e3;
+    unit = "ka";
+  }
+  const _value = value.toLocaleString("en-US", { maximumFractionDigits: 2 });
+  return `${_value} ${unit}`;
+}
+function asNumber(value) {
+  if (value == null || value === "") return null;
+  const n = Number(value);
+  if (Number.isNaN(n)) return null;
+  return n;
+}
+exports.INTERVAL_BOUND_TOLERANCE = INTERVAL_BOUND_TOLERANCE;
+exports.TIE_POINT_STATUSES = TIE_POINT_STATUSES;
+exports.boundaryToSurface = boundaryToSurface;
+exports.compareAlongAxis = compareAlongAxis;
+exports.formatAge = formatAge;
+exports.formatProportion = formatProportion;
+exports.inferTiePointStatuses = inferTiePointStatuses;
+exports.isPositionAxis = isPositionAxis;
+exports.nullifyUnitID = nullifyUnitID;
+exports.surfaceLabel = surfaceLabel;
+exports.surfacePosition = surfacePosition;
+exports.surfaceStatusDescriptions = surfaceStatusDescriptions;
+exports.surfaceStatusLabels = surfaceStatusLabels;
+exports.surfaceStatuses = surfaceStatuses;
+exports.surfaceToken = surfaceToken;
+exports.surfacesFromBoundaries = surfacesFromBoundaries;
+exports.surfacesFromUnits = surfacesFromUnits;
+//# sourceMappingURL=types.cjs.map

package/dist/age-model/types.cjs.map

@@ -0,0 +1 @@
+{"version":3,"file":"types.cjs","sources":["../../src/age-model/types.ts"],"sourcesContent":["/** Surfaces: the view model for a column's age model.\n *\n * A `ColumnSurface` is one horizontal boundary in a column. It is what the\n * surfaces view renders, and it is deliberately data-agnostic: a surface can\n * come from an `/age_model` boundary (`boundaryToSurface`), be derived from\n * unit tops and bottoms when a column has no age-model rows\n * (`surfacesFromUnits`), or be built by an editor from its own state.\n */\nimport type {\n  AgeModelBoundary,\n  BoundaryStatus,\n  BoundaryType,\n  UnitLong,\n} from \"@macrostrat/api-types\";\nimport { ColumnAxisType } from \"@macrostrat/column-components\";\n\nexport type { AgeModelBoundary, BoundaryStatus, BoundaryType };\n\n/** A boundary status, plus `derived` for surfaces computed from units rather\n * than stored in the age model. */\nexport type SurfaceStatus = BoundaryStatus | \"derived\";\n\n/** The interval a surface is calibrated against */\nexport interface SurfaceCalibration {\n  id: number;\n  name: string;\n  b_age: number;\n  t_age: number;\n}\n\nexport interface ColumnSurface {\n  id: number | string;\n  /** Modeled age of the surface (Ma) */\n  age: number;\n  /** Measured position (height or depth) along the column, when known */\n  position?: number | null;\n  status: SurfaceStatus;\n  type: BoundaryType;\n  /** The calibration interval, and the surface's position within it\n   * (0 at the interval's base, 1 at its top) */\n  calibration?: SurfaceCalibration | null;\n  proportion?: number | null;\n  /** True when `status` was inferred rather than read from the data — see\n   * `inferTiePointStatuses`. */\n  statusInferred?: boolean;\n  unitsAbove: number[];\n  unitsBelow: number[];\n  section_id?: number | null;\n  ref_id?: number | null;\n  /** The API record, when the surface came from `/age_model` */\n  boundary?: AgeModelBoundary;\n}\n\nexport const surfaceStatuses: SurfaceStatus[] = [\n  \"absolute\",\n  \"relative\",\n  \"spike\",\n  \"imposed\",\n  \"modeled\",\n  \"\",\n  \"derived\",\n];\n\n/** The statuses at which the age model was actually constrained — the tie\n * points — as opposed to surfaces whose age falls out of the model\n * (`modeled`), carries no status, or is merely derived from units. These are\n * the surfaces worth labeling by default. */\nexport const TIE_POINT_STATUSES: SurfaceStatus[] = [\n  \"absolute\",\n  \"relative\",\n  \"spike\",\n  \"imposed\",\n];\n\nexport const surfaceStatusLabels: Record<SurfaceStatus, string> = {\n  absolute: \"Absolute\",\n  relative: \"Relative\",\n  spike: \"Spike\",\n  imposed: \"Imposed\",\n  modeled: \"Modeled\",\n  \"\": \"Unspecified\",\n  derived: \"Derived\",\n};\n\nexport const surfaceStatusDescriptions: Record<SurfaceStatus, string> = {\n  absolute:\n    \"Independently dated; the age is authoritative and the interval position is derived from it\",\n  relative:\n    \"Calibrated to a position within a geologic interval; the age follows the timescale\",\n  spike: \"Tied to a formally defined boundary point\",\n  imposed:\n    \"Interval-relative representation imposed on a boundary with no chronostratigraphic tie of its own (for instance, thickness interpolation)\",\n  modeled: \"Interpolated by the age model between calibrated surfaces\",\n  \"\": \"No status recorded\",\n  derived:\n    \"Computed from unit tops and bottoms; not stored in the column's age model\",\n};\n\n/** A CSS-safe token for a status or type (`angular unconformity` →\n * `angular-unconformity`, `\"\"` → `unspecified`). */\nexport function surfaceToken(value: string | null | undefined): string {\n  if (value == null || value === \"\") return \"unspecified\";\n  return value.replace(/[^a-z0-9]+/gi, \"-\").toLowerCase();\n}\n\n/** Legacy `unit_boundaries` rows use `0` as the \"no unit\" sentinel. */\nexport function nullifyUnitID(id: number | null | undefined): number | null {\n  if (id == null || id === 0) return null;\n  return id;\n}\n\n/** How close to an interval's base or top a surface must sit to count as\n * bounding it (as a fraction of the interval). */\nexport const INTERVAL_BOUND_TOLERANCE = 0.001;\n\n/** Promote the `modeled` surfaces that cannot in fact have been interpolated\n * to `relative`.\n *\n * This is a client-side affordance over a known gap in the data. An age model\n * interpolates *between* tie points, so a surface with no constraint to\n * interpolate from has to be one itself, whatever `boundary_status` says. Two\n * cases are recognizable from the boundaries alone:\n *\n * - **It sits on the base or top of its calibration interval** (proportion 0\n *   or 1). A height referred directly to an interval bounds the model rather\n *   than falling out of it.\n * - **It is the edge of a gap-bound package** — the youngest surface in its\n *   section with no unit above, or the oldest with none below. There is\n *   nothing beyond it to interpolate against.\n *\n * Both are common in columns whose age model was built by assigning positions\n * within intervals, where the importer recorded everything as `modeled`;\n * filtering to the tie-point statuses would otherwise hide real ones. Promoted\n * surfaces are marked `statusInferred`, so views can say where the status came\n * from; nothing about the underlying record changes.\n *\n * Two cases are deliberately left alone. A *mid-interval* position is\n * ambiguous: an assigned position and an interpolated one are the same record\n * there, and only the data can settle it. And an open side in the *interior*\n * of a section is a lateral relationship, not a package edge — a unit that\n * pinches out while others span past it — which the model can interpolate\n * normally.\n */\nexport function inferTiePointStatuses(\n  surfaces: ColumnSurface[],\n  tolerance: number = INTERVAL_BOUND_TOLERANCE,\n): ColumnSurface[] {\n  const extents = sectionAgeExtents(surfaces);\n\n  return surfaces.map((surface) => {\n    if (surface.status !== \"modeled\") return surface;\n    if (surface.calibration == null) return surface;\n\n    const { proportion } = surface;\n    const onIntervalBound =\n      proportion != null &&\n      (Math.abs(proportion) <= tolerance ||\n        Math.abs(proportion - 1) <= tolerance);\n\n    if (!onIntervalBound && !boundsPackage(surface, extents)) {\n      return surface;\n    }\n    return { ...surface, status: \"relative\", statusInferred: true };\n  });\n}\n\n/** An open side only means a package edge at the ends of the section. In its\n * interior it means a unit pinching out while others span past it, which the\n * age model interpolates like any other surface. */\nfunction boundsPackage(\n  surface: ColumnSurface,\n  extents: Map<ColumnSurface[\"section_id\"], [number, number]>,\n): boolean {\n  const extent = extents.get(surface.section_id ?? null);\n  if (extent == null) return false;\n  const [youngest, oldest] = extent;\n  if (surface.unitsAbove.length === 0 && surface.age === youngest) return true;\n  return surface.unitsBelow.length === 0 && surface.age === oldest;\n}\n\n/** The youngest and oldest surface age in each section. */\nfunction sectionAgeExtents(\n  surfaces: ColumnSurface[],\n): Map<ColumnSurface[\"section_id\"], [number, number]> {\n  const extents = new Map<ColumnSurface[\"section_id\"], [number, number]>();\n  for (const surface of surfaces) {\n    const key = surface.section_id ?? null;\n    const extent = extents.get(key);\n    if (extent == null) {\n      extents.set(key, [surface.age, surface.age]);\n      continue;\n    }\n    extent[0] = Math.min(extent[0], surface.age);\n    extent[1] = Math.max(extent[1], surface.age);\n  }\n  return extents;\n}\n\nexport function boundaryToSurface(boundary: AgeModelBoundary): ColumnSurface {\n  const above = nullifyUnitID(boundary.unit_above);\n  const below = nullifyUnitID(boundary.unit_below);\n  let calibration: SurfaceCalibration | null = null;\n  if (boundary.interval_id != null) {\n    calibration = {\n      id: boundary.interval_id,\n      name: boundary.interval_name,\n      b_age: boundary.age_bottom,\n      t_age: boundary.age_top,\n    };\n  }\n  return {\n    id: boundary.boundary_id,\n    age: boundary.model_age,\n    position: boundary.boundary_position,\n    status: boundary.boundary_status ?? \"\",\n    type: boundary.boundary_type ?? \"\",\n    calibration,\n    proportion: boundary.rel_position,\n    unitsAbove: above == null ? [] : [above],\n    unitsBelow: below == null ? [] : [below],\n    section_id: boundary.section_id,\n    ref_id: boundary.ref_id,\n    boundary,\n  };\n}\n\nexport interface SurfacesFromBoundariesOptions {\n  /** Promote `modeled` surfaces that must in fact constrain the age model to\n   * `relative` (default `true`) — see `inferTiePointStatuses`. */\n  inferTiePoints?: boolean;\n}\n\n/** Surfaces for a column's age-model boundaries, youngest first. */\nexport function surfacesFromBoundaries(\n  boundaries: AgeModelBoundary[] | null | undefined,\n  options: SurfacesFromBoundariesOptions = {},\n): ColumnSurface[] {\n  const { inferTiePoints = true } = options;\n  if (boundaries == null) return [];\n  let surfaces = boundaries.map(boundaryToSurface);\n  if (inferTiePoints) {\n    surfaces = inferTiePointStatuses(surfaces);\n  }\n  return surfaces.sort((a, b) => a.age - b.age);\n}\n\nexport interface SurfacesFromUnitsOptions {\n  /** Surfaces closer than this (in the axis's units) are merged */\n  tolerance?: number;\n  /** Which unit field the surfaces are keyed on: ages for age columns,\n   * measured positions for height and depth columns */\n  axisType?: ColumnAxisType;\n}\n\n/** Derive surfaces from unit tops and bottoms. Every distinct top or bottom\n * becomes one surface with status `derived`; units sharing a boundary share\n * the surface. This is the fallback when a column has no `unit_boundaries`. */\nexport function surfacesFromUnits<T extends UnitLong>(\n  units: T[] | null | undefined,\n  options: SurfacesFromUnitsOptions = {},\n): ColumnSurface[] {\n  const { tolerance = 0.001, axisType = ColumnAxisType.AGE } = options;\n  if (units == null) return [];\n  const byPosition = isPositionAxis(axisType);\n\n  const surfaces: ColumnSurface[] = [];\n\n  function keyOf(surface: ColumnSurface): number | null {\n    if (byPosition) return surface.position ?? null;\n    return surface.age;\n  }\n\n  function addBound(unit: T, top: boolean) {\n    const age = top ? unit.t_age : unit.b_age;\n    const position = asNumber(top ? unit.t_pos : unit.b_pos);\n    const key = byPosition ? position : age;\n    if (key == null) return;\n\n    let surface = surfaces.find((d) => {\n      const k = keyOf(d);\n      return k != null && Math.abs(k - key) < tolerance;\n    });\n    if (surface == null) {\n      surface = {\n        id: \"\",\n        age,\n        position,\n        status: \"derived\",\n        type: \"\",\n        calibration: null,\n        proportion: null,\n        unitsAbove: [],\n        unitsBelow: [],\n        section_id: unit.section_id,\n      };\n      surfaces.push(surface);\n    }\n    // The unit whose top this is lies *below* the surface, and vice versa\n    if (top) {\n      surface.unitsBelow.push(unit.unit_id);\n    } else {\n      surface.unitsAbove.push(unit.unit_id);\n    }\n  }\n\n  for (const unit of units) {\n    addBound(unit, true);\n    addBound(unit, false);\n  }\n\n  surfaces.sort((a, b) => compareAlongAxis(a, b, axisType));\n  return surfaces.map((s, i) => ({ ...s, id: `derived-${i}` }));\n}\n\n/** The coordinate a surface is drawn at for a given axis type: its age for\n * age columns, its measured position for height and depth columns. `null`\n * when the surface has no coordinate on that axis. */\nexport function surfacePosition(\n  surface: ColumnSurface,\n  axisType: ColumnAxisType | null | undefined,\n): number | null {\n  if (isPositionAxis(axisType)) {\n    return surface.position ?? null;\n  }\n  return surface.age ?? null;\n}\n\nexport function isPositionAxis(\n  axisType: ColumnAxisType | null | undefined,\n): boolean {\n  return (\n    axisType === ColumnAxisType.HEIGHT || axisType === ColumnAxisType.DEPTH\n  );\n}\n\n/** Order surfaces from the top of the column downwards. */\nexport function compareAlongAxis(\n  a: ColumnSurface,\n  b: ColumnSurface,\n  axisType: ColumnAxisType | null | undefined,\n): number {\n  const pa = surfacePosition(a, axisType);\n  const pb = surfacePosition(b, axisType);\n  if (pa == null || pb == null) return 0;\n  if (axisType === ColumnAxisType.HEIGHT) {\n    // Heights increase upwards\n    return pb - pa;\n  }\n  // Ages and depths increase downwards\n  return pa - pb;\n}\n\n/** Percent of the way up its calibration interval, e.g. `85%`. */\nexport function formatProportion(\n  proportion: number | null | undefined,\n): string | null {\n  if (proportion == null) return null;\n  return `${Math.round(proportion * 100)}%`;\n}\n\n/** A short text label for a surface: its calibration (\"Calymmian · 85%\")\n * or, failing that, its age. */\nexport function surfaceLabel(surface: ColumnSurface): string {\n  const { calibration } = surface;\n  if (calibration != null) {\n    const prop = formatProportion(surface.proportion);\n    if (prop != null) {\n      return `${calibration.name} · ${prop}`;\n    }\n    return calibration.name;\n  }\n  return formatAge(surface.age);\n}\n\nexport function formatAge(age: number | null | undefined): string {\n  if (age == null) return \"\";\n  let value = age;\n  let unit = \"Ma\";\n  if (Math.abs(age) < 1) {\n    value = age * 1000;\n    unit = \"ka\";\n  }\n  const _value = value.toLocaleString(\"en-US\", { maximumFractionDigits: 2 });\n  return `${_value} ${unit}`;\n}\n\nfunction asNumber(value: unknown): number | null {\n  if (value == null || value === \"\") return null;\n  const n = Number(value);\n  if (Number.isNaN(n)) return null;\n  return n;\n}\n"],"names":["ColumnAxisType"],"mappings":";;;AAqDO,MAAM,kBAAmC;AAAA,EAC9C;AAAA,EACA;AAAA,EACA;AAAA,EACA;AAAA,EACA;AAAA,EACA;AAAA,EACA;AACF;AAMO,MAAM,qBAAsC;AAAA,EACjD;AAAA,EACA;AAAA,EACA;AAAA,EACA;AACF;AAEO,MAAM,sBAAqD;AAAA,EAChE,UAAU;AAAA,EACV,UAAU;AAAA,EACV,OAAO;AAAA,EACP,SAAS;AAAA,EACT,SAAS;AAAA,EACT,IAAI;AAAA,EACJ,SAAS;AACX;AAEO,MAAM,4BAA2D;AAAA,EACtE,UACE;AAAA,EACF,UACE;AAAA,EACF,OAAO;AAAA,EACP,SACE;AAAA,EACF,SAAS;AAAA,EACT,IAAI;AAAA,EACJ,SACE;AACJ;AAIO,SAAS,aAAa,OAA0C;AACrE,MAAI,SAAS,QAAQ,UAAU,GAAI,QAAO;AAC1C,SAAO,MAAM,QAAQ,gBAAgB,GAAG,EAAE,YAAA;AAC5C;AAGO,SAAS,cAAc,IAA8C;AAC1E,MAAI,MAAM,QAAQ,OAAO,EAAG,QAAO;AACnC,SAAO;AACT;AAIO,MAAM,2BAA2B;AA8BjC,SAAS,sBACd,UACA,YAAoB,0BACH;AACjB,QAAM,UAAU,kBAAkB,QAAQ;AAE1C,SAAO,SAAS,IAAI,CAAC,YAAY;AAC/B,QAAI,QAAQ,WAAW,UAAW,QAAO;AACzC,QAAI,QAAQ,eAAe,KAAM,QAAO;AAExC,UAAM,EAAE,eAAe;AACvB,UAAM,kBACJ,cAAc,SACb,KAAK,IAAI,UAAU,KAAK,aACvB,KAAK,IAAI,aAAa,CAAC,KAAK;AAEhC,QAAI,CAAC,mBAAmB,CAAC,cAAc,SAAS,OAAO,GAAG;AACxD,aAAO;AAAA,IACT;AACA,WAAO,EAAE,GAAG,SAAS,QAAQ,YAAY,gBAAgB,KAAA;AAAA,EAC3D,CAAC;AACH;AAKA,SAAS,cACP,SACA,SACS;AACT,QAAM,SAAS,QAAQ,IAAI,QAAQ,cAAc,IAAI;AACrD,MAAI,UAAU,KAAM,QAAO;AAC3B,QAAM,CAAC,UAAU,MAAM,IAAI;AAC3B,MAAI,QAAQ,WAAW,WAAW,KAAK,QAAQ,QAAQ,SAAU,QAAO;AACxE,SAAO,QAAQ,WAAW,WAAW,KAAK,QAAQ,QAAQ;AAC5D;AAGA,SAAS,kBACP,UACoD;AACpD,QAAM,8BAAc,IAAA;AACpB,aAAW,WAAW,UAAU;AAC9B,UAAM,MAAM,QAAQ,cAAc;AAClC,UAAM,SAAS,QAAQ,IAAI,GAAG;AAC9B,QAAI,UAAU,MAAM;AAClB,cAAQ,IAAI,KAAK,CAAC,QAAQ,KAAK,QAAQ,GAAG,CAAC;AAC3C;AAAA,IACF;AACA,WAAO,CAAC,IAAI,KAAK,IAAI,OAAO,CAAC,GAAG,QAAQ,GAAG;AAC3C,WAAO,CAAC,IAAI,KAAK,IAAI,OAAO,CAAC,GAAG,QAAQ,GAAG;AAAA,EAC7C;AACA,SAAO;AACT;AAEO,SAAS,kBAAkB,UAA2C;AAC3E,QAAM,QAAQ,cAAc,SAAS,UAAU;AAC/C,QAAM,QAAQ,cAAc,SAAS,UAAU;AAC/C,MAAI,cAAyC;AAC7C,MAAI,SAAS,eAAe,MAAM;AAChC,kBAAc;AAAA,MACZ,IAAI,SAAS;AAAA,MACb,MAAM,SAAS;AAAA,MACf,OAAO,SAAS;AAAA,MAChB,OAAO,SAAS;AAAA,IAAA;AAAA,EAEpB;AACA,SAAO;AAAA,IACL,IAAI,SAAS;AAAA,IACb,KAAK,SAAS;AAAA,IACd,UAAU,SAAS;AAAA,IACnB,QAAQ,SAAS,mBAAmB;AAAA,IACpC,MAAM,SAAS,iBAAiB;AAAA,IAChC;AAAA,IACA,YAAY,SAAS;AAAA,IACrB,YAAY,SAAS,OAAO,CAAA,IAAK,CAAC,KAAK;AAAA,IACvC,YAAY,SAAS,OAAO,CAAA,IAAK,CAAC,KAAK;AAAA,IACvC,YAAY,SAAS;AAAA,IACrB,QAAQ,SAAS;AAAA,IACjB;AAAA,EAAA;AAEJ;AASO,SAAS,uBACd,YACA,UAAyC,IACxB;AACjB,QAAM,EAAE,iBAAiB,KAAA,IAAS;AAClC,MAAI,cAAc,KAAM,QAAO,CAAA;AAC/B,MAAI,WAAW,WAAW,IAAI,iBAAiB;AAC/C,MAAI,gBAAgB;AAClB,eAAW,sBAAsB,QAAQ;AAAA,EAC3C;AACA,SAAO,SAAS,KAAK,CAAC,GAAG,MAAM,EAAE,MAAM,EAAE,GAAG;AAC9C;AAaO,SAAS,kBACd,OACA,UAAoC,IACnB;AACjB,QAAM,EAAE,YAAY,MAAO,WAAWA,iBAAAA,eAAe,QAAQ;AAC7D,MAAI,SAAS,KAAM,QAAO,CAAA;AAC1B,QAAM,aAAa,eAAe,QAAQ;AAE1C,QAAM,WAA4B,CAAA;AAElC,WAAS,MAAM,SAAuC;AACpD,QAAI,WAAY,QAAO,QAAQ,YAAY;AAC3C,WAAO,QAAQ;AAAA,EACjB;AAEA,WAAS,SAAS,MAAS,KAAc;AACvC,UAAM,MAAM,MAAM,KAAK,QAAQ,KAAK;AACpC,UAAM,WAAW,SAAS,MAAM,KAAK,QAAQ,KAAK,KAAK;AACvD,UAAM,MAAM,aAAa,WAAW;AACpC,QAAI,OAAO,KAAM;AAEjB,QAAI,UAAU,SAAS,KAAK,CAAC,MAAM;AACjC,YAAM,IAAI,MAAM,CAAC;AACjB,aAAO,KAAK,QAAQ,KAAK,IAAI,IAAI,GAAG,IAAI;AAAA,IAC1C,CAAC;AACD,QAAI,WAAW,MAAM;AACnB,gBAAU;AAAA,QACR,IAAI;AAAA,QACJ;AAAA,QACA;AAAA,QACA,QAAQ;AAAA,QACR,MAAM;AAAA,QACN,aAAa;AAAA,QACb,YAAY;AAAA,QACZ,YAAY,CAAA;AAAA,QACZ,YAAY,CAAA;AAAA,QACZ,YAAY,KAAK;AAAA,MAAA;AAEnB,eAAS,KAAK,OAAO;AAAA,IACvB;AAEA,QAAI,KAAK;AACP,cAAQ,WAAW,KAAK,KAAK,OAAO;AAAA,IACtC,OAAO;AACL,cAAQ,WAAW,KAAK,KAAK,OAAO;AAAA,IACtC;AAAA,EACF;AAEA,aAAW,QAAQ,OAAO;AACxB,aAAS,MAAM,IAAI;AACnB,aAAS,MAAM,KAAK;AAAA,EACtB;AAEA,WAAS,KAAK,CAAC,GAAG,MAAM,iBAAiB,GAAG,GAAG,QAAQ,CAAC;AACxD,SAAO,SAAS,IAAI,CAAC,GAAG,OAAO,EAAE,GAAG,GAAG,IAAI,WAAW,CAAC,GAAA,EAAK;AAC9D;AAKO,SAAS,gBACd,SACA,UACe;AACf,MAAI,eAAe,QAAQ,GAAG;AAC5B,WAAO,QAAQ,YAAY;AAAA,EAC7B;AACA,SAAO,QAAQ,OAAO;AACxB;AAEO,SAAS,eACd,UACS;AACT,SACE,aAAaA,iBAAAA,eAAe,UAAU,aAAaA,iBAAAA,eAAe;AAEtE;AAGO,SAAS,iBACd,GACA,GACA,UACQ;AACR,QAAM,KAAK,gBAAgB,GAAG,QAAQ;AACtC,QAAM,KAAK,gBAAgB,GAAG,QAAQ;AACtC,MAAI,MAAM,QAAQ,MAAM,KAAM,QAAO;AACrC,MAAI,aAAaA,iBAAAA,eAAe,QAAQ;AAEtC,WAAO,KAAK;AAAA,EACd;AAEA,SAAO,KAAK;AACd;AAGO,SAAS,iBACd,YACe;AACf,MAAI,cAAc,KAAM,QAAO;AAC/B,SAAO,GAAG,KAAK,MAAM,aAAa,GAAG,CAAC;AACxC;AAIO,SAAS,aAAa,SAAgC;AAC3D,QAAM,EAAE,gBAAgB;AACxB,MAAI,eAAe,MAAM;AACvB,UAAM,OAAO,iBAAiB,QAAQ,UAAU;AAChD,QAAI,QAAQ,MAAM;AAChB,aAAO,GAAG,YAAY,IAAI,MAAM,IAAI;AAAA,IACtC;AACA,WAAO,YAAY;AAAA,EACrB;AACA,SAAO,UAAU,QAAQ,GAAG;AAC9B;AAEO,SAAS,UAAU,KAAwC;AAChE,MAAI,OAAO,KAAM,QAAO;AACxB,MAAI,QAAQ;AACZ,MAAI,OAAO;AACX,MAAI,KAAK,IAAI,GAAG,IAAI,GAAG;AACrB,YAAQ,MAAM;AACd,WAAO;AAAA,EACT;AACA,QAAM,SAAS,MAAM,eAAe,SAAS,EAAE,uBAAuB,GAAG;AACzE,SAAO,GAAG,MAAM,IAAI,IAAI;AAC1B;AAEA,SAAS,SAAS,OAA+B;AAC/C,MAAI,SAAS,QAAQ,UAAU,GAAI,QAAO;AAC1C,QAAM,IAAI,OAAO,KAAK;AACtB,MAAI,OAAO,MAAM,CAAC,EAAG,QAAO;AAC5B,SAAO;AACT;;;;;;;;;;;;;;;;;;"}

package/dist/age-model/types.d.ts

@@ -0,0 +1,112 @@
+import { AgeModelBoundary, BoundaryStatus, BoundaryType, UnitLong } from '@macrostrat/api-types';
+import { ColumnAxisType } from '@macrostrat/column-components';
+export type { AgeModelBoundary, BoundaryStatus, BoundaryType };
+/** A boundary status, plus `derived` for surfaces computed from units rather
+ * than stored in the age model. */
+export type SurfaceStatus = BoundaryStatus | "derived";
+/** The interval a surface is calibrated against */
+export interface SurfaceCalibration {
+    id: number;
+    name: string;
+    b_age: number;
+    t_age: number;
+}
+export interface ColumnSurface {
+    id: number | string;
+    /** Modeled age of the surface (Ma) */
+    age: number;
+    /** Measured position (height or depth) along the column, when known */
+    position?: number | null;
+    status: SurfaceStatus;
+    type: BoundaryType;
+    /** The calibration interval, and the surface's position within it
+     * (0 at the interval's base, 1 at its top) */
+    calibration?: SurfaceCalibration | null;
+    proportion?: number | null;
+    /** True when `status` was inferred rather than read from the data — see
+     * `inferTiePointStatuses`. */
+    statusInferred?: boolean;
+    unitsAbove: number[];
+    unitsBelow: number[];
+    section_id?: number | null;
+    ref_id?: number | null;
+    /** The API record, when the surface came from `/age_model` */
+    boundary?: AgeModelBoundary;
+}
+export declare const surfaceStatuses: SurfaceStatus[];
+/** The statuses at which the age model was actually constrained — the tie
+ * points — as opposed to surfaces whose age falls out of the model
+ * (`modeled`), carries no status, or is merely derived from units. These are
+ * the surfaces worth labeling by default. */
+export declare const TIE_POINT_STATUSES: SurfaceStatus[];
+export declare const surfaceStatusLabels: Record<SurfaceStatus, string>;
+export declare const surfaceStatusDescriptions: Record<SurfaceStatus, string>;
+/** A CSS-safe token for a status or type (`angular unconformity` →
+ * `angular-unconformity`, `""` → `unspecified`). */
+export declare function surfaceToken(value: string | null | undefined): string;
+/** Legacy `unit_boundaries` rows use `0` as the "no unit" sentinel. */
+export declare function nullifyUnitID(id: number | null | undefined): number | null;
+/** How close to an interval's base or top a surface must sit to count as
+ * bounding it (as a fraction of the interval). */
+export declare const INTERVAL_BOUND_TOLERANCE = 0.001;
+/** Promote the `modeled` surfaces that cannot in fact have been interpolated
+ * to `relative`.
+ *
+ * This is a client-side affordance over a known gap in the data. An age model
+ * interpolates *between* tie points, so a surface with no constraint to
+ * interpolate from has to be one itself, whatever `boundary_status` says. Two
+ * cases are recognizable from the boundaries alone:
+ *
+ * - **It sits on the base or top of its calibration interval** (proportion 0
+ *   or 1). A height referred directly to an interval bounds the model rather
+ *   than falling out of it.
+ * - **It is the edge of a gap-bound package** — the youngest surface in its
+ *   section with no unit above, or the oldest with none below. There is
+ *   nothing beyond it to interpolate against.
+ *
+ * Both are common in columns whose age model was built by assigning positions
+ * within intervals, where the importer recorded everything as `modeled`;
+ * filtering to the tie-point statuses would otherwise hide real ones. Promoted
+ * surfaces are marked `statusInferred`, so views can say where the status came
+ * from; nothing about the underlying record changes.
+ *
+ * Two cases are deliberately left alone. A *mid-interval* position is
+ * ambiguous: an assigned position and an interpolated one are the same record
+ * there, and only the data can settle it. And an open side in the *interior*
+ * of a section is a lateral relationship, not a package edge — a unit that
+ * pinches out while others span past it — which the model can interpolate
+ * normally.
+ */
+export declare function inferTiePointStatuses(surfaces: ColumnSurface[], tolerance?: number): ColumnSurface[];
+export declare function boundaryToSurface(boundary: AgeModelBoundary): ColumnSurface;
+export interface SurfacesFromBoundariesOptions {
+    /** Promote `modeled` surfaces that must in fact constrain the age model to
+     * `relative` (default `true`) — see `inferTiePointStatuses`. */
+    inferTiePoints?: boolean;
+}
+/** Surfaces for a column's age-model boundaries, youngest first. */
+export declare function surfacesFromBoundaries(boundaries: AgeModelBoundary[] | null | undefined, options?: SurfacesFromBoundariesOptions): ColumnSurface[];
+export interface SurfacesFromUnitsOptions {
+    /** Surfaces closer than this (in the axis's units) are merged */
+    tolerance?: number;
+    /** Which unit field the surfaces are keyed on: ages for age columns,
+     * measured positions for height and depth columns */
+    axisType?: ColumnAxisType;
+}
+/** Derive surfaces from unit tops and bottoms. Every distinct top or bottom
+ * becomes one surface with status `derived`; units sharing a boundary share
+ * the surface. This is the fallback when a column has no `unit_boundaries`. */
+export declare function surfacesFromUnits<T extends UnitLong>(units: T[] | null | undefined, options?: SurfacesFromUnitsOptions): ColumnSurface[];
+/** The coordinate a surface is drawn at for a given axis type: its age for
+ * age columns, its measured position for height and depth columns. `null`
+ * when the surface has no coordinate on that axis. */
+export declare function surfacePosition(surface: ColumnSurface, axisType: ColumnAxisType | null | undefined): number | null;
+export declare function isPositionAxis(axisType: ColumnAxisType | null | undefined): boolean;
+/** Order surfaces from the top of the column downwards. */
+export declare function compareAlongAxis(a: ColumnSurface, b: ColumnSurface, axisType: ColumnAxisType | null | undefined): number;
+/** Percent of the way up its calibration interval, e.g. `85%`. */
+export declare function formatProportion(proportion: number | null | undefined): string | null;
+/** A short text label for a surface: its calibration ("Calymmian · 85%")
+ * or, failing that, its age. */
+export declare function surfaceLabel(surface: ColumnSurface): string;
+export declare function formatAge(age: number | null | undefined): string;

package/dist/age-model/types.js

@@ -0,0 +1,229 @@
+import { ColumnAxisType } from "@macrostrat/column-components";
+const surfaceStatuses = [
+  "absolute",
+  "relative",
+  "spike",
+  "imposed",
+  "modeled",
+  "",
+  "derived"
+];
+const TIE_POINT_STATUSES = [
+  "absolute",
+  "relative",
+  "spike",
+  "imposed"
+];
+const surfaceStatusLabels = {
+  absolute: "Absolute",
+  relative: "Relative",
+  spike: "Spike",
+  imposed: "Imposed",
+  modeled: "Modeled",
+  "": "Unspecified",
+  derived: "Derived"
+};
+const surfaceStatusDescriptions = {
+  absolute: "Independently dated; the age is authoritative and the interval position is derived from it",
+  relative: "Calibrated to a position within a geologic interval; the age follows the timescale",
+  spike: "Tied to a formally defined boundary point",
+  imposed: "Interval-relative representation imposed on a boundary with no chronostratigraphic tie of its own (for instance, thickness interpolation)",
+  modeled: "Interpolated by the age model between calibrated surfaces",
+  "": "No status recorded",
+  derived: "Computed from unit tops and bottoms; not stored in the column's age model"
+};
+function surfaceToken(value) {
+  if (value == null || value === "") return "unspecified";
+  return value.replace(/[^a-z0-9]+/gi, "-").toLowerCase();
+}
+function nullifyUnitID(id) {
+  if (id == null || id === 0) return null;
+  return id;
+}
+const INTERVAL_BOUND_TOLERANCE = 1e-3;
+function inferTiePointStatuses(surfaces, tolerance = INTERVAL_BOUND_TOLERANCE) {
+  const extents = sectionAgeExtents(surfaces);
+  return surfaces.map((surface) => {
+    if (surface.status !== "modeled") return surface;
+    if (surface.calibration == null) return surface;
+    const { proportion } = surface;
+    const onIntervalBound = proportion != null && (Math.abs(proportion) <= tolerance || Math.abs(proportion - 1) <= tolerance);
+    if (!onIntervalBound && !boundsPackage(surface, extents)) {
+      return surface;
+    }
+    return { ...surface, status: "relative", statusInferred: true };
+  });
+}
+function boundsPackage(surface, extents) {
+  const extent = extents.get(surface.section_id ?? null);
+  if (extent == null) return false;
+  const [youngest, oldest] = extent;
+  if (surface.unitsAbove.length === 0 && surface.age === youngest) return true;
+  return surface.unitsBelow.length === 0 && surface.age === oldest;
+}
+function sectionAgeExtents(surfaces) {
+  const extents = /* @__PURE__ */ new Map();
+  for (const surface of surfaces) {
+    const key = surface.section_id ?? null;
+    const extent = extents.get(key);
+    if (extent == null) {
+      extents.set(key, [surface.age, surface.age]);
+      continue;
+    }
+    extent[0] = Math.min(extent[0], surface.age);
+    extent[1] = Math.max(extent[1], surface.age);
+  }
+  return extents;
+}
+function boundaryToSurface(boundary) {
+  const above = nullifyUnitID(boundary.unit_above);
+  const below = nullifyUnitID(boundary.unit_below);
+  let calibration = null;
+  if (boundary.interval_id != null) {
+    calibration = {
+      id: boundary.interval_id,
+      name: boundary.interval_name,
+      b_age: boundary.age_bottom,
+      t_age: boundary.age_top
+    };
+  }
+  return {
+    id: boundary.boundary_id,
+    age: boundary.model_age,
+    position: boundary.boundary_position,
+    status: boundary.boundary_status ?? "",
+    type: boundary.boundary_type ?? "",
+    calibration,
+    proportion: boundary.rel_position,
+    unitsAbove: above == null ? [] : [above],
+    unitsBelow: below == null ? [] : [below],
+    section_id: boundary.section_id,
+    ref_id: boundary.ref_id,
+    boundary
+  };
+}
+function surfacesFromBoundaries(boundaries, options = {}) {
+  const { inferTiePoints = true } = options;
+  if (boundaries == null) return [];
+  let surfaces = boundaries.map(boundaryToSurface);
+  if (inferTiePoints) {
+    surfaces = inferTiePointStatuses(surfaces);
+  }
+  return surfaces.sort((a, b) => a.age - b.age);
+}
+function surfacesFromUnits(units, options = {}) {
+  const { tolerance = 1e-3, axisType = ColumnAxisType.AGE } = options;
+  if (units == null) return [];
+  const byPosition = isPositionAxis(axisType);
+  const surfaces = [];
+  function keyOf(surface) {
+    if (byPosition) return surface.position ?? null;
+    return surface.age;
+  }
+  function addBound(unit, top) {
+    const age = top ? unit.t_age : unit.b_age;
+    const position = asNumber(top ? unit.t_pos : unit.b_pos);
+    const key = byPosition ? position : age;
+    if (key == null) return;
+    let surface = surfaces.find((d) => {
+      const k = keyOf(d);
+      return k != null && Math.abs(k - key) < tolerance;
+    });
+    if (surface == null) {
+      surface = {
+        id: "",
+        age,
+        position,
+        status: "derived",
+        type: "",
+        calibration: null,
+        proportion: null,
+        unitsAbove: [],
+        unitsBelow: [],
+        section_id: unit.section_id
+      };
+      surfaces.push(surface);
+    }
+    if (top) {
+      surface.unitsBelow.push(unit.unit_id);
+    } else {
+      surface.unitsAbove.push(unit.unit_id);
+    }
+  }
+  for (const unit of units) {
+    addBound(unit, true);
+    addBound(unit, false);
+  }
+  surfaces.sort((a, b) => compareAlongAxis(a, b, axisType));
+  return surfaces.map((s, i) => ({ ...s, id: `derived-${i}` }));
+}
+function surfacePosition(surface, axisType) {
+  if (isPositionAxis(axisType)) {
+    return surface.position ?? null;
+  }
+  return surface.age ?? null;
+}
+function isPositionAxis(axisType) {
+  return axisType === ColumnAxisType.HEIGHT || axisType === ColumnAxisType.DEPTH;
+}
+function compareAlongAxis(a, b, axisType) {
+  const pa = surfacePosition(a, axisType);
+  const pb = surfacePosition(b, axisType);
+  if (pa == null || pb == null) return 0;
+  if (axisType === ColumnAxisType.HEIGHT) {
+    return pb - pa;
+  }
+  return pa - pb;
+}
+function formatProportion(proportion) {
+  if (proportion == null) return null;
+  return `${Math.round(proportion * 100)}%`;
+}
+function surfaceLabel(surface) {
+  const { calibration } = surface;
+  if (calibration != null) {
+    const prop = formatProportion(surface.proportion);
+    if (prop != null) {
+      return `${calibration.name} · ${prop}`;
+    }
+    return calibration.name;
+  }
+  return formatAge(surface.age);
+}
+function formatAge(age) {
+  if (age == null) return "";
+  let value = age;
+  let unit = "Ma";
+  if (Math.abs(age) < 1) {
+    value = age * 1e3;
+    unit = "ka";
+  }
+  const _value = value.toLocaleString("en-US", { maximumFractionDigits: 2 });
+  return `${_value} ${unit}`;
+}
+function asNumber(value) {
+  if (value == null || value === "") return null;
+  const n = Number(value);
+  if (Number.isNaN(n)) return null;
+  return n;
+}
+export {
+  INTERVAL_BOUND_TOLERANCE,
+  TIE_POINT_STATUSES,
+  boundaryToSurface,
+  compareAlongAxis,
+  formatAge,
+  formatProportion,
+  inferTiePointStatuses,
+  isPositionAxis,
+  nullifyUnitID,
+  surfaceLabel,
+  surfacePosition,
+  surfaceStatusDescriptions,
+  surfaceStatusLabels,
+  surfaceStatuses,
+  surfaceToken,
+  surfacesFromBoundaries,
+  surfacesFromUnits
+};
+//# sourceMappingURL=types.js.map

package/dist/age-model/types.js.map

@@ -0,0 +1 @@
+{"version":3,"file":"types.js","sources":["../../src/age-model/types.ts"],"sourcesContent":["/** Surfaces: the view model for a column's age model.\n *\n * A `ColumnSurface` is one horizontal boundary in a column. It is what the\n * surfaces view renders, and it is deliberately data-agnostic: a surface can\n * come from an `/age_model` boundary (`boundaryToSurface`), be derived from\n * unit tops and bottoms when a column has no age-model rows\n * (`surfacesFromUnits`), or be built by an editor from its own state.\n */\nimport type {\n  AgeModelBoundary,\n  BoundaryStatus,\n  BoundaryType,\n  UnitLong,\n} from \"@macrostrat/api-types\";\nimport { ColumnAxisType } from \"@macrostrat/column-components\";\n\nexport type { AgeModelBoundary, BoundaryStatus, BoundaryType };\n\n/** A boundary status, plus `derived` for surfaces computed from units rather\n * than stored in the age model. */\nexport type SurfaceStatus = BoundaryStatus | \"derived\";\n\n/** The interval a surface is calibrated against */\nexport interface SurfaceCalibration {\n  id: number;\n  name: string;\n  b_age: number;\n  t_age: number;\n}\n\nexport interface ColumnSurface {\n  id: number | string;\n  /** Modeled age of the surface (Ma) */\n  age: number;\n  /** Measured position (height or depth) along the column, when known */\n  position?: number | null;\n  status: SurfaceStatus;\n  type: BoundaryType;\n  /** The calibration interval, and the surface's position within it\n   * (0 at the interval's base, 1 at its top) */\n  calibration?: SurfaceCalibration | null;\n  proportion?: number | null;\n  /** True when `status` was inferred rather than read from the data — see\n   * `inferTiePointStatuses`. */\n  statusInferred?: boolean;\n  unitsAbove: number[];\n  unitsBelow: number[];\n  section_id?: number | null;\n  ref_id?: number | null;\n  /** The API record, when the surface came from `/age_model` */\n  boundary?: AgeModelBoundary;\n}\n\nexport const surfaceStatuses: SurfaceStatus[] = [\n  \"absolute\",\n  \"relative\",\n  \"spike\",\n  \"imposed\",\n  \"modeled\",\n  \"\",\n  \"derived\",\n];\n\n/** The statuses at which the age model was actually constrained — the tie\n * points — as opposed to surfaces whose age falls out of the model\n * (`modeled`), carries no status, or is merely derived from units. These are\n * the surfaces worth labeling by default. */\nexport const TIE_POINT_STATUSES: SurfaceStatus[] = [\n  \"absolute\",\n  \"relative\",\n  \"spike\",\n  \"imposed\",\n];\n\nexport const surfaceStatusLabels: Record<SurfaceStatus, string> = {\n  absolute: \"Absolute\",\n  relative: \"Relative\",\n  spike: \"Spike\",\n  imposed: \"Imposed\",\n  modeled: \"Modeled\",\n  \"\": \"Unspecified\",\n  derived: \"Derived\",\n};\n\nexport const surfaceStatusDescriptions: Record<SurfaceStatus, string> = {\n  absolute:\n    \"Independently dated; the age is authoritative and the interval position is derived from it\",\n  relative:\n    \"Calibrated to a position within a geologic interval; the age follows the timescale\",\n  spike: \"Tied to a formally defined boundary point\",\n  imposed:\n    \"Interval-relative representation imposed on a boundary with no chronostratigraphic tie of its own (for instance, thickness interpolation)\",\n  modeled: \"Interpolated by the age model between calibrated surfaces\",\n  \"\": \"No status recorded\",\n  derived:\n    \"Computed from unit tops and bottoms; not stored in the column's age model\",\n};\n\n/** A CSS-safe token for a status or type (`angular unconformity` →\n * `angular-unconformity`, `\"\"` → `unspecified`). */\nexport function surfaceToken(value: string | null | undefined): string {\n  if (value == null || value === \"\") return \"unspecified\";\n  return value.replace(/[^a-z0-9]+/gi, \"-\").toLowerCase();\n}\n\n/** Legacy `unit_boundaries` rows use `0` as the \"no unit\" sentinel. */\nexport function nullifyUnitID(id: number | null | undefined): number | null {\n  if (id == null || id === 0) return null;\n  return id;\n}\n\n/** How close to an interval's base or top a surface must sit to count as\n * bounding it (as a fraction of the interval). */\nexport const INTERVAL_BOUND_TOLERANCE = 0.001;\n\n/** Promote the `modeled` surfaces that cannot in fact have been interpolated\n * to `relative`.\n *\n * This is a client-side affordance over a known gap in the data. An age model\n * interpolates *between* tie points, so a surface with no constraint to\n * interpolate from has to be one itself, whatever `boundary_status` says. Two\n * cases are recognizable from the boundaries alone:\n *\n * - **It sits on the base or top of its calibration interval** (proportion 0\n *   or 1). A height referred directly to an interval bounds the model rather\n *   than falling out of it.\n * - **It is the edge of a gap-bound package** — the youngest surface in its\n *   section with no unit above, or the oldest with none below. There is\n *   nothing beyond it to interpolate against.\n *\n * Both are common in columns whose age model was built by assigning positions\n * within intervals, where the importer recorded everything as `modeled`;\n * filtering to the tie-point statuses would otherwise hide real ones. Promoted\n * surfaces are marked `statusInferred`, so views can say where the status came\n * from; nothing about the underlying record changes.\n *\n * Two cases are deliberately left alone. A *mid-interval* position is\n * ambiguous: an assigned position and an interpolated one are the same record\n * there, and only the data can settle it. And an open side in the *interior*\n * of a section is a lateral relationship, not a package edge — a unit that\n * pinches out while others span past it — which the model can interpolate\n * normally.\n */\nexport function inferTiePointStatuses(\n  surfaces: ColumnSurface[],\n  tolerance: number = INTERVAL_BOUND_TOLERANCE,\n): ColumnSurface[] {\n  const extents = sectionAgeExtents(surfaces);\n\n  return surfaces.map((surface) => {\n    if (surface.status !== \"modeled\") return surface;\n    if (surface.calibration == null) return surface;\n\n    const { proportion } = surface;\n    const onIntervalBound =\n      proportion != null &&\n      (Math.abs(proportion) <= tolerance ||\n        Math.abs(proportion - 1) <= tolerance);\n\n    if (!onIntervalBound && !boundsPackage(surface, extents)) {\n      return surface;\n    }\n    return { ...surface, status: \"relative\", statusInferred: true };\n  });\n}\n\n/** An open side only means a package edge at the ends of the section. In its\n * interior it means a unit pinching out while others span past it, which the\n * age model interpolates like any other surface. */\nfunction boundsPackage(\n  surface: ColumnSurface,\n  extents: Map<ColumnSurface[\"section_id\"], [number, number]>,\n): boolean {\n  const extent = extents.get(surface.section_id ?? null);\n  if (extent == null) return false;\n  const [youngest, oldest] = extent;\n  if (surface.unitsAbove.length === 0 && surface.age === youngest) return true;\n  return surface.unitsBelow.length === 0 && surface.age === oldest;\n}\n\n/** The youngest and oldest surface age in each section. */\nfunction sectionAgeExtents(\n  surfaces: ColumnSurface[],\n): Map<ColumnSurface[\"section_id\"], [number, number]> {\n  const extents = new Map<ColumnSurface[\"section_id\"], [number, number]>();\n  for (const surface of surfaces) {\n    const key = surface.section_id ?? null;\n    const extent = extents.get(key);\n    if (extent == null) {\n      extents.set(key, [surface.age, surface.age]);\n      continue;\n    }\n    extent[0] = Math.min(extent[0], surface.age);\n    extent[1] = Math.max(extent[1], surface.age);\n  }\n  return extents;\n}\n\nexport function boundaryToSurface(boundary: AgeModelBoundary): ColumnSurface {\n  const above = nullifyUnitID(boundary.unit_above);\n  const below = nullifyUnitID(boundary.unit_below);\n  let calibration: SurfaceCalibration | null = null;\n  if (boundary.interval_id != null) {\n    calibration = {\n      id: boundary.interval_id,\n      name: boundary.interval_name,\n      b_age: boundary.age_bottom,\n      t_age: boundary.age_top,\n    };\n  }\n  return {\n    id: boundary.boundary_id,\n    age: boundary.model_age,\n    position: boundary.boundary_position,\n    status: boundary.boundary_status ?? \"\",\n    type: boundary.boundary_type ?? \"\",\n    calibration,\n    proportion: boundary.rel_position,\n    unitsAbove: above == null ? [] : [above],\n    unitsBelow: below == null ? [] : [below],\n    section_id: boundary.section_id,\n    ref_id: boundary.ref_id,\n    boundary,\n  };\n}\n\nexport interface SurfacesFromBoundariesOptions {\n  /** Promote `modeled` surfaces that must in fact constrain the age model to\n   * `relative` (default `true`) — see `inferTiePointStatuses`. */\n  inferTiePoints?: boolean;\n}\n\n/** Surfaces for a column's age-model boundaries, youngest first. */\nexport function surfacesFromBoundaries(\n  boundaries: AgeModelBoundary[] | null | undefined,\n  options: SurfacesFromBoundariesOptions = {},\n): ColumnSurface[] {\n  const { inferTiePoints = true } = options;\n  if (boundaries == null) return [];\n  let surfaces = boundaries.map(boundaryToSurface);\n  if (inferTiePoints) {\n    surfaces = inferTiePointStatuses(surfaces);\n  }\n  return surfaces.sort((a, b) => a.age - b.age);\n}\n\nexport interface SurfacesFromUnitsOptions {\n  /** Surfaces closer than this (in the axis's units) are merged */\n  tolerance?: number;\n  /** Which unit field the surfaces are keyed on: ages for age columns,\n   * measured positions for height and depth columns */\n  axisType?: ColumnAxisType;\n}\n\n/** Derive surfaces from unit tops and bottoms. Every distinct top or bottom\n * becomes one surface with status `derived`; units sharing a boundary share\n * the surface. This is the fallback when a column has no `unit_boundaries`. */\nexport function surfacesFromUnits<T extends UnitLong>(\n  units: T[] | null | undefined,\n  options: SurfacesFromUnitsOptions = {},\n): ColumnSurface[] {\n  const { tolerance = 0.001, axisType = ColumnAxisType.AGE } = options;\n  if (units == null) return [];\n  const byPosition = isPositionAxis(axisType);\n\n  const surfaces: ColumnSurface[] = [];\n\n  function keyOf(surface: ColumnSurface): number | null {\n    if (byPosition) return surface.position ?? null;\n    return surface.age;\n  }\n\n  function addBound(unit: T, top: boolean) {\n    const age = top ? unit.t_age : unit.b_age;\n    const position = asNumber(top ? unit.t_pos : unit.b_pos);\n    const key = byPosition ? position : age;\n    if (key == null) return;\n\n    let surface = surfaces.find((d) => {\n      const k = keyOf(d);\n      return k != null && Math.abs(k - key) < tolerance;\n    });\n    if (surface == null) {\n      surface = {\n        id: \"\",\n        age,\n        position,\n        status: \"derived\",\n        type: \"\",\n        calibration: null,\n        proportion: null,\n        unitsAbove: [],\n        unitsBelow: [],\n        section_id: unit.section_id,\n      };\n      surfaces.push(surface);\n    }\n    // The unit whose top this is lies *below* the surface, and vice versa\n    if (top) {\n      surface.unitsBelow.push(unit.unit_id);\n    } else {\n      surface.unitsAbove.push(unit.unit_id);\n    }\n  }\n\n  for (const unit of units) {\n    addBound(unit, true);\n    addBound(unit, false);\n  }\n\n  surfaces.sort((a, b) => compareAlongAxis(a, b, axisType));\n  return surfaces.map((s, i) => ({ ...s, id: `derived-${i}` }));\n}\n\n/** The coordinate a surface is drawn at for a given axis type: its age for\n * age columns, its measured position for height and depth columns. `null`\n * when the surface has no coordinate on that axis. */\nexport function surfacePosition(\n  surface: ColumnSurface,\n  axisType: ColumnAxisType | null | undefined,\n): number | null {\n  if (isPositionAxis(axisType)) {\n    return surface.position ?? null;\n  }\n  return surface.age ?? null;\n}\n\nexport function isPositionAxis(\n  axisType: ColumnAxisType | null | undefined,\n): boolean {\n  return (\n    axisType === ColumnAxisType.HEIGHT || axisType === ColumnAxisType.DEPTH\n  );\n}\n\n/** Order surfaces from the top of the column downwards. */\nexport function compareAlongAxis(\n  a: ColumnSurface,\n  b: ColumnSurface,\n  axisType: ColumnAxisType | null | undefined,\n): number {\n  const pa = surfacePosition(a, axisType);\n  const pb = surfacePosition(b, axisType);\n  if (pa == null || pb == null) return 0;\n  if (axisType === ColumnAxisType.HEIGHT) {\n    // Heights increase upwards\n    return pb - pa;\n  }\n  // Ages and depths increase downwards\n  return pa - pb;\n}\n\n/** Percent of the way up its calibration interval, e.g. `85%`. */\nexport function formatProportion(\n  proportion: number | null | undefined,\n): string | null {\n  if (proportion == null) return null;\n  return `${Math.round(proportion * 100)}%`;\n}\n\n/** A short text label for a surface: its calibration (\"Calymmian · 85%\")\n * or, failing that, its age. */\nexport function surfaceLabel(surface: ColumnSurface): string {\n  const { calibration } = surface;\n  if (calibration != null) {\n    const prop = formatProportion(surface.proportion);\n    if (prop != null) {\n      return `${calibration.name} · ${prop}`;\n    }\n    return calibration.name;\n  }\n  return formatAge(surface.age);\n}\n\nexport function formatAge(age: number | null | undefined): string {\n  if (age == null) return \"\";\n  let value = age;\n  let unit = \"Ma\";\n  if (Math.abs(age) < 1) {\n    value = age * 1000;\n    unit = \"ka\";\n  }\n  const _value = value.toLocaleString(\"en-US\", { maximumFractionDigits: 2 });\n  return `${_value} ${unit}`;\n}\n\nfunction asNumber(value: unknown): number | null {\n  if (value == null || value === \"\") return null;\n  const n = Number(value);\n  if (Number.isNaN(n)) return null;\n  return n;\n}\n"],"names":[],"mappings":";AAqDO,MAAM,kBAAmC;AAAA,EAC9C;AAAA,EACA;AAAA,EACA;AAAA,EACA;AAAA,EACA;AAAA,EACA;AAAA,EACA;AACF;AAMO,MAAM,qBAAsC;AAAA,EACjD;AAAA,EACA;AAAA,EACA;AAAA,EACA;AACF;AAEO,MAAM,sBAAqD;AAAA,EAChE,UAAU;AAAA,EACV,UAAU;AAAA,EACV,OAAO;AAAA,EACP,SAAS;AAAA,EACT,SAAS;AAAA,EACT,IAAI;AAAA,EACJ,SAAS;AACX;AAEO,MAAM,4BAA2D;AAAA,EACtE,UACE;AAAA,EACF,UACE;AAAA,EACF,OAAO;AAAA,EACP,SACE;AAAA,EACF,SAAS;AAAA,EACT,IAAI;AAAA,EACJ,SACE;AACJ;AAIO,SAAS,aAAa,OAA0C;AACrE,MAAI,SAAS,QAAQ,UAAU,GAAI,QAAO;AAC1C,SAAO,MAAM,QAAQ,gBAAgB,GAAG,EAAE,YAAA;AAC5C;AAGO,SAAS,cAAc,IAA8C;AAC1E,MAAI,MAAM,QAAQ,OAAO,EAAG,QAAO;AACnC,SAAO;AACT;AAIO,MAAM,2BAA2B;AA8BjC,SAAS,sBACd,UACA,YAAoB,0BACH;AACjB,QAAM,UAAU,kBAAkB,QAAQ;AAE1C,SAAO,SAAS,IAAI,CAAC,YAAY;AAC/B,QAAI,QAAQ,WAAW,UAAW,QAAO;AACzC,QAAI,QAAQ,eAAe,KAAM,QAAO;AAExC,UAAM,EAAE,eAAe;AACvB,UAAM,kBACJ,cAAc,SACb,KAAK,IAAI,UAAU,KAAK,aACvB,KAAK,IAAI,aAAa,CAAC,KAAK;AAEhC,QAAI,CAAC,mBAAmB,CAAC,cAAc,SAAS,OAAO,GAAG;AACxD,aAAO;AAAA,IACT;AACA,WAAO,EAAE,GAAG,SAAS,QAAQ,YAAY,gBAAgB,KAAA;AAAA,EAC3D,CAAC;AACH;AAKA,SAAS,cACP,SACA,SACS;AACT,QAAM,SAAS,QAAQ,IAAI,QAAQ,cAAc,IAAI;AACrD,MAAI,UAAU,KAAM,QAAO;AAC3B,QAAM,CAAC,UAAU,MAAM,IAAI;AAC3B,MAAI,QAAQ,WAAW,WAAW,KAAK,QAAQ,QAAQ,SAAU,QAAO;AACxE,SAAO,QAAQ,WAAW,WAAW,KAAK,QAAQ,QAAQ;AAC5D;AAGA,SAAS,kBACP,UACoD;AACpD,QAAM,8BAAc,IAAA;AACpB,aAAW,WAAW,UAAU;AAC9B,UAAM,MAAM,QAAQ,cAAc;AAClC,UAAM,SAAS,QAAQ,IAAI,GAAG;AAC9B,QAAI,UAAU,MAAM;AAClB,cAAQ,IAAI,KAAK,CAAC,QAAQ,KAAK,QAAQ,GAAG,CAAC;AAC3C;AAAA,IACF;AACA,WAAO,CAAC,IAAI,KAAK,IAAI,OAAO,CAAC,GAAG,QAAQ,GAAG;AAC3C,WAAO,CAAC,IAAI,KAAK,IAAI,OAAO,CAAC,GAAG,QAAQ,GAAG;AAAA,EAC7C;AACA,SAAO;AACT;AAEO,SAAS,kBAAkB,UAA2C;AAC3E,QAAM,QAAQ,cAAc,SAAS,UAAU;AAC/C,QAAM,QAAQ,cAAc,SAAS,UAAU;AAC/C,MAAI,cAAyC;AAC7C,MAAI,SAAS,eAAe,MAAM;AAChC,kBAAc;AAAA,MACZ,IAAI,SAAS;AAAA,MACb,MAAM,SAAS;AAAA,MACf,OAAO,SAAS;AAAA,MAChB,OAAO,SAAS;AAAA,IAAA;AAAA,EAEpB;AACA,SAAO;AAAA,IACL,IAAI,SAAS;AAAA,IACb,KAAK,SAAS;AAAA,IACd,UAAU,SAAS;AAAA,IACnB,QAAQ,SAAS,mBAAmB;AAAA,IACpC,MAAM,SAAS,iBAAiB;AAAA,IAChC;AAAA,IACA,YAAY,SAAS;AAAA,IACrB,YAAY,SAAS,OAAO,CAAA,IAAK,CAAC,KAAK;AAAA,IACvC,YAAY,SAAS,OAAO,CAAA,IAAK,CAAC,KAAK;AAAA,IACvC,YAAY,SAAS;AAAA,IACrB,QAAQ,SAAS;AAAA,IACjB;AAAA,EAAA;AAEJ;AASO,SAAS,uBACd,YACA,UAAyC,IACxB;AACjB,QAAM,EAAE,iBAAiB,KAAA,IAAS;AAClC,MAAI,cAAc,KAAM,QAAO,CAAA;AAC/B,MAAI,WAAW,WAAW,IAAI,iBAAiB;AAC/C,MAAI,gBAAgB;AAClB,eAAW,sBAAsB,QAAQ;AAAA,EAC3C;AACA,SAAO,SAAS,KAAK,CAAC,GAAG,MAAM,EAAE,MAAM,EAAE,GAAG;AAC9C;AAaO,SAAS,kBACd,OACA,UAAoC,IACnB;AACjB,QAAM,EAAE,YAAY,MAAO,WAAW,eAAe,QAAQ;AAC7D,MAAI,SAAS,KAAM,QAAO,CAAA;AAC1B,QAAM,aAAa,eAAe,QAAQ;AAE1C,QAAM,WAA4B,CAAA;AAElC,WAAS,MAAM,SAAuC;AACpD,QAAI,WAAY,QAAO,QAAQ,YAAY;AAC3C,WAAO,QAAQ;AAAA,EACjB;AAEA,WAAS,SAAS,MAAS,KAAc;AACvC,UAAM,MAAM,MAAM,KAAK,QAAQ,KAAK;AACpC,UAAM,WAAW,SAAS,MAAM,KAAK,QAAQ,KAAK,KAAK;AACvD,UAAM,MAAM,aAAa,WAAW;AACpC,QAAI,OAAO,KAAM;AAEjB,QAAI,UAAU,SAAS,KAAK,CAAC,MAAM;AACjC,YAAM,IAAI,MAAM,CAAC;AACjB,aAAO,KAAK,QAAQ,KAAK,IAAI,IAAI,GAAG,IAAI;AAAA,IAC1C,CAAC;AACD,QAAI,WAAW,MAAM;AACnB,gBAAU;AAAA,QACR,IAAI;AAAA,QACJ;AAAA,QACA;AAAA,QACA,QAAQ;AAAA,QACR,MAAM;AAAA,QACN,aAAa;AAAA,QACb,YAAY;AAAA,QACZ,YAAY,CAAA;AAAA,QACZ,YAAY,CAAA;AAAA,QACZ,YAAY,KAAK;AAAA,MAAA;AAEnB,eAAS,KAAK,OAAO;AAAA,IACvB;AAEA,QAAI,KAAK;AACP,cAAQ,WAAW,KAAK,KAAK,OAAO;AAAA,IACtC,OAAO;AACL,cAAQ,WAAW,KAAK,KAAK,OAAO;AAAA,IACtC;AAAA,EACF;AAEA,aAAW,QAAQ,OAAO;AACxB,aAAS,MAAM,IAAI;AACnB,aAAS,MAAM,KAAK;AAAA,EACtB;AAEA,WAAS,KAAK,CAAC,GAAG,MAAM,iBAAiB,GAAG,GAAG,QAAQ,CAAC;AACxD,SAAO,SAAS,IAAI,CAAC,GAAG,OAAO,EAAE,GAAG,GAAG,IAAI,WAAW,CAAC,GAAA,EAAK;AAC9D;AAKO,SAAS,gBACd,SACA,UACe;AACf,MAAI,eAAe,QAAQ,GAAG;AAC5B,WAAO,QAAQ,YAAY;AAAA,EAC7B;AACA,SAAO,QAAQ,OAAO;AACxB;AAEO,SAAS,eACd,UACS;AACT,SACE,aAAa,eAAe,UAAU,aAAa,eAAe;AAEtE;AAGO,SAAS,iBACd,GACA,GACA,UACQ;AACR,QAAM,KAAK,gBAAgB,GAAG,QAAQ;AACtC,QAAM,KAAK,gBAAgB,GAAG,QAAQ;AACtC,MAAI,MAAM,QAAQ,MAAM,KAAM,QAAO;AACrC,MAAI,aAAa,eAAe,QAAQ;AAEtC,WAAO,KAAK;AAAA,EACd;AAEA,SAAO,KAAK;AACd;AAGO,SAAS,iBACd,YACe;AACf,MAAI,cAAc,KAAM,QAAO;AAC/B,SAAO,GAAG,KAAK,MAAM,aAAa,GAAG,CAAC;AACxC;AAIO,SAAS,aAAa,SAAgC;AAC3D,QAAM,EAAE,gBAAgB;AACxB,MAAI,eAAe,MAAM;AACvB,UAAM,OAAO,iBAAiB,QAAQ,UAAU;AAChD,QAAI,QAAQ,MAAM;AAChB,aAAO,GAAG,YAAY,IAAI,MAAM,IAAI;AAAA,IACtC;AACA,WAAO,YAAY;AAAA,EACrB;AACA,SAAO,UAAU,QAAQ,GAAG;AAC9B;AAEO,SAAS,UAAU,KAAwC;AAChE,MAAI,OAAO,KAAM,QAAO;AACxB,MAAI,QAAQ;AACZ,MAAI,OAAO;AACX,MAAI,KAAK,IAAI,GAAG,IAAI,GAAG;AACrB,YAAQ,MAAM;AACd,WAAO;AAAA,EACT;AACA,QAAM,SAAS,MAAM,eAAe,SAAS,EAAE,uBAAuB,GAAG;AACzE,SAAO,GAAG,MAAM,IAAI,IAAI;AAC1B;AAEA,SAAS,SAAS,OAA+B;AAC/C,MAAI,SAAS,QAAQ,UAAU,GAAI,QAAO;AAC1C,QAAM,IAAI,OAAO,KAAK;AACtB,MAAI,OAAO,MAAM,CAAC,EAAG,QAAO;AAC5B,SAAO;AACT;"}