@onerjs/core 8.55.9 → 8.56.2

package/Maths/colorTemperature.functions.js

@@ -0,0 +1,196 @@
+import { Clamp } from "./math.scalar.functions.js";
+import { Matrix, Vector2, Vector3 } from "./math.vector.pure.js";
+/**
+ * Tabulated Planckian locus in CIE 1960 UCS (u, v) space, with the local isotherm slope `t` used to offset
+ * a chromaticity along the tint axis. Source: Wyszecki & Stiles, "Color Science", 2nd ed., 1982, pp. 227-228.
+ * Each row is [mired, u, v, t].
+ */
+const PlanckianLocusTable = [
+    [0, 0.18006, 0.26352, -0.24341],
+    [10, 0.18066, 0.26589, -0.25479],
+    [20, 0.18133, 0.26846, -0.26876],
+    [30, 0.18208, 0.27119, -0.28539],
+    [40, 0.18293, 0.27407, -0.3047],
+    [50, 0.18388, 0.27709, -0.32675],
+    [60, 0.18494, 0.28021, -0.35156],
+    [70, 0.18611, 0.28342, -0.37915],
+    [80, 0.1874, 0.28668, -0.40955],
+    [90, 0.1888, 0.28997, -0.44278],
+    [100, 0.19032, 0.29326, -0.47888],
+    [125, 0.19462, 0.30141, -0.58204],
+    [150, 0.19962, 0.30921, -0.70471],
+    [175, 0.20525, 0.31647, -0.84901],
+    [200, 0.21142, 0.32312, -1.0182],
+    [225, 0.21807, 0.32909, -1.2168],
+    [250, 0.22511, 0.33439, -1.4512],
+    [275, 0.23247, 0.33904, -1.7298],
+    [300, 0.2401, 0.34308, -2.0637],
+    [325, 0.24792, 0.34655, -2.4681],
+    [350, 0.25591, 0.34951, -2.9641],
+    [375, 0.264, 0.352, -3.5814],
+    [400, 0.27218, 0.35407, -4.3633],
+    [425, 0.28039, 0.35577, -5.3762],
+    [450, 0.28863, 0.35714, -6.7262],
+    [475, 0.29685, 0.35823, -8.5955],
+    [500, 0.30505, 0.35907, -11.324],
+    [525, 0.3132, 0.35968, -15.628],
+    [550, 0.32129, 0.36011, -23.325],
+    [575, 0.32931, 0.36038, -40.77],
+    [600, 0.33724, 0.36051, -116.45],
+];
+/**
+ * The minimum correlated color temperature, in Kelvin, representable by {@link TemperatureTintToXyz} (there is no
+ * upper bound: increasingly high temperatures approach mired 0, already within the tabulated range).
+ */
+export const MinTemperatureKelvin = 1e6 / PlanckianLocusTable[PlanckianLocusTable.length - 1][0];
+/**
+ * The maximum magnitude of the tint offset accepted by {@link TemperatureTintToXyz}, in either direction.
+ */
+export const MaxTintMagnitude = 150;
+/**
+ * Converts a correlated color temperature and tint offset into the CIE XYZ (Y = 1) coordinates of the
+ * corresponding illuminant white point, using a tabulated approximation of the Planckian locus in CIE 1960 UCS
+ * (u, v) space.
+ * @param temperatureKelvin The correlated color temperature of the illuminant, in Kelvin
+ * @param tint An offset perpendicular to the Planckian locus (the green/magenta axis), in the range
+ * [-{@link MaxTintMagnitude}, {@link MaxTintMagnitude}]
+ * @returns The CIE XYZ (Y = 1) coordinates of the illuminant white point
+ */
+export function TemperatureTintToXyz(temperatureKelvin, tint) {
+    // Clamp before taking the reciprocal: a non-positive or NaN Kelvin value has no valid reciprocal-mired
+    // representation, and taking 1e6 / temperatureKelvin first (before clamping the result) would map negative
+    // inputs to the *hottest* end of the table instead of the coldest/minimum one.
+    const safeTemperatureKelvin = Number.isNaN(temperatureKelvin) ? MinTemperatureKelvin : Math.max(temperatureKelvin, MinTemperatureKelvin);
+    const mired = Clamp(1e6 / safeTemperatureKelvin, PlanckianLocusTable[0][0], PlanckianLocusTable[PlanckianLocusTable.length - 1][0]);
+    // Binary search for the table row pair bracketing `mired` (the table is sorted ascending by mired).
+    let low = 0;
+    let high = PlanckianLocusTable.length - 1;
+    while (high - low > 1) {
+        const mid = (low + high) >> 1;
+        if (PlanckianLocusTable[mid][0] <= mired) {
+            low = mid;
+        }
+        else {
+            high = mid;
+        }
+    }
+    const index = low;
+    const [lowMired, lowU, lowV, lowT] = PlanckianLocusTable[index];
+    const [highMired, highU, highV, highT] = PlanckianLocusTable[index + 1];
+    const f = (mired - lowMired) / (highMired - lowMired);
+    const uv = Vector2.Lerp(new Vector2(lowU, lowV), new Vector2(highU, highV), f);
+    const isothermLow = Vector2.Normalize(new Vector2(1, lowT));
+    const isothermHigh = Vector2.Normalize(new Vector2(1, highT));
+    const isotherm = Vector2.Normalize(Vector2.Lerp(isothermLow, isothermHigh, f));
+    const clampedTint = Clamp(tint, -MaxTintMagnitude, MaxTintMagnitude);
+    uv.subtractInPlace(isotherm.scale(clampedTint / 3000));
+    // CIE 1960 (u, v) -> CIE xyY (Y = 1) -> CIE XYZ.
+    const denom = 2 * uv.x - 8 * uv.y + 4;
+    const x = (3 * uv.x) / denom;
+    const y = (2 * uv.y) / denom;
+    return new Vector3(x / y, 1, (1 - x - y) / y);
+}
+// `Matrix.FromArray` reads values column-major, so each 3x3 matrix below is written column-by-column, not row-by-row.
+/**
+ * The Bradford cone-response matrix used for chromatic adaptation.
+ */
+// prettier-ignore
+const BradfordMatrix = Matrix.FromArray([
+    0.8951, -0.7502, 0.0389, 0,
+    0.2664, 1.7135, -0.0685, 0,
+    -0.1614, 0.0367, 1.0296, 0,
+    0, 0, 0, 1,
+]);
+const InverseBradfordMatrix = Matrix.Invert(BradfordMatrix);
+/**
+ * Computes a per-channel Bradford adaptation ratio, clamped to avoid extreme temperature/tint combinations
+ * blowing up the resulting matrix when a cone-response component approaches zero.
+ * @param destComponent The destination white's cone-response component
+ * @param sourceComponent The source white's cone-response component
+ * @returns The ratio to use for this channel, clamped to [0.1, 10]
+ */
+function ClampedAdaptationRatio(destComponent, sourceComponent) {
+    const ratio = sourceComponent > 0 ? destComponent / sourceComponent : 10;
+    return Clamp(ratio, 0.1, 10);
+}
+/**
+ * Computes the Bradford chromatic-adaptation matrix that maps CIE XYZ values from `sourceWhiteXYZ` to `destWhiteXYZ`.
+ * @param sourceWhiteXYZ The CIE XYZ (Y = 1) coordinates of the source reference white
+ * @param destWhiteXYZ The CIE XYZ (Y = 1) coordinates of the destination reference white
+ * @returns A matrix mapping source-referenced XYZ to destination-referenced XYZ
+ */
+function BradfordAdapt(sourceWhiteXYZ, destWhiteXYZ) {
+    const sourceLMS = Vector3.TransformNormal(sourceWhiteXYZ, BradfordMatrix);
+    const destLMS = Vector3.TransformNormal(destWhiteXYZ, BradfordMatrix);
+    // prettier-ignore
+    const scale = Matrix.FromArray([
+        ClampedAdaptationRatio(destLMS.x, sourceLMS.x), 0, 0, 0,
+        0, ClampedAdaptationRatio(destLMS.y, sourceLMS.y), 0, 0,
+        0, 0, ClampedAdaptationRatio(destLMS.z, sourceLMS.z), 0,
+        0, 0, 0, 1,
+    ]);
+    // adapt = InverseBradfordMatrix * scale * BradfordMatrix, written back-to-front since `A.multiply(B)` computes `B x A`.
+    const scaledByBradford = BradfordMatrix.multiply(scale);
+    return scaledByBradford.multiply(InverseBradfordMatrix);
+}
+/**
+ * The standard CIE Illuminant D65 white point, in CIE XYZ (Y = 1).
+ */
+const StandardD65Xyz = new Vector3(0.95047, 1.0, 1.08883);
+/**
+ * The standard sRGB / Rec.709 linear RGB to CIE XYZ matrix, calibrated to the externally-standardized CIE
+ * Illuminant D65 chromaticity - i.e. `StandardRgbToXyzD65 * (1, 1, 1) = StandardD65Xyz` exactly.
+ */
+// prettier-ignore
+const StandardRgbToXyzD65 = Matrix.FromArray([
+    0.4124564, 0.2126729, 0.0193339, 0,
+    0.3575761, 0.7151522, 0.119192, 0,
+    0.1804375, 0.072175, 0.9503041, 0,
+    0, 0, 0, 1,
+]);
+/**
+ * The linear RGB to CIE XYZ matrix actually used by {@link GetWhiteBalanceMatrix}, recalibrated from
+ * {@link StandardRgbToXyzD65} via a single Bradford adaptation step so that neutral (1, 1, 1) maps exactly onto
+ * this module's own Planckian-locus approximation of 6500 K (see {@link TemperatureTintToXyz}) instead of the
+ * externally-standardized CIE D65 chromaticity. Those two are close but never exactly equal, by definition: D65
+ * is defined via the CIE daylight locus, a model fit to measured daylight spectra, which deliberately deviates
+ * from the Planckian locus (real daylight isn't quite blackbody radiation) - no table or precision fix can close
+ * that gap without distorting temperature/tint for every other illuminant.
+ *
+ * This is what lets `GetWhiteBalanceMatrix(6500, 0)` be exactly the identity matrix - the default "no adjustment"
+ * temperature/tint are a true no-op - without needing to special-case that one input: the reverse adaptation is
+ * applied once, here, to the matrix itself, rather than to every call's result.
+ */
+const RgbToXyzD65 = StandardRgbToXyzD65.multiply(BradfordAdapt(StandardD65Xyz, TemperatureTintToXyz(6500, 0)));
+/**
+ * The CIE XYZ to linear RGB matrix, the inverse of {@link RgbToXyzD65}.
+ */
+const XyzToRgbD65 = Matrix.Invert(RgbToXyzD65);
+/**
+ * The working color space's reference white, in CIE XYZ (Y = 1) - what neutral (1, 1, 1) linear RGB maps to
+ * through {@link RgbToXyzD65}. Since that matrix is itself calibrated to this module's own 6500 K point (see
+ * above), this is exactly `TemperatureTintToXyz(6500, 0)`. Deriving it as "whatever (1, 1, 1) maps to" (rather
+ * than calling `TemperatureTintToXyz` a second time) keeps the invariant self-documenting: a fully white-balanced
+ * neutral patch round-trips back to exactly (1, 1, 1) for any input illuminant, because the adaptation target is
+ * always tied directly to whatever the working color space itself already treats as neutral.
+ */
+const ReferenceWhiteXyz = Vector3.TransformNormal(new Vector3(1, 1, 1), RgbToXyzD65);
+/**
+ * Computes the linear RGB (sRGB / Rec.709 primaries) color-correction matrix that white-balances the given
+ * illuminant, by chromatically adapting its white point (see {@link TemperatureTintToXyz}) to the working
+ * color space's reference white using the Bradford transform.
+ * @param temperatureKelvin The correlated color temperature of the illuminant to neutralize, in Kelvin
+ * @param tint An offset perpendicular to the Planckian locus (the green/magenta axis), in the range [-150, 150]
+ * @returns A column-major 3x3 matrix (9 values), ready to be bound as a `mat3` shader uniform
+ */
+export function GetWhiteBalanceMatrix(temperatureKelvin, tint) {
+    const targetWhiteXYZ = TemperatureTintToXyz(temperatureKelvin, tint);
+    const adapt = BradfordAdapt(targetWhiteXYZ, ReferenceWhiteXyz);
+    // finalMatrix = XyzToRgbD65 * adapt * RgbToXyzD65, written back-to-front (see BradfordAdapt).
+    const adaptedFromRgb = RgbToXyzD65.multiply(adapt);
+    const finalMatrix = adaptedFromRgb.multiply(XyzToRgbD65);
+    // `Matrix` stores its values column-major, matching the layout GLSL/WGSL `mat3` uniforms expect,
+    // so the upper-left 3x3 block can be extracted directly with no transpose.
+    return Matrix.GetAsMatrix3x3(finalMatrix);
+}
+//# sourceMappingURL=colorTemperature.functions.js.map

