image-js 0.0.0-next-e30968a2d8

package/lib/operations/threshold.js

@@ -0,0 +1,114 @@
+"use strict";
+var __importDefault = (this && this.__importDefault) || function (mod) {
+    return (mod && mod.__esModule) ? mod : { "default": mod };
+};
+Object.defineProperty(exports, "__esModule", { value: true });
+exports.threshold = exports.computeThreshold = exports.ThresholdAlgorithm = void 0;
+const getOutputImage_1 = require("../utils/getOutputImage");
+const validators_1 = require("../utils/validators");
+const huang_1 = __importDefault(require("./thresholds/huang"));
+const intermodes_1 = __importDefault(require("./thresholds/intermodes"));
+const isodata_1 = __importDefault(require("./thresholds/isodata"));
+const li_1 = __importDefault(require("./thresholds/li"));
+const maxEntropy_1 = __importDefault(require("./thresholds/maxEntropy"));
+const mean_1 = __importDefault(require("./thresholds/mean"));
+const minError_1 = __importDefault(require("./thresholds/minError"));
+const minimum_1 = __importDefault(require("./thresholds/minimum"));
+const moments_1 = __importDefault(require("./thresholds/moments"));
+const otsu_1 = require("./thresholds/otsu");
+const percentile_1 = __importDefault(require("./thresholds/percentile"));
+const renyiEntropy_1 = __importDefault(require("./thresholds/renyiEntropy"));
+const shanbhag_1 = __importDefault(require("./thresholds/shanbhag"));
+const triangle_1 = require("./thresholds/triangle");
+const yen_1 = __importDefault(require("./thresholds/yen"));
+var ThresholdAlgorithm;
+(function (ThresholdAlgorithm) {
+    ThresholdAlgorithm["HUANG"] = "HUANG";
+    ThresholdAlgorithm["INTERMODES"] = "INTERMODES";
+    ThresholdAlgorithm["ISODATA"] = "ISODATA";
+    ThresholdAlgorithm["LI"] = "LI";
+    ThresholdAlgorithm["MAX_ENTROPY"] = "MAX_ENTROPY";
+    ThresholdAlgorithm["MEAN"] = "MEAN";
+    ThresholdAlgorithm["MIN_ERROR"] = "MIN_ERROR";
+    ThresholdAlgorithm["MINIMUM"] = "MINIMUM";
+    ThresholdAlgorithm["MOMENTS"] = "MOMENTS";
+    ThresholdAlgorithm["OTSU"] = "OTSU";
+    ThresholdAlgorithm["PERCENTILE"] = "PERCENTILE";
+    ThresholdAlgorithm["RENYI_ENTROPY"] = "RENYI_ENTROPY";
+    ThresholdAlgorithm["SHANBHAG"] = "SHANBHAG";
+    ThresholdAlgorithm["TRIANGLE"] = "TRIANGLE";
+    ThresholdAlgorithm["YEN"] = "YEN";
+})(ThresholdAlgorithm = exports.ThresholdAlgorithm || (exports.ThresholdAlgorithm = {}));
+/**
+ * Compute threshold value for an image using the specified algorithm.
+ *
+ * @param image - The grey image.
+ * @param algorithm - Algorithm that defines the threshold.
+ * @returns The threshold value for the image.
+ */
+function computeThreshold(image, algorithm = ThresholdAlgorithm.OTSU) {
+    if (image.channels !== 1) {
+        throw new Error('threshold can only be computed on images with one channel');
+    }
+    const histogram = image.histogram();
+    switch (algorithm) {
+        case ThresholdAlgorithm.HUANG:
+            return (0, huang_1.default)(histogram);
+        case ThresholdAlgorithm.INTERMODES:
+            return (0, intermodes_1.default)(histogram);
+        case ThresholdAlgorithm.ISODATA:
+            return (0, isodata_1.default)(histogram);
+        case ThresholdAlgorithm.LI:
+            return (0, li_1.default)(histogram, image.size);
+        case ThresholdAlgorithm.MAX_ENTROPY:
+            return (0, maxEntropy_1.default)(histogram, image.size);
+        case ThresholdAlgorithm.MEAN:
+            return (0, mean_1.default)(histogram, image.size);
+        case ThresholdAlgorithm.MINIMUM:
+            return (0, minimum_1.default)(histogram);
+        case ThresholdAlgorithm.MIN_ERROR:
+            return (0, minError_1.default)(histogram, image.size);
+        case ThresholdAlgorithm.MOMENTS:
+            return (0, moments_1.default)(histogram, image.size);
+        case ThresholdAlgorithm.OTSU:
+            return (0, otsu_1.otsu)(histogram, image.size);
+        case ThresholdAlgorithm.PERCENTILE:
+            return (0, percentile_1.default)(histogram);
+        case ThresholdAlgorithm.RENYI_ENTROPY:
+            return (0, renyiEntropy_1.default)(histogram, image.size);
+        case ThresholdAlgorithm.SHANBHAG:
+            return (0, shanbhag_1.default)(histogram, image.size);
+        case ThresholdAlgorithm.TRIANGLE:
+            return (0, triangle_1.triangle)(histogram);
+        case ThresholdAlgorithm.YEN:
+            return (0, yen_1.default)(histogram, image.size);
+        default:
+            throw new RangeError(`unsupported threshold algorithm: ${algorithm}`);
+    }
+}
+exports.computeThreshold = computeThreshold;
+// See: https://docs.opencv.org/4.0.1/d7/d1b/group__imgproc__misc.html#gaa9e58d2860d4afa658ef70a9b1115576
+/**
+ * Create a black and white image based on a threshold value.
+ *
+ * @param image - The grey image to convert.
+ * @param options - Threshold options.
+ * @returns The resulting mask.
+ */
+function threshold(image, options = {}) {
+    let thresholdValue;
+    if ('threshold' in options) {
+        thresholdValue = options.threshold;
+    }
+    else {
+        thresholdValue = computeThreshold(image, options.algorithm);
+    }
+    (0, validators_1.validateValue)(thresholdValue, image);
+    const result = (0, getOutputImage_1.imageToOutputMask)(image, options);
+    for (let i = 0; i < image.size; i++) {
+        result.setBitByIndex(i, image.getValueByIndex(i, 0) > thresholdValue ? 1 : 0);
+    }
+    return result;
+}
+exports.threshold = threshold;
+//# sourceMappingURL=threshold.js.map

package/lib/operations/threshold.js.map

@@ -0,0 +1 @@
+{"version":3,"file":"threshold.js","sourceRoot":"","sources":["../../src/operations/threshold.ts"],"names":[],"mappings":";;;;;;AAEA,4DAA4D;AAC5D,oDAAoD;AAEpD,+DAAuC;AACvC,yEAAiD;AACjD,mEAA2C;AAC3C,yDAAiC;AACjC,yEAAiD;AACjD,6DAAqC;AACrC,qEAA6C;AAC7C,mEAA2C;AAC3C,mEAA2C;AAC3C,4CAAyC;AACzC,yEAAiD;AACjD,6EAAqD;AACrD,qEAA6C;AAC7C,oDAAiD;AACjD,2DAAmC;AAEnC,IAAY,kBAgBX;AAhBD,WAAY,kBAAkB;IAC5B,qCAAe,CAAA;IACf,+CAAyB,CAAA;IACzB,yCAAmB,CAAA;IACnB,+BAAS,CAAA;IACT,iDAA2B,CAAA;IAC3B,mCAAa,CAAA;IACb,6CAAuB,CAAA;IACvB,yCAAmB,CAAA;IACnB,yCAAmB,CAAA;IACnB,mCAAa,CAAA;IACb,+CAAyB,CAAA;IACzB,qDAA+B,CAAA;IAC/B,2CAAqB,CAAA;IACrB,2CAAqB,CAAA;IACrB,iCAAW,CAAA;AACb,CAAC,EAhBW,kBAAkB,GAAlB,0BAAkB,KAAlB,0BAAkB,QAgB7B;AA6BD;;;;;;GAMG;AACH,SAAgB,gBAAgB,CAC9B,KAAU,EACV,YAAgC,kBAAkB,CAAC,IAAI;IAEvD,IAAI,KAAK,CAAC,QAAQ,KAAK,CAAC,EAAE;QACxB,MAAM,IAAI,KAAK,CACb,2DAA2D,CAC5D,CAAC;KACH;IACD,MAAM,SAAS,GAAG,KAAK,CAAC,SAAS,EAAE,CAAC;IAEpC,QAAQ,SAAS,EAAE;QACjB,KAAK,kBAAkB,CAAC,KAAK;YAC3B,OAAO,IAAA,eAAK,EAAC,SAAS,CAAC,CAAC;QAC1B,KAAK,kBAAkB,CAAC,UAAU;YAChC,OAAO,IAAA,oBAAU,EAAC,SAAS,CAAC,CAAC;QAC/B,KAAK,kBAAkB,CAAC,OAAO;YAC7B,OAAO,IAAA,iBAAO,EAAC,SAAS,CAAC,CAAC;QAC5B,KAAK,kBAAkB,CAAC,EAAE;YACxB,OAAO,IAAA,YAAE,EAAC,SAAS,EAAE,KAAK,CAAC,IAAI,CAAC,CAAC;QACnC,KAAK,kBAAkB,CAAC,WAAW;YACjC,OAAO,IAAA,oBAAU,EAAC,SAAS,EAAE,KAAK,CAAC,IAAI,CAAC,CAAC;QAC3C,KAAK,kBAAkB,CAAC,IAAI;YAC1B,OAAO,IAAA,cAAI,EAAC,SAAS,EAAE,KAAK,CAAC,IAAI,CAAC,CAAC;QACrC,KAAK,kBAAkB,CAAC,OAAO;YAC7B,OAAO,IAAA,iBAAO,EAAC,SAAS,CAAC,CAAC;QAC5B,KAAK,kBAAkB,CAAC,SAAS;YAC/B,OAAO,IAAA,kBAAQ,EAAC,SAAS,EAAE,KAAK,CAAC,IAAI,CAAC,CAAC;QACzC,KAAK,kBAAkB,CAAC,OAAO;YAC7B,OAAO,IAAA,iBAAO,EAAC,SAAS,EAAE,KAAK,CAAC,IAAI,CAAC,CAAC;QACxC,KAAK,kBAAkB,CAAC,IAAI;YAC1B,OAAO,IAAA,WAAI,EAAC,SAAS,EAAE,KAAK,CAAC,IAAI,CAAC,CAAC;QACrC,KAAK,kBAAkB,CAAC,UAAU;YAChC,OAAO,IAAA,oBAAU,EAAC,SAAS,CAAC,CAAC;QAC/B,KAAK,kBAAkB,CAAC,aAAa;YACnC,OAAO,IAAA,sBAAY,EAAC,SAAS,EAAE,KAAK,CAAC,IAAI,CAAC,CAAC;QAC7C,KAAK,kBAAkB,CAAC,QAAQ;YAC9B,OAAO,IAAA,kBAAQ,EAAC,SAAS,EAAE,KAAK,CAAC,IAAI,CAAC,CAAC;QACzC,KAAK,kBAAkB,CAAC,QAAQ;YAC9B,OAAO,IAAA,mBAAQ,EAAC,SAAS,CAAC,CAAC;QAC7B,KAAK,kBAAkB,CAAC,GAAG;YACzB,OAAO,IAAA,aAAG,EAAC,SAAS,EAAE,KAAK,CAAC,IAAI,CAAC,CAAC;QACpC;YACE,MAAM,IAAI,UAAU,CAAC,oCAAoC,SAAS,EAAE,CAAC,CAAC;KACzE;AACH,CAAC;AA7CD,4CA6CC;AAED,yGAAyG;AACzG;;;;;;GAMG;AACH,SAAgB,SAAS,CAAC,KAAU,EAAE,UAA4B,EAAE;IAClE,IAAI,cAAsB,CAAC;IAC3B,IAAI,WAAW,IAAI,OAAO,EAAE;QAC1B,cAAc,GAAG,OAAO,CAAC,SAAS,CAAC;KACpC;SAAM;QACL,cAAc,GAAG,gBAAgB,CAAC,KAAK,EAAE,OAAO,CAAC,SAAS,CAAC,CAAC;KAC7D;IACD,IAAA,0BAAa,EAAC,cAAc,EAAE,KAAK,CAAC,CAAC;IACrC,MAAM,MAAM,GAAG,IAAA,kCAAiB,EAAC,KAAK,EAAE,OAAO,CAAC,CAAC;IACjD,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,KAAK,CAAC,IAAI,EAAE,CAAC,EAAE,EAAE;QACnC,MAAM,CAAC,aAAa,CAClB,CAAC,EACD,KAAK,CAAC,eAAe,CAAC,CAAC,EAAE,CAAC,CAAC,GAAG,cAAc,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CACrD,CAAC;KACH;IACD,OAAO,MAAM,CAAC;AAChB,CAAC;AAhBD,8BAgBC"}

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

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

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

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

package/lib/operations/thresholds/huang.js

@@ -0,0 +1,81 @@
+"use strict";
+/*
+ * * see http://rsb.info.nih.gov/ij/developer/source/ij/process/AutoThresholder.java.html.
+ * Huang: Implements Huang's fuzzy thresholding method: Huang, L-K & Wang, M-J J (1995),
+ * "Image thresholding by minimizing the measure of fuzziness", Pattern Recognition 28(1): 41-51
+ */
+Object.defineProperty(exports, "__esModule", { value: true });
+/**
+ * Return a threshold for a histogram using Huang algorithm.
+ *
+ * @param histogram - Image histogram.
+ * @returns The threshold.
+ */
+function huang(histogram) {
+    /* Determine the first non-zero bin */
+    let firstBin = 0;
+    for (let ih = 0; ih < histogram.length; ih++) {
+        if (histogram[ih] !== 0) {
+            firstBin = ih;
+            break;
+        }
+    }
+    /* Determine the last non-zero bin */
+    let lastBin = histogram.length - 1;
+    for (let ih = histogram.length - 1; ih >= firstBin; ih--) {
+        if (histogram[ih] !== 0) {
+            lastBin = ih;
+            break;
+        }
+    }
+    const term = 1.0 / (lastBin - firstBin);
+    const mu0 = new Array(histogram.length);
+    let sumPix = 0;
+    let numPix = 0;
+    for (let ih = firstBin; ih < histogram.length; ih++) {
+        sumPix += ih * histogram[ih];
+        numPix += histogram[ih];
+        mu0[ih] = sumPix / numPix;
+    }
+    const mu1 = new Array(histogram.length);
+    sumPix = numPix = 0;
+    for (let ih = lastBin; ih > 0; ih--) {
+        sumPix += ih * histogram[ih];
+        numPix += histogram[ih];
+        mu1[ih - 1] = sumPix / numPix;
+    }
+    /* Determine the threshold that minimizes the fuzzy entropy*/
+    let threshold = -1;
+    let minEnt = Number.MAX_VALUE;
+    for (let it = 0; it < histogram.length; it++) {
+        let ent = 0;
+        let muX;
+        for (let ih = 0; ih <= it; ih++) {
+            /* Equation (4) in Ref. 1 */
+            muX = 1 / (1 + term * Math.abs(ih - mu0[it]));
+            if (!(muX < 1e-6 || muX > 0.999999)) {
+                /* Equation (6) & (8) in Ref. 1 */
+                ent +=
+                    histogram[ih] *
+                        (-muX * Math.log(muX) - (1 - muX) * Math.log(1 - muX));
+            }
+        }
+        for (let ih = it + 1; ih < histogram.length; ih++) {
+            /* Equation (4) in Ref. 1 */
+            muX = 1 / (1 + term * Math.abs(ih - mu1[it]));
+            if (!(muX < 1e-6 || muX > 0.999999)) {
+                /* Equation (6) & (8) in Ref. 1 */
+                ent +=
+                    histogram[ih] *
+                        (-muX * Math.log(muX) - (1 - muX) * Math.log(1 - muX));
+            }
+        }
+        if (ent < minEnt) {
+            minEnt = ent;
+            threshold = it;
+        }
+    }
+    return threshold;
+}
+exports.default = huang;
+//# sourceMappingURL=huang.js.map

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

@@ -0,0 +1 @@
+{"version":3,"file":"huang.js","sourceRoot":"","sources":["../../../src/operations/thresholds/huang.ts"],"names":[],"mappings":";AAAA;;;;GAIG;;AAEH;;;;;GAKG;AACH,SAAwB,KAAK,CAAC,SAAsB;IAClD,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,SAAS,CAAC,EAAE,CAAC,KAAK,CAAC,EAAE;YACvB,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,SAAS,CAAC,EAAE,CAAC,KAAK,CAAC,EAAE;YACvB,OAAO,GAAG,EAAE,CAAC;YACb,MAAM;SACP;KACF;IAED,MAAM,IAAI,GAAG,GAAG,GAAG,CAAC,OAAO,GAAG,QAAQ,CAAC,CAAC;IACxC,MAAM,GAAG,GAAG,IAAI,KAAK,CAAC,SAAS,CAAC,MAAM,CAAC,CAAC;IACxC,IAAI,MAAM,GAAG,CAAC,CAAC;IACf,IAAI,MAAM,GAAG,CAAC,CAAC;IACf,KAAK,IAAI,EAAE,GAAG,QAAQ,EAAE,EAAE,GAAG,SAAS,CAAC,MAAM,EAAE,EAAE,EAAE,EAAE;QACnD,MAAM,IAAI,EAAE,GAAG,SAAS,CAAC,EAAE,CAAC,CAAC;QAC7B,MAAM,IAAI,SAAS,CAAC,EAAE,CAAC,CAAC;QACxB,GAAG,CAAC,EAAE,CAAC,GAAG,MAAM,GAAG,MAAM,CAAC;KAC3B;IAED,MAAM,GAAG,GAAG,IAAI,KAAK,CAAC,SAAS,CAAC,MAAM,CAAC,CAAC;IACxC,MAAM,GAAG,MAAM,GAAG,CAAC,CAAC;IACpB,KAAK,IAAI,EAAE,GAAG,OAAO,EAAE,EAAE,GAAG,CAAC,EAAE,EAAE,EAAE,EAAE;QACnC,MAAM,IAAI,EAAE,GAAG,SAAS,CAAC,EAAE,CAAC,CAAC;QAC7B,MAAM,IAAI,SAAS,CAAC,EAAE,CAAC,CAAC;QACxB,GAAG,CAAC,EAAE,GAAG,CAAC,CAAC,GAAG,MAAM,GAAG,MAAM,CAAC;KAC/B;IAED,6DAA6D;IAC7D,IAAI,SAAS,GAAG,CAAC,CAAC,CAAC;IACnB,IAAI,MAAM,GAAG,MAAM,CAAC,SAAS,CAAC;IAC9B,KAAK,IAAI,EAAE,GAAG,CAAC,EAAE,EAAE,GAAG,SAAS,CAAC,MAAM,EAAE,EAAE,EAAE,EAAE;QAC5C,IAAI,GAAG,GAAG,CAAC,CAAC;QACZ,IAAI,GAAG,CAAC;QACR,KAAK,IAAI,EAAE,GAAG,CAAC,EAAE,EAAE,IAAI,EAAE,EAAE,EAAE,EAAE,EAAE;YAC/B,4BAA4B;YAC5B,GAAG,GAAG,CAAC,GAAG,CAAC,CAAC,GAAG,IAAI,GAAG,IAAI,CAAC,GAAG,CAAC,EAAE,GAAG,GAAG,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC;YAC9C,IAAI,CAAC,CAAC,GAAG,GAAG,IAAI,IAAI,GAAG,GAAG,QAAQ,CAAC,EAAE;gBACnC,kCAAkC;gBAClC,GAAG;oBACD,SAAS,CAAC,EAAE,CAAC;wBACb,CAAC,CAAC,GAAG,GAAG,IAAI,CAAC,GAAG,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC,GAAG,GAAG,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,CAAC,GAAG,GAAG,CAAC,CAAC,CAAC;aAC1D;SACF;QAED,KAAK,IAAI,EAAE,GAAG,EAAE,GAAG,CAAC,EAAE,EAAE,GAAG,SAAS,CAAC,MAAM,EAAE,EAAE,EAAE,EAAE;YACjD,4BAA4B;YAC5B,GAAG,GAAG,CAAC,GAAG,CAAC,CAAC,GAAG,IAAI,GAAG,IAAI,CAAC,GAAG,CAAC,EAAE,GAAG,GAAG,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC;YAC9C,IAAI,CAAC,CAAC,GAAG,GAAG,IAAI,IAAI,GAAG,GAAG,QAAQ,CAAC,EAAE;gBACnC,kCAAkC;gBAClC,GAAG;oBACD,SAAS,CAAC,EAAE,CAAC;wBACb,CAAC,CAAC,GAAG,GAAG,IAAI,CAAC,GAAG,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC,GAAG,GAAG,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,CAAC,GAAG,GAAG,CAAC,CAAC,CAAC;aAC1D;SACF;QAED,IAAI,GAAG,GAAG,MAAM,EAAE;YAChB,MAAM,GAAG,GAAG,CAAC;YACb,SAAS,GAAG,EAAE,CAAC;SAChB;KACF;IACD,OAAO,SAAS,CAAC;AACnB,CAAC;AAvED,wBAuEC"}

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

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

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

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

package/lib/operations/thresholds/intermodes.js

@@ -0,0 +1,62 @@
+"use strict";
+/*
+ *
+ * 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
+ *
+ */
+Object.defineProperty(exports, "__esModule", { value: true });
+/**
+ * Return a threshold for a histogram using Intermodes algorithm.
+ *
+ * @param histogram - Image histogram.
+ * @returns The threshold.
+ */
+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);
+}
+exports.default = intermodes;
+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/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,SAAwB,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;AA7BD,6BA6BC;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/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/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/operations/thresholds/isodata.js

@@ -0,0 +1,55 @@
+"use strict";
+Object.defineProperty(exports, "__esModule", { value: true });
+/*
+ * 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.
+ */
+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;
+}
+exports.default = isodata;
+//# sourceMappingURL=isodata.js.map

package/lib/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,SAAwB,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;AAxCD,0BAwCC"}

package/lib/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/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/operations/thresholds/li.js

@@ -0,0 +1,77 @@
+"use strict";
+/*
+ * 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
+ */
+Object.defineProperty(exports, "__esModule", { value: true });
+/**
+ * Return a threshold for a histogram using Li algorithm.
+ *
+ * @param histogram - Image histogram.
+ * @param total - Number of pixels in the image.
+ * @returns The threshold.
+ */
+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;
+}
+exports.default = li;
+//# sourceMappingURL=li.js.map

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;;;;;;GAMG;AACH,SAAwB,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;AA3DD,qBA2DC"}

package/lib/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/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/operations/thresholds/maxEntropy.js

@@ -0,0 +1,80 @@
+"use strict";
+/*
+ * 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
+ */
+Object.defineProperty(exports, "__esModule", { value: true });
+/**
+ * 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.
+ */
+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;
+}
+exports.default = maxEntropy;
+//# 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;;;;;;GAMG;AACH,SAAwB,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;AAvED,6BAuEC"}

package/lib/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/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/operations/thresholds/mean.js

@@ -0,0 +1,26 @@
+"use strict";
+/*
+ * 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
+ */
+Object.defineProperty(exports, "__esModule", { value: true });
+/**
+ * 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.
+ */
+function mean(histogram, total) {
+    let sum = 0;
+    for (let i = 0; i < histogram.length; i++) {
+        sum += i * histogram[i];
+    }
+    return Math.floor(sum / total);
+}
+exports.default = mean;
+//# 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;;;;;;GAMG;AACH,SAAwB,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;AAND,uBAMC"}

package/lib/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/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"}