image-js 1.0.0-alpha.3 → 1.1.0

package/lib/compare/subtract.js.map

@@ -0,0 +1 @@
+{"version":3,"file":"subtract.js","sourceRoot":"","sources":["../../src/compare/subtract.ts"],"names":[],"mappings":"AAAA,OAAO,EAAE,KAAK,EAAE,MAAM,aAAa,CAAC;AACpC,OAAO,EAAE,IAAI,EAAE,MAAM,YAAY,CAAC;AAClC,OAAO,gBAAgB,MAAM,yCAAyC,CAAC;AACvE,OAAO,EAAE,qBAAqB,EAAE,MAAM,mCAAmC,CAAC;AAyB1E;;;;;;GAMG;AACH,MAAM,UAAU,QAAQ,CACtB,KAAmB,EACnB,UAAwB,EACxB,UAAgC,EAAE;IAElC,MAAM,EAAE,QAAQ,GAAG,KAAK,EAAE,GAAG,OAAO,CAAC;IAErC,IAAI,KAAK,YAAY,KAAK,EAAE,CAAC;QAC3B,gBAAgB,CAAC,KAAK,EAAE;YACtB,QAAQ,EAAE,CAAC,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC;YACpB,UAAU,EAAE,CAAC,CAAC,EAAE,CAAC,CAAC;YAClB,KAAK,EAAE,KAAK;SACb,CAAC,CAAC;IACL,CAAC;IAED,qBAAqB,CAAC,KAAK,EAAE,UAAU,CAAC,CAAC;IAEzC,MAAM,QAAQ,GAAG,KAAK,CAAC,KAAK,EAAE,CAAC;IAC/B,IAAI,QAAQ,YAAY,KAAK,EAAE,CAAC;QAC9B,KAAK,IAAI,KAAK,GAAG,CAAC,EAAE,KAAK,GAAG,KAAK,CAAC,IAAI,EAAE,KAAK,EAAE,EAAE,CAAC;YAChD,KAAK,IAAI,OAAO,GAAG,CAAC,EAAE,OAAO,GAAG,KAAK,CAAC,QAAQ,EAAE,OAAO,EAAE,EAAE,CAAC;gBAC1D,MAAM,KAAK,GACT,KAAK,CAAC,eAAe,CAAC,KAAK,EAAE,OAAO,CAAC;oBACrC,UAAU,CAAC,eAAe,CAAC,KAAK,EAAE,OAAO,CAAC,CAAC;gBAC7C,IAAI,QAAQ,EAAE,CAAC;oBACb,IAAI,KAAK,GAAG,CAAC,EAAE,CAAC;wBACd,QAAQ,CAAC,eAAe,CAAC,KAAK,EAAE,OAAO,EAAE,CAAC,KAAK,CAAC,CAAC;oBACnD,CAAC;yBAAM,CAAC;wBACN,QAAQ,CAAC,eAAe,CAAC,KAAK,EAAE,OAAO,EAAE,KAAK,CAAC,CAAC;oBAClD,CAAC;gBACH,CAAC;qBAAM,IAAI,KAAK,GAAG,CAAC,EAAE,CAAC;oBACrB,QAAQ,CAAC,eAAe,CAAC,KAAK,EAAE,OAAO,EAAE,CAAC,CAAC,CAAC;gBAC9C,CAAC;qBAAM,CAAC;oBACN,QAAQ,CAAC,eAAe,CAAC,KAAK,EAAE,OAAO,EAAE,KAAK,CAAC,CAAC;gBAClD,CAAC;YACH,CAAC;QACH,CAAC;IACH,CAAC;SAAM,IAAI,KAAK,YAAY,IAAI,IAAI,UAAU,YAAY,IAAI,EAAE,CAAC;QAC/D,KAAK,IAAI,KAAK,GAAG,CAAC,EAAE,KAAK,GAAG,KAAK,CAAC,IAAI,EAAE,KAAK,EAAE,EAAE,CAAC;YAChD,MAAM,KAAK,GACT,KAAK,CAAC,aAAa,CAAC,KAAK,CAAC,GAAG,UAAU,CAAC,aAAa,CAAC,KAAK,CAAC,CAAC;YAC/D,IAAI,QAAQ,EAAE,CAAC;gBACb,QAAQ,CAAC,aAAa,CAAC,KAAK,EAAE,KAAK,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;YAC/C,CAAC;iBAAM,CAAC;gBACN,QAAQ,CAAC,aAAa,CAAC,KAAK,EAAE,KAAK,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;YACnD,CAAC;QACH,CAAC;IACH,CAAC;IAED,OAAO,QAAQ,CAAC;AAClB,CAAC"}

package/lib/compute/getExtrema.d.ts

@@ -0,0 +1,33 @@
+import type { Image } from '../Image.js';
+import type { Mask } from '../Mask.js';
+import type { Point } from '../geometry/index.js';
+export interface ExtremaOptions {
+    /**
+     * Chooses what kind of extremum to compute.
+     * @default `'maximum'`
+     */
+    kind?: 'minimum' | 'maximum';
+    /**
+     * Uses mask to check if a point belongs to a ROI or not
+     * @default `undefined`
+     */
+    mask?: Mask;
+    /**
+     * Chooses what kind of coverage algorithm to use to compute the extremum.
+     * @default `'star'`
+     */
+    algorithm?: 'cross' | 'square' | 'star';
+    /**
+     * Maximum number of points that can be equal to the extremum
+     * @default `2`
+     */
+    maxEquals?: number;
+}
+/**
+ * Checks the surrounding values of a point. If they are all bigger or smaller than the pixel in question then this point is considered an extremum.
+ * @param image - Image to find extrema from.
+ * @param options - ExtremaOptions
+ * @returns Array of Points.
+ */
+export declare function getExtrema(image: Image, options: ExtremaOptions): Point[];
+//# sourceMappingURL=getExtrema.d.ts.map

package/lib/compute/getExtrema.d.ts.map

@@ -0,0 +1 @@
+{"version":3,"file":"getExtrema.d.ts","sourceRoot":"","sources":["../../src/compute/getExtrema.ts"],"names":[],"mappings":"AAAA,OAAO,KAAK,EAAE,KAAK,EAAE,MAAM,aAAa,CAAC;AACzC,OAAO,KAAK,EAAE,IAAI,EAAE,MAAM,YAAY,CAAC;AACvC,OAAO,KAAK,EAAE,KAAK,EAAE,MAAM,sBAAsB,CAAC;AAIlD,MAAM,WAAW,cAAc;IAC7B;;;OAGG;IACH,IAAI,CAAC,EAAE,SAAS,GAAG,SAAS,CAAC;IAC7B;;;OAGG;IACH,IAAI,CAAC,EAAE,IAAI,CAAC;IACZ;;;OAGG;IACH,SAAS,CAAC,EAAE,OAAO,GAAG,QAAQ,GAAG,MAAM,CAAC;IACxC;;;OAGG;IACH,SAAS,CAAC,EAAE,MAAM,CAAC;CACpB;AACD;;;;;GAKG;AACH,wBAAgB,UAAU,CAAC,KAAK,EAAE,KAAK,EAAE,OAAO,EAAE,cAAc,GAAG,KAAK,EAAE,CAoEzE"}

package/lib/compute/getExtrema.js