package/Maths/colorTemperature.functions.js.map

@@ -0,0 +1 @@
+{"version":3,"file":"colorTemperature.functions.js","sourceRoot":"","sources":["../../../../dev/core/src/Maths/colorTemperature.functions.ts"],"names":[],"mappings":"AAAA,OAAO,EAAE,KAAK,EAAE,MAAM,yBAAyB,CAAC;AAChD,OAAO,EAAE,MAAM,EAAE,OAAO,EAAE,OAAO,EAAE,MAAM,oBAAoB,CAAC;AAE9D;;;;GAIG;AACH,MAAM,mBAAmB,GAA6D;IAClF,CAAC,CAAC,EAAE,OAAO,EAAE,OAAO,EAAE,CAAC,OAAO,CAAC;IAC/B,CAAC,EAAE,EAAE,OAAO,EAAE,OAAO,EAAE,CAAC,OAAO,CAAC;IAChC,CAAC,EAAE,EAAE,OAAO,EAAE,OAAO,EAAE,CAAC,OAAO,CAAC;IAChC,CAAC,EAAE,EAAE,OAAO,EAAE,OAAO,EAAE,CAAC,OAAO,CAAC;IAChC,CAAC,EAAE,EAAE,OAAO,EAAE,OAAO,EAAE,CAAC,MAAM,CAAC;IAC/B,CAAC,EAAE,EAAE,OAAO,EAAE,OAAO,EAAE,CAAC,OAAO,CAAC;IAChC,CAAC,EAAE,EAAE,OAAO,EAAE,OAAO,EAAE,CAAC,OAAO,CAAC;IAChC,CAAC,EAAE,EAAE,OAAO,EAAE,OAAO,EAAE,CAAC,OAAO,CAAC;IAChC,CAAC,EAAE,EAAE,MAAM,EAAE,OAAO,EAAE,CAAC,OAAO,CAAC;IAC/B,CAAC,EAAE,EAAE,MAAM,EAAE,OAAO,EAAE,CAAC,OAAO,CAAC;IAC/B,CAAC,GAAG,EAAE,OAAO,EAAE,OAAO,EAAE,CAAC,OAAO,CAAC;IACjC,CAAC,GAAG,EAAE,OAAO,EAAE,OAAO,EAAE,CAAC,OAAO,CAAC;IACjC,CAAC,GAAG,EAAE,OAAO,EAAE,OAAO,EAAE,CAAC,OAAO,CAAC;IACjC,CAAC,GAAG,EAAE,OAAO,EAAE,OAAO,EAAE,CAAC,OAAO,CAAC;IACjC,CAAC,GAAG,EAAE,OAAO,EAAE,OAAO,EAAE,CAAC,MAAM,CAAC;IAChC,CAAC,GAAG,EAAE,OAAO,EAAE,OAAO,EAAE,CAAC,MAAM,CAAC;IAChC,CAAC,GAAG,EAAE,OAAO,EAAE,OAAO,EAAE,CAAC,MAAM,CAAC;IAChC,CAAC,GAAG,EAAE,OAAO,EAAE,OAAO,EAAE,CAAC,MAAM,CAAC;IAChC,CAAC,GAAG,EAAE,MAAM,EAAE,OAAO,EAAE,CAAC,MAAM,CAAC;IAC/B,CAAC,GAAG,EAAE,OAAO,EAAE,OAAO,EAAE,CAAC,MAAM,CAAC;IAChC,CAAC,GAAG,EAAE,OAAO,EAAE,OAAO,EAAE,CAAC,MAAM,CAAC;IAChC,CAAC,GAAG,EAAE,KAAK,EAAE,KAAK,EAAE,CAAC,MAAM,CAAC;IAC5B,CAAC,GAAG,EAAE,OAAO,EAAE,OAAO,EAAE,CAAC,MAAM,CAAC;IAChC,CAAC,GAAG,EAAE,OAAO,EAAE,OAAO,EAAE,CAAC,MAAM,CAAC;IAChC,CAAC,GAAG,EAAE,OAAO,EAAE,OAAO,EAAE,CAAC,MAAM,CAAC;IAChC,CAAC,GAAG,EAAE,OAAO,EAAE,OAAO,EAAE,CAAC,MAAM,CAAC;IAChC,CAAC,GAAG,EAAE,OAAO,EAAE,OAAO,EAAE,CAAC,MAAM,CAAC;IAChC,CAAC,GAAG,EAAE,MAAM,EAAE,OAAO,EAAE,CAAC,MAAM,CAAC;IAC/B,CAAC,GAAG,EAAE,OAAO,EAAE,OAAO,EAAE,CAAC,MAAM,CAAC;IAChC,CAAC,GAAG,EAAE,OAAO,EAAE,OAAO,EAAE,CAAC,KAAK,CAAC;IAC/B,CAAC,GAAG,EAAE,OAAO,EAAE,OAAO,EAAE,CAAC,MAAM,CAAC;CACnC,CAAC;AAEF;;;GAGG;AACH,MAAM,CAAC,MAAM,oBAAoB,GAAG,GAAG,GAAG,mBAAmB,CAAC,mBAAmB,CAAC,MAAM,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;AAEjG;;GAEG;AACH,MAAM,CAAC,MAAM,gBAAgB,GAAG,GAAG,CAAC;AAEpC;;;;;;;;GAQG;AACH,MAAM,UAAU,oBAAoB,CAAC,iBAAyB,EAAE,IAAY;IACxE,uGAAuG;IACvG,2GAA2G;IAC3G,+EAA+E;IAC/E,MAAM,qBAAqB,GAAG,MAAM,CAAC,KAAK,CAAC,iBAAiB,CAAC,CAAC,CAAC,CAAC,oBAAoB,CAAC,CAAC,CAAC,IAAI,CAAC,GAAG,CAAC,iBAAiB,EAAE,oBAAoB,CAAC,CAAC;IACzI,MAAM,KAAK,GAAG,KAAK,CAAC,GAAG,GAAG,qBAAqB,EAAE,mBAAmB,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,EAAE,mBAAmB,CAAC,mBAAmB,CAAC,MAAM,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;IAEpI,oGAAoG;IACpG,IAAI,GAAG,GAAG,CAAC,CAAC;IACZ,IAAI,IAAI,GAAG,mBAAmB,CAAC,MAAM,GAAG,CAAC,CAAC;IAC1C,OAAO,IAAI,GAAG,GAAG,GAAG,CAAC,EAAE,CAAC;QACpB,MAAM,GAAG,GAAG,CAAC,GAAG,GAAG,IAAI,CAAC,IAAI,CAAC,CAAC;QAC9B,IAAI,mBAAmB,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,IAAI,KAAK,EAAE,CAAC;YACvC,GAAG,GAAG,GAAG,CAAC;QACd,CAAC;aAAM,CAAC;YACJ,IAAI,GAAG,GAAG,CAAC;QACf,CAAC;IACL,CAAC;IACD,MAAM,KAAK,GAAG,GAAG,CAAC;IAElB,MAAM,CAAC,QAAQ,EAAE,IAAI,EAAE,IAAI,EAAE,IAAI,CAAC,GAAG,mBAAmB,CAAC,KAAK,CAAC,CAAC;IAChE,MAAM,CAAC,SAAS,EAAE,KAAK,EAAE,KAAK,EAAE,KAAK,CAAC,GAAG,mBAAmB,CAAC,KAAK,GAAG,CAAC,CAAC,CAAC;IACxE,MAAM,CAAC,GAAG,CAAC,KAAK,GAAG,QAAQ,CAAC,GAAG,CAAC,SAAS,GAAG,QAAQ,CAAC,CAAC;IAEtD,MAAM,EAAE,GAAG,OAAO,CAAC,IAAI,CAAC,IAAI,OAAO,CAAC,IAAI,EAAE,IAAI,CAAC,EAAE,IAAI,OAAO,CAAC,KAAK,EAAE,KAAK,CAAC,EAAE,CAAC,CAAC,CAAC;IAE/E,MAAM,WAAW,GAAG,OAAO,CAAC,SAAS,CAAC,IAAI,OAAO,CAAC,CAAC,EAAE,IAAI,CAAC,CAAC,CAAC;IAC5D,MAAM,YAAY,GAAG,OAAO,CAAC,SAAS,CAAC,IAAI,OAAO,CAAC,CAAC,EAAE,KAAK,CAAC,CAAC,CAAC;IAC9D,MAAM,QAAQ,GAAG,OAAO,CAAC,SAAS,CAAC,OAAO,CAAC,IAAI,CAAC,WAAW,EAAE,YAAY,EAAE,CAAC,CAAC,CAAC,CAAC;IAE/E,MAAM,WAAW,GAAG,KAAK,CAAC,IAAI,EAAE,CAAC,gBAAgB,EAAE,gBAAgB,CAAC,CAAC;IACrE,EAAE,CAAC,eAAe,CAAC,QAAQ,CAAC,KAAK,CAAC,WAAW,GAAG,IAAI,CAAC,CAAC,CAAC;IAEvD,iDAAiD;IACjD,MAAM,KAAK,GAAG,CAAC,GAAG,EAAE,CAAC,CAAC,GAAG,CAAC,GAAG,EAAE,CAAC,CAAC,GAAG,CAAC,CAAC;IACtC,MAAM,CAAC,GAAG,CAAC,CAAC,GAAG,EAAE,CAAC,CAAC,CAAC,GAAG,KAAK,CAAC;IAC7B,MAAM,CAAC,GAAG,CAAC,CAAC,GAAG,EAAE,CAAC,CAAC,CAAC,GAAG,KAAK,CAAC;IAE7B,OAAO,IAAI,OAAO,CAAC,CAAC,GAAG,CAAC,EAAE,CAAC,EAAE,CAAC,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC;AAClD,CAAC;AAED,sHAAsH;AAEtH;;GAEG;AACH,kBAAkB;AAClB,MAAM,cAAc,GAAG,MAAM,CAAC,SAAS,CAAC;IACpC,MAAM,EAAE,CAAC,MAAM,EAAE,MAAM,EAAE,CAAC;IAC1B,MAAM,EAAE,MAAM,EAAE,CAAC,MAAM,EAAE,CAAC;IAC1B,CAAC,MAAM,EAAE,MAAM,EAAE,MAAM,EAAE,CAAC;IAC1B,CAAC,EAAE,CAAC,EAAE,CAAC,EAAE,CAAC;CACb,CAAC,CAAC;AAEH,MAAM,qBAAqB,GAAG,MAAM,CAAC,MAAM,CAAC,cAAc,CAAC,CAAC;AAE5D;;;;;;GAMG;AACH,SAAS,sBAAsB,CAAC,aAAqB,EAAE,eAAuB;IAC1E,MAAM,KAAK,GAAG,eAAe,GAAG,CAAC,CAAC,CAAC,CAAC,aAAa,GAAG,eAAe,CAAC,CAAC,CAAC,EAAE,CAAC;IACzE,OAAO,KAAK,CAAC,KAAK,EAAE,GAAG,EAAE,EAAE,CAAC,CAAC;AACjC,CAAC;AAED;;;;;GAKG;AACH,SAAS,aAAa,CAAC,cAAuB,EAAE,YAAqB;IACjE,MAAM,SAAS,GAAG,OAAO,CAAC,eAAe,CAAC,cAAc,EAAE,cAAc,CAAC,CAAC;IAC1E,MAAM,OAAO,GAAG,OAAO,CAAC,eAAe,CAAC,YAAY,EAAE,cAAc,CAAC,CAAC;IAEtE,kBAAkB;IAClB,MAAM,KAAK,GAAG,MAAM,CAAC,SAAS,CAAC;QAC3B,sBAAsB,CAAC,OAAO,CAAC,CAAC,EAAE,SAAS,CAAC,CAAC,CAAC,EAAE,CAAC,EAAE,CAAC,EAAE,CAAC;QACvD,CAAC,EAAE,sBAAsB,CAAC,OAAO,CAAC,CAAC,EAAE,SAAS,CAAC,CAAC,CAAC,EAAE,CAAC,EAAE,CAAC;QACvD,CAAC,EAAE,CAAC,EAAE,sBAAsB,CAAC,OAAO,CAAC,CAAC,EAAE,SAAS,CAAC,CAAC,CAAC,EAAE,CAAC;QACvD,CAAC,EAAE,CAAC,EAAE,CAAC,EAAE,CAAC;KACb,CAAC,CAAC;IAEH,wHAAwH;IACxH,MAAM,gBAAgB,GAAG,cAAc,CAAC,QAAQ,CAAC,KAAK,CAAC,CAAC;IACxD,OAAO,gBAAgB,CAAC,QAAQ,CAAC,qBAAqB,CAAC,CAAC;AAC5D,CAAC;AAED;;GAEG;AACH,MAAM,cAAc,GAAG,IAAI,OAAO,CAAC,OAAO,EAAE,GAAG,EAAE,OAAO,CAAC,CAAC;AAE1D;;;GAGG;AACH,kBAAkB;AAClB,MAAM,mBAAmB,GAAG,MAAM,CAAC,SAAS,CAAC;IACzC,SAAS,EAAE,SAAS,EAAE,SAAS,EAAE,CAAC;IAClC,SAAS,EAAE,SAAS,EAAE,QAAQ,EAAE,CAAC;IACjC,SAAS,EAAE,QAAQ,EAAE,SAAS,EAAE,CAAC;IACjC,CAAC,EAAE,CAAC,EAAE,CAAC,EAAE,CAAC;CACb,CAAC,CAAC;AAEH;;;;;;;;;;;;GAYG;AACH,MAAM,WAAW,GAAG,mBAAmB,CAAC,QAAQ,CAAC,aAAa,CAAC,cAAc,EAAE,oBAAoB,CAAC,IAAI,EAAE,CAAC,CAAC,CAAC,CAAC,CAAC;AAE/G;;GAEG;AACH,MAAM,WAAW,GAAG,MAAM,CAAC,MAAM,CAAC,WAAW,CAAC,CAAC;AAE/C;;;;;;;GAOG;AACH,MAAM,iBAAiB,GAAG,OAAO,CAAC,eAAe,CAAC,IAAI,OAAO,CAAC,CAAC,EAAE,CAAC,EAAE,CAAC,CAAC,EAAE,WAAW,CAAC,CAAC;AAErF;;;;;;;GAOG;AACH,MAAM,UAAU,qBAAqB,CAAC,iBAAyB,EAAE,IAAY;IACzE,MAAM,cAAc,GAAG,oBAAoB,CAAC,iBAAiB,EAAE,IAAI,CAAC,CAAC;IACrE,MAAM,KAAK,GAAG,aAAa,CAAC,cAAc,EAAE,iBAAiB,CAAC,CAAC;IAE/D,8FAA8F;IAC9F,MAAM,cAAc,GAAG,WAAW,CAAC,QAAQ,CAAC,KAAK,CAAC,CAAC;IACnD,MAAM,WAAW,GAAG,cAAc,CAAC,QAAQ,CAAC,WAAW,CAAC,CAAC;IAEzD,iGAAiG;IACjG,2EAA2E;IAC3E,OAAO,MAAM,CAAC,cAAc,CAAC,WAAW,CAAC,CAAC;AAC9C,CAAC","sourcesContent":["import { Clamp } from \"./math.scalar.functions\";\nimport { Matrix, Vector2, Vector3 } from \"./math.vector.pure\";\n\n/**\n * Tabulated Planckian locus in CIE 1960 UCS (u, v) space, with the local isotherm slope `t` used to offset\n * a chromaticity along the tint axis. Source: Wyszecki &amp; Stiles, \"Color Science\", 2nd ed., 1982, pp. 227-228.\n * Each row is [mired, u, v, t].\n */\nconst PlanckianLocusTable: ReadonlyArray<readonly [number, number, number, number]> = [\n    [0, 0.18006, 0.26352, -0.24341],\n    [10, 0.18066, 0.26589, -0.25479],\n    [20, 0.18133, 0.26846, -0.26876],\n    [30, 0.18208, 0.27119, -0.28539],\n    [40, 0.18293, 0.27407, -0.3047],\n    [50, 0.18388, 0.27709, -0.32675],\n    [60, 0.18494, 0.28021, -0.35156],\n    [70, 0.18611, 0.28342, -0.37915],\n    [80, 0.1874, 0.28668, -0.40955],\n    [90, 0.1888, 0.28997, -0.44278],\n    [100, 0.19032, 0.29326, -0.47888],\n    [125, 0.19462, 0.30141, -0.58204],\n    [150, 0.19962, 0.30921, -0.70471],\n    [175, 0.20525, 0.31647, -0.84901],\n    [200, 0.21142, 0.32312, -1.0182],\n    [225, 0.21807, 0.32909, -1.2168],\n    [250, 0.22511, 0.33439, -1.4512],\n    [275, 0.23247, 0.33904, -1.7298],\n    [300, 0.2401, 0.34308, -2.0637],\n    [325, 0.24792, 0.34655, -2.4681],\n    [350, 0.25591, 0.34951, -2.9641],\n    [375, 0.264, 0.352, -3.5814],\n    [400, 0.27218, 0.35407, -4.3633],\n    [425, 0.28039, 0.35577, -5.3762],\n    [450, 0.28863, 0.35714, -6.7262],\n    [475, 0.29685, 0.35823, -8.5955],\n    [500, 0.30505, 0.35907, -11.324],\n    [525, 0.3132, 0.35968, -15.628],\n    [550, 0.32129, 0.36011, -23.325],\n    [575, 0.32931, 0.36038, -40.77],\n    [600, 0.33724, 0.36051, -116.45],\n];\n\n/**\n * The minimum correlated color temperature, in Kelvin, representable by {@link TemperatureTintToXyz} (there is no\n * upper bound: increasingly high temperatures approach mired 0, already within the tabulated range).\n */\nexport const MinTemperatureKelvin = 1e6 / PlanckianLocusTable[PlanckianLocusTable.length - 1][0];\n\n/**\n * The maximum magnitude of the tint offset accepted by {@link TemperatureTintToXyz}, in either direction.\n */\nexport const MaxTintMagnitude = 150;\n\n/**\n * Converts a correlated color temperature and tint offset into the CIE XYZ (Y = 1) coordinates of the\n * corresponding illuminant white point, using a tabulated approximation of the Planckian locus in CIE 1960 UCS\n * (u, v) space.\n * @param temperatureKelvin The correlated color temperature of the illuminant, in Kelvin\n * @param tint An offset perpendicular to the Planckian locus (the green/magenta axis), in the range\n * [-{@link MaxTintMagnitude}, {@link MaxTintMagnitude}]\n * @returns The CIE XYZ (Y = 1) coordinates of the illuminant white point\n */\nexport function TemperatureTintToXyz(temperatureKelvin: number, tint: number): Vector3 {\n    // Clamp before taking the reciprocal: a non-positive or NaN Kelvin value has no valid reciprocal-mired\n    // representation, and taking 1e6 / temperatureKelvin first (before clamping the result) would map negative\n    // inputs to the *hottest* end of the table instead of the coldest/minimum one.\n    const safeTemperatureKelvin = Number.isNaN(temperatureKelvin) ? MinTemperatureKelvin : Math.max(temperatureKelvin, MinTemperatureKelvin);\n    const mired = Clamp(1e6 / safeTemperatureKelvin, PlanckianLocusTable[0][0], PlanckianLocusTable[PlanckianLocusTable.length - 1][0]);\n\n    // Binary search for the table row pair bracketing `mired` (the table is sorted ascending by mired).\n    let low = 0;\n    let high = PlanckianLocusTable.length - 1;\n    while (high - low > 1) {\n        const mid = (low + high) >> 1;\n        if (PlanckianLocusTable[mid][0] <= mired) {\n            low = mid;\n        } else {\n            high = mid;\n        }\n    }\n    const index = low;\n\n    const [lowMired, lowU, lowV, lowT] = PlanckianLocusTable[index];\n    const [highMired, highU, highV, highT] = PlanckianLocusTable[index + 1];\n    const f = (mired - lowMired) / (highMired - lowMired);\n\n    const uv = Vector2.Lerp(new Vector2(lowU, lowV), new Vector2(highU, highV), f);\n\n    const isothermLow = Vector2.Normalize(new Vector2(1, lowT));\n    const isothermHigh = Vector2.Normalize(new Vector2(1, highT));\n    const isotherm = Vector2.Normalize(Vector2.Lerp(isothermLow, isothermHigh, f));\n\n    const clampedTint = Clamp(tint, -MaxTintMagnitude, MaxTintMagnitude);\n    uv.subtractInPlace(isotherm.scale(clampedTint / 3000));\n\n    // CIE 1960 (u, v) -> CIE xyY (Y = 1) -> CIE XYZ.\n    const denom = 2 * uv.x - 8 * uv.y + 4;\n    const x = (3 * uv.x) / denom;\n    const y = (2 * uv.y) / denom;\n\n    return new Vector3(x / y, 1, (1 - x - y) / y);\n}\n\n// `Matrix.FromArray` reads values column-major, so each 3x3 matrix below is written column-by-column, not row-by-row.\n\n/**\n * The Bradford cone-response matrix used for chromatic adaptation.\n */\n// prettier-ignore\nconst BradfordMatrix = Matrix.FromArray([\n    0.8951, -0.7502, 0.0389, 0,\n    0.2664, 1.7135, -0.0685, 0,\n    -0.1614, 0.0367, 1.0296, 0,\n    0, 0, 0, 1,\n]);\n\nconst InverseBradfordMatrix = Matrix.Invert(BradfordMatrix);\n\n/**\n * Computes a per-channel Bradford adaptation ratio, clamped to avoid extreme temperature/tint combinations\n * blowing up the resulting matrix when a cone-response component approaches zero.\n * @param destComponent The destination white's cone-response component\n * @param sourceComponent The source white's cone-response component\n * @returns The ratio to use for this channel, clamped to [0.1, 10]\n */\nfunction ClampedAdaptationRatio(destComponent: number, sourceComponent: number): number {\n    const ratio = sourceComponent > 0 ? destComponent / sourceComponent : 10;\n    return Clamp(ratio, 0.1, 10);\n}\n\n/**\n * Computes the Bradford chromatic-adaptation matrix that maps CIE XYZ values from `sourceWhiteXYZ` to `destWhiteXYZ`.\n * @param sourceWhiteXYZ The CIE XYZ (Y = 1) coordinates of the source reference white\n * @param destWhiteXYZ The CIE XYZ (Y = 1) coordinates of the destination reference white\n * @returns A matrix mapping source-referenced XYZ to destination-referenced XYZ\n */\nfunction BradfordAdapt(sourceWhiteXYZ: Vector3, destWhiteXYZ: Vector3): Matrix {\n    const sourceLMS = Vector3.TransformNormal(sourceWhiteXYZ, BradfordMatrix);\n    const destLMS = Vector3.TransformNormal(destWhiteXYZ, BradfordMatrix);\n\n    // prettier-ignore\n    const scale = Matrix.FromArray([\n        ClampedAdaptationRatio(destLMS.x, sourceLMS.x), 0, 0, 0,\n        0, ClampedAdaptationRatio(destLMS.y, sourceLMS.y), 0, 0,\n        0, 0, ClampedAdaptationRatio(destLMS.z, sourceLMS.z), 0,\n        0, 0, 0, 1,\n    ]);\n\n    // adapt = InverseBradfordMatrix * scale * BradfordMatrix, written back-to-front since `A.multiply(B)` computes `B x A`.\n    const scaledByBradford = BradfordMatrix.multiply(scale);\n    return scaledByBradford.multiply(InverseBradfordMatrix);\n}\n\n/**\n * The standard CIE Illuminant D65 white point, in CIE XYZ (Y = 1).\n */\nconst StandardD65Xyz = new Vector3(0.95047, 1.0, 1.08883);\n\n/**\n * The standard sRGB / Rec.709 linear RGB to CIE XYZ matrix, calibrated to the externally-standardized CIE\n * Illuminant D65 chromaticity - i.e. `StandardRgbToXyzD65 * (1, 1, 1) = StandardD65Xyz` exactly.\n */\n// prettier-ignore\nconst StandardRgbToXyzD65 = Matrix.FromArray([\n    0.4124564, 0.2126729, 0.0193339, 0,\n    0.3575761, 0.7151522, 0.119192, 0,\n    0.1804375, 0.072175, 0.9503041, 0,\n    0, 0, 0, 1,\n]);\n\n/**\n * The linear RGB to CIE XYZ matrix actually used by {@link GetWhiteBalanceMatrix}, recalibrated from\n * {@link StandardRgbToXyzD65} via a single Bradford adaptation step so that neutral (1, 1, 1) maps exactly onto\n * this module's own Planckian-locus approximation of 6500 K (see {@link TemperatureTintToXyz}) instead of the\n * externally-standardized CIE D65 chromaticity. Those two are close but never exactly equal, by definition: D65\n * is defined via the CIE daylight locus, a model fit to measured daylight spectra, which deliberately deviates\n * from the Planckian locus (real daylight isn't quite blackbody radiation) - no table or precision fix can close\n * that gap without distorting temperature/tint for every other illuminant.\n *\n * This is what lets `GetWhiteBalanceMatrix(6500, 0)` be exactly the identity matrix - the default \"no adjustment\"\n * temperature/tint are a true no-op - without needing to special-case that one input: the reverse adaptation is\n * applied once, here, to the matrix itself, rather than to every call's result.\n */\nconst RgbToXyzD65 = StandardRgbToXyzD65.multiply(BradfordAdapt(StandardD65Xyz, TemperatureTintToXyz(6500, 0)));\n\n/**\n * The CIE XYZ to linear RGB matrix, the inverse of {@link RgbToXyzD65}.\n */\nconst XyzToRgbD65 = Matrix.Invert(RgbToXyzD65);\n\n/**\n * The working color space's reference white, in CIE XYZ (Y = 1) - what neutral (1, 1, 1) linear RGB maps to\n * through {@link RgbToXyzD65}. Since that matrix is itself calibrated to this module's own 6500 K point (see\n * above), this is exactly `TemperatureTintToXyz(6500, 0)`. Deriving it as \"whatever (1, 1, 1) maps to\" (rather\n * than calling `TemperatureTintToXyz` a second time) keeps the invariant self-documenting: a fully white-balanced\n * neutral patch round-trips back to exactly (1, 1, 1) for any input illuminant, because the adaptation target is\n * always tied directly to whatever the working color space itself already treats as neutral.\n */\nconst ReferenceWhiteXyz = Vector3.TransformNormal(new Vector3(1, 1, 1), RgbToXyzD65);\n\n/**\n * Computes the linear RGB (sRGB / Rec.709 primaries) color-correction matrix that white-balances the given\n * illuminant, by chromatically adapting its white point (see {@link TemperatureTintToXyz}) to the working\n * color space's reference white using the Bradford transform.\n * @param temperatureKelvin The correlated color temperature of the illuminant to neutralize, in Kelvin\n * @param tint An offset perpendicular to the Planckian locus (the green/magenta axis), in the range [-150, 150]\n * @returns A column-major 3x3 matrix (9 values), ready to be bound as a `mat3` shader uniform\n */\nexport function GetWhiteBalanceMatrix(temperatureKelvin: number, tint: number): Float32Array | Array<number> {\n    const targetWhiteXYZ = TemperatureTintToXyz(temperatureKelvin, tint);\n    const adapt = BradfordAdapt(targetWhiteXYZ, ReferenceWhiteXyz);\n\n    // finalMatrix = XyzToRgbD65 * adapt * RgbToXyzD65, written back-to-front (see BradfordAdapt).\n    const adaptedFromRgb = RgbToXyzD65.multiply(adapt);\n    const finalMatrix = adaptedFromRgb.multiply(XyzToRgbD65);\n\n    // `Matrix` stores its values column-major, matching the layout GLSL/WGSL `mat3` uniforms expect,\n    // so the upper-left 3x3 block can be extracted directly with no transpose.\n    return Matrix.GetAsMatrix3x3(finalMatrix);\n}\n"]}

