image-js 1.0.0-alpha.3 → 1.1.0

package/lib/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;;;;;GAKG;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,CAAC,CAAC;IACT,KAAK,IAAI,EAAE,GAAG,CAAC,EAAE,EAAE,GAAG,SAAS,CAAC,MAAM,EAAE,EAAE,EAAE,EAAE,CAAC;QAC7C,IAAI,IAAI,EAAE,GAAG,SAAS,CAAC,EAAE,CAAC,CAAC;IAC7B,CAAC;IAED,IAAI,IAAI,KAAK,CAAC;IACd,sBAAsB;IACtB,SAAS,GAAG,IAAI,CAAC;IAEjB,GAAG,CAAC;QACF,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,CAAC;YACvC,OAAO,IAAI,EAAE,GAAG,SAAS,CAAC,EAAE,CAAC,CAAC;YAC9B,OAAO,IAAI,SAAS,CAAC,EAAE,CAAC,CAAC;QAC3B,CAAC;QACD,QAAQ,GAAG,OAAO,KAAK,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,OAAO,GAAG,OAAO,CAAC;QAEjD,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,CAAC;YACzD,MAAM,IAAI,EAAE,GAAG,SAAS,CAAC,EAAE,CAAC,CAAC;YAC7B,MAAM,IAAI,SAAS,CAAC,EAAE,CAAC,CAAC;QAC1B,CAAC;QACD,OAAO,GAAG,MAAM,KAAK,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,MAAM,GAAG,MAAM,CAAC;QAC7C,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,CAAC;YAC3B,SAAS,GAAG,CAAC,IAAI,GAAG,GAAG,CAAC,GAAG,CAAC,CAAC;QAC/B,CAAC;aAAM,CAAC;YACN,SAAS,GAAG,CAAC,IAAI,GAAG,GAAG,CAAC,GAAG,CAAC,CAAC;QAC/B,CAAC;QACD;uEAC+D;IACjE,CAAC,QAAQ,IAAI,CAAC,GAAG,CAAC,SAAS,GAAG,SAAS,CAAC,GAAG,SAAS,EAAE;IAEtD,OAAO,SAAS,CAAC;AACnB,CAAC"}

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