@@ -0,0 +1,72 @@
+import { assertUnreachable } from '../utils/validators/assert.js';
+import checkProcessable from '../utils/validators/checkProcessable.js';
+/**
+ * Checks the surrounding values of a point. If they are all bigger or smaller than the pixel in question then this point is considered an extremum.
+ * @param image - Image to find extrema from.
+ * @param options - ExtremaOptions
+ * @returns Array of Points.
+ */
+export function getExtrema(image, options) {
+    const { kind = 'maximum', mask, algorithm = 'star', maxEquals = 2 } = options;
+    checkProcessable(image, {
+        bitDepth: [8, 16],
+    });
+    const searchingMinimum = kind === 'minimum';
+    const maskExpectedValue = searchingMinimum ? 0 : 1;
+    const dx = [1, 0, -1, 0, 1, 1, -1, -1, 2, 0, -2, 0, 2, 2, -2, -2];
+    const dy = [0, 1, 0, -1, 1, -1, 1, -1, 0, 2, 0, -2, 2, -2, 2, -2];
+    switch (algorithm) {
+        case 'cross':
+            dx.length = 4;
+            dy.length = 4;
+            break;
+        case 'square':
+            dx.length = 8;
+            dy.length = 8;
+            break;
+        case 'star':
+            break;
+        default:
+            assertUnreachable(algorithm);
+    }
+    const shift = dx.length <= 8 ? 1 : 2; // deal with borders
+    const points = [];
+    for (let channel = 0; channel < image.channels; channel++) {
+        for (let currentY = shift; currentY < image.height - shift; currentY++) {
+            for (let currentX = shift; currentX < image.width - shift; currentX++) {
+                if (mask && mask.getBit(currentX, currentY) !== maskExpectedValue) {
+                    continue;
+                }
+                let counter = 0;
+                let nbEquals = 0;
+                const currentValue = image.getValue(currentX, currentY, channel);
+                for (let dir = 0; dir < dx.length; dir++) {
+                    const currentAroundValue = image.getValue(currentX + dx[dir], currentY + dy[dir], channel);
+                    if (searchingMinimum) {
+                        // we search for minima
+                        if (currentAroundValue > currentValue) {
+                            counter++;
+                        }
+                    }
+                    else if (currentAroundValue < currentValue) {
+                        counter++;
+                    }
+                    if (currentAroundValue === currentValue) {
+                        nbEquals++;
+                    }
+                }
+                if (counter + nbEquals === dx.length && nbEquals <= maxEquals) {
+                    points.push({ column: currentX, row: currentY });
+                }
+            }
+        }
+    }
+    // TODO How to make a more performant and general way
+    // we don't deal correctly here with groups of points that should be grouped if at the
+    // beginning one of them is closer to another
+    // Seems that we would ened to calculate a matrix and then split this matrix in 'independant matrices'
+    // Or to assign a cluster to each point and regroup them if 2 clusters are close to each other
+    // later approach seems much better
+    return points;
+}
+//# sourceMappingURL=getExtrema.js.map

package/lib/compute/getExtrema.js.map

@@ -0,0 +1 @@
+{"version":3,"file":"getExtrema.js","sourceRoot":"","sources":["../../src/compute/getExtrema.ts"],"names":[],"mappings":"AAGA,OAAO,EAAE,iBAAiB,EAAE,MAAM,+BAA+B,CAAC;AAClE,OAAO,gBAAgB,MAAM,yCAAyC,CAAC;AAwBvE;;;;;GAKG;AACH,MAAM,UAAU,UAAU,CAAC,KAAY,EAAE,OAAuB;IAC9D,MAAM,EAAE,IAAI,GAAG,SAAS,EAAE,IAAI,EAAE,SAAS,GAAG,MAAM,EAAE,SAAS,GAAG,CAAC,EAAE,GAAG,OAAO,CAAC;IAC9E,gBAAgB,CAAC,KAAK,EAAE;QACtB,QAAQ,EAAE,CAAC,CAAC,EAAE,EAAE,CAAC;KAClB,CAAC,CAAC;IACH,MAAM,gBAAgB,GAAG,IAAI,KAAK,SAAS,CAAC;IAE5C,MAAM,iBAAiB,GAAG,gBAAgB,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;IAEnD,MAAM,EAAE,GAAG,CAAC,CAAC,EAAE,CAAC,EAAE,CAAC,CAAC,EAAE,CAAC,EAAE,CAAC,EAAE,CAAC,EAAE,CAAC,CAAC,EAAE,CAAC,CAAC,EAAE,CAAC,EAAE,CAAC,EAAE,CAAC,CAAC,EAAE,CAAC,EAAE,CAAC,EAAE,CAAC,EAAE,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC;IAClE,MAAM,EAAE,GAAG,CAAC,CAAC,EAAE,CAAC,EAAE,CAAC,EAAE,CAAC,CAAC,EAAE,CAAC,EAAE,CAAC,CAAC,EAAE,CAAC,EAAE,CAAC,CAAC,EAAE,CAAC,EAAE,CAAC,EAAE,CAAC,EAAE,CAAC,CAAC,EAAE,CAAC,EAAE,CAAC,CAAC,EAAE,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC;IAClE,QAAQ,SAAS,EAAE,CAAC;QAClB,KAAK,OAAO;YACV,EAAE,CAAC,MAAM,GAAG,CAAC,CAAC;YACd,EAAE,CAAC,MAAM,GAAG,CAAC,CAAC;YACd,MAAM;QACR,KAAK,QAAQ;YACX,EAAE,CAAC,MAAM,GAAG,CAAC,CAAC;YACd,EAAE,CAAC,MAAM,GAAG,CAAC,CAAC;YACd,MAAM;QACR,KAAK,MAAM;YACT,MAAM;QACR;YACE,iBAAiB,CAAC,SAAS,CAAC,CAAC;IACjC,CAAC;IACD,MAAM,KAAK,GAAG,EAAE,CAAC,MAAM,IAAI,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,oBAAoB;IAC1D,MAAM,MAAM,GAAY,EAAE,CAAC;IAC3B,KAAK,IAAI,OAAO,GAAG,CAAC,EAAE,OAAO,GAAG,KAAK,CAAC,QAAQ,EAAE,OAAO,EAAE,EAAE,CAAC;QAC1D,KAAK,IAAI,QAAQ,GAAG,KAAK,EAAE,QAAQ,GAAG,KAAK,CAAC,MAAM,GAAG,KAAK,EAAE,QAAQ,EAAE,EAAE,CAAC;YACvE,KAAK,IAAI,QAAQ,GAAG,KAAK,EAAE,QAAQ,GAAG,KAAK,CAAC,KAAK,GAAG,KAAK,EAAE,QAAQ,EAAE,EAAE,CAAC;gBACtE,IAAI,IAAI,IAAI,IAAI,CAAC,MAAM,CAAC,QAAQ,EAAE,QAAQ,CAAC,KAAK,iBAAiB,EAAE,CAAC;oBAClE,SAAS;gBACX,CAAC;gBACD,IAAI,OAAO,GAAG,CAAC,CAAC;gBAChB,IAAI,QAAQ,GAAG,CAAC,CAAC;gBACjB,MAAM,YAAY,GAAG,KAAK,CAAC,QAAQ,CAAC,QAAQ,EAAE,QAAQ,EAAE,OAAO,CAAC,CAAC;gBACjE,KAAK,IAAI,GAAG,GAAG,CAAC,EAAE,GAAG,GAAG,EAAE,CAAC,MAAM,EAAE,GAAG,EAAE,EAAE,CAAC;oBACzC,MAAM,kBAAkB,GAAG,KAAK,CAAC,QAAQ,CACvC,QAAQ,GAAG,EAAE,CAAC,GAAG,CAAC,EAClB,QAAQ,GAAG,EAAE,CAAC,GAAG,CAAC,EAClB,OAAO,CACR,CAAC;oBACF,IAAI,gBAAgB,EAAE,CAAC;wBACrB,uBAAuB;wBACvB,IAAI,kBAAkB,GAAG,YAAY,EAAE,CAAC;4BACtC,OAAO,EAAE,CAAC;wBACZ,CAAC;oBACH,CAAC;yBAAM,IAAI,kBAAkB,GAAG,YAAY,EAAE,CAAC;wBAC7C,OAAO,EAAE,CAAC;oBACZ,CAAC;oBACD,IAAI,kBAAkB,KAAK,YAAY,EAAE,CAAC;wBACxC,QAAQ,EAAE,CAAC;oBACb,CAAC;gBACH,CAAC;gBACD,IAAI,OAAO,GAAG,QAAQ,KAAK,EAAE,CAAC,MAAM,IAAI,QAAQ,IAAI,SAAS,EAAE,CAAC;oBAC9D,MAAM,CAAC,IAAI,CAAC,EAAE,MAAM,EAAE,QAAQ,EAAE,GAAG,EAAE,QAAQ,EAAE,CAAC,CAAC;gBACnD,CAAC;YACH,CAAC;QACH,CAAC;IACH,CAAC;IAED,qDAAqD;IACrD,sFAAsF;IACtF,6CAA6C;IAC7C,sGAAsG;IACtG,8FAA8F;IAC9F,mCAAmC;IACnC,OAAO,MAAM,CAAC;AAChB,CAAC"}

package/lib/compute/histogram.d.ts

@@ -1,15 +1,23 @@
-import { Image } from '../Image';
-export interface IHistogramOptions {
-    /**
-     * The channel for which to compute the histogram.
-     * If it is unspecified, the image must have one channel or the method will
-     * throw an error.
-     */
-    channel?: number;
-}
-/**
- * Returns a histogram of pixel intensities.
- * @param image
- */
-export declare function histogram(image: Image, options?: IHistogramOptions): number[];
+import type { Image } from '../Image.js';
+export interface HistogramOptions {
+    /**
+     * The channel for which to compute the histogram.
+     * If it is unspecified, the image must have one channel or the method will
+     * throw an error.
+     * @default `0`
+     */
+    channel?: number;
+    /**
+     * The number of slots that histogram can have.
+     * @default 2 ** image.bitDepth
+     */
+    slots?: number;
+}
+/**
+ * Returns a histogram of pixel intensities.
+ * @param image - The original image.
+ * @param options - Histogram options.
+ * @returns - The histogram.
+ */
+export declare function histogram(image: Image, options?: HistogramOptions): Uint32Array;
 //# sourceMappingURL=histogram.d.ts.map

package/lib/compute/histogram.d.ts.map

@@ -1 +1 @@
-{"version":3,"file":"histogram.d.ts","sourceRoot":"","sources":["../../src/compute/histogram.ts"],"names":[],"mappings":"AAAA,OAAO,EAAE,KAAK,EAAE,MAAM,UAAU,CAAC;AAGjC,MAAM,WAAW,iBAAiB;IAChC;;;;OAIG;IACH,OAAO,CAAC,EAAE,MAAM,CAAC;CAClB;AAED;;;GAGG;AACH,wBAAgB,SAAS,CACvB,KAAK,EAAE,KAAK,EACZ,OAAO,GAAE,iBAAsB,GAC9B,MAAM,EAAE,CAcV"}
+{"version":3,"file":"histogram.d.ts","sourceRoot":"","sources":["../../src/compute/histogram.ts"],"names":[],"mappings":"AAAA,OAAO,KAAK,EAAE,KAAK,EAAE,MAAM,aAAa,CAAC;AAGzC,MAAM,WAAW,gBAAgB;IAC/B;;;;;OAKG;IACH,OAAO,CAAC,EAAE,MAAM,CAAC;IACjB;;;OAGG;IACH,KAAK,CAAC,EAAE,MAAM,CAAC;CAChB;AAED;;;;;GAKG;AACH,wBAAgB,SAAS,CACvB,KAAK,EAAE,KAAK,EACZ,OAAO,GAAE,gBAAqB,GAC7B,WAAW,CA4Bb"}

package/lib/compute/histogram.js

@@ -1,24 +1,30 @@
-"use strict";
-Object.defineProperty(exports, "__esModule", { value: true });
-const validators_1 = require("../utils/validators");
-/**
- * Returns a histogram of pixel intensities.
- * @param image
- */
-function histogram(image, options = {}) {
-    let { channel } = options;
-    if (typeof channel !== 'number') {
-        if (image.channels !== 1) {
-            throw new Error('channel option is mandatory for multi-channel images');
-        }
-        channel = 0;
-    }
-    validators_1.validateChannel(channel, image);
-    const hist = new Array(image.maxValue + 1).fill(0);
-    for (let i = channel; i < image.data.length; i += image.channels) {
-        hist[image.data[i]]++;
-    }
-    return hist;
-}
-exports.histogram = histogram;
+import { validateChannel } from '../utils/validators/validators.js';
+/**
+ * Returns a histogram of pixel intensities.
+ * @param image - The original image.
+ * @param options - Histogram options.
+ * @returns - The histogram.
+ */
+export function histogram(image, options = {}) {
+    let { channel } = options;
+    const { slots = 2 ** image.bitDepth } = options;
+    if (!(slots !== 0 && (slots & (slots - 1)) === 0)) {
+        throw new RangeError('slots must be a power of 2, for example: 64, 256, 1024');
+    }
+    if (typeof channel !== 'number') {
+        if (image.channels !== 1) {
+            throw new TypeError('channel option is mandatory for multi-channel images');
+        }
+        channel = 0;
+    }
+    validateChannel(channel, image);
+    const hist = new Uint32Array(slots);
+    let bitShift = 0;
+    const bitSlots = Math.log2(slots);
+    bitShift = image.bitDepth - bitSlots;
+    for (let i = 0; i < image.size; i++) {
+        hist[image.getValueByIndex(i, channel) >> bitShift]++;
+    }
+    return hist;
+}
 //# sourceMappingURL=histogram.js.map

package/lib/compute/histogram.js.map

@@ -1 +1 @@
-{"version":3,"file":"histogram.js","sourceRoot":"","sources":["../../src/compute/histogram.ts"],"names":[],"mappings":";;AACA,oDAAsD;AAWtD;;;GAGG;AACH,SAAgB,SAAS,CACvB,KAAY,EACZ,UAA6B,EAAE;IAE/B,IAAI,EAAE,OAAO,EAAE,GAAG,OAAO,CAAC;IAC1B,IAAI,OAAO,OAAO,KAAK,QAAQ,EAAE;QAC/B,IAAI,KAAK,CAAC,QAAQ,KAAK,CAAC,EAAE;YACxB,MAAM,IAAI,KAAK,CAAC,sDAAsD,CAAC,CAAC;SACzE;QACD,OAAO,GAAG,CAAC,CAAC;KACb;IACD,4BAAe,CAAC,OAAO,EAAE,KAAK,CAAC,CAAC;IAChC,MAAM,IAAI,GAAG,IAAI,KAAK,CAAC,KAAK,CAAC,QAAQ,GAAG,CAAC,CAAC,CAAC,IAAI,CAAC,CAAC,CAAC,CAAC;IACnD,KAAK,IAAI,CAAC,GAAG,OAAO,EAAE,CAAC,GAAG,KAAK,CAAC,IAAI,CAAC,MAAM,EAAE,CAAC,IAAI,KAAK,CAAC,QAAQ,EAAE;QAChE,IAAI,CAAC,KAAK,CAAC,IAAI,CAAC,CAAC,CAAC,CAAC,EAAE,CAAC;KACvB;IACD,OAAO,IAAI,CAAC;AACd,CAAC;AAjBD,8BAiBC"}
+{"version":3,"file":"histogram.js","sourceRoot":"","sources":["../../src/compute/histogram.ts"],"names":[],"mappings":"AACA,OAAO,EAAE,eAAe,EAAE,MAAM,mCAAmC,CAAC;AAiBpE;;;;;GAKG;AACH,MAAM,UAAU,SAAS,CACvB,KAAY,EACZ,UAA4B,EAAE;IAE9B,IAAI,EAAE,OAAO,EAAE,GAAG,OAAO,CAAC;IAC1B,MAAM,EAAE,KAAK,GAAG,CAAC,IAAI,KAAK,CAAC,QAAQ,EAAE,GAAG,OAAO,CAAC;IAChD,IAAI,CAAC,CAAC,KAAK,KAAK,CAAC,IAAI,CAAC,KAAK,GAAG,CAAC,KAAK,GAAG,CAAC,CAAC,CAAC,KAAK,CAAC,CAAC,EAAE,CAAC;QAClD,MAAM,IAAI,UAAU,CAClB,wDAAwD,CACzD,CAAC;IACJ,CAAC;IACD,IAAI,OAAO,OAAO,KAAK,QAAQ,EAAE,CAAC;QAChC,IAAI,KAAK,CAAC,QAAQ,KAAK,CAAC,EAAE,CAAC;YACzB,MAAM,IAAI,SAAS,CACjB,sDAAsD,CACvD,CAAC;QACJ,CAAC;QACD,OAAO,GAAG,CAAC,CAAC;IACd,CAAC;IACD,eAAe,CAAC,OAAO,EAAE,KAAK,CAAC,CAAC;IAEhC,MAAM,IAAI,GAAG,IAAI,WAAW,CAAC,KAAK,CAAC,CAAC;IAEpC,IAAI,QAAQ,GAAG,CAAC,CAAC;IACjB,MAAM,QAAQ,GAAG,IAAI,CAAC,IAAI,CAAC,KAAK,CAAC,CAAC;IAClC,QAAQ,GAAG,KAAK,CAAC,QAAQ,GAAG,QAAQ,CAAC;IACrC,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,KAAK,CAAC,IAAI,EAAE,CAAC,EAAE,EAAE,CAAC;QACpC,IAAI,CAAC,KAAK,CAAC,eAAe,CAAC,CAAC,EAAE,OAAO,CAAC,IAAI,QAAQ,CAAC,EAAE,CAAC;IACxD,CAAC;IAED,OAAO,IAAI,CAAC;AACd,CAAC"}

package/lib/compute/index.d.ts

@@ -1,2 +1,6 @@
-export * from './histogram';
+export * from './mean.js';
+export * from './histogram.js';
+export * from './median.js';
+export * from './getExtrema.js';
+export * from './variance.js';
 //# sourceMappingURL=index.d.ts.map

package/lib/compute/index.d.ts.map

@@ -1 +1 @@
-{"version":3,"file":"index.d.ts","sourceRoot":"","sources":["../../src/compute/index.ts"],"names":[],"mappings":"AAAA,cAAc,aAAa,CAAC"}
+{"version":3,"file":"index.d.ts","sourceRoot":"","sources":["../../src/compute/index.ts"],"names":[],"mappings":"AAAA,cAAc,WAAW,CAAC;AAC1B,cAAc,gBAAgB,CAAC;AAC/B,cAAc,aAAa,CAAC;AAC5B,cAAc,iBAAiB,CAAC;AAChC,cAAc,eAAe,CAAC"}

package/lib/compute/index.js

@@ -1,7 +1,6 @@
-"use strict";
-function __export(m) {
-    for (var p in m) if (!exports.hasOwnProperty(p)) exports[p] = m[p];
-}
-Object.defineProperty(exports, "__esModule", { value: true });
-__export(require("./histogram"));
+export * from './mean.js';
+export * from './histogram.js';
+export * from './median.js';
+export * from './getExtrema.js';
+export * from './variance.js';
 //# sourceMappingURL=index.js.map

package/lib/compute/index.js.map

@@ -1 +1 @@
-{"version":3,"file":"index.js","sourceRoot":"","sources":["../../src/compute/index.ts"],"names":[],"mappings":";;;;;AAAA,iCAA4B"}
+{"version":3,"file":"index.js","sourceRoot":"","sources":["../../src/compute/index.ts"],"names":[],"mappings":"AAAA,cAAc,WAAW,CAAC;AAC1B,cAAc,gBAAgB,CAAC;AAC/B,cAAc,aAAa,CAAC;AAC5B,cAAc,iBAAiB,CAAC;AAChC,cAAc,eAAe,CAAC"}

package/lib/compute/mean.d.ts

@@ -0,0 +1,17 @@
+import type { Image } from '../Image.js';
+import type { Point } from '../geometry/index.js';
+export interface MeanOptions {
+    /**
+     * Points to calculate mean from.
+     */
+    points: Point[];
+}
+/**
+ * Compute the mean of an image. The mean can be either computed on each channel
+ * individually or on the whole image.
+ * @param image - Image to process.
+ * @param options - Mean options.
+ * @returns The mean pixel.
+ */
+export declare function mean(image: Image, options?: MeanOptions): number[];
+//# sourceMappingURL=mean.d.ts.map

package/lib/compute/mean.d.ts.map

@@ -0,0 +1 @@
+{"version":3,"file":"mean.d.ts","sourceRoot":"","sources":["../../src/compute/mean.ts"],"names":[],"mappings":"AAAA,OAAO,KAAK,EAAE,KAAK,EAAE,MAAM,aAAa,CAAC;AACzC,OAAO,KAAK,EAAE,KAAK,EAAE,MAAM,sBAAsB,CAAC;AAElD,MAAM,WAAW,WAAW;IAC1B;;OAEG;IACH,MAAM,EAAE,KAAK,EAAE,CAAC;CACjB;AAED;;;;;;GAMG;AACH,wBAAgB,IAAI,CAAC,KAAK,EAAE,KAAK,EAAE,OAAO,CAAC,EAAE,WAAW,GAAG,MAAM,EAAE,CA8BlE"}

package/lib/compute/mean.js

@@ -0,0 +1,37 @@
+/**
+ * Compute the mean of an image. The mean can be either computed on each channel
+ * individually or on the whole image.
+ * @param image - Image to process.
+ * @param options - Mean options.
+ * @returns The mean pixel.
+ */
+export function mean(image, options) {
+    const pixelSum = new Array(image.channels).fill(0);
+    const nbValues = options ? options.points.length : image.size;
+    if (nbValues === 0)
+        throw new RangeError('Array of coordinates is empty.');
+    if (options) {
+        for (const point of options.points) {
+            for (let channel = 0; channel < image.channels; channel++) {
+                if (point.column < 0 ||
+                    point.column >= image.width ||
+                    point.row < 0 ||
+                    point.row >= image.height) {
+                    throw new RangeError(`Invalid coordinate: {column: ${point.column}, row: ${point.row}}.`);
+                }
+                pixelSum[channel] += image.getValueByPoint(point, channel);
+            }
+        }
+    }
+    else {
+        for (let row = 0; row < image.height; row++) {
+            for (let column = 0; column < image.width; column++) {
+                for (let channel = 0; channel < image.channels; channel++) {
+                    pixelSum[channel] += image.getValue(column, row, channel);
+                }
+            }
+        }
+    }
+    return pixelSum.map((channelSum) => channelSum / nbValues);
+}
+//# sourceMappingURL=mean.js.map

package/lib/compute/mean.js.map

@@ -0,0 +1 @@
+{"version":3,"file":"mean.js","sourceRoot":"","sources":["../../src/compute/mean.ts"],"names":[],"mappings":"AAUA;;;;;;GAMG;AACH,MAAM,UAAU,IAAI,CAAC,KAAY,EAAE,OAAqB;IACtD,MAAM,QAAQ,GAAG,IAAI,KAAK,CAAS,KAAK,CAAC,QAAQ,CAAC,CAAC,IAAI,CAAC,CAAC,CAAC,CAAC;IAC3D,MAAM,QAAQ,GAAG,OAAO,CAAC,CAAC,CAAC,OAAO,CAAC,MAAM,CAAC,MAAM,CAAC,CAAC,CAAC,KAAK,CAAC,IAAI,CAAC;IAC9D,IAAI,QAAQ,KAAK,CAAC;QAAE,MAAM,IAAI,UAAU,CAAC,gCAAgC,CAAC,CAAC;IAC3E,IAAI,OAAO,EAAE,CAAC;QACZ,KAAK,MAAM,KAAK,IAAI,OAAO,CAAC,MAAM,EAAE,CAAC;YACnC,KAAK,IAAI,OAAO,GAAG,CAAC,EAAE,OAAO,GAAG,KAAK,CAAC,QAAQ,EAAE,OAAO,EAAE,EAAE,CAAC;gBAC1D,IACE,KAAK,CAAC,MAAM,GAAG,CAAC;oBAChB,KAAK,CAAC,MAAM,IAAI,KAAK,CAAC,KAAK;oBAC3B,KAAK,CAAC,GAAG,GAAG,CAAC;oBACb,KAAK,CAAC,GAAG,IAAI,KAAK,CAAC,MAAM,EACzB,CAAC;oBACD,MAAM,IAAI,UAAU,CAClB,gCAAgC,KAAK,CAAC,MAAM,UAAU,KAAK,CAAC,GAAG,IAAI,CACpE,CAAC;gBACJ,CAAC;gBACD,QAAQ,CAAC,OAAO,CAAC,IAAI,KAAK,CAAC,eAAe,CAAC,KAAK,EAAE,OAAO,CAAC,CAAC;YAC7D,CAAC;QACH,CAAC;IACH,CAAC;SAAM,CAAC;QACN,KAAK,IAAI,GAAG,GAAG,CAAC,EAAE,GAAG,GAAG,KAAK,CAAC,MAAM,EAAE,GAAG,EAAE,EAAE,CAAC;YAC5C,KAAK,IAAI,MAAM,GAAG,CAAC,EAAE,MAAM,GAAG,KAAK,CAAC,KAAK,EAAE,MAAM,EAAE,EAAE,CAAC;gBACpD,KAAK,IAAI,OAAO,GAAG,CAAC,EAAE,OAAO,GAAG,KAAK,CAAC,QAAQ,EAAE,OAAO,EAAE,EAAE,CAAC;oBAC1D,QAAQ,CAAC,OAAO,CAAC,IAAI,KAAK,CAAC,QAAQ,CAAC,MAAM,EAAE,GAAG,EAAE,OAAO,CAAC,CAAC;gBAC5D,CAAC;YACH,CAAC;QACH,CAAC;IACH,CAAC;IACD,OAAO,QAAQ,CAAC,GAAG,CAAC,CAAC,UAAU,EAAE,EAAE,CAAC,UAAU,GAAG,QAAQ,CAAC,CAAC;AAC7D,CAAC"}

package/lib/compute/median.d.ts

@@ -0,0 +1,16 @@
+import type { Image } from '../Image.js';
+import type { Point } from '../geometry/index.js';
+export interface MedianOptions {
+    /**
+     * Points to calculate median from.
+     */
+    points: Point[];
+}
+/**
+ * Returns the median pixel of the image. The median is computed on each channel individually.
+ * @param image - Image to process.
+ * @param options - Median options.
+ * @returns Median pixel.
+ */
+export declare function median(image: Image, options?: MedianOptions): number[];
+//# sourceMappingURL=median.d.ts.map

package/lib/compute/median.d.ts.map

@@ -0,0 +1 @@
+{"version":3,"file":"median.d.ts","sourceRoot":"","sources":["../../src/compute/median.ts"],"names":[],"mappings":"AAGA,OAAO,KAAK,EAAE,KAAK,EAAE,MAAM,aAAa,CAAC;AACzC,OAAO,KAAK,EAAE,KAAK,EAAE,MAAM,sBAAsB,CAAC;AAElD,MAAM,WAAW,aAAa;IAC5B;;OAEG;IACH,MAAM,EAAE,KAAK,EAAE,CAAC;CACjB;AACD;;;;;GAKG;AACH,wBAAgB,MAAM,CAAC,KAAK,EAAE,KAAK,EAAE,OAAO,CAAC,EAAE,aAAa,GAAG,MAAM,EAAE,CA+BtE"}

package/lib/compute/median.js

@@ -0,0 +1,37 @@
+// @ts-expect-error: median-quisckselect has no types
+import quickMedian from 'median-quickselect';
+/**
+ * Returns the median pixel of the image. The median is computed on each channel individually.
+ * @param image - Image to process.
+ * @param options - Median options.
+ * @returns Median pixel.
+ */
+export function median(image, options) {
+    const pixel = new Array(image.channels).fill(0);
+    if (options) {
+        if (options.points.length === 0) {
+            throw new RangeError('Array of coordinates is empty.');
+        }
+        for (let i = 0; i < image.channels; i++) {
+            const channel = [];
+            for (const point of options.points) {
+                if (point.column < 0 ||
+                    point.column >= image.width ||
+                    point.row < 0 ||
+                    point.row >= image.height) {
+                    throw new RangeError(`Invalid coordinate: {column: ${point.column}, row: ${point.row}}.`);
+                }
+                channel.push(image.getValueByPoint(point, i));
+            }
+            pixel[i] = quickMedian(channel);
+        }
+    }
+    else {
+        for (let i = 0; i < image.channels; i++) {
+            const channel = image.getChannel(i);
+            pixel[i] = quickMedian(channel);
+        }
+    }
+    return pixel;
+}
+//# sourceMappingURL=median.js.map

package/lib/compute/median.js.map

@@ -0,0 +1 @@
+{"version":3,"file":"median.js","sourceRoot":"","sources":["../../src/compute/median.ts"],"names":[],"mappings":"AAAA,qDAAqD;AACrD,OAAO,WAAW,MAAM,oBAAoB,CAAC;AAW7C;;;;;GAKG;AACH,MAAM,UAAU,MAAM,CAAC,KAAY,EAAE,OAAuB;IAC1D,MAAM,KAAK,GAAG,IAAI,KAAK,CAAS,KAAK,CAAC,QAAQ,CAAC,CAAC,IAAI,CAAC,CAAC,CAAC,CAAC;IACxD,IAAI,OAAO,EAAE,CAAC;QACZ,IAAI,OAAO,CAAC,MAAM,CAAC,MAAM,KAAK,CAAC,EAAE,CAAC;YAChC,MAAM,IAAI,UAAU,CAAC,gCAAgC,CAAC,CAAC;QACzD,CAAC;QACD,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,KAAK,CAAC,QAAQ,EAAE,CAAC,EAAE,EAAE,CAAC;YACxC,MAAM,OAAO,GAAa,EAAE,CAAC;YAC7B,KAAK,MAAM,KAAK,IAAI,OAAO,CAAC,MAAM,EAAE,CAAC;gBACnC,IACE,KAAK,CAAC,MAAM,GAAG,CAAC;oBAChB,KAAK,CAAC,MAAM,IAAI,KAAK,CAAC,KAAK;oBAC3B,KAAK,CAAC,GAAG,GAAG,CAAC;oBACb,KAAK,CAAC,GAAG,IAAI,KAAK,CAAC,MAAM,EACzB,CAAC;oBACD,MAAM,IAAI,UAAU,CAClB,gCAAgC,KAAK,CAAC,MAAM,UAAU,KAAK,CAAC,GAAG,IAAI,CACpE,CAAC;gBACJ,CAAC;gBACD,OAAO,CAAC,IAAI,CAAC,KAAK,CAAC,eAAe,CAAC,KAAK,EAAE,CAAC,CAAC,CAAC,CAAC;YAChD,CAAC;YACD,KAAK,CAAC,CAAC,CAAC,GAAG,WAAW,CAAC,OAAO,CAAC,CAAC;QAClC,CAAC;IACH,CAAC;SAAM,CAAC;QACN,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,KAAK,CAAC,QAAQ,EAAE,CAAC,EAAE,EAAE,CAAC;YACxC,MAAM,OAAO,GAAG,KAAK,CAAC,UAAU,CAAC,CAAC,CAAC,CAAC;YACpC,KAAK,CAAC,CAAC,CAAC,GAAG,WAAW,CAAC,OAAO,CAAC,CAAC;QAClC,CAAC;IACH,CAAC;IAED,OAAO,KAAK,CAAC;AACf,CAAC"}

package/lib/compute/variance.d.ts

@@ -0,0 +1,17 @@
+import type { Image } from '../Image.js';
+import type { Point } from '../geometry/index.js';
+export interface VarianceOptions {
+    /**
+     * Points to calculate variance from.
+     */
+    points: Point[];
+}
+/**
+ * Compute the variance of each channel of an image.
+ * @see {@link https://en.wikipedia.org/wiki/Variance}
+ * @param image - Image to process.
+ * @param options - Variance options.
+ * @returns The variance of the channels of the image.
+ */
+export declare function variance(image: Image, options?: VarianceOptions): number[];
+//# sourceMappingURL=variance.d.ts.map

package/lib/compute/variance.d.ts.map

@@ -0,0 +1 @@
+{"version":3,"file":"variance.d.ts","sourceRoot":"","sources":["../../src/compute/variance.ts"],"names":[],"mappings":"AAAA,OAAO,KAAK,EAAE,KAAK,EAAE,MAAM,aAAa,CAAC;AACzC,OAAO,KAAK,EAAE,KAAK,EAAE,MAAM,sBAAsB,CAAC;AAElD,MAAM,WAAW,eAAe;IAC9B;;OAEG;IACH,MAAM,EAAE,KAAK,EAAE,CAAC;CACjB;AAED;;;;;;GAMG;AACH,wBAAgB,QAAQ,CAAC,KAAK,EAAE,KAAK,EAAE,OAAO,CAAC,EAAE,eAAe,GAAG,MAAM,EAAE,CAuB1E"}

package/lib/compute/variance.js

@@ -0,0 +1,33 @@
+/**
+ * Compute the variance of each channel of an image.
+ * @see {@link https://en.wikipedia.org/wiki/Variance}
+ * @param image - Image to process.
+ * @param options - Variance options.
+ * @returns The variance of the channels of the image.
+ */
+export function variance(image, options) {
+    const mean = image.mean(options);
+    const sum = new Array(image.channels).fill(0);
+    if (options) {
+        for (const point of options.points) {
+            for (let channel = 0; channel < image.channels; channel++) {
+                sum[channel] +=
+                    (image.getValue(point.column, point.row, channel) - mean[channel]) **
+                        2;
+            }
+        }
+    }
+    else {
+        for (let i = 0; i < image.size; i++) {
+            for (let channel = 0; channel < image.channels; channel++) {
+                sum[channel] +=
+                    (image.getValueByIndex(i, channel) - mean[channel]) ** 2;
+            }
+        }
+    }
+    const nbValues = options
+        ? options.points.length * image.channels
+        : image.size * image.channels;
+    return sum.map((channel) => channel / nbValues);
+}
+//# sourceMappingURL=variance.js.map

package/lib/compute/variance.js.map

@@ -0,0 +1 @@
+{"version":3,"file":"variance.js","sourceRoot":"","sources":["../../src/compute/variance.ts"],"names":[],"mappings":"AAUA;;;;;;GAMG;AACH,MAAM,UAAU,QAAQ,CAAC,KAAY,EAAE,OAAyB;IAC9D,MAAM,IAAI,GAAG,KAAK,CAAC,IAAI,CAAC,OAAO,CAAC,CAAC;IACjC,MAAM,GAAG,GAAG,IAAI,KAAK,CAAS,KAAK,CAAC,QAAQ,CAAC,CAAC,IAAI,CAAC,CAAC,CAAC,CAAC;IACtD,IAAI,OAAO,EAAE,CAAC;QACZ,KAAK,MAAM,KAAK,IAAI,OAAO,CAAC,MAAM,EAAE,CAAC;YACnC,KAAK,IAAI,OAAO,GAAG,CAAC,EAAE,OAAO,GAAG,KAAK,CAAC,QAAQ,EAAE,OAAO,EAAE,EAAE,CAAC;gBAC1D,GAAG,CAAC,OAAO,CAAC;oBACV,CAAC,KAAK,CAAC,QAAQ,CAAC,KAAK,CAAC,MAAM,EAAE,KAAK,CAAC,GAAG,EAAE,OAAO,CAAC,GAAG,IAAI,CAAC,OAAO,CAAC,CAAC;wBAClE,CAAC,CAAC;YACN,CAAC;QACH,CAAC;IACH,CAAC;SAAM,CAAC;QACN,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,KAAK,CAAC,IAAI,EAAE,CAAC,EAAE,EAAE,CAAC;YACpC,KAAK,IAAI,OAAO,GAAG,CAAC,EAAE,OAAO,GAAG,KAAK,CAAC,QAAQ,EAAE,OAAO,EAAE,EAAE,CAAC;gBAC1D,GAAG,CAAC,OAAO,CAAC;oBACV,CAAC,KAAK,CAAC,eAAe,CAAC,CAAC,EAAE,OAAO,CAAC,GAAG,IAAI,CAAC,OAAO,CAAC,CAAC,IAAI,CAAC,CAAC;YAC7D,CAAC;QACH,CAAC;IACH,CAAC;IACD,MAAM,QAAQ,GAAG,OAAO;QACtB,CAAC,CAAC,OAAO,CAAC,MAAM,CAAC,MAAM,GAAG,KAAK,CAAC,QAAQ;QACxC,CAAC,CAAC,KAAK,CAAC,IAAI,GAAG,KAAK,CAAC,QAAQ,CAAC;IAChC,OAAO,GAAG,CAAC,GAAG,CAAC,CAAC,OAAO,EAAE,EAAE,CAAC,OAAO,GAAG,QAAQ,CAAC,CAAC;AAClD,CAAC"}

package/lib/correctColor/correctColor.d.ts

@@ -0,0 +1,22 @@
+import type { RgbColor } from 'colord';
+import { Image } from '../Image.js';
+/**
+ * Correct the colors in an image using the reference colors.
+ * Algorithm is based on the paper "Color correction using improved linear regression algorithm".
+ * DOI: 10.1109/ICTS.2015.7379874.
+ * @param image - Image to process.
+ * @param measuredColors - Colors from the image, which will be compared to the reference.
+ * @param referenceColors - Reference colors.
+ * @returns Image with the colors corrected.
+ */
+export declare function correctColor(image: Image, measuredColors: RgbColor[], referenceColors: RgbColor[]): Image;
+/**
+ * Compute the third order variables for the regression from an RGB color.
+ * @param r - Red component.
+ * @param g - Green component.
+ * @param b - Blue component.
+ * @param maxValue - Maximal acceptable value for the image to process.
+ * @returns The variables for the multivariate linear regression.
+ */
+export declare function getRegressionVariables(r: number, g: number, b: number, maxValue: number): number[];
+//# sourceMappingURL=correctColor.d.ts.map

package/lib/correctColor/correctColor.d.ts.map

@@ -0,0 +1 @@
+{"version":3,"file":"correctColor.d.ts","sourceRoot":"","sources":["../../src/correctColor/correctColor.ts"],"names":[],"mappings":"AAAA,OAAO,KAAK,EAAE,QAAQ,EAAE,MAAM,QAAQ,CAAC;AAGvC,OAAO,EAAE,KAAK,EAAE,MAAM,aAAa,CAAC;AASpC;;;;;;;;GAQG;AACH,wBAAgB,YAAY,CAC1B,KAAK,EAAE,KAAK,EACZ,cAAc,EAAE,QAAQ,EAAE,EAC1B,eAAe,EAAE,QAAQ,EAAE,GAC1B,KAAK,CAkDP;AAED;;;;;;;GAOG;AACH,wBAAgB,sBAAsB,CACpC,CAAC,EAAE,MAAM,EACT,CAAC,EAAE,MAAM,EACT,CAAC,EAAE,MAAM,EACT,QAAQ,EAAE,MAAM,GACf,MAAM,EAAE,CAkBV"}

package/lib/correctColor/correctColor.js

@@ -0,0 +1,75 @@
+import MLR from 'ml-regression-multivariate-linear';
+import { Image } from '../Image.js';
+import { getClamp } from '../utils/clamp.js';
+import checkProcessable from '../utils/validators/checkProcessable.js';
+import { formatInputForMlr, formatReferenceForMlr, } from './utils/formatData.js';
+/**
+ * Correct the colors in an image using the reference colors.
+ * Algorithm is based on the paper "Color correction using improved linear regression algorithm".
+ * DOI: 10.1109/ICTS.2015.7379874.
+ * @param image - Image to process.
+ * @param measuredColors - Colors from the image, which will be compared to the reference.
+ * @param referenceColors - Reference colors.
+ * @returns Image with the colors corrected.
+ */
+export function correctColor(image, measuredColors, referenceColors) {
+    checkProcessable(image, {
+        colorModel: ['RGB', 'RGBA'],
+    });
+    if (measuredColors.length !== referenceColors.length) {
+        throw new RangeError('number of measured colors and reference colors must be the same');
+    }
+    const inputData = formatInputForMlr(measuredColors, image.maxValue);
+    const referenceData = formatReferenceForMlr(referenceColors, image.maxValue);
+    const mlrRed = new MLR(inputData, referenceData.r);
+    const mlrGreen = new MLR(inputData, referenceData.g);
+    const mlrBlue = new MLR(inputData, referenceData.b);
+    const result = Image.createFrom(image);
+    for (let row = 0; row < image.height; row++) {
+        for (let column = 0; column < image.width; column++) {
+            const pixel = image.getPixel(column, row);
+            const variables = getRegressionVariables(pixel[0], pixel[1], pixel[2], image.maxValue);
+            const clamp = getClamp(image);
+            const newPixel = [0, 0, 0];
+            const red = mlrRed.predict(variables)[0] * image.maxValue;
+            const green = mlrGreen.predict(variables)[0] * image.maxValue;
+            const blue = mlrBlue.predict(variables)[0] * image.maxValue;
+            newPixel[0] = clamp(red);
+            newPixel[1] = clamp(green);
+            newPixel[2] = clamp(blue);
+            if (image.alpha) {
+                newPixel[3] = image.getValue(column, row, 3);
+            }
+            result.setPixel(column, row, newPixel);
+        }
+    }
+    return result;
+}
+/**
+ * Compute the third order variables for the regression from an RGB color.
+ * @param r - Red component.
+ * @param g - Green component.
+ * @param b - Blue component.
+ * @param maxValue - Maximal acceptable value for the image to process.
+ * @returns The variables for the multivariate linear regression.
+ */
+export function getRegressionVariables(r, g, b, maxValue) {
+    r /= maxValue;
+    g /= maxValue;
+    b /= maxValue;
+    return [
+        r,
+        g,
+        b,
+        r ** 2,
+        g ** 2,
+        b ** 2,
+        r ** 3,
+        g ** 3,
+        b ** 3,
+        r * g,
+        r * b,
+        b * g,
+    ];
+}
+//# sourceMappingURL=correctColor.js.map

package/lib/correctColor/correctColor.js.map

@@ -0,0 +1 @@
+{"version":3,"file":"correctColor.js","sourceRoot":"","sources":["../../src/correctColor/correctColor.ts"],"names":[],"mappings":"AACA,OAAO,GAAG,MAAM,mCAAmC,CAAC;AAEpD,OAAO,EAAE,KAAK,EAAE,MAAM,aAAa,CAAC;AACpC,OAAO,EAAE,QAAQ,EAAE,MAAM,mBAAmB,CAAC;AAC7C,OAAO,gBAAgB,MAAM,yCAAyC,CAAC;AAEvE,OAAO,EACL,iBAAiB,EACjB,qBAAqB,GACtB,MAAM,uBAAuB,CAAC;AAE/B;;;;;;;;GAQG;AACH,MAAM,UAAU,YAAY,CAC1B,KAAY,EACZ,cAA0B,EAC1B,eAA2B;IAE3B,gBAAgB,CAAC,KAAK,EAAE;QACtB,UAAU,EAAE,CAAC,KAAK,EAAE,MAAM,CAAC;KAC5B,CAAC,CAAC;IAEH,IAAI,cAAc,CAAC,MAAM,KAAK,eAAe,CAAC,MAAM,EAAE,CAAC;QACrD,MAAM,IAAI,UAAU,CAClB,iEAAiE,CAClE,CAAC;IACJ,CAAC;IAED,MAAM,SAAS,GAAG,iBAAiB,CAAC,cAAc,EAAE,KAAK,CAAC,QAAQ,CAAC,CAAC;IACpE,MAAM,aAAa,GAAG,qBAAqB,CAAC,eAAe,EAAE,KAAK,CAAC,QAAQ,CAAC,CAAC;IAE7E,MAAM,MAAM,GAAG,IAAI,GAAG,CAAC,SAAS,EAAE,aAAa,CAAC,CAAC,CAAC,CAAC;IACnD,MAAM,QAAQ,GAAG,IAAI,GAAG,CAAC,SAAS,EAAE,aAAa,CAAC,CAAC,CAAC,CAAC;IACrD,MAAM,OAAO,GAAG,IAAI,GAAG,CAAC,SAAS,EAAE,aAAa,CAAC,CAAC,CAAC,CAAC;IAEpD,MAAM,MAAM,GAAG,KAAK,CAAC,UAAU,CAAC,KAAK,CAAC,CAAC;IAEvC,KAAK,IAAI,GAAG,GAAG,CAAC,EAAE,GAAG,GAAG,KAAK,CAAC,MAAM,EAAE,GAAG,EAAE,EAAE,CAAC;QAC5C,KAAK,IAAI,MAAM,GAAG,CAAC,EAAE,MAAM,GAAG,KAAK,CAAC,KAAK,EAAE,MAAM,EAAE,EAAE,CAAC;YACpD,MAAM,KAAK,GAAG,KAAK,CAAC,QAAQ,CAAC,MAAM,EAAE,GAAG,CAAC,CAAC;YAC1C,MAAM,SAAS,GAAG,sBAAsB,CACtC,KAAK,CAAC,CAAC,CAAC,EACR,KAAK,CAAC,CAAC,CAAC,EACR,KAAK,CAAC,CAAC,CAAC,EACR,KAAK,CAAC,QAAQ,CACf,CAAC;YAEF,MAAM,KAAK,GAAG,QAAQ,CAAC,KAAK,CAAC,CAAC;YAE9B,MAAM,QAAQ,GAAG,CAAC,CAAC,EAAE,CAAC,EAAE,CAAC,CAAC,CAAC;YAE3B,MAAM,GAAG,GAAG,MAAM,CAAC,OAAO,CAAC,SAAS,CAAC,CAAC,CAAC,CAAC,GAAG,KAAK,CAAC,QAAQ,CAAC;YAC1D,MAAM,KAAK,GAAG,QAAQ,CAAC,OAAO,CAAC,SAAS,CAAC,CAAC,CAAC,CAAC,GAAG,KAAK,CAAC,QAAQ,CAAC;YAC9D,MAAM,IAAI,GAAG,OAAO,CAAC,OAAO,CAAC,SAAS,CAAC,CAAC,CAAC,CAAC,GAAG,KAAK,CAAC,QAAQ,CAAC;YAE5D,QAAQ,CAAC,CAAC,CAAC,GAAG,KAAK,CAAC,GAAG,CAAC,CAAC;YACzB,QAAQ,CAAC,CAAC,CAAC,GAAG,KAAK,CAAC,KAAK,CAAC,CAAC;YAC3B,QAAQ,CAAC,CAAC,CAAC,GAAG,KAAK,CAAC,IAAI,CAAC,CAAC;YAC1B,IAAI,KAAK,CAAC,KAAK,EAAE,CAAC;gBAChB,QAAQ,CAAC,CAAC,CAAC,GAAG,KAAK,CAAC,QAAQ,CAAC,MAAM,EAAE,GAAG,EAAE,CAAC,CAAC,CAAC;YAC/C,CAAC;YAED,MAAM,CAAC,QAAQ,CAAC,MAAM,EAAE,GAAG,EAAE,QAAQ,CAAC,CAAC;QACzC,CAAC;IACH,CAAC;IAED,OAAO,MAAM,CAAC;AAChB,CAAC;AAED;;;;;;;GAOG;AACH,MAAM,UAAU,sBAAsB,CACpC,CAAS,EACT,CAAS,EACT,CAAS,EACT,QAAgB;IAEhB,CAAC,IAAI,QAAQ,CAAC;IACd,CAAC,IAAI,QAAQ,CAAC;IACd,CAAC,IAAI,QAAQ,CAAC;IACd,OAAO;QACL,CAAC;QACD,CAAC;QACD,CAAC;QACD,CAAC,IAAI,CAAC;QACN,CAAC,IAAI,CAAC;QACN,CAAC,IAAI,CAAC;QACN,CAAC,IAAI,CAAC;QACN,CAAC,IAAI,CAAC;QACN,CAAC,IAAI,CAAC;QACN,CAAC,GAAG,CAAC;QACL,CAAC,GAAG,CAAC;QACL,CAAC,GAAG,CAAC;KACN,CAAC;AACJ,CAAC"}

package/lib/correctColor/index.d.ts

@@ -0,0 +1,2 @@
+export * from './correctColor.js';
+//# sourceMappingURL=index.d.ts.map

package/lib/correctColor/index.d.ts.map

@@ -0,0 +1 @@
+{"version":3,"file":"index.d.ts","sourceRoot":"","sources":["../../src/correctColor/index.ts"],"names":[],"mappings":"AAAA,cAAc,mBAAmB,CAAC"}

package/lib/correctColor/index.js

@@ -0,0 +1,2 @@
+export * from './correctColor.js';
+//# sourceMappingURL=index.js.map

package/lib/correctColor/index.js.map

@@ -0,0 +1 @@
+{"version":3,"file":"index.js","sourceRoot":"","sources":["../../src/correctColor/index.ts"],"names":[],"mappings":"AAAA,cAAc,mBAAmB,CAAC"}

package/lib/correctColor/utils/formatData.d.ts

@@ -0,0 +1,37 @@
+import type { RgbColor } from 'colord';
+import type { ColorCard } from './referenceColorCard.ts';
+/**
+ * Reference QP card data formatted to be easier to use for the multivariate linear regression.
+ */
+export interface ReferenceDataForMlr {
+    r: number[][];
+    g: number[][];
+    b: number[][];
+}
+/**
+ * Convert RGB array colors to RGB object colors. Used to get the properly formatted measured colors.
+ * @param arrayColors - Array of RGB colors as 3 elements array.
+ * @returns Array of RGB objects.
+ */
+export declare function getMeasuredColors(arrayColors: number[][]): RgbColor[];
+/**
+ * Extract the colors from a QP card and convert them to RGB.
+ * @param qpCard - QP card containing the color reference values in L*a*b*.
+ * @returns Array of reference RGB colors.
+ */
+export declare function getReferenceColors(qpCard: ColorCard): RgbColor[];
+/**
+ * Format and normalize data from a QP card to use as a reference in a multivariate linear regression.
+ * @param referenceColors - Array of RGB colors used as a reference.
+ * @param maxValue - Maximal acceptable value for the image to process.
+ * @returns The formatted data.
+ */
+export declare function formatReferenceForMlr(referenceColors: RgbColor[], maxValue: number): ReferenceDataForMlr;
+/**
+ * Compute the variables for the multivariate linear regression based on the the input colors. Values are normalized between 0 and 1.
+ * @param inputColors - The input colors as an array of rgb objects.
+ * @param maxValue - Maximal acceptable value for the image to process.
+ * @returns The formatted input data for the regression.
+ */
+export declare function formatInputForMlr(inputColors: RgbColor[], maxValue: number): number[][];
+//# sourceMappingURL=formatData.d.ts.map

package/lib/correctColor/utils/formatData.d.ts.map

@@ -0,0 +1 @@
+{"version":3,"file":"formatData.d.ts","sourceRoot":"","sources":["../../../src/correctColor/utils/formatData.ts"],"names":[],"mappings":"AAAA,OAAO,KAAK,EAAE,QAAQ,EAAE,MAAM,QAAQ,CAAC;AAMvC,OAAO,KAAK,EAAE,SAAS,EAAE,MAAM,yBAAyB,CAAC;AAOzD;;GAEG;AACH,MAAM,WAAW,mBAAmB;IAClC,CAAC,EAAE,MAAM,EAAE,EAAE,CAAC;IACd,CAAC,EAAE,MAAM,EAAE,EAAE,CAAC;IACd,CAAC,EAAE,MAAM,EAAE,EAAE,CAAC;CACf;AAED;;;;GAIG;AACH,wBAAgB,iBAAiB,CAAC,WAAW,EAAE,MAAM,EAAE,EAAE,GAAG,QAAQ,EAAE,CAMrE;AAED;;;;GAIG;AACH,wBAAgB,kBAAkB,CAAC,MAAM,EAAE,SAAS,GAAG,QAAQ,EAAE,CAOhE;AAED;;;;;GAKG;AACH,wBAAgB,qBAAqB,CACnC,eAAe,EAAE,QAAQ,EAAE,EAC3B,QAAQ,EAAE,MAAM,GACf,mBAAmB,CAUrB;AAED;;;;;GAKG;AACH,wBAAgB,iBAAiB,CAC/B,WAAW,EAAE,QAAQ,EAAE,EACvB,QAAQ,EAAE,MAAM,GACf,MAAM,EAAE,EAAE,CAMZ"}