image-js 0.0.0-next-e30968a2d8

package/lib-esm/operations/thresholds/intermodes.js

@@ -0,0 +1,59 @@
+/*
+ *
+ * see https://github.com/fiji/Auto_Threshold/blob/master/src/main/java/fiji/threshold/Auto_Threshold.java
+ * Intermodes: This assumes a bimodal histogram. Implements the thresholding Prewitt, JMS & Mendelsohn, ML (1966),
+ * "The analysis of cell images", Annals of the NewYork Academy of Sciences 128: 1035-1053
+ *
+ */
+/**
+ * Return a threshold for a histogram using Intermodes algorithm.
+ *
+ * @param histogram - Image histogram.
+ * @returns The threshold.
+ */
+export default function intermodes(histogram) {
+    const iHisto = histogram.slice();
+    let iter = 0;
+    while (!bimodalTest(iHisto)) {
+        // smooth with a 3 point running mean filter
+        let previous = 0;
+        let current = 0;
+        let next = iHisto[0];
+        for (let i = 0; i < histogram.length - 1; i++) {
+            previous = current;
+            current = next;
+            next = iHisto[i + 1];
+            iHisto[i] = (previous + current + next) / 3;
+        }
+        iHisto[histogram.length - 1] = (current + next) / 3;
+        iter++;
+        if (iter > 10000) {
+            throw new Error('Intermodes Threshold not found after 10000 iterations');
+        }
+    }
+    // The threshold is the mean between the two peaks.
+    let tt = 0;
+    for (let i = 1; i < histogram.length - 1; i++) {
+        if (iHisto[i - 1] < iHisto[i] && iHisto[i + 1] < iHisto[i]) {
+            tt += i;
+        }
+    }
+    return Math.floor(tt / 2.0);
+}
+function bimodalTest(iHisto) {
+    let b = false;
+    let modes = 0;
+    for (let k = 1; k < iHisto.length - 1; k++) {
+        if (iHisto[k - 1] < iHisto[k] && iHisto[k + 1] < iHisto[k]) {
+            modes++;
+            if (modes > 2) {
+                return false;
+            }
+        }
+    }
+    if (modes === 2) {
+        b = true;
+    }
+    return b;
+}
+//# sourceMappingURL=intermodes.js.map

package/lib-esm/operations/thresholds/intermodes.js.map

@@ -0,0 +1 @@
+{"version":3,"file":"intermodes.js","sourceRoot":"","sources":["../../../src/operations/thresholds/intermodes.ts"],"names":[],"mappings":"AAAA;;;;;;GAMG;AAEH;;;;;GAKG;AACH,MAAM,CAAC,OAAO,UAAU,UAAU,CAAC,SAAsB;IACvD,MAAM,MAAM,GAAG,SAAS,CAAC,KAAK,EAAE,CAAC;IACjC,IAAI,IAAI,GAAG,CAAC,CAAC;IACb,OAAO,CAAC,WAAW,CAAC,MAAM,CAAC,EAAE;QAC3B,4CAA4C;QAC5C,IAAI,QAAQ,GAAG,CAAC,CAAC;QACjB,IAAI,OAAO,GAAG,CAAC,CAAC;QAChB,IAAI,IAAI,GAAG,MAAM,CAAC,CAAC,CAAC,CAAC;QACrB,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,SAAS,CAAC,MAAM,GAAG,CAAC,EAAE,CAAC,EAAE,EAAE;YAC7C,QAAQ,GAAG,OAAO,CAAC;YACnB,OAAO,GAAG,IAAI,CAAC;YACf,IAAI,GAAG,MAAM,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC;YACrB,MAAM,CAAC,CAAC,CAAC,GAAG,CAAC,QAAQ,GAAG,OAAO,GAAG,IAAI,CAAC,GAAG,CAAC,CAAC;SAC7C;QACD,MAAM,CAAC,SAAS,CAAC,MAAM,GAAG,CAAC,CAAC,GAAG,CAAC,OAAO,GAAG,IAAI,CAAC,GAAG,CAAC,CAAC;QACpD,IAAI,EAAE,CAAC;QACP,IAAI,IAAI,GAAG,KAAK,EAAE;YAChB,MAAM,IAAI,KAAK,CAAC,uDAAuD,CAAC,CAAC;SAC1E;KACF;IAED,mDAAmD;IACnD,IAAI,EAAE,GAAG,CAAC,CAAC;IACX,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,SAAS,CAAC,MAAM,GAAG,CAAC,EAAE,CAAC,EAAE,EAAE;QAC7C,IAAI,MAAM,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,MAAM,CAAC,CAAC,CAAC,IAAI,MAAM,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,MAAM,CAAC,CAAC,CAAC,EAAE;YAC1D,EAAE,IAAI,CAAC,CAAC;SACT;KACF;IACD,OAAO,IAAI,CAAC,KAAK,CAAC,EAAE,GAAG,GAAG,CAAC,CAAC;AAC9B,CAAC;AAED,SAAS,WAAW,CAAC,MAAmB;IACtC,IAAI,CAAC,GAAG,KAAK,CAAC;IACd,IAAI,KAAK,GAAG,CAAC,CAAC;IAEd,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,MAAM,CAAC,MAAM,GAAG,CAAC,EAAE,CAAC,EAAE,EAAE;QAC1C,IAAI,MAAM,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,MAAM,CAAC,CAAC,CAAC,IAAI,MAAM,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,MAAM,CAAC,CAAC,CAAC,EAAE;YAC1D,KAAK,EAAE,CAAC;YACR,IAAI,KAAK,GAAG,CAAC,EAAE;gBACb,OAAO,KAAK,CAAC;aACd;SACF;KACF;IACD,IAAI,KAAK,KAAK,CAAC,EAAE;QACf,CAAC,GAAG,IAAI,CAAC;KACV;IACD,OAAO,CAAC,CAAC;AACX,CAAC"}

package/lib-esm/operations/thresholds/isodata.d.ts

@@ -0,0 +1,8 @@
+/**
+ * Return a threshold for a histogram using Isodata algorithm.
+ *
+ * @param histogram - Image histogram.
+ * @returns The threshold.
+ */
+export default function isodata(histogram: Uint32Array): number;
+//# sourceMappingURL=isodata.d.ts.map

package/lib-esm/operations/thresholds/isodata.d.ts.map

@@ -0,0 +1 @@
+{"version":3,"file":"isodata.d.ts","sourceRoot":"","sources":["../../../src/operations/thresholds/isodata.ts"],"names":[],"mappings":"AAMA;;;;;GAKG;AACH,MAAM,CAAC,OAAO,UAAU,OAAO,CAAC,SAAS,EAAE,WAAW,GAAG,MAAM,CAwC9D"}

package/lib-esm/operations/thresholds/isodata.js

@@ -0,0 +1,52 @@
+/*
+ * see https://github.com/fiji/Auto_Threshold/blob/master/src/main/java/fiji/threshold/Auto_Threshold.java
+ * Isodata: Ridler, TW & Calvard, S (1978), "Picture thresholding using an iterative selection method"
+ * IEEE Transactions on Systems, Man and Cybernetics 8: 630-632.
+ *
+ */
+/**
+ * Return a threshold for a histogram using Isodata algorithm.
+ *
+ * @param histogram - Image histogram.
+ * @returns The threshold.
+ */
+export default function isodata(histogram) {
+    let l; // the average grey value of pixels with intensities < g
+    let toth; // the the average grey value of pixels with intensities > g
+    let totl; // the total the average grey value of pixels with intensities < g
+    let h; // the average grey value of pixels with intensities > g
+    let g = 0; // threshold value
+    for (let i = 1; i < histogram.length; i++) {
+        if (histogram[i] > 0) {
+            g = i + 1;
+            break;
+        }
+    }
+    while (true) {
+        l = 0;
+        totl = 0;
+        for (let i = 0; i < g; i++) {
+            totl = totl + histogram[i];
+            l = l + histogram[i] * i;
+        }
+        h = 0;
+        toth = 0;
+        for (let i = g + 1; i < histogram.length; i++) {
+            toth += histogram[i];
+            h += histogram[i] * i;
+        }
+        if (totl > 0 && toth > 0) {
+            l /= totl;
+            h /= toth;
+            if (g === Math.round((l + h) / 2.0)) {
+                break;
+            }
+        }
+        g++;
+        if (g > histogram.length - 2) {
+            throw new Error('Threshold not found');
+        }
+    }
+    return g;
+}
+//# sourceMappingURL=isodata.js.map

package/lib-esm/operations/thresholds/isodata.js.map

@@ -0,0 +1 @@
+{"version":3,"file":"isodata.js","sourceRoot":"","sources":["../../../src/operations/thresholds/isodata.ts"],"names":[],"mappings":"AAAA;;;;;GAKG;AACH;;;;;GAKG;AACH,MAAM,CAAC,OAAO,UAAU,OAAO,CAAC,SAAsB;IACpD,IAAI,CAAC,CAAC,CAAC,wDAAwD;IAC/D,IAAI,IAAI,CAAC,CAAC,4DAA4D;IACtE,IAAI,IAAI,CAAC,CAAC,kEAAkE;IAC5E,IAAI,CAAC,CAAC,CAAC,wDAAwD;IAC/D,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC,kBAAkB;IAE7B,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,SAAS,CAAC,MAAM,EAAE,CAAC,EAAE,EAAE;QACzC,IAAI,SAAS,CAAC,CAAC,CAAC,GAAG,CAAC,EAAE;YACpB,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC;YACV,MAAM;SACP;KACF;IAED,OAAO,IAAI,EAAE;QACX,CAAC,GAAG,CAAC,CAAC;QACN,IAAI,GAAG,CAAC,CAAC;QACT,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,CAAC,EAAE,CAAC,EAAE,EAAE;YAC1B,IAAI,GAAG,IAAI,GAAG,SAAS,CAAC,CAAC,CAAC,CAAC;YAC3B,CAAC,GAAG,CAAC,GAAG,SAAS,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC;SAC1B;QACD,CAAC,GAAG,CAAC,CAAC;QACN,IAAI,GAAG,CAAC,CAAC;QACT,KAAK,IAAI,CAAC,GAAG,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,SAAS,CAAC,MAAM,EAAE,CAAC,EAAE,EAAE;YAC7C,IAAI,IAAI,SAAS,CAAC,CAAC,CAAC,CAAC;YACrB,CAAC,IAAI,SAAS,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC;SACvB;QACD,IAAI,IAAI,GAAG,CAAC,IAAI,IAAI,GAAG,CAAC,EAAE;YACxB,CAAC,IAAI,IAAI,CAAC;YACV,CAAC,IAAI,IAAI,CAAC;YACV,IAAI,CAAC,KAAK,IAAI,CAAC,KAAK,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,GAAG,CAAC,EAAE;gBACnC,MAAM;aACP;SACF;QACD,CAAC,EAAE,CAAC;QACJ,IAAI,CAAC,GAAG,SAAS,CAAC,MAAM,GAAG,CAAC,EAAE;YAC5B,MAAM,IAAI,KAAK,CAAC,qBAAqB,CAAC,CAAC;SACxC;KACF;IACD,OAAO,CAAC,CAAC;AACX,CAAC"}

package/lib-esm/operations/thresholds/li.d.ts

@@ -0,0 +1,9 @@
+/**
+ * Return a threshold for a histogram using Li algorithm.
+ *
+ * @param histogram - Image histogram.
+ * @param total - Number of pixels in the image.
+ * @returns The threshold.
+ */
+export default function li(histogram: Uint32Array, total: number): number;
+//# sourceMappingURL=li.d.ts.map

package/lib-esm/operations/thresholds/li.d.ts.map

@@ -0,0 +1 @@
+{"version":3,"file":"li.d.ts","sourceRoot":"","sources":["../../../src/operations/thresholds/li.ts"],"names":[],"mappings":"AAcA;;;;;;GAMG;AACH,MAAM,CAAC,OAAO,UAAU,EAAE,CAAC,SAAS,EAAE,WAAW,EAAE,KAAK,EAAE,MAAM,GAAG,MAAM,CA2DxE"}

package/lib-esm/operations/thresholds/li.js

@@ -0,0 +1,74 @@
+/*
+ * see http://rsb.info.nih.gov/ij/developer/source/ij/process/AutoThresholder.java.html
+ * The method is present in: Implements Li's Minimum Cross Entropy thresholding method
+ * This implementation is based on the iterative version (Ref. 2nd reference below) of the algorithm.
+ *  1) Li, CH & Lee, CK (1993), "Minimum Cross Entropy Thresholding", Pattern Recognition 26(4): 61 625
+ *  2) Li, CH & Tam, PKS (1998), "An Iterative Algorithm for Minimum Cross Entropy Thresholding",
+ *     Pattern Recognition Letters 18(8): 771-776
+ *  3) Sezgin, M & Sankur, B (2004), "Survey over Image Thresholding Techniques and Quantitative Performance
+ *     Evaluation",Journal of Electronic Imaging 13(1): 146-165
+ * @param histogram - the histogram of the image
+ * @param total - the number of pixels in the image
+ * @returns {number} - the threshold
+ */
+/**
+ * Return a threshold for a histogram using Li algorithm.
+ *
+ * @param histogram - Image histogram.
+ * @param total - Number of pixels in the image.
+ * @returns The threshold.
+ */
+export default function li(histogram, total) {
+    let threshold;
+    let sumBack; /* sum of the background pixels at a given threshold */
+    let sumObj; /* sum of the object pixels at a given threshold */
+    let numBack; /* number of background pixels at a given threshold */
+    let numObj; /* number of object pixels at a given threshold */
+    let oldThresh;
+    let newThresh;
+    let meanBack; /* mean of the background pixels at a given threshold */
+    let meanObj; /* mean of the object pixels at a given threshold */
+    let mean; /* mean gray-level in the image */
+    const tolerance = 0.5; /* threshold tolerance */
+    let temp;
+    /* Calculate the mean gray-level */
+    mean = 0.0;
+    for (let ih = 0; ih < histogram.length; ih++) {
+        mean += ih * histogram[ih];
+    }
+    mean /= total;
+    /* Initial estimate */
+    newThresh = mean;
+    do {
+        oldThresh = newThresh;
+        threshold = (oldThresh + 0.5) | 0; /* range */
+        /* Calculate the means of background and object pixels */
+        /* Background */
+        sumBack = 0;
+        numBack = 0;
+        for (let ih = 0; ih <= threshold; ih++) {
+            sumBack += ih * histogram[ih];
+            numBack += histogram[ih];
+        }
+        meanBack = numBack === 0 ? 0.0 : sumBack / numBack;
+        /* Object */
+        sumObj = 0;
+        numObj = 0;
+        for (let ih = threshold + 1; ih < histogram.length; ih++) {
+            sumObj += ih * histogram[ih];
+            numObj += histogram[ih];
+        }
+        meanObj = numObj === 0 ? 0.0 : sumObj / numObj;
+        temp = (meanBack - meanObj) / (Math.log(meanBack) - Math.log(meanObj));
+        if (temp < -Number.EPSILON) {
+            newThresh = (temp - 0.5) | 0;
+        }
+        else {
+            newThresh = (temp + 0.5) | 0;
+        }
+        /*  Stop the iterations when the difference between the
+             new and old threshold values is less than the tolerance */
+    } while (Math.abs(newThresh - oldThresh) > tolerance);
+    return threshold;
+}
+//# sourceMappingURL=li.js.map

package/lib-esm/operations/thresholds/li.js.map

@@ -0,0 +1 @@
+{"version":3,"file":"li.js","sourceRoot":"","sources":["../../../src/operations/thresholds/li.ts"],"names":[],"mappings":"AAAA;;;;;;;;;;;;GAYG;AAEH;;;;;;GAMG;AACH,MAAM,CAAC,OAAO,UAAU,EAAE,CAAC,SAAsB,EAAE,KAAa;IAC9D,IAAI,SAAS,CAAC;IACd,IAAI,OAAO,CAAC,CAAC,uDAAuD;IACpE,IAAI,MAAM,CAAC,CAAC,mDAAmD;IAC/D,IAAI,OAAO,CAAC,CAAC,sDAAsD;IACnE,IAAI,MAAM,CAAC,CAAC,kDAAkD;IAC9D,IAAI,SAAS,CAAC;IACd,IAAI,SAAS,CAAC;IACd,IAAI,QAAQ,CAAC,CAAC,wDAAwD;IACtE,IAAI,OAAO,CAAC,CAAC,oDAAoD;IACjE,IAAI,IAAI,CAAC,CAAC,kCAAkC;IAC5C,MAAM,SAAS,GAAG,GAAG,CAAC,CAAC,yBAAyB;IAChD,IAAI,IAAI,CAAC;IAET,mCAAmC;IACnC,IAAI,GAAG,GAAG,CAAC;IACX,KAAK,IAAI,EAAE,GAAG,CAAC,EAAE,EAAE,GAAG,SAAS,CAAC,MAAM,EAAE,EAAE,EAAE,EAAE;QAC5C,IAAI,IAAI,EAAE,GAAG,SAAS,CAAC,EAAE,CAAC,CAAC;KAC5B;IAED,IAAI,IAAI,KAAK,CAAC;IACd,sBAAsB;IACtB,SAAS,GAAG,IAAI,CAAC;IAEjB,GAAG;QACD,SAAS,GAAG,SAAS,CAAC;QACtB,SAAS,GAAG,CAAC,SAAS,GAAG,GAAG,CAAC,GAAG,CAAC,CAAC,CAAC,WAAW;QAE9C,yDAAyD;QACzD,gBAAgB;QAChB,OAAO,GAAG,CAAC,CAAC;QACZ,OAAO,GAAG,CAAC,CAAC;QAEZ,KAAK,IAAI,EAAE,GAAG,CAAC,EAAE,EAAE,IAAI,SAAS,EAAE,EAAE,EAAE,EAAE;YACtC,OAAO,IAAI,EAAE,GAAG,SAAS,CAAC,EAAE,CAAC,CAAC;YAC9B,OAAO,IAAI,SAAS,CAAC,EAAE,CAAC,CAAC;SAC1B;QACD,QAAQ,GAAG,OAAO,KAAK,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,OAAO,GAAG,OAAO,CAAC;QAEnD,YAAY;QACZ,MAAM,GAAG,CAAC,CAAC;QACX,MAAM,GAAG,CAAC,CAAC;QACX,KAAK,IAAI,EAAE,GAAG,SAAS,GAAG,CAAC,EAAE,EAAE,GAAG,SAAS,CAAC,MAAM,EAAE,EAAE,EAAE,EAAE;YACxD,MAAM,IAAI,EAAE,GAAG,SAAS,CAAC,EAAE,CAAC,CAAC;YAC7B,MAAM,IAAI,SAAS,CAAC,EAAE,CAAC,CAAC;SACzB;QACD,OAAO,GAAG,MAAM,KAAK,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,MAAM,GAAG,MAAM,CAAC;QAC/C,IAAI,GAAG,CAAC,QAAQ,GAAG,OAAO,CAAC,GAAG,CAAC,IAAI,CAAC,GAAG,CAAC,QAAQ,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,OAAO,CAAC,CAAC,CAAC;QAEvE,IAAI,IAAI,GAAG,CAAC,MAAM,CAAC,OAAO,EAAE;YAC1B,SAAS,GAAG,CAAC,IAAI,GAAG,GAAG,CAAC,GAAG,CAAC,CAAC;SAC9B;aAAM;YACL,SAAS,GAAG,CAAC,IAAI,GAAG,GAAG,CAAC,GAAG,CAAC,CAAC;SAC9B;QACD;uEAC+D;KAChE,QAAQ,IAAI,CAAC,GAAG,CAAC,SAAS,GAAG,SAAS,CAAC,GAAG,SAAS,EAAE;IAEtD,OAAO,SAAS,CAAC;AACnB,CAAC"}

package/lib-esm/operations/thresholds/maxEntropy.d.ts

@@ -0,0 +1,9 @@
+/**
+ * Return a threshold for a histogram using maximal enthropy algorithm.
+ *
+ * @param histogram - Image histogram.
+ * @param total - Number of pixels in the image.
+ * @returns The threshold.
+ */
+export default function maxEntropy(histogram: Uint32Array, total: number): number;
+//# sourceMappingURL=maxEntropy.d.ts.map

package/lib-esm/operations/thresholds/maxEntropy.d.ts.map

@@ -0,0 +1 @@
+{"version":3,"file":"maxEntropy.d.ts","sourceRoot":"","sources":["../../../src/operations/thresholds/maxEntropy.ts"],"names":[],"mappings":"AAUA;;;;;;GAMG;AACH,MAAM,CAAC,OAAO,UAAU,UAAU,CAChC,SAAS,EAAE,WAAW,EACtB,KAAK,EAAE,MAAM,GACZ,MAAM,CAoER"}

package/lib-esm/operations/thresholds/maxEntropy.js

@@ -0,0 +1,77 @@
+/*
+ * see http://rsb.info.nih.gov/ij/developer/source/ij/process/AutoThresholder.java.html
+ * The method is present in: Implements Kapur-Sahoo-Wong (Maximum Entropy) thresholding method:
+ * Kapur, JN; Sahoo, PK & Wong, ACK (1985), "A New Method for Gray-Level Picture Thresholding Using the Entropy of the Histogram",
+ * Graphical Models and Image Processing 29(3): 273-285
+ * @param histogram - the histogram of the image
+ *        total - the number of pixels in the image
+ * @returns {number} - the threshold
+ */
+/**
+ * Return a threshold for a histogram using maximal enthropy algorithm.
+ *
+ * @param histogram - Image histogram.
+ * @param total - Number of pixels in the image.
+ * @returns The threshold.
+ */
+export default function maxEntropy(histogram, total) {
+    const normHisto = new Array(histogram.length); // normalized histogram
+    for (let ih = 0; ih < histogram.length; ih++) {
+        normHisto[ih] = histogram[ih] / total;
+    }
+    const P1 = new Array(histogram.length); // cumulative normalized histogram
+    const P2 = new Array(histogram.length);
+    P1[0] = normHisto[0];
+    P2[0] = 1.0 - P1[0];
+    for (let ih = 1; ih < histogram.length; ih++) {
+        P1[ih] = P1[ih - 1] + normHisto[ih];
+        P2[ih] = 1.0 - P1[ih];
+    }
+    /* Determine the first non-zero bin */
+    let firstBin = 0;
+    for (let ih = 0; ih < histogram.length; ih++) {
+        if (Math.abs(P1[ih]) >= Number.EPSILON) {
+            firstBin = ih;
+            break;
+        }
+    }
+    /* Determine the last non-zero bin */
+    let lastBin = histogram.length - 1;
+    for (let ih = histogram.length - 1; ih >= firstBin; ih--) {
+        if (Math.abs(P2[ih]) >= Number.EPSILON) {
+            lastBin = ih;
+            break;
+        }
+    }
+    // Calculate the total entropy each gray-level
+    // and find the threshold that maximizes it
+    let threshold = -1;
+    let totEnt; // total entropy
+    let maxEnt = Number.MIN_VALUE; // max entropy
+    let entBack; // entropy of the background pixels at a given threshold
+    let entObj; // entropy of the object pixels at a given threshold
+    for (let it = firstBin; it <= lastBin; it++) {
+        /* Entropy of the background pixels */
+        entBack = 0.0;
+        for (let ih = 0; ih <= it; ih++) {
+            if (histogram[ih] !== 0) {
+                entBack -= (normHisto[ih] / P1[it]) * Math.log(normHisto[ih] / P1[it]);
+            }
+        }
+        /* Entropy of the object pixels */
+        entObj = 0.0;
+        for (let ih = it + 1; ih < histogram.length; ih++) {
+            if (histogram[ih] !== 0) {
+                entObj -= (normHisto[ih] / P2[it]) * Math.log(normHisto[ih] / P2[it]);
+            }
+        }
+        /* Total entropy */
+        totEnt = entBack + entObj;
+        if (maxEnt < totEnt) {
+            maxEnt = totEnt;
+            threshold = it;
+        }
+    }
+    return threshold;
+}
+//# sourceMappingURL=maxEntropy.js.map

package/lib-esm/operations/thresholds/maxEntropy.js.map

@@ -0,0 +1 @@
+{"version":3,"file":"maxEntropy.js","sourceRoot":"","sources":["../../../src/operations/thresholds/maxEntropy.ts"],"names":[],"mappings":"AAAA;;;;;;;;GAQG;AAEH;;;;;;GAMG;AACH,MAAM,CAAC,OAAO,UAAU,UAAU,CAChC,SAAsB,EACtB,KAAa;IAEb,MAAM,SAAS,GAAG,IAAI,KAAK,CAAC,SAAS,CAAC,MAAM,CAAC,CAAC,CAAC,uBAAuB;IACtE,KAAK,IAAI,EAAE,GAAG,CAAC,EAAE,EAAE,GAAG,SAAS,CAAC,MAAM,EAAE,EAAE,EAAE,EAAE;QAC5C,SAAS,CAAC,EAAE,CAAC,GAAG,SAAS,CAAC,EAAE,CAAC,GAAG,KAAK,CAAC;KACvC;IAED,MAAM,EAAE,GAAG,IAAI,KAAK,CAAC,SAAS,CAAC,MAAM,CAAC,CAAC,CAAC,kCAAkC;IAC1E,MAAM,EAAE,GAAG,IAAI,KAAK,CAAC,SAAS,CAAC,MAAM,CAAC,CAAC;IACvC,EAAE,CAAC,CAAC,CAAC,GAAG,SAAS,CAAC,CAAC,CAAC,CAAC;IACrB,EAAE,CAAC,CAAC,CAAC,GAAG,GAAG,GAAG,EAAE,CAAC,CAAC,CAAC,CAAC;IAEpB,KAAK,IAAI,EAAE,GAAG,CAAC,EAAE,EAAE,GAAG,SAAS,CAAC,MAAM,EAAE,EAAE,EAAE,EAAE;QAC5C,EAAE,CAAC,EAAE,CAAC,GAAG,EAAE,CAAC,EAAE,GAAG,CAAC,CAAC,GAAG,SAAS,CAAC,EAAE,CAAC,CAAC;QACpC,EAAE,CAAC,EAAE,CAAC,GAAG,GAAG,GAAG,EAAE,CAAC,EAAE,CAAC,CAAC;KACvB;IAED,sCAAsC;IACtC,IAAI,QAAQ,GAAG,CAAC,CAAC;IACjB,KAAK,IAAI,EAAE,GAAG,CAAC,EAAE,EAAE,GAAG,SAAS,CAAC,MAAM,EAAE,EAAE,EAAE,EAAE;QAC5C,IAAI,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,EAAE,CAAC,CAAC,IAAI,MAAM,CAAC,OAAO,EAAE;YACtC,QAAQ,GAAG,EAAE,CAAC;YACd,MAAM;SACP;KACF;IAED,qCAAqC;IACrC,IAAI,OAAO,GAAG,SAAS,CAAC,MAAM,GAAG,CAAC,CAAC;IACnC,KAAK,IAAI,EAAE,GAAG,SAAS,CAAC,MAAM,GAAG,CAAC,EAAE,EAAE,IAAI,QAAQ,EAAE,EAAE,EAAE,EAAE;QACxD,IAAI,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,EAAE,CAAC,CAAC,IAAI,MAAM,CAAC,OAAO,EAAE;YACtC,OAAO,GAAG,EAAE,CAAC;YACb,MAAM;SACP;KACF;IAED,8CAA8C;IAC9C,2CAA2C;IAC3C,IAAI,SAAS,GAAG,CAAC,CAAC,CAAC;IACnB,IAAI,MAAM,CAAC,CAAC,gBAAgB;IAC5B,IAAI,MAAM,GAAG,MAAM,CAAC,SAAS,CAAC,CAAC,cAAc;IAC7C,IAAI,OAAO,CAAC,CAAC,wDAAwD;IACrE,IAAI,MAAM,CAAC,CAAC,oDAAoD;IAEhE,KAAK,IAAI,EAAE,GAAG,QAAQ,EAAE,EAAE,IAAI,OAAO,EAAE,EAAE,EAAE,EAAE;QAC3C,sCAAsC;QACtC,OAAO,GAAG,GAAG,CAAC;QACd,KAAK,IAAI,EAAE,GAAG,CAAC,EAAE,EAAE,IAAI,EAAE,EAAE,EAAE,EAAE,EAAE;YAC/B,IAAI,SAAS,CAAC,EAAE,CAAC,KAAK,CAAC,EAAE;gBACvB,OAAO,IAAI,CAAC,SAAS,CAAC,EAAE,CAAC,GAAG,EAAE,CAAC,EAAE,CAAC,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,SAAS,CAAC,EAAE,CAAC,GAAG,EAAE,CAAC,EAAE,CAAC,CAAC,CAAC;aACxE;SACF;QAED,kCAAkC;QAClC,MAAM,GAAG,GAAG,CAAC;QACb,KAAK,IAAI,EAAE,GAAG,EAAE,GAAG,CAAC,EAAE,EAAE,GAAG,SAAS,CAAC,MAAM,EAAE,EAAE,EAAE,EAAE;YACjD,IAAI,SAAS,CAAC,EAAE,CAAC,KAAK,CAAC,EAAE;gBACvB,MAAM,IAAI,CAAC,SAAS,CAAC,EAAE,CAAC,GAAG,EAAE,CAAC,EAAE,CAAC,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,SAAS,CAAC,EAAE,CAAC,GAAG,EAAE,CAAC,EAAE,CAAC,CAAC,CAAC;aACvE;SACF;QAED,mBAAmB;QACnB,MAAM,GAAG,OAAO,GAAG,MAAM,CAAC;QAE1B,IAAI,MAAM,GAAG,MAAM,EAAE;YACnB,MAAM,GAAG,MAAM,CAAC;YAChB,SAAS,GAAG,EAAE,CAAC;SAChB;KACF;IACD,OAAO,SAAS,CAAC;AACnB,CAAC"}

package/lib-esm/operations/thresholds/mean.d.ts

@@ -0,0 +1,9 @@
+/**
+ * Return a threshold for a histogram by making its average.
+ *
+ * @param histogram - Image histogram.
+ * @param total - Number of pixels in the image.
+ * @returns The threshold.
+ */
+export default function mean(histogram: Uint32Array, total: number): number;
+//# sourceMappingURL=mean.d.ts.map

package/lib-esm/operations/thresholds/mean.d.ts.map

@@ -0,0 +1 @@
+{"version":3,"file":"mean.d.ts","sourceRoot":"","sources":["../../../src/operations/thresholds/mean.ts"],"names":[],"mappings":"AASA;;;;;;GAMG;AACH,MAAM,CAAC,OAAO,UAAU,IAAI,CAAC,SAAS,EAAE,WAAW,EAAE,KAAK,EAAE,MAAM,GAAG,MAAM,CAM1E"}

package/lib-esm/operations/thresholds/mean.js

@@ -0,0 +1,23 @@
+/*
+ * The method is present in: Uses the mean of grey levels as the threshold. It is described in:
+ * Glasbey, CA (1993), "An analysis of histogram-based thresholding algorithms",
+ * CVGIP: Graphical Models and Image Processing 55: 532-537
+ * @param histogram - the histogram of the image
+ * @param total - the number of pixels in the image
+ * @returns {number} - the threshold
+ */
+/**
+ * Return a threshold for a histogram by making its average.
+ *
+ * @param histogram - Image histogram.
+ * @param total - Number of pixels in the image.
+ * @returns The threshold.
+ */
+export default function mean(histogram, total) {
+    let sum = 0;
+    for (let i = 0; i < histogram.length; i++) {
+        sum += i * histogram[i];
+    }
+    return Math.floor(sum / total);
+}
+//# sourceMappingURL=mean.js.map

package/lib-esm/operations/thresholds/mean.js.map

@@ -0,0 +1 @@
+{"version":3,"file":"mean.js","sourceRoot":"","sources":["../../../src/operations/thresholds/mean.ts"],"names":[],"mappings":"AAAA;;;;;;;GAOG;AAEH;;;;;;GAMG;AACH,MAAM,CAAC,OAAO,UAAU,IAAI,CAAC,SAAsB,EAAE,KAAa;IAChE,IAAI,GAAG,GAAG,CAAC,CAAC;IACZ,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,SAAS,CAAC,MAAM,EAAE,CAAC,EAAE,EAAE;QACzC,GAAG,IAAI,CAAC,GAAG,SAAS,CAAC,CAAC,CAAC,CAAC;KACzB;IACD,OAAO,IAAI,CAAC,KAAK,CAAC,GAAG,GAAG,KAAK,CAAC,CAAC;AACjC,CAAC"}

package/lib-esm/operations/thresholds/minError.d.ts

@@ -0,0 +1,9 @@
+/**
+ * Return a threshold for a histogram using Li algorithm.
+ *
+ * @param histogram - Image histogram.
+ * @param total - Number of pixels in the image.
+ * @returns The threshold.
+ */
+export default function minError(histogram: Uint32Array, total: number): number;
+//# sourceMappingURL=minError.d.ts.map

package/lib-esm/operations/thresholds/minError.d.ts.map

@@ -0,0 +1 @@
+{"version":3,"file":"minError.d.ts","sourceRoot":"","sources":["../../../src/operations/thresholds/minError.ts"],"names":[],"mappings":"AASA;;;;;;GAMG;AACH,MAAM,CAAC,OAAO,UAAU,QAAQ,CAC9B,SAAS,EAAE,WAAW,EACtB,KAAK,EAAE,MAAM,GACZ,MAAM,CAwDR"}

package/lib-esm/operations/thresholds/minError.js

@@ -0,0 +1,87 @@
+/*
+ * see http://rsb.info.nih.gov/ij/developer/source/ij/process/AutoThresholder.java.html
+ * The method is present in: An iterative implementation of Kittler and Illingworth's Minimum Error
+ * thresholding:Kittler, J & Illingworth, J (1986), "Minimum error thresholding", Pattern Recognition 19: 41-47
+ * @param histogram - the histogram of the image
+ * @param total - the number of pixels in the image
+ * @returns {number} - the threshold
+ */
+/**
+ * Return a threshold for a histogram using Li algorithm.
+ *
+ * @param histogram - Image histogram.
+ * @param total - Number of pixels in the image.
+ * @returns The threshold.
+ */
+export default function minError(histogram, total) {
+    let threshold;
+    let Tprev = -2;
+    let mu, nu, p, q, sigma2, tau2, w0, w1, w2, sqterm, temp;
+    /* Calculate the mean gray-level */
+    let mean = 0.0;
+    for (let ih = 0; ih < histogram.length; ih++) {
+        mean += ih * histogram[ih];
+    }
+    mean /= total;
+    threshold = mean;
+    while (threshold !== Tprev) {
+        // Calculate some statistics.
+        const sumA1 = sumA(histogram, threshold);
+        const sumA2 = sumA(histogram, histogram.length - 1);
+        const sumB1 = sumB(histogram, threshold);
+        const sumB2 = sumB(histogram, histogram.length - 1);
+        const sumC1 = sumC(histogram, threshold);
+        const sumC2 = sumC(histogram, histogram.length - 1);
+        mu = sumB1 / sumA1;
+        nu = (sumB2 - sumB1) / (sumA2 - sumA1);
+        p = sumA1 / sumA2;
+        q = (sumA2 - sumA1) / sumA2;
+        sigma2 = sumC1 / sumA1 - mu * mu;
+        tau2 = (sumC2 - sumC1) / (sumA2 - sumA1) - nu * nu;
+        // The terms of the quadratic equation to be solved.
+        w0 = 1.0 / sigma2 - 1.0 / tau2;
+        w1 = mu / sigma2 - nu / tau2;
+        w2 =
+            (mu * mu) / sigma2 -
+                (nu * nu) / tau2 +
+                Math.log10((sigma2 * (q * q)) / (tau2 * (p * p)));
+        // If the next threshold would be imaginary, return with the current one.
+        sqterm = w1 * w1 - w0 * w2;
+        if (sqterm < 0) {
+            return threshold;
+        }
+        // The updated threshold is the integer part of the solution of the quadratic equation.
+        Tprev = threshold;
+        temp = (w1 + Math.sqrt(sqterm)) / w0;
+        if (isNaN(temp)) {
+            threshold = Tprev;
+        }
+        else {
+            threshold = Math.floor(temp);
+        }
+    }
+    return threshold;
+}
+// aux func
+function sumA(y, j) {
+    let x = 0;
+    for (let i = 0; i <= j; i++) {
+        x += y[i];
+    }
+    return x;
+}
+function sumB(y, j) {
+    let x = 0;
+    for (let i = 0; i <= j; i++) {
+        x += i * y[i];
+    }
+    return x;
+}
+function sumC(y, j) {
+    let x = 0;
+    for (let i = 0; i <= j; i++) {
+        x += i * i * y[i];
+    }
+    return x;
+}
+//# sourceMappingURL=minError.js.map

package/lib-esm/operations/thresholds/minError.js.map

@@ -0,0 +1 @@
+{"version":3,"file":"minError.js","sourceRoot":"","sources":["../../../src/operations/thresholds/minError.ts"],"names":[],"mappings":"AAAA;;;;;;;GAOG;AAEH;;;;;;GAMG;AACH,MAAM,CAAC,OAAO,UAAU,QAAQ,CAC9B,SAAsB,EACtB,KAAa;IAEb,IAAI,SAAS,CAAC;IACd,IAAI,KAAK,GAAG,CAAC,CAAC,CAAC;IACf,IAAI,EAAE,EAAE,EAAE,EAAE,CAAC,EAAE,CAAC,EAAE,MAAM,EAAE,IAAI,EAAE,EAAE,EAAE,EAAE,EAAE,EAAE,EAAE,MAAM,EAAE,IAAI,CAAC;IAEzD,mCAAmC;IACnC,IAAI,IAAI,GAAG,GAAG,CAAC;IACf,KAAK,IAAI,EAAE,GAAG,CAAC,EAAE,EAAE,GAAG,SAAS,CAAC,MAAM,EAAE,EAAE,EAAE,EAAE;QAC5C,IAAI,IAAI,EAAE,GAAG,SAAS,CAAC,EAAE,CAAC,CAAC;KAC5B;IAED,IAAI,IAAI,KAAK,CAAC;IAEd,SAAS,GAAG,IAAI,CAAC;IAEjB,OAAO,SAAS,KAAK,KAAK,EAAE;QAC1B,6BAA6B;QAC7B,MAAM,KAAK,GAAG,IAAI,CAAC,SAAS,EAAE,SAAS,CAAC,CAAC;QACzC,MAAM,KAAK,GAAG,IAAI,CAAC,SAAS,EAAE,SAAS,CAAC,MAAM,GAAG,CAAC,CAAC,CAAC;QACpD,MAAM,KAAK,GAAG,IAAI,CAAC,SAAS,EAAE,SAAS,CAAC,CAAC;QACzC,MAAM,KAAK,GAAG,IAAI,CAAC,SAAS,EAAE,SAAS,CAAC,MAAM,GAAG,CAAC,CAAC,CAAC;QACpD,MAAM,KAAK,GAAG,IAAI,CAAC,SAAS,EAAE,SAAS,CAAC,CAAC;QACzC,MAAM,KAAK,GAAG,IAAI,CAAC,SAAS,EAAE,SAAS,CAAC,MAAM,GAAG,CAAC,CAAC,CAAC;QAEpD,EAAE,GAAG,KAAK,GAAG,KAAK,CAAC;QACnB,EAAE,GAAG,CAAC,KAAK,GAAG,KAAK,CAAC,GAAG,CAAC,KAAK,GAAG,KAAK,CAAC,CAAC;QACvC,CAAC,GAAG,KAAK,GAAG,KAAK,CAAC;QAClB,CAAC,GAAG,CAAC,KAAK,GAAG,KAAK,CAAC,GAAG,KAAK,CAAC;QAC5B,MAAM,GAAG,KAAK,GAAG,KAAK,GAAG,EAAE,GAAG,EAAE,CAAC;QACjC,IAAI,GAAG,CAAC,KAAK,GAAG,KAAK,CAAC,GAAG,CAAC,KAAK,GAAG,KAAK,CAAC,GAAG,EAAE,GAAG,EAAE,CAAC;QAEnD,oDAAoD;QACpD,EAAE,GAAG,GAAG,GAAG,MAAM,GAAG,GAAG,GAAG,IAAI,CAAC;QAC/B,EAAE,GAAG,EAAE,GAAG,MAAM,GAAG,EAAE,GAAG,IAAI,CAAC;QAC7B,EAAE;YACA,CAAC,EAAE,GAAG,EAAE,CAAC,GAAG,MAAM;gBAClB,CAAC,EAAE,GAAG,EAAE,CAAC,GAAG,IAAI;gBAChB,IAAI,CAAC,KAAK,CAAC,CAAC,MAAM,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,GAAG,CAAC,IAAI,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC;QAEpD,yEAAyE;QACzE,MAAM,GAAG,EAAE,GAAG,EAAE,GAAG,EAAE,GAAG,EAAE,CAAC;QAC3B,IAAI,MAAM,GAAG,CAAC,EAAE;YACd,OAAO,SAAS,CAAC;SAClB;QAED,uFAAuF;QACvF,KAAK,GAAG,SAAS,CAAC;QAClB,IAAI,GAAG,CAAC,EAAE,GAAG,IAAI,CAAC,IAAI,CAAC,MAAM,CAAC,CAAC,GAAG,EAAE,CAAC;QAErC,IAAI,KAAK,CAAC,IAAI,CAAC,EAAE;YACf,SAAS,GAAG,KAAK,CAAC;SACnB;aAAM;YACL,SAAS,GAAG,IAAI,CAAC,KAAK,CAAC,IAAI,CAAC,CAAC;SAC9B;KACF;IACD,OAAO,SAAS,CAAC;AACnB,CAAC;AAED,WAAW;AAEX,SAAS,IAAI,CAAC,CAAc,EAAE,CAAS;IACrC,IAAI,CAAC,GAAG,CAAC,CAAC;IACV,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,IAAI,CAAC,EAAE,CAAC,EAAE,EAAE;QAC3B,CAAC,IAAI,CAAC,CAAC,CAAC,CAAC,CAAC;KACX;IACD,OAAO,CAAC,CAAC;AACX,CAAC;AAED,SAAS,IAAI,CAAC,CAAc,EAAE,CAAS;IACrC,IAAI,CAAC,GAAG,CAAC,CAAC;IACV,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,IAAI,CAAC,EAAE,CAAC,EAAE,EAAE;QAC3B,CAAC,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC;KACf;IACD,OAAO,CAAC,CAAC;AACX,CAAC;AAED,SAAS,IAAI,CAAC,CAAc,EAAE,CAAS;IACrC,IAAI,CAAC,GAAG,CAAC,CAAC;IACV,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,IAAI,CAAC,EAAE,CAAC,EAAE,EAAE;QAC3B,CAAC,IAAI,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC;KACnB;IACD,OAAO,CAAC,CAAC;AACX,CAAC"}

package/lib-esm/operations/thresholds/minimum.d.ts

@@ -0,0 +1,8 @@
+/**
+ * Return a threshold for a histogram.
+ *
+ * @param histogram - Image histogram.
+ * @returns The threshold.
+ */
+export default function minimum(histogram: Uint32Array): number;
+//# sourceMappingURL=minimum.d.ts.map

package/lib-esm/operations/thresholds/minimum.d.ts.map

@@ -0,0 +1 @@
+{"version":3,"file":"minimum.d.ts","sourceRoot":"","sources":["../../../src/operations/thresholds/minimum.ts"],"names":[],"mappings":"AAQA;;;;;GAKG;AACH,MAAM,CAAC,OAAO,UAAU,OAAO,CAAC,SAAS,EAAE,WAAW,GAAG,MAAM,CAiC9D"}

package/lib-esm/operations/thresholds/minimum.js

@@ -0,0 +1,75 @@
+// see https://github.com/fiji/Auto_Threshold/blob/master/src/main/java/fiji/threshold/Auto_Threshold.java
+// J. M. S. Prewitt and M. L. Mendelsohn, "The analysis of cell images," in
+// Annals of the New York Academy of Sciences, vol. 128, pp. 1035-1053, 1966.
+// ported to ImageJ plugin by G.Landini from Antti Niemisto's Matlab code (GPL)
+// Original Matlab code Copyright (C) 2004 Antti Niemisto
+// See http://www.cs.tut.fi/~ant/histthresh/ for an excellent slide presentation
+// and the original Matlab code
+/**
+ * Return a threshold for a histogram.
+ *
+ * @param histogram - Image histogram.
+ * @returns The threshold.
+ */
+export default function minimum(histogram) {
+    if (histogram.length < 2) {
+        // validate that the histogram has at least two color values
+        return 0;
+    }
+    let iterations = 0; // number of iterations of the smoothing process
+    let threshold = -1;
+    let max = -1; // maximum color value with a greater number of pixels to 0
+    let histogramCopy = new Array(histogram.length); // a copy of the histogram
+    for (let i = 0; i < histogram.length; i++) {
+        histogramCopy[i] = histogram[i];
+        if (histogram[i] > 0) {
+            max = i;
+        }
+    }
+    while (!bimodalTest(histogramCopy)) {
+        histogramCopy = smoothed(histogramCopy);
+        iterations++;
+        if (iterations > 10000) {
+            // if they occur more than 10000 iterations it returns -1
+            return threshold;
+        }
+    }
+    for (let i = 1; i < max; i++) {
+        if (histogramCopy[i - 1] > histogramCopy[i] &&
+            histogramCopy[i + 1] >= histogramCopy[i]) {
+            threshold = i;
+            break;
+        }
+    }
+    return threshold;
+}
+function smoothed(histogram) {
+    // Smooth with a 3 point running mean filter
+    const auHistogram = new Array(histogram.length); // a copy of the histograma for the smoothing process
+    for (let i = 1; i < histogram.length - 1; i++) {
+        auHistogram[i] = (histogram[i - 1] + histogram[i] + histogram[i + 1]) / 3;
+    }
+    auHistogram[0] = (histogram[0] + histogram[1]) / 3;
+    auHistogram[histogram.length - 1] =
+        (histogram[histogram.length - 2] + histogram[histogram.length - 1]) / 3;
+    return auHistogram;
+}
+function bimodalTest(histogram) {
+    // It is responsible for determining if a histogram is bimodal
+    const len = histogram.length;
+    let isBimodal = false;
+    let peaks = 0;
+    for (let k = 1; k < len - 1; k++) {
+        if (histogram[k - 1] < histogram[k] && histogram[k + 1] < histogram[k]) {
+            peaks++;
+            if (peaks > 2) {
+                return false;
+            }
+        }
+    }
+    if (peaks === 2) {
+        isBimodal = true;
+    }
+    return isBimodal;
+}
+//# sourceMappingURL=minimum.js.map

package/lib-esm/operations/thresholds/minimum.js.map

@@ -0,0 +1 @@
+{"version":3,"file":"minimum.js","sourceRoot":"","sources":["../../../src/operations/thresholds/minimum.ts"],"names":[],"mappings":"AAAA,0GAA0G;AAC1G,2EAA2E;AAC3E,6EAA6E;AAC7E,+EAA+E;AAC/E,yDAAyD;AACzD,gFAAgF;AAChF,+BAA+B;AAE/B;;;;;GAKG;AACH,MAAM,CAAC,OAAO,UAAU,OAAO,CAAC,SAAsB;IACpD,IAAI,SAAS,CAAC,MAAM,GAAG,CAAC,EAAE;QACxB,4DAA4D;QAC5D,OAAO,CAAC,CAAC;KACV;IACD,IAAI,UAAU,GAAG,CAAC,CAAC,CAAC,gDAAgD;IACpE,IAAI,SAAS,GAAG,CAAC,CAAC,CAAC;IACnB,IAAI,GAAG,GAAG,CAAC,CAAC,CAAC,CAAC,2DAA2D;IACzE,IAAI,aAAa,GAAG,IAAI,KAAK,CAAC,SAAS,CAAC,MAAM,CAAC,CAAC,CAAC,0BAA0B;IAC3E,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,SAAS,CAAC,MAAM,EAAE,CAAC,EAAE,EAAE;QACzC,aAAa,CAAC,CAAC,CAAC,GAAG,SAAS,CAAC,CAAC,CAAC,CAAC;QAChC,IAAI,SAAS,CAAC,CAAC,CAAC,GAAG,CAAC,EAAE;YACpB,GAAG,GAAG,CAAC,CAAC;SACT;KACF;IACD,OAAO,CAAC,WAAW,CAAC,aAAa,CAAC,EAAE;QAClC,aAAa,GAAG,QAAQ,CAAC,aAAa,CAAC,CAAC;QACxC,UAAU,EAAE,CAAC;QACb,IAAI,UAAU,GAAG,KAAK,EAAE;YACtB,yDAAyD;YACzD,OAAO,SAAS,CAAC;SAClB;KACF;IACD,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,GAAG,EAAE,CAAC,EAAE,EAAE;QAC5B,IACE,aAAa,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,aAAa,CAAC,CAAC,CAAC;YACvC,aAAa,CAAC,CAAC,GAAG,CAAC,CAAC,IAAI,aAAa,CAAC,CAAC,CAAC,EACxC;YACA,SAAS,GAAG,CAAC,CAAC;YACd,MAAM;SACP;KACF;IACD,OAAO,SAAS,CAAC;AACnB,CAAC;AACD,SAAS,QAAQ,CAAC,SAAmB;IACnC,4CAA4C;IAC5C,MAAM,WAAW,GAAG,IAAI,KAAK,CAAC,SAAS,CAAC,MAAM,CAAC,CAAC,CAAC,qDAAqD;IACtG,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,SAAS,CAAC,MAAM,GAAG,CAAC,EAAE,CAAC,EAAE,EAAE;QAC7C,WAAW,CAAC,CAAC,CAAC,GAAG,CAAC,SAAS,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,SAAS,CAAC,CAAC,CAAC,GAAG,SAAS,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC;KAC3E;IACD,WAAW,CAAC,CAAC,CAAC,GAAG,CAAC,SAAS,CAAC,CAAC,CAAC,GAAG,SAAS,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC;IACnD,WAAW,CAAC,SAAS,CAAC,MAAM,GAAG,CAAC,CAAC;QAC/B,CAAC,SAAS,CAAC,SAAS,CAAC,MAAM,GAAG,CAAC,CAAC,GAAG,SAAS,CAAC,SAAS,CAAC,MAAM,GAAG,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC;IAC1E,OAAO,WAAW,CAAC;AACrB,CAAC;AAED,SAAS,WAAW,CAAC,SAAmB;IACtC,8DAA8D;IAC9D,MAAM,GAAG,GAAG,SAAS,CAAC,MAAM,CAAC;IAC7B,IAAI,SAAS,GAAG,KAAK,CAAC;IACtB,IAAI,KAAK,GAAG,CAAC,CAAC;IACd,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,GAAG,GAAG,CAAC,EAAE,CAAC,EAAE,EAAE;QAChC,IAAI,SAAS,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,SAAS,CAAC,CAAC,CAAC,IAAI,SAAS,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,SAAS,CAAC,CAAC,CAAC,EAAE;YACtE,KAAK,EAAE,CAAC;YACR,IAAI,KAAK,GAAG,CAAC,EAAE;gBACb,OAAO,KAAK,CAAC;aACd;SACF;KACF;IACD,IAAI,KAAK,KAAK,CAAC,EAAE;QACf,SAAS,GAAG,IAAI,CAAC;KAClB;IACD,OAAO,SAAS,CAAC;AACnB,CAAC"}

package/lib-esm/operations/thresholds/moments.d.ts

@@ -0,0 +1,9 @@
+/**
+ * Return a threshold for a histogram.
+ *
+ * @param histogram - Image histogram.
+ * @param total - Number of pixels in the image.
+ * @returns The threshold.
+ */
+export default function moments(histogram: Uint32Array, total: number): number;
+//# sourceMappingURL=moments.d.ts.map

package/lib-esm/operations/thresholds/moments.d.ts.map

@@ -0,0 +1 @@
+{"version":3,"file":"moments.d.ts","sourceRoot":"","sources":["../../../src/operations/thresholds/moments.ts"],"names":[],"mappings":"AASA;;;;;;GAMG;AACH,MAAM,CAAC,OAAO,UAAU,OAAO,CAAC,SAAS,EAAE,WAAW,EAAE,KAAK,EAAE,MAAM,GAAG,MAAM,CA2C7E"}