package/Maths/index.d.ts

@@ -18,3 +18,4 @@ export * from "./ThinMaths/index.js";
 export * from "./math.scalar.functions.js";
 export * from "./math.geospatial.js";
 export * from "./math.geospatial.functions.js";
+export * from "./colorTemperature.functions.js";

package/Maths/index.js

@@ -19,4 +19,5 @@ export * from "./ThinMaths/index.js";
 export * from "./math.scalar.functions.js";
 export * from "./math.geospatial.js";
 export * from "./math.geospatial.functions.js";
+export * from "./colorTemperature.functions.js";
 //# sourceMappingURL=index.js.map

package/Maths/index.js.map

@@ -1 +1 @@
-{"version":3,"file":"index.js","sourceRoot":"","sources":["../../../../dev/core/src/Maths/index.ts"],"names":[],"mappings":"AAAA,6DAA6D;AAC7D,cAAc,eAAe,CAAC;AAC9B,cAAc,kBAAkB,CAAC;AACjC,cAAc,cAAc,CAAC;AAC7B,cAAc,aAAa,CAAC;AAC5B,cAAc,aAAa,CAAC;AAC5B,cAAc,kBAAkB,CAAC;AACjC,cAAc,gBAAgB,CAAC;AAC/B,cAAc,aAAa,CAAC;AAC5B,cAAc,cAAc,CAAC;AAC7B,cAAc,aAAa,CAAC;AAC5B,cAAc,qBAAqB,CAAC;AACpC,cAAc,iBAAiB,CAAC;AAChC,cAAc,QAAQ,CAAC;AACvB,cAAc,uBAAuB,CAAC;AACtC,cAAc,oBAAoB,CAAC;AACnC,cAAc,yBAAyB,CAAC;AACxC,cAAc,mBAAmB,CAAC;AAClC,cAAc,yBAAyB,CAAC;AACxC,cAAc,mBAAmB,CAAC;AAClC,cAAc,6BAA6B,CAAC","sourcesContent":["/* eslint-disable @typescript-eslint/no-restricted-imports */\r\nexport * from \"./math.scalar\";\r\nexport * from \"./math.functions\";\r\nexport * from \"./math.polar\";\r\nexport * from \"./math.like\";\r\nexport * from \"./math.axis\";\r\nexport * from \"./math.constants\";\r\nexport * from \"./math.frustum\";\r\nexport * from \"./math.path\";\r\nexport * from \"./math.plane\";\r\nexport * from \"./math.size\";\r\nexport * from \"./math.vertexFormat\";\r\nexport * from \"./math.viewport\";\r\nexport * from \"./math\";\r\nexport * from \"./sphericalPolynomial\";\r\nexport * from \"./halton2DSequence\";\r\nexport * from \"./math.vector.functions\";\r\nexport * from \"./ThinMaths/index\";\r\nexport * from \"./math.scalar.functions\";\r\nexport * from \"./math.geospatial\";\r\nexport * from \"./math.geospatial.functions\";\r\n"]}
+{"version":3,"file":"index.js","sourceRoot":"","sources":["../../../../dev/core/src/Maths/index.ts"],"names":[],"mappings":"AAAA,6DAA6D;AAC7D,cAAc,eAAe,CAAC;AAC9B,cAAc,kBAAkB,CAAC;AACjC,cAAc,cAAc,CAAC;AAC7B,cAAc,aAAa,CAAC;AAC5B,cAAc,aAAa,CAAC;AAC5B,cAAc,kBAAkB,CAAC;AACjC,cAAc,gBAAgB,CAAC;AAC/B,cAAc,aAAa,CAAC;AAC5B,cAAc,cAAc,CAAC;AAC7B,cAAc,aAAa,CAAC;AAC5B,cAAc,qBAAqB,CAAC;AACpC,cAAc,iBAAiB,CAAC;AAChC,cAAc,QAAQ,CAAC;AACvB,cAAc,uBAAuB,CAAC;AACtC,cAAc,oBAAoB,CAAC;AACnC,cAAc,yBAAyB,CAAC;AACxC,cAAc,mBAAmB,CAAC;AAClC,cAAc,yBAAyB,CAAC;AACxC,cAAc,mBAAmB,CAAC;AAClC,cAAc,6BAA6B,CAAC;AAC5C,cAAc,8BAA8B,CAAC","sourcesContent":["/* eslint-disable @typescript-eslint/no-restricted-imports */\r\nexport * from \"./math.scalar\";\r\nexport * from \"./math.functions\";\r\nexport * from \"./math.polar\";\r\nexport * from \"./math.like\";\r\nexport * from \"./math.axis\";\r\nexport * from \"./math.constants\";\r\nexport * from \"./math.frustum\";\r\nexport * from \"./math.path\";\r\nexport * from \"./math.plane\";\r\nexport * from \"./math.size\";\r\nexport * from \"./math.vertexFormat\";\r\nexport * from \"./math.viewport\";\r\nexport * from \"./math\";\r\nexport * from \"./sphericalPolynomial\";\r\nexport * from \"./halton2DSequence\";\r\nexport * from \"./math.vector.functions\";\r\nexport * from \"./ThinMaths/index\";\r\nexport * from \"./math.scalar.functions\";\r\nexport * from \"./math.geospatial\";\r\nexport * from \"./math.geospatial.functions\";\r\nexport * from \"./colorTemperature.functions\";\r\n"]}

package/Maths/math.vector.pure.d.ts

@@ -842,11 +842,7 @@ export declare class Vector3 implements Vector<Tuple<number, 3>, IVector3LikeInt
      * Intercept the x,y,z value change
      * @param interceptor
      */
-    setInterceptor(interceptor?: (newValue: {
-        x?: number;
-        y?: number;
-        z?: number;
-    }) => boolean): void;
+    setInterceptor(interceptor?: VectorInterceptor): void;
     /**
      * Get the current interceptor
      * @returns the current vector interceptor