@@ -0,0 +1,8 @@
+/**
+ * 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/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;;;;;GAKG;AACH,MAAM,CAAC,OAAO,UAAU,UAAU,CAChC,SAAS,EAAE,WAAW,EACtB,KAAK,EAAE,MAAM,GACZ,MAAM,CAoER"}

package/lib/operations/thresholds/maxEntropy.js

@@ -0,0 +1,76 @@
+/*
+ * 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 - P1[0];
+    for (let ih = 1; ih < histogram.length; ih++) {
+        P1[ih] = P1[ih - 1] + normHisto[ih];
+        P2[ih] = 1 - 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;
+        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;
+        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/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;;;;;GAKG;AACH,MAAM,CAAC,OAAO,UAAU,UAAU,CAChC,SAAsB,EACtB,KAAa;IAEb,MAAM,SAAS,GAAG,IAAI,KAAK,CAAS,SAAS,CAAC,MAAM,CAAC,CAAC,CAAC,uBAAuB;IAC9E,KAAK,IAAI,EAAE,GAAG,CAAC,EAAE,EAAE,GAAG,SAAS,CAAC,MAAM,EAAE,EAAE,EAAE,EAAE,CAAC;QAC7C,SAAS,CAAC,EAAE,CAAC,GAAG,SAAS,CAAC,EAAE,CAAC,GAAG,KAAK,CAAC;IACxC,CAAC;IAED,MAAM,EAAE,GAAG,IAAI,KAAK,CAAS,SAAS,CAAC,MAAM,CAAC,CAAC,CAAC,kCAAkC;IAClF,MAAM,EAAE,GAAG,IAAI,KAAK,CAAS,SAAS,CAAC,MAAM,CAAC,CAAC;IAC/C,EAAE,CAAC,CAAC,CAAC,GAAG,SAAS,CAAC,CAAC,CAAC,CAAC;IACrB,EAAE,CAAC,CAAC,CAAC,GAAG,CAAC,GAAG,EAAE,CAAC,CAAC,CAAC,CAAC;IAElB,KAAK,IAAI,EAAE,GAAG,CAAC,EAAE,EAAE,GAAG,SAAS,CAAC,MAAM,EAAE,EAAE,EAAE,EAAE,CAAC;QAC7C,EAAE,CAAC,EAAE,CAAC,GAAG,EAAE,CAAC,EAAE,GAAG,CAAC,CAAC,GAAG,SAAS,CAAC,EAAE,CAAC,CAAC;QACpC,EAAE,CAAC,EAAE,CAAC,GAAG,CAAC,GAAG,EAAE,CAAC,EAAE,CAAC,CAAC;IACtB,CAAC;IAED,sCAAsC;IACtC,IAAI,QAAQ,GAAG,CAAC,CAAC;IACjB,KAAK,IAAI,EAAE,GAAG,CAAC,EAAE,EAAE,GAAG,SAAS,CAAC,MAAM,EAAE,EAAE,EAAE,EAAE,CAAC;QAC7C,IAAI,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,EAAE,CAAC,CAAC,IAAI,MAAM,CAAC,OAAO,EAAE,CAAC;YACvC,QAAQ,GAAG,EAAE,CAAC;YACd,MAAM;QACR,CAAC;IACH,CAAC;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,CAAC;QACzD,IAAI,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,EAAE,CAAC,CAAC,IAAI,MAAM,CAAC,OAAO,EAAE,CAAC;YACvC,OAAO,GAAG,EAAE,CAAC;YACb,MAAM;QACR,CAAC;IACH,CAAC;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,CAAC;QAC5C,sCAAsC;QACtC,OAAO,GAAG,CAAC,CAAC;QACZ,KAAK,IAAI,EAAE,GAAG,CAAC,EAAE,EAAE,IAAI,EAAE,EAAE,EAAE,EAAE,EAAE,CAAC;YAChC,IAAI,SAAS,CAAC,EAAE,CAAC,KAAK,CAAC,EAAE,CAAC;gBACxB,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;YACzE,CAAC;QACH,CAAC;QAED,kCAAkC;QAClC,MAAM,GAAG,CAAC,CAAC;QACX,KAAK,IAAI,EAAE,GAAG,EAAE,GAAG,CAAC,EAAE,EAAE,GAAG,SAAS,CAAC,MAAM,EAAE,EAAE,EAAE,EAAE,CAAC;YAClD,IAAI,SAAS,CAAC,EAAE,CAAC,KAAK,CAAC,EAAE,CAAC;gBACxB,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;YACxE,CAAC;QACH,CAAC;QAED,mBAAmB;QACnB,MAAM,GAAG,OAAO,GAAG,MAAM,CAAC;QAE1B,IAAI,MAAM,GAAG,MAAM,EAAE,CAAC;YACpB,MAAM,GAAG,MAAM,CAAC;YAChB,SAAS,GAAG,EAAE,CAAC;QACjB,CAAC;IACH,CAAC;IACD,OAAO,SAAS,CAAC;AACnB,CAAC"}

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

@@ -0,0 +1,8 @@
+/**
+ * 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/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;;;;;GAKG;AACH,MAAM,CAAC,OAAO,UAAU,IAAI,CAAC,SAAS,EAAE,WAAW,EAAE,KAAK,EAAE,MAAM,GAAG,MAAM,CAM1E"}

package/lib/operations/thresholds/mean.js

@@ -0,0 +1,22 @@
+/*
+ * 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/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;;;;;GAKG;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,CAAC;QAC1C,GAAG,IAAI,CAAC,GAAG,SAAS,CAAC,CAAC,CAAC,CAAC;IAC1B,CAAC;IACD,OAAO,IAAI,CAAC,KAAK,CAAC,GAAG,GAAG,KAAK,CAAC,CAAC;AACjC,CAAC"}

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

@@ -0,0 +1,8 @@
+/**
+ * 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/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;;;;;GAKG;AACH,MAAM,CAAC,OAAO,UAAU,QAAQ,CAC9B,SAAS,EAAE,WAAW,EACtB,KAAK,EAAE,MAAM,GACZ,MAAM,CAwDR"}

package/lib/operations/thresholds/minError.js

@@ -0,0 +1,86 @@
+/*
+ * 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;
+    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 / sigma2 - 1 / 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 (Number.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/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;;;;;GAKG;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,CAAC,CAAC;IACb,KAAK,IAAI,EAAE,GAAG,CAAC,EAAE,EAAE,GAAG,SAAS,CAAC,MAAM,EAAE,EAAE,EAAE,EAAE,CAAC;QAC7C,IAAI,IAAI,EAAE,GAAG,SAAS,CAAC,EAAE,CAAC,CAAC;IAC7B,CAAC;IAED,IAAI,IAAI,KAAK,CAAC;IAEd,SAAS,GAAG,IAAI,CAAC;IAEjB,OAAO,SAAS,KAAK,KAAK,EAAE,CAAC;QAC3B,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,CAAC,GAAG,MAAM,GAAG,CAAC,GAAG,IAAI,CAAC;QAC3B,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,CAAC;YACf,OAAO,SAAS,CAAC;QACnB,CAAC;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,MAAM,CAAC,KAAK,CAAC,IAAI,CAAC,EAAE,CAAC;YACvB,SAAS,GAAG,KAAK,CAAC;QACpB,CAAC;aAAM,CAAC;YACN,SAAS,GAAG,IAAI,CAAC,KAAK,CAAC,IAAI,CAAC,CAAC;QAC/B,CAAC;IACH,CAAC;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,CAAC;QAC5B,CAAC,IAAI,CAAC,CAAC,CAAC,CAAC,CAAC;IACZ,CAAC;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,CAAC;QAC5B,CAAC,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC;IAChB,CAAC;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,CAAC;QAC5B,CAAC,IAAI,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC;IACpB,CAAC;IACD,OAAO,CAAC,CAAC;AACX,CAAC"}

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

@@ -0,0 +1,7 @@
+/**
+ * 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/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;;;;GAIG;AACH,MAAM,CAAC,OAAO,UAAU,OAAO,CAAC,SAAS,EAAE,WAAW,GAAG,MAAM,CAiC9D"}

package/lib/operations/thresholds/minimum.js

@@ -0,0 +1,74 @@
+// 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 histogram 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.at(-2) + histogram.at(-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/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;;;;GAIG;AACH,MAAM,CAAC,OAAO,UAAU,OAAO,CAAC,SAAsB;IACpD,IAAI,SAAS,CAAC,MAAM,GAAG,CAAC,EAAE,CAAC;QACzB,4DAA4D;QAC5D,OAAO,CAAC,CAAC;IACX,CAAC;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,CAAC;QAC1C,aAAa,CAAC,CAAC,CAAC,GAAG,SAAS,CAAC,CAAC,CAAC,CAAC;QAChC,IAAI,SAAS,CAAC,CAAC,CAAC,GAAG,CAAC,EAAE,CAAC;YACrB,GAAG,GAAG,CAAC,CAAC;QACV,CAAC;IACH,CAAC;IACD,OAAO,CAAC,WAAW,CAAC,aAAa,CAAC,EAAE,CAAC;QACnC,aAAa,GAAG,QAAQ,CAAC,aAAa,CAAC,CAAC;QACxC,UAAU,EAAE,CAAC;QACb,IAAI,UAAU,GAAG,KAAK,EAAE,CAAC;YACvB,yDAAyD;YACzD,OAAO,SAAS,CAAC;QACnB,CAAC;IACH,CAAC;IACD,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,GAAG,EAAE,CAAC,EAAE,EAAE,CAAC;QAC7B,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,CAAC;YACD,SAAS,GAAG,CAAC,CAAC;YACd,MAAM;QACR,CAAC;IACH,CAAC;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,oDAAoD;IACrG,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,SAAS,CAAC,MAAM,GAAG,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC;QAC9C,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;IAC5E,CAAC;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,CAAE,SAAS,CAAC,EAAE,CAAC,CAAC,CAAC,CAAY,GAAI,SAAS,CAAC,EAAE,CAAC,CAAC,CAAC,CAAY,CAAC,GAAG,CAAC,CAAC;IACpE,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,CAAC;QACjC,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,CAAC;YACvE,KAAK,EAAE,CAAC;YACR,IAAI,KAAK,GAAG,CAAC,EAAE,CAAC;gBACd,OAAO,KAAK,CAAC;YACf,CAAC;QACH,CAAC;IACH,CAAC;IACD,IAAI,KAAK,KAAK,CAAC,EAAE,CAAC;QAChB,SAAS,GAAG,IAAI,CAAC;IACnB,CAAC;IACD,OAAO,SAAS,CAAC;AACnB,CAAC"}

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

@@ -0,0 +1,8 @@
+/**
+ * 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/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;;;;;GAKG;AACH,MAAM,CAAC,OAAO,UAAU,OAAO,CAAC,SAAS,EAAE,WAAW,EAAE,KAAK,EAAE,MAAM,GAAG,MAAM,CA2C7E"}

package/lib/operations/thresholds/moments.js

@@ -0,0 +1,58 @@
+// see https://github.com/fiji/Auto_Threshold/blob/master/src/main/java/fiji/threshold/Auto_Threshold.java
+// W. Tsai, "Moment-preserving thresholding: a new approach," Computer Vision,
+// Graphics, and Image Processing, vol. 29, pp. 377-393, 1985.
+// Ported to ImageJ plugin by G.Landini from the the open source project FOURIER 0.8
+// by M. Emre Celebi , Department of Computer Science, Louisiana State University in Shreveport
+// Shreveport, LA 71115, USA
+// http://sourceforge.net/projects/fourier-ipal
+// http://www.lsus.edu/faculty/~ecelebi/fourier.htm
+/**
+ * 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, total) {
+    // moments
+    const m0 = 1;
+    let m1 = 0;
+    let m2 = 0;
+    let m3 = 0;
+    let sum = 0;
+    let threshold = -1;
+    const histogramLength = histogram.length;
+    const normalizedHistogram = [];
+    for (let i = 0; i < histogramLength; i++) {
+        normalizedHistogram.push(histogram[i] / total);
+    }
+    /* Calculate the first, second, and third order moments */
+    for (let i = 0; i < histogramLength; i++) {
+        m1 += i * normalizedHistogram[i];
+        m2 += i * i * normalizedHistogram[i];
+        m3 += i * i * i * normalizedHistogram[i];
+    }
+    /*
+       First 4 moments of the gray-level image should match the first 4 moments
+       of the target binary image. This leads to 4 equalities whose solutions
+       are given in the Appendix of Ref. 1
+       */
+    const cd = m0 * m2 - m1 * m1; // determinant of the matriz of hankel for moments 2x2
+    const c0 = (-m2 * m2 + m1 * m3) / cd;
+    const c1 = (m0 * -m3 + m2 * m1) / cd;
+    // new two gray values where z0<z1
+    const z0 = 0.5 * (-c1 - Math.sqrt(c1 * c1 - 4 * c0));
+    const z1 = 0.5 * (-c1 + Math.sqrt(c1 * c1 - 4 * c0));
+    const p0 = (z1 - m1) /
+        (z1 -
+            z0); /* Fraction of the object pixels in the target binary image (p0z0+p1z1=m1) */
+    // The threshold is the gray-level closest to the p0-tile of the normalized histogram
+    for (let i = 0; i < histogramLength; i++) {
+        sum += normalizedHistogram[i];
+        if (sum > p0) {
+            threshold = i;
+            break;
+        }
+    }
+    return threshold;
+}
+//# sourceMappingURL=moments.js.map

package/lib/operations/thresholds/moments.js.map

@@ -0,0 +1 @@
+{"version":3,"file":"moments.js","sourceRoot":"","sources":["../../../src/operations/thresholds/moments.ts"],"names":[],"mappings":"AAAA,0GAA0G;AAC1G,8EAA8E;AAC9E,8DAA8D;AAC9D,oFAAoF;AACpF,+FAA+F;AAC/F,4BAA4B;AAC5B,+CAA+C;AAC/C,mDAAmD;AAEnD;;;;;GAKG;AACH,MAAM,CAAC,OAAO,UAAU,OAAO,CAAC,SAAsB,EAAE,KAAa;IACnE,UAAU;IACV,MAAM,EAAE,GAAG,CAAC,CAAC;IACb,IAAI,EAAE,GAAG,CAAC,CAAC;IACX,IAAI,EAAE,GAAG,CAAC,CAAC;IACX,IAAI,EAAE,GAAG,CAAC,CAAC;IACX,IAAI,GAAG,GAAG,CAAC,CAAC;IACZ,IAAI,SAAS,GAAG,CAAC,CAAC,CAAC;IACnB,MAAM,eAAe,GAAG,SAAS,CAAC,MAAM,CAAC;IACzC,MAAM,mBAAmB,GAAG,EAAE,CAAC;IAC/B,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,eAAe,EAAE,CAAC,EAAE,EAAE,CAAC;QACzC,mBAAmB,CAAC,IAAI,CAAC,SAAS,CAAC,CAAC,CAAC,GAAG,KAAK,CAAC,CAAC;IACjD,CAAC;IACD,0DAA0D;IAC1D,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,eAAe,EAAE,CAAC,EAAE,EAAE,CAAC;QACzC,EAAE,IAAI,CAAC,GAAG,mBAAmB,CAAC,CAAC,CAAC,CAAC;QACjC,EAAE,IAAI,CAAC,GAAG,CAAC,GAAG,mBAAmB,CAAC,CAAC,CAAC,CAAC;QACrC,EAAE,IAAI,CAAC,GAAG,CAAC,GAAG,CAAC,GAAG,mBAAmB,CAAC,CAAC,CAAC,CAAC;IAC3C,CAAC;IACD;;;;SAIK;IACL,MAAM,EAAE,GAAG,EAAE,GAAG,EAAE,GAAG,EAAE,GAAG,EAAE,CAAC,CAAC,sDAAsD;IACpF,MAAM,EAAE,GAAG,CAAC,CAAC,EAAE,GAAG,EAAE,GAAG,EAAE,GAAG,EAAE,CAAC,GAAG,EAAE,CAAC;IACrC,MAAM,EAAE,GAAG,CAAC,EAAE,GAAG,CAAC,EAAE,GAAG,EAAE,GAAG,EAAE,CAAC,GAAG,EAAE,CAAC;IACrC,kCAAkC;IAClC,MAAM,EAAE,GAAG,GAAG,GAAG,CAAC,CAAC,EAAE,GAAG,IAAI,CAAC,IAAI,CAAC,EAAE,GAAG,EAAE,GAAG,CAAC,GAAG,EAAE,CAAC,CAAC,CAAC;IACrD,MAAM,EAAE,GAAG,GAAG,GAAG,CAAC,CAAC,EAAE,GAAG,IAAI,CAAC,IAAI,CAAC,EAAE,GAAG,EAAE,GAAG,CAAC,GAAG,EAAE,CAAC,CAAC,CAAC;IACrD,MAAM,EAAE,GACN,CAAC,EAAE,GAAG,EAAE,CAAC;QACT,CAAC,EAAE;YACD,EAAE,CAAC,CAAC,CAAC,6EAA6E;IACtF,qFAAqF;IACrF,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,eAAe,EAAE,CAAC,EAAE,EAAE,CAAC;QACzC,GAAG,IAAI,mBAAmB,CAAC,CAAC,CAAC,CAAC;QAC9B,IAAI,GAAG,GAAG,EAAE,EAAE,CAAC;YACb,SAAS,GAAG,CAAC,CAAC;YACd,MAAM;QACR,CAAC;IACH,CAAC;IACD,OAAO,SAAS,CAAC;AACnB,CAAC"}

package/lib/operations/thresholds/otsu.d.ts

@@ -1,2 +1,8 @@
-export declare function otsu(histogramCounts: number[], total: number): number;
+/**
+ * Return a threshold for a histogram using Otsu algorithm.
+ * @param histogramCounts - The image histogram.
+ * @param total - Total number of pixels of the image.
+ * @returns The threshold.
+ */
+export declare function otsu(histogramCounts: Uint32Array, total: number): number;
 //# sourceMappingURL=otsu.d.ts.map

package/lib/operations/thresholds/otsu.d.ts.map

@@ -1 +1 @@
-{"version":3,"file":"otsu.d.ts","sourceRoot":"","sources":["../../../src/operations/thresholds/otsu.ts"],"names":[],"mappings":"AAOA,wBAAgB,IAAI,CAAC,eAAe,EAAE,MAAM,EAAE,EAAE,KAAK,EAAE,MAAM,GAAG,MAAM,CA2BrE"}
+{"version":3,"file":"otsu.d.ts","sourceRoot":"","sources":["../../../src/operations/thresholds/otsu.ts"],"names":[],"mappings":"AAOA;;;;;GAKG;AACH,wBAAgB,IAAI,CAAC,eAAe,EAAE,WAAW,EAAE,KAAK,EAAE,MAAM,GAAG,MAAM,CA2BxE"}

package/lib/operations/thresholds/otsu.js

@@ -1,35 +1,38 @@
-"use strict";
-/*
- * The method is present in: Otsu, N (1979), "A threshold selection method from gray-level histograms", IEEE Trans. Sys., Man., Cyber. 9: 62-66
- * The Otsu implementation is based on: https://en.wikipedia.org/wiki/Otsu's_method
- * @param histogram - the histogram of the image
- * @returns {number} - the threshold
- */
-Object.defineProperty(exports, "__esModule", { value: true });
-function otsu(histogramCounts, total) {
-    let sumB = 0;
-    let wB = 0;
-    let maximum = 0;
-    let level = 0;
-    let sum1 = 0;
-    for (let i = 0; i < histogramCounts.length; i++) {
-        sum1 += i * histogramCounts[i];
-    }
-    for (let ii = 0; ii < histogramCounts.length; ii++) {
-        wB = wB + histogramCounts[ii];
-        const wF = total - wB;
-        if (wB === 0 || wF === 0) {
-            continue;
-        }
-        sumB = sumB + ii * histogramCounts[ii];
-        const mF = (sum1 - sumB) / wF;
-        const between = wB * wF * (sumB / wB - mF) * (sumB / wB - mF);
-        if (between >= maximum) {
-            level = ii;
-            maximum = between;
-        }
-    }
-    return level;
-}
-exports.otsu = otsu;
+/*
+ * The method is present in: Otsu, N (1979), "A threshold selection method from gray-level histograms", IEEE Trans. Sys., Man., Cyber. 9: 62-66
+ * The Otsu implementation is based on: https://en.wikipedia.org/wiki/Otsu's_method
+ * @param histogram - the histogram of the image
+ * @returns {number} - the threshold
+ */
+/**
+ * Return a threshold for a histogram using Otsu algorithm.
+ * @param histogramCounts - The image histogram.
+ * @param total - Total number of pixels of the image.
+ * @returns The threshold.
+ */
+export function otsu(histogramCounts, total) {
+    let sumB = 0;
+    let wB = 0;
+    let maximum = 0;
+    let level = 0;
+    let sum1 = 0;
+    for (let i = 0; i < histogramCounts.length; i++) {
+        sum1 += i * histogramCounts[i];
+    }
+    for (let ii = 0; ii < histogramCounts.length; ii++) {
+        wB = wB + histogramCounts[ii];
+        const wF = total - wB;
+        if (wB === 0 || wF === 0) {
+            continue;
+        }
+        sumB = sumB + ii * histogramCounts[ii];
+        const mF = (sum1 - sumB) / wF;
+        const between = wB * wF * (sumB / wB - mF) * (sumB / wB - mF);
+        if (between >= maximum) {
+            level = ii;
+            maximum = between;
+        }
+    }
+    return level;
+}
 //# sourceMappingURL=otsu.js.map

package/lib/operations/thresholds/otsu.js.map

@@ -1 +1 @@
-{"version":3,"file":"otsu.js","sourceRoot":"","sources":["../../../src/operations/thresholds/otsu.ts"],"names":[],"mappings":";AAAA;;;;;GAKG;;AAEH,SAAgB,IAAI,CAAC,eAAyB,EAAE,KAAa;IAC3D,IAAI,IAAI,GAAG,CAAC,CAAC;IACb,IAAI,EAAE,GAAG,CAAC,CAAC;IACX,IAAI,OAAO,GAAG,CAAC,CAAC;IAChB,IAAI,KAAK,GAAG,CAAC,CAAC;IAEd,IAAI,IAAI,GAAG,CAAC,CAAC;IACb,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,eAAe,CAAC,MAAM,EAAE,CAAC,EAAE,EAAE;QAC/C,IAAI,IAAI,CAAC,GAAG,eAAe,CAAC,CAAC,CAAC,CAAC;KAChC;IAED,KAAK,IAAI,EAAE,GAAG,CAAC,EAAE,EAAE,GAAG,eAAe,CAAC,MAAM,EAAE,EAAE,EAAE,EAAE;QAClD,EAAE,GAAG,EAAE,GAAG,eAAe,CAAC,EAAE,CAAC,CAAC;QAC9B,MAAM,EAAE,GAAG,KAAK,GAAG,EAAE,CAAC;QACtB,IAAI,EAAE,KAAK,CAAC,IAAI,EAAE,KAAK,CAAC,EAAE;YACxB,SAAS;SACV;QACD,IAAI,GAAG,IAAI,GAAG,EAAE,GAAG,eAAe,CAAC,EAAE,CAAC,CAAC;QACvC,MAAM,EAAE,GAAG,CAAC,IAAI,GAAG,IAAI,CAAC,GAAG,EAAE,CAAC;QAC9B,MAAM,OAAO,GAAG,EAAE,GAAG,EAAE,GAAG,CAAC,IAAI,GAAG,EAAE,GAAG,EAAE,CAAC,GAAG,CAAC,IAAI,GAAG,EAAE,GAAG,EAAE,CAAC,CAAC;QAC9D,IAAI,OAAO,IAAI,OAAO,EAAE;YACtB,KAAK,GAAG,EAAE,CAAC;YACX,OAAO,GAAG,OAAO,CAAC;SACnB;KACF;IAED,OAAO,KAAK,CAAC;AACf,CAAC;AA3BD,oBA2BC"}
+{"version":3,"file":"otsu.js","sourceRoot":"","sources":["../../../src/operations/thresholds/otsu.ts"],"names":[],"mappings":"AAAA;;;;;GAKG;AAEH;;;;;GAKG;AACH,MAAM,UAAU,IAAI,CAAC,eAA4B,EAAE,KAAa;IAC9D,IAAI,IAAI,GAAG,CAAC,CAAC;IACb,IAAI,EAAE,GAAG,CAAC,CAAC;IACX,IAAI,OAAO,GAAG,CAAC,CAAC;IAChB,IAAI,KAAK,GAAG,CAAC,CAAC;IAEd,IAAI,IAAI,GAAG,CAAC,CAAC;IACb,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,eAAe,CAAC,MAAM,EAAE,CAAC,EAAE,EAAE,CAAC;QAChD,IAAI,IAAI,CAAC,GAAG,eAAe,CAAC,CAAC,CAAC,CAAC;IACjC,CAAC;IAED,KAAK,IAAI,EAAE,GAAG,CAAC,EAAE,EAAE,GAAG,eAAe,CAAC,MAAM,EAAE,EAAE,EAAE,EAAE,CAAC;QACnD,EAAE,GAAG,EAAE,GAAG,eAAe,CAAC,EAAE,CAAC,CAAC;QAC9B,MAAM,EAAE,GAAG,KAAK,GAAG,EAAE,CAAC;QACtB,IAAI,EAAE,KAAK,CAAC,IAAI,EAAE,KAAK,CAAC,EAAE,CAAC;YACzB,SAAS;QACX,CAAC;QACD,IAAI,GAAG,IAAI,GAAG,EAAE,GAAG,eAAe,CAAC,EAAE,CAAC,CAAC;QACvC,MAAM,EAAE,GAAG,CAAC,IAAI,GAAG,IAAI,CAAC,GAAG,EAAE,CAAC;QAC9B,MAAM,OAAO,GAAG,EAAE,GAAG,EAAE,GAAG,CAAC,IAAI,GAAG,EAAE,GAAG,EAAE,CAAC,GAAG,CAAC,IAAI,GAAG,EAAE,GAAG,EAAE,CAAC,CAAC;QAC9D,IAAI,OAAO,IAAI,OAAO,EAAE,CAAC;YACvB,KAAK,GAAG,EAAE,CAAC;YACX,OAAO,GAAG,OAAO,CAAC;QACpB,CAAC;IACH,CAAC;IAED,OAAO,KAAK,CAAC;AACf,CAAC"}

package/lib/operations/thresholds/percentile.d.ts

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

package/lib/operations/thresholds/percentile.d.ts.map

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

package/lib/operations/thresholds/percentile.js

@@ -0,0 +1,35 @@
+// See http://imagej.nih.gov/ij/download/tools/source/ij/process/AutoThresholder.java
+// W. Doyle, "Operation useful for similarity-invariant pattern recognition,"
+// Journal of the Association for Computing Machinery, vol. 9,pp. 259-267, 1962.
+// 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 using percentiles.
+ * @param histogram - Image histogram.
+ * @returns The threshold.
+ */
+export default function percentile(histogram) {
+    let threshold = -1;
+    const percentile = 0.5; // default fraction of foreground pixels
+    const avec = new Array(histogram.length);
+    const total = partialSum(histogram, histogram.length - 1);
+    let temp = 1;
+    for (let i = 0; i < histogram.length; i++) {
+        avec[i] = Math.abs(partialSum(histogram, i) / total - percentile);
+        if (avec[i] < temp) {
+            temp = avec[i];
+            threshold = i;
+        }
+    }
+    return threshold;
+}
+function partialSum(histogram, endIndex) {
+    let x = 0;
+    for (let i = 0; i <= endIndex; i++) {
+        x += histogram[i];
+    }
+    return x;
+}
+//# sourceMappingURL=percentile.js.map

package/lib/operations/thresholds/percentile.js.map

@@ -0,0 +1 @@
+{"version":3,"file":"percentile.js","sourceRoot":"","sources":["../../../src/operations/thresholds/percentile.ts"],"names":[],"mappings":"AAAA,qFAAqF;AACrF,6EAA6E;AAC7E,gFAAgF;AAChF,+EAA+E;AAC/E,yDAAyD;AACzD,gFAAgF;AAChF,gCAAgC;AAEhC;;;;GAIG;AACH,MAAM,CAAC,OAAO,UAAU,UAAU,CAAC,SAAsB;IACvD,IAAI,SAAS,GAAG,CAAC,CAAC,CAAC;IACnB,MAAM,UAAU,GAAG,GAAG,CAAC,CAAC,wCAAwC;IAChE,MAAM,IAAI,GAAG,IAAI,KAAK,CAAC,SAAS,CAAC,MAAM,CAAC,CAAC;IAEzC,MAAM,KAAK,GAAG,UAAU,CAAC,SAAS,EAAE,SAAS,CAAC,MAAM,GAAG,CAAC,CAAC,CAAC;IAC1D,IAAI,IAAI,GAAG,CAAC,CAAC;IAEb,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,SAAS,CAAC,MAAM,EAAE,CAAC,EAAE,EAAE,CAAC;QAC1C,IAAI,CAAC,CAAC,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,UAAU,CAAC,SAAS,EAAE,CAAC,CAAC,GAAG,KAAK,GAAG,UAAU,CAAC,CAAC;QAClE,IAAI,IAAI,CAAC,CAAC,CAAC,GAAG,IAAI,EAAE,CAAC;YACnB,IAAI,GAAG,IAAI,CAAC,CAAC,CAAC,CAAC;YACf,SAAS,GAAG,CAAC,CAAC;QAChB,CAAC;IACH,CAAC;IAED,OAAO,SAAS,CAAC;AACnB,CAAC;AAED,SAAS,UAAU,CAAC,SAAsB,EAAE,QAAgB;IAC1D,IAAI,CAAC,GAAG,CAAC,CAAC;IACV,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,IAAI,QAAQ,EAAE,CAAC,EAAE,EAAE,CAAC;QACnC,CAAC,IAAI,SAAS,CAAC,CAAC,CAAC,CAAC;IACpB,CAAC;IACD,OAAO,CAAC,CAAC;AACX,CAAC"}

package/lib/operations/thresholds/renyiEntropy.d.ts

@@ -0,0 +1,8 @@
+/**
+ * Return a threshold for a histogram using its entropy.
+ * @param histogram - The image histogram.
+ * @param total - Total number of pixels of the image.
+ * @returns The threshold.
+ */
+export default function renyiEntropy(histogram: Uint32Array, total: number): number;
+//# sourceMappingURL=renyiEntropy.d.ts.map

package/lib/operations/thresholds/renyiEntropy.d.ts.map

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