image-js 1.0.0-alpha.3 → 1.0.0

package/lib/maskAnalysis/utils/getMbrFromPoints.js

@@ -0,0 +1,102 @@
+import { rotate } from '../../utils/geometry/points.js';
+import { getAngle } from './getAngle.js';
+import { getMbrAngle } from './getMbrAngle.js';
+/**
+ * Get the four corners of the minimum bounding rectangle from a set of points defining a simple convex polygon.
+ * @see {@link https://www.researchgate.net/profile/Lennert_Den_Boer2/publication/303783472_A_Fast_Algorithm_for_Generating_a_Minimal_Bounding_Rectangle/links/5751a14108ae6807fafb2aa5.pdf}
+ * @param points - Points from which to compute the MBR.
+ * @returns The array of corners.
+ */
+export function getMbrFromPoints(points) {
+    if (points.length === 0) {
+        return {
+            points: [],
+            angle: 0,
+            width: 0,
+            height: 0,
+            surface: 0,
+            perimeter: 0,
+            aspectRatio: 0,
+        };
+    }
+    if (points.length === 1) {
+        return {
+            points: [points[0], points[0], points[0], points[0]],
+            perimeter: 0,
+            surface: 0,
+            angle: 0,
+            width: 0,
+            height: 0,
+            aspectRatio: 1,
+        };
+    }
+    let rotatedVertices = [];
+    let minSurface = Number.POSITIVE_INFINITY;
+    let minSurfaceAngle = 0;
+    let mbr = [];
+    for (let i = 0; i < points.length; i++) {
+        const angle = getAngle(points[i], points[(i + 1) % points.length]);
+        rotatedVertices = rotate(-angle, points);
+        // Rotate and translate so that this segment is at the bottom.
+        const aX = rotatedVertices[i].column;
+        const aY = rotatedVertices[i].row;
+        const bX = rotatedVertices[(i + 1) % rotatedVertices.length].column;
+        const bY = rotatedVertices[(i + 1) % rotatedVertices.length].row;
+        let tUndefined = true;
+        let tMin = 0;
+        let tMax = 0;
+        let maxWidth = 0;
+        for (const point of rotatedVertices) {
+            const cX = point.column;
+            const cY = point.row;
+            const t = (cX - aX) / (bX - aX);
+            if (tUndefined) {
+                tUndefined = false;
+                tMin = t;
+                tMax = t;
+            }
+            else {
+                if (t < tMin)
+                    tMin = t;
+                if (t > tMax)
+                    tMax = t;
+            }
+            const width = (-(bX - aX) * cY + bX * aY - bY * aX) / (bX - aX);
+            if (Math.abs(width) > Math.abs(maxWidth)) {
+                maxWidth = width;
+            }
+        }
+        const minPoint = { column: aX + tMin * (bX - aX), row: aY };
+        const maxPoint = { column: aX + tMax * (bX - aX), row: aY };
+        const currentSurface = Math.abs(maxWidth * (tMin - tMax) * (bX - aX));
+        if (currentSurface < minSurface) {
+            minSurfaceAngle = angle;
+            minSurface = currentSurface;
+            mbr = [
+                maxPoint,
+                minPoint,
+                { column: minPoint.column, row: minPoint.row - maxWidth },
+                { column: maxPoint.column, row: maxPoint.row - maxWidth },
+            ];
+        }
+    }
+    const mbrRotated = rotate(minSurfaceAngle, mbr);
+    const sides = [
+        Math.hypot(mbr[0].column - mbr[1].column, mbr[0].row - mbr[1].row),
+        Math.hypot(mbr[0].column - mbr[3].column, mbr[0].row - mbr[3].row),
+    ];
+    const maxSide = Math.max(...sides);
+    const minSide = Math.min(...sides);
+    const mbrAngle = getMbrAngle(mbrRotated);
+    const ratio = minSide / maxSide;
+    return {
+        points: mbrRotated,
+        surface: minSurface,
+        angle: mbrAngle,
+        width: maxSide,
+        height: minSide,
+        perimeter: 2 * maxSide + 2 * minSide,
+        aspectRatio: ratio,
+    };
+}
+//# sourceMappingURL=getMbrFromPoints.js.map

package/lib/maskAnalysis/utils/getMbrFromPoints.js.map

@@ -0,0 +1 @@
+{"version":3,"file":"getMbrFromPoints.js","sourceRoot":"","sources":["../../../src/maskAnalysis/utils/getMbrFromPoints.ts"],"names":[],"mappings":"AACA,OAAO,EAAE,MAAM,EAAE,MAAM,gCAAgC,CAAC;AAGxD,OAAO,EAAE,QAAQ,EAAE,MAAM,eAAe,CAAC;AACzC,OAAO,EAAE,WAAW,EAAE,MAAM,kBAAkB,CAAC;AAE/C;;;;;GAKG;AACH,MAAM,UAAU,gBAAgB,CAAC,MAAwB;IACvD,IAAI,MAAM,CAAC,MAAM,KAAK,CAAC,EAAE,CAAC;QACxB,OAAO;YACL,MAAM,EAAE,EAAE;YACV,KAAK,EAAE,CAAC;YACR,KAAK,EAAE,CAAC;YACR,MAAM,EAAE,CAAC;YACT,OAAO,EAAE,CAAC;YACV,SAAS,EAAE,CAAC;YACZ,WAAW,EAAE,CAAC;SACf,CAAC;IACJ,CAAC;IACD,IAAI,MAAM,CAAC,MAAM,KAAK,CAAC,EAAE,CAAC;QACxB,OAAO;YACL,MAAM,EAAE,CAAC,MAAM,CAAC,CAAC,CAAC,EAAE,MAAM,CAAC,CAAC,CAAC,EAAE,MAAM,CAAC,CAAC,CAAC,EAAE,MAAM,CAAC,CAAC,CAAC,CAAC;YACpD,SAAS,EAAE,CAAC;YACZ,OAAO,EAAE,CAAC;YACV,KAAK,EAAE,CAAC;YACR,KAAK,EAAE,CAAC;YACR,MAAM,EAAE,CAAC;YACT,WAAW,EAAE,CAAC;SACf,CAAC;IACJ,CAAC;IAED,IAAI,eAAe,GAAY,EAAE,CAAC;IAClC,IAAI,UAAU,GAAG,MAAM,CAAC,iBAAiB,CAAC;IAC1C,IAAI,eAAe,GAAG,CAAC,CAAC;IACxB,IAAI,GAAG,GAAY,EAAE,CAAC;IAEtB,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,MAAM,CAAC,MAAM,EAAE,CAAC,EAAE,EAAE,CAAC;QACvC,MAAM,KAAK,GAAG,QAAQ,CAAC,MAAM,CAAC,CAAC,CAAC,EAAE,MAAM,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,MAAM,CAAC,MAAM,CAAC,CAAC,CAAC;QAEnE,eAAe,GAAG,MAAM,CAAC,CAAC,KAAK,EAAE,MAAM,CAAC,CAAC;QAEzC,8DAA8D;QAC9D,MAAM,EAAE,GAAG,eAAe,CAAC,CAAC,CAAC,CAAC,MAAM,CAAC;QACrC,MAAM,EAAE,GAAG,eAAe,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC;QAClC,MAAM,EAAE,GAAG,eAAe,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,eAAe,CAAC,MAAM,CAAC,CAAC,MAAM,CAAC;QACpE,MAAM,EAAE,GAAG,eAAe,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,eAAe,CAAC,MAAM,CAAC,CAAC,GAAG,CAAC;QAEjE,IAAI,UAAU,GAAG,IAAI,CAAC;QACtB,IAAI,IAAI,GAAG,CAAC,CAAC;QACb,IAAI,IAAI,GAAG,CAAC,CAAC;QACb,IAAI,QAAQ,GAAG,CAAC,CAAC;QACjB,KAAK,MAAM,KAAK,IAAI,eAAe,EAAE,CAAC;YACpC,MAAM,EAAE,GAAG,KAAK,CAAC,MAAM,CAAC;YACxB,MAAM,EAAE,GAAG,KAAK,CAAC,GAAG,CAAC;YACrB,MAAM,CAAC,GAAG,CAAC,EAAE,GAAG,EAAE,CAAC,GAAG,CAAC,EAAE,GAAG,EAAE,CAAC,CAAC;YAChC,IAAI,UAAU,EAAE,CAAC;gBACf,UAAU,GAAG,KAAK,CAAC;gBACnB,IAAI,GAAG,CAAC,CAAC;gBACT,IAAI,GAAG,CAAC,CAAC;YACX,CAAC;iBAAM,CAAC;gBACN,IAAI,CAAC,GAAG,IAAI;oBAAE,IAAI,GAAG,CAAC,CAAC;gBACvB,IAAI,CAAC,GAAG,IAAI;oBAAE,IAAI,GAAG,CAAC,CAAC;YACzB,CAAC;YACD,MAAM,KAAK,GAAG,CAAC,CAAC,CAAC,EAAE,GAAG,EAAE,CAAC,GAAG,EAAE,GAAG,EAAE,GAAG,EAAE,GAAG,EAAE,GAAG,EAAE,CAAC,GAAG,CAAC,EAAE,GAAG,EAAE,CAAC,CAAC;YAEhE,IAAI,IAAI,CAAC,GAAG,CAAC,KAAK,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,QAAQ,CAAC,EAAE,CAAC;gBACzC,QAAQ,GAAG,KAAK,CAAC;YACnB,CAAC;QACH,CAAC;QACD,MAAM,QAAQ,GAAG,EAAE,MAAM,EAAE,EAAE,GAAG,IAAI,GAAG,CAAC,EAAE,GAAG,EAAE,CAAC,EAAE,GAAG,EAAE,EAAE,EAAE,CAAC;QAC5D,MAAM,QAAQ,GAAG,EAAE,MAAM,EAAE,EAAE,GAAG,IAAI,GAAG,CAAC,EAAE,GAAG,EAAE,CAAC,EAAE,GAAG,EAAE,EAAE,EAAE,CAAC;QAE5D,MAAM,cAAc,GAAG,IAAI,CAAC,GAAG,CAAC,QAAQ,GAAG,CAAC,IAAI,GAAG,IAAI,CAAC,GAAG,CAAC,EAAE,GAAG,EAAE,CAAC,CAAC,CAAC;QAEtE,IAAI,cAAc,GAAG,UAAU,EAAE,CAAC;YAChC,eAAe,GAAG,KAAK,CAAC;YACxB,UAAU,GAAG,cAAc,CAAC;YAC5B,GAAG,GAAG;gBACJ,QAAQ;gBACR,QAAQ;gBACR,EAAE,MAAM,EAAE,QAAQ,CAAC,MAAM,EAAE,GAAG,EAAE,QAAQ,CAAC,GAAG,GAAG,QAAQ,EAAE;gBACzD,EAAE,MAAM,EAAE,QAAQ,CAAC,MAAM,EAAE,GAAG,EAAE,QAAQ,CAAC,GAAG,GAAG,QAAQ,EAAE;aAC1D,CAAC;QACJ,CAAC;IACH,CAAC;IAED,MAAM,UAAU,GAAG,MAAM,CAAC,eAAe,EAAE,GAAG,CAAC,CAAC;IAChD,MAAM,KAAK,GAAG;QACZ,IAAI,CAAC,KAAK,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,MAAM,GAAG,GAAG,CAAC,CAAC,CAAC,CAAC,MAAM,EAAE,GAAG,CAAC,CAAC,CAAC,CAAC,GAAG,GAAG,GAAG,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC;QAClE,IAAI,CAAC,KAAK,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,MAAM,GAAG,GAAG,CAAC,CAAC,CAAC,CAAC,MAAM,EAAE,GAAG,CAAC,CAAC,CAAC,CAAC,GAAG,GAAG,GAAG,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC;KACnE,CAAC;IACF,MAAM,OAAO,GAAG,IAAI,CAAC,GAAG,CAAC,GAAG,KAAK,CAAC,CAAC;IACnC,MAAM,OAAO,GAAG,IAAI,CAAC,GAAG,CAAC,GAAG,KAAK,CAAC,CAAC;IACnC,MAAM,QAAQ,GAAG,WAAW,CAAC,UAAU,CAAC,CAAC;IACzC,MAAM,KAAK,GAAG,OAAO,GAAG,OAAO,CAAC;IAEhC,OAAO;QACL,MAAM,EAAE,UAAU;QAClB,OAAO,EAAE,UAAU;QACnB,KAAK,EAAE,QAAQ;QACf,KAAK,EAAE,OAAO;QACd,MAAM,EAAE,OAAO;QACf,SAAS,EAAE,CAAC,GAAG,OAAO,GAAG,CAAC,GAAG,OAAO;QACpC,WAAW,EAAE,KAAK;KACnB,CAAC;AACJ,CAAC"}

package/lib/maskAnalysis/utils/monotoneChainConvexHull.d.ts

@@ -0,0 +1,18 @@
+import type { Point } from '../../utils/geometry/points.js';
+interface McchOptions {
+    /**
+     * Whether the input points are already sorted or not (sort by column and then row).
+     * @default `false`
+     */
+    sorted?: boolean;
+}
+/**
+ * Computes the convex hull of a binary image using Andrew's Monotone Chain Algorithm
+ * @see {@link http://www.algorithmist.com/index.php/Monotone_Chain_Convex_Hull}
+ * @param points - An array of points.
+ * @param options - MCCH Algorithm options.
+ * @returns Coordinates of the convex hull in clockwise order.
+ */
+export declare function monotoneChainConvexHull(points: Point[], options?: McchOptions): Point[];
+export {};
+//# sourceMappingURL=monotoneChainConvexHull.d.ts.map

package/lib/maskAnalysis/utils/monotoneChainConvexHull.d.ts.map

@@ -0,0 +1 @@
+{"version":3,"file":"monotoneChainConvexHull.d.ts","sourceRoot":"","sources":["../../../src/maskAnalysis/utils/monotoneChainConvexHull.ts"],"names":[],"mappings":"AAAA,OAAO,KAAK,EAAE,KAAK,EAAE,MAAM,gCAAgC,CAAC;AAE5D,UAAU,WAAW;IACnB;;;OAGG;IACH,MAAM,CAAC,EAAE,OAAO,CAAC;CAClB;AAED;;;;;;GAMG;AACH,wBAAgB,uBAAuB,CACrC,MAAM,EAAE,KAAK,EAAE,EACf,OAAO,GAAE,WAAgB,GACxB,KAAK,EAAE,CA4BT"}

package/lib/maskAnalysis/utils/monotoneChainConvexHull.js

@@ -0,0 +1,44 @@
+/**
+ * Computes the convex hull of a binary image using Andrew's Monotone Chain Algorithm
+ * @see {@link http://www.algorithmist.com/index.php/Monotone_Chain_Convex_Hull}
+ * @param points - An array of points.
+ * @param options - MCCH Algorithm options.
+ * @returns Coordinates of the convex hull in clockwise order.
+ */
+export function monotoneChainConvexHull(points, options = {}) {
+    const { sorted = false } = options;
+    if (!sorted) {
+        points = points.slice().sort(byXThenY);
+    }
+    const n = points.length;
+    const result = new Array(n * 2);
+    let k = 0;
+    for (let i = 0; i < n; i++) {
+        const point = points[i];
+        while (k >= 2 && cw(result[k - 2], result[k - 1], point) <= 0) {
+            k--;
+        }
+        result[k++] = point;
+    }
+    const t = k + 1;
+    for (let i = n - 2; i >= 0; i--) {
+        const point = points[i];
+        while (k >= t && cw(result[k - 2], result[k - 1], point) <= 0) {
+            k--;
+        }
+        result[k++] = point;
+    }
+    return result.slice(0, k - 1);
+}
+function cw(p1, p2, p3) {
+    return ((p2.row - p1.row) * (p3.column - p1.column) -
+        (p2.column - p1.column) * (p3.row - p1.row));
+}
+function byXThenY(point1, point2) {
+    if (point1.column === point2.column) {
+        return point1.row - point2.row;
+    }
+    return point1.column - point2.column;
+}
+// 0 -> column, 1 -> row
+//# sourceMappingURL=monotoneChainConvexHull.js.map

package/lib/maskAnalysis/utils/monotoneChainConvexHull.js.map

@@ -0,0 +1 @@
+{"version":3,"file":"monotoneChainConvexHull.js","sourceRoot":"","sources":["../../../src/maskAnalysis/utils/monotoneChainConvexHull.ts"],"names":[],"mappings":"AAUA;;;;;;GAMG;AACH,MAAM,UAAU,uBAAuB,CACrC,MAAe,EACf,UAAuB,EAAE;IAEzB,MAAM,EAAE,MAAM,GAAG,KAAK,EAAE,GAAG,OAAO,CAAC;IACnC,IAAI,CAAC,MAAM,EAAE,CAAC;QACZ,MAAM,GAAG,MAAM,CAAC,KAAK,EAAE,CAAC,IAAI,CAAC,QAAQ,CAAC,CAAC;IACzC,CAAC;IAED,MAAM,CAAC,GAAG,MAAM,CAAC,MAAM,CAAC;IACxB,MAAM,MAAM,GAAG,IAAI,KAAK,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC;IAChC,IAAI,CAAC,GAAG,CAAC,CAAC;IAEV,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC;QAC3B,MAAM,KAAK,GAAG,MAAM,CAAC,CAAC,CAAC,CAAC;QACxB,OAAO,CAAC,IAAI,CAAC,IAAI,EAAE,CAAC,MAAM,CAAC,CAAC,GAAG,CAAC,CAAC,EAAE,MAAM,CAAC,CAAC,GAAG,CAAC,CAAC,EAAE,KAAK,CAAC,IAAI,CAAC,EAAE,CAAC;YAC9D,CAAC,EAAE,CAAC;QACN,CAAC;QACD,MAAM,CAAC,CAAC,EAAE,CAAC,GAAG,KAAK,CAAC;IACtB,CAAC;IAED,MAAM,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC;IAChB,KAAK,IAAI,CAAC,GAAG,CAAC,GAAG,CAAC,EAAE,CAAC,IAAI,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC;QAChC,MAAM,KAAK,GAAG,MAAM,CAAC,CAAC,CAAC,CAAC;QACxB,OAAO,CAAC,IAAI,CAAC,IAAI,EAAE,CAAC,MAAM,CAAC,CAAC,GAAG,CAAC,CAAC,EAAE,MAAM,CAAC,CAAC,GAAG,CAAC,CAAC,EAAE,KAAK,CAAC,IAAI,CAAC,EAAE,CAAC;YAC9D,CAAC,EAAE,CAAC;QACN,CAAC;QACD,MAAM,CAAC,CAAC,EAAE,CAAC,GAAG,KAAK,CAAC;IACtB,CAAC;IAED,OAAO,MAAM,CAAC,KAAK,CAAC,CAAC,EAAE,CAAC,GAAG,CAAC,CAAC,CAAC;AAChC,CAAC;AAED,SAAS,EAAE,CAAC,EAAS,EAAE,EAAS,EAAE,EAAS;IACzC,OAAO,CACL,CAAC,EAAE,CAAC,GAAG,GAAG,EAAE,CAAC,GAAG,CAAC,GAAG,CAAC,EAAE,CAAC,MAAM,GAAG,EAAE,CAAC,MAAM,CAAC;QAC3C,CAAC,EAAE,CAAC,MAAM,GAAG,EAAE,CAAC,MAAM,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,GAAG,EAAE,CAAC,GAAG,CAAC,CAC5C,CAAC;AACJ,CAAC;AAED,SAAS,QAAQ,CAAC,MAAa,EAAE,MAAa;IAC5C,IAAI,MAAM,CAAC,MAAM,KAAK,MAAM,CAAC,MAAM,EAAE,CAAC;QACpC,OAAO,MAAM,CAAC,GAAG,GAAG,MAAM,CAAC,GAAG,CAAC;IACjC,CAAC;IACD,OAAO,MAAM,CAAC,MAAM,GAAG,MAAM,CAAC,MAAM,CAAC;AACvC,CAAC;AAED,wBAAwB"}

package/lib/morphology/bottomHat.d.ts

@@ -0,0 +1,18 @@
+import { Image } from '../Image.js';
+import type { Mask } from '../Mask.js';
+export interface BottomHatOptions {
+    /**
+     * 3x3 matrix. The kernel can only have ones and zeros.
+     * Accessing a value: kernel[row][column].
+     * @default `[[1, 1, 1], [1, 1, 1], [1, 1, 1]]`
+     */
+    kernel?: number[][];
+    /**
+     * Number of iterations of the algorithm.
+     * @default `1`
+     */
+    iterations?: number;
+}
+export declare function bottomHat(image: Image, options?: BottomHatOptions): Image;
+export declare function bottomHat(image: Mask, options?: BottomHatOptions): Mask;
+//# sourceMappingURL=bottomHat.d.ts.map

package/lib/morphology/bottomHat.d.ts.map

@@ -0,0 +1 @@
+{"version":3,"file":"bottomHat.d.ts","sourceRoot":"","sources":["../../src/morphology/bottomHat.ts"],"names":[],"mappings":"AAAA,OAAO,EAAE,KAAK,EAAE,MAAM,aAAa,CAAC;AACpC,OAAO,KAAK,EAAE,IAAI,EAAE,MAAM,YAAY,CAAC;AAKvC,MAAM,WAAW,gBAAgB;IAC/B;;;;OAIG;IACH,MAAM,CAAC,EAAE,MAAM,EAAE,EAAE,CAAC;IACpB;;;OAGG;IACH,UAAU,CAAC,EAAE,MAAM,CAAC;CACrB;AAED,wBAAgB,SAAS,CAAC,KAAK,EAAE,KAAK,EAAE,OAAO,CAAC,EAAE,gBAAgB,GAAG,KAAK,CAAC;AAC3E,wBAAgB,SAAS,CAAC,KAAK,EAAE,IAAI,EAAE,OAAO,CAAC,EAAE,gBAAgB,GAAG,IAAI,CAAC"}

package/lib/morphology/bottomHat.js

@@ -0,0 +1,38 @@
+import { Image } from '../Image.js';
+import { subtract } from '../compare/index.js';
+import { checkKernel } from '../utils/validators/checkKernel.js';
+import checkProcessable from '../utils/validators/checkProcessable.js';
+/**
+ * This function is the black top hat (also called bottom hat).
+ * In mathematical morphology and digital image processing,
+ * top-hat transform is an operation that extracts small elements and details from given images.
+ * The black top-hat transform is defined dually as the difference between the closed and the input image.
+ * Top-hat transforms are used for various image processing tasks, such as feature extraction, background equalization,
+ * image enhancement, and others. (Wikipedia)
+ * @see {@link http://docs.opencv.org/2.4/doc/tutorials/imgproc/opening_closing_hats/opening_closing_hats.html}
+ * @param image - Image to process.
+ * @param options - Bottom hat options.
+ * @returns The bottom-hatted image.
+ */
+export function bottomHat(image, options = {}) {
+    const { kernel = [
+        [1, 1, 1],
+        [1, 1, 1],
+        [1, 1, 1],
+    ], iterations = 1, } = options;
+    if (image instanceof Image) {
+        checkProcessable(image, {
+            bitDepth: [1, 8, 16],
+            components: 1,
+            alpha: false,
+        });
+    }
+    checkKernel(kernel);
+    let newImage = image;
+    for (let i = 0; i < iterations; i++) {
+        const openImage = newImage.close({ kernel });
+        newImage = subtract(openImage, newImage, { absolute: true });
+    }
+    return newImage;
+}
+//# sourceMappingURL=bottomHat.js.map

package/lib/morphology/bottomHat.js.map

@@ -0,0 +1 @@
+{"version":3,"file":"bottomHat.js","sourceRoot":"","sources":["../../src/morphology/bottomHat.ts"],"names":[],"mappings":"AAAA,OAAO,EAAE,KAAK,EAAE,MAAM,aAAa,CAAC;AAEpC,OAAO,EAAE,QAAQ,EAAE,MAAM,qBAAqB,CAAC;AAC/C,OAAO,EAAE,WAAW,EAAE,MAAM,oCAAoC,CAAC;AACjE,OAAO,gBAAgB,MAAM,yCAAyC,CAAC;AAkBvE;;;;;;;;;;;GAWG;AACH,MAAM,UAAU,SAAS,CACvB,KAAmB,EACnB,UAA4B,EAAE;IAE9B,MAAM,EACJ,MAAM,GAAG;QACP,CAAC,CAAC,EAAE,CAAC,EAAE,CAAC,CAAC;QACT,CAAC,CAAC,EAAE,CAAC,EAAE,CAAC,CAAC;QACT,CAAC,CAAC,EAAE,CAAC,EAAE,CAAC,CAAC;KACV,EACD,UAAU,GAAG,CAAC,GACf,GAAG,OAAO,CAAC;IAEZ,IAAI,KAAK,YAAY,KAAK,EAAE,CAAC;QAC3B,gBAAgB,CAAC,KAAK,EAAE;YACtB,QAAQ,EAAE,CAAC,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC;YACpB,UAAU,EAAE,CAAC;YACb,KAAK,EAAE,KAAK;SACb,CAAC,CAAC;IACL,CAAC;IAED,WAAW,CAAC,MAAM,CAAC,CAAC;IAEpB,IAAI,QAAQ,GAAG,KAAK,CAAC;IACrB,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,UAAU,EAAE,CAAC,EAAE,EAAE,CAAC;QACpC,MAAM,SAAS,GAAG,QAAQ,CAAC,KAAK,CAAC,EAAE,MAAM,EAAE,CAAC,CAAC;QAC7C,QAAQ,GAAG,QAAQ,CAAC,SAAS,EAAE,QAAQ,EAAE,EAAE,QAAQ,EAAE,IAAI,EAAE,CAAC,CAAC;IAC/D,CAAC;IACD,OAAO,QAAQ,CAAC;AAClB,CAAC"}

package/lib/morphology/cannyEdgeDetector.d.ts

@@ -0,0 +1,44 @@
+import type { Image } from '../Image.js';
+import type { Mask } from '../Mask.js';
+import type { GaussianBlurOptions } from '../filters/index.js';
+export interface CannyEdgeOptions {
+    /**
+     * Lower threshold of the gaussian blur (indicates the weak edges to discard).
+     * @default `0.04`
+     */
+    lowThreshold?: number;
+    /**
+     * Higher threshold of the gaussian blur (indicates the strong edges to keep). Value must be between 0 and 1.
+     * @default `0.1`
+     */
+    highThreshold?: number;
+    /**
+     * Standard deviation of the gaussian blur (sigma). Value must be between 0 and 1.
+     * @default `{ sigma: 1 }`
+     */
+    gaussianBlurOptions?: GaussianBlurOptions;
+    /**
+     * Enable/ disable hysteresis steps.
+     * @default `true`
+     */
+    hysteresis?: boolean;
+    /**
+     * Image to which the resulting image has to be put.
+     */
+    out?: Mask;
+}
+/**
+ * Apply Canny edge detection to an image.
+ * @param image - Image to process.
+ * @param options - Canny edge detection options.
+ * @returns The processed image.
+ */
+export declare function cannyEdgeDetector(image: Image, options?: CannyEdgeOptions): Mask;
+/**
+ * Return a 0 to 3 value indicating the four main directions (horizontal, upward diagonal, vertical, downward diagonal).
+ * @param x - The x coordinate.
+ * @param y - The y coordinate.
+ * @returns The direction as a 0 to 4 value.
+ */
+export declare function getDirection(x: number, y: number): number;
+//# sourceMappingURL=cannyEdgeDetector.d.ts.map

package/lib/morphology/cannyEdgeDetector.d.ts.map

@@ -0,0 +1 @@
+{"version":3,"file":"cannyEdgeDetector.d.ts","sourceRoot":"","sources":["../../src/morphology/cannyEdgeDetector.ts"],"names":[],"mappings":"AAAA,OAAO,KAAK,EAAE,KAAK,EAAE,MAAM,aAAa,CAAC;AACzC,OAAO,KAAK,EAAE,IAAI,EAAE,MAAM,YAAY,CAAC;AACvC,OAAO,KAAK,EAAE,mBAAmB,EAAE,MAAM,qBAAqB,CAAC;AAK/D,MAAM,WAAW,gBAAgB;IAC/B;;;OAGG;IACH,YAAY,CAAC,EAAE,MAAM,CAAC;IACtB;;;OAGG;IACH,aAAa,CAAC,EAAE,MAAM,CAAC;IACvB;;;OAGG;IACH,mBAAmB,CAAC,EAAE,mBAAmB,CAAC;IAC1C;;;OAGG;IACH,UAAU,CAAC,EAAE,OAAO,CAAC;IACrB;;OAEG;IACH,GAAG,CAAC,EAAE,IAAI,CAAC;CACZ;AAcD;;;;;GAKG;AACH,wBAAgB,iBAAiB,CAC/B,KAAK,EAAE,KAAK,EACZ,OAAO,GAAE,gBAAqB,GAC7B,IAAI,CA+IN;AAED;;;;;GAKG;AACH,wBAAgB,YAAY,CAAC,CAAC,EAAE,MAAM,EAAE,CAAC,EAAE,MAAM,GAAG,MAAM,CAEzD"}

package/lib/morphology/cannyEdgeDetector.js

@@ -0,0 +1,146 @@
+import { getIndex } from '../utils/getIndex.js';
+import { imageToOutputMask } from '../utils/getOutputImage.js';
+import checkProcessable from '../utils/validators/checkProcessable.js';
+const kernelX = [
+    [-1, 0, 1],
+    [-2, 0, 2],
+    [-1, 0, 1],
+];
+const kernelY = [
+    [-1, -2, -1],
+    [0, 0, 0],
+    [1, 2, 1],
+];
+/**
+ * Apply Canny edge detection to an image.
+ * @param image - Image to process.
+ * @param options - Canny edge detection options.
+ * @returns The processed image.
+ */
+export function cannyEdgeDetector(image, options = {}) {
+    const { hysteresis = true, lowThreshold = 0.04, highThreshold = 0.1, gaussianBlurOptions = { sigma: 1 }, } = options;
+    const minValue = lowThreshold * image.maxValue;
+    const maxValue = highThreshold * image.maxValue;
+    checkProcessable(image, {
+        colorModel: 'GREY',
+    });
+    const width = image.width;
+    const height = image.height;
+    const blurred = image.gaussianBlur(gaussianBlurOptions);
+    const gradientX = blurred.rawDirectConvolution(kernelY);
+    const gradientY = blurred.rawDirectConvolution(kernelX);
+    const gradient = new Float64Array(image.size);
+    for (let i = 0; i < image.size; i++) {
+        gradient[i] = Math.hypot(gradientX[i], gradientY[i]);
+    }
+    const nonMaxSuppression = new Float64Array(image.size);
+    const edges = new Float64Array(image.size);
+    const finalImage = imageToOutputMask(image, options);
+    // Non-Maximum suppression
+    for (let column = 1; column < width - 1; column++) {
+        for (let row = 1; row < height - 1; row++) {
+            const currentGradientX = gradientX[getIndex(column, row, image, 0)];
+            const currentGradientY = gradientY[getIndex(column, row, image, 0)];
+            const direction = getDirection(currentGradientX, currentGradientY);
+            const currentGradient = gradient[getIndex(column, row, image, 0)];
+            if (
+            // horizontal
+            (direction === 0 &&
+                currentGradient >= gradient[getIndex(column, row - 1, image, 0)] &&
+                currentGradient >= gradient[getIndex(column, row + 1, image, 0)]) ||
+                // upward slope
+                (direction === 1 &&
+                    currentGradient >=
+                        gradient[getIndex(column - 1, row - 1, image, 0)] &&
+                    currentGradient >=
+                        gradient[getIndex(column + 1, row + 1, image, 0)]) ||
+                // vertical
+                (direction === 2 &&
+                    currentGradient >= gradient[getIndex(column - 1, row, image, 0)] &&
+                    currentGradient >= gradient[getIndex(column + 1, row, image, 0)]) ||
+                // downward slope
+                (direction === 3 &&
+                    currentGradient >=
+                        gradient[getIndex(column - 1, row + 1, image, 0)] &&
+                    currentGradient >= gradient[getIndex(column + 1, row - 1, image, 0)])) {
+                // pixels to remove from the final image
+                nonMaxSuppression[getIndex(column, row, image, 0)] = currentGradient;
+            }
+        }
+    }
+    for (let i = 0; i < width * height; ++i) {
+        const currentNms = nonMaxSuppression[i];
+        let currentEdge = 0;
+        if (currentNms > maxValue) {
+            currentEdge++;
+            finalImage.setBitByIndex(i, 1);
+        }
+        if (currentNms > minValue) {
+            currentEdge++;
+        }
+        edges[i] = currentEdge;
+    }
+    // Hysteresis: first pass
+    if (hysteresis) {
+        let currentPixels = [];
+        for (let column = 1; column < width - 1; ++column) {
+            for (let row = 1; row < height - 1; ++row) {
+                if (edges[getIndex(column, row, image, 0)] !== 1) {
+                    continue;
+                }
+                outer: for (let hystColumn = column - 1; hystColumn < column + 2; ++hystColumn) {
+                    for (let hystRow = row - 1; hystRow < row + 2; ++hystRow) {
+                        if (edges[getIndex(hystColumn, hystRow, image, 0)] === 2) {
+                            currentPixels.push([column, row]);
+                            finalImage.setValue(column, row, 0, 1);
+                            break outer;
+                        }
+                    }
+                }
+            }
+        }
+        // Hysteresis: second pass
+        while (currentPixels.length > 0) {
+            const newPixels = [];
+            for (const currentPixel of currentPixels) {
+                for (let j = -1; j < 2; ++j) {
+                    for (let k = -1; k < 2; ++k) {
+                        if (j === 0 && k === 0) {
+                            continue;
+                        }
+                        const row = currentPixel[0] + j;
+                        const column = currentPixel[1] + k;
+                        if (
+                        // there could be an error here
+                        edges[getIndex(column, row, image, 0)] === 1 &&
+                            finalImage.getValue(column, row, 0) === 0) {
+                            newPixels.push([row, column]);
+                            finalImage.setValue(column, row, 0, 1);
+                        }
+                    }
+                }
+            }
+            currentPixels = newPixels;
+        }
+    }
+    return finalImage;
+    /* Function for debug
+    import { Matrix } from 'ml-matrix';
+  
+    function printArray(array: Float64Array): void {
+      // @ts-expect-error: only for debug
+      const matrix = Matrix.from1DArray(height, width, array);
+      console.log(matrix);
+    }
+    */
+}
+/**
+ * Return a 0 to 3 value indicating the four main directions (horizontal, upward diagonal, vertical, downward diagonal).
+ * @param x - The x coordinate.
+ * @param y - The y coordinate.
+ * @returns The direction as a 0 to 4 value.
+ */
+export function getDirection(x, y) {
+    return (Math.round(Math.atan2(y, x) * (4 / Math.PI)) + 4) % 4;
+}
+//# sourceMappingURL=cannyEdgeDetector.js.map

package/lib/morphology/cannyEdgeDetector.js.map

@@ -0,0 +1 @@
+{"version":3,"file":"cannyEdgeDetector.js","sourceRoot":"","sources":["../../src/morphology/cannyEdgeDetector.ts"],"names":[],"mappings":"AAGA,OAAO,EAAE,QAAQ,EAAE,MAAM,sBAAsB,CAAC;AAChD,OAAO,EAAE,iBAAiB,EAAE,MAAM,4BAA4B,CAAC;AAC/D,OAAO,gBAAgB,MAAM,yCAAyC,CAAC;AA6BvE,MAAM,OAAO,GAAG;IACd,CAAC,CAAC,CAAC,EAAE,CAAC,EAAE,CAAC,CAAC;IACV,CAAC,CAAC,CAAC,EAAE,CAAC,EAAE,CAAC,CAAC;IACV,CAAC,CAAC,CAAC,EAAE,CAAC,EAAE,CAAC,CAAC;CACX,CAAC;AAEF,MAAM,OAAO,GAAG;IACd,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC;IACZ,CAAC,CAAC,EAAE,CAAC,EAAE,CAAC,CAAC;IACT,CAAC,CAAC,EAAE,CAAC,EAAE,CAAC,CAAC;CACV,CAAC;AAEF;;;;;GAKG;AACH,MAAM,UAAU,iBAAiB,CAC/B,KAAY,EACZ,UAA4B,EAAE;IAE9B,MAAM,EACJ,UAAU,GAAG,IAAI,EACjB,YAAY,GAAG,IAAI,EACnB,aAAa,GAAG,GAAG,EACnB,mBAAmB,GAAG,EAAE,KAAK,EAAE,CAAC,EAAE,GACnC,GAAG,OAAO,CAAC;IAEZ,MAAM,QAAQ,GAAG,YAAY,GAAG,KAAK,CAAC,QAAQ,CAAC;IAC/C,MAAM,QAAQ,GAAG,aAAa,GAAG,KAAK,CAAC,QAAQ,CAAC;IAEhD,gBAAgB,CAAC,KAAK,EAAE;QACtB,UAAU,EAAE,MAAM;KACnB,CAAC,CAAC;IAEH,MAAM,KAAK,GAAG,KAAK,CAAC,KAAK,CAAC;IAC1B,MAAM,MAAM,GAAG,KAAK,CAAC,MAAM,CAAC;IAE5B,MAAM,OAAO,GAAG,KAAK,CAAC,YAAY,CAAC,mBAAmB,CAAC,CAAC;IAExD,MAAM,SAAS,GAAG,OAAO,CAAC,oBAAoB,CAAC,OAAO,CAAC,CAAC;IACxD,MAAM,SAAS,GAAG,OAAO,CAAC,oBAAoB,CAAC,OAAO,CAAC,CAAC;IAExD,MAAM,QAAQ,GAAG,IAAI,YAAY,CAAC,KAAK,CAAC,IAAI,CAAC,CAAC;IAC9C,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,KAAK,CAAC,IAAI,EAAE,CAAC,EAAE,EAAE,CAAC;QACpC,QAAQ,CAAC,CAAC,CAAC,GAAG,IAAI,CAAC,KAAK,CAAC,SAAS,CAAC,CAAC,CAAC,EAAE,SAAS,CAAC,CAAC,CAAC,CAAC,CAAC;IACvD,CAAC;IAED,MAAM,iBAAiB,GAAG,IAAI,YAAY,CAAC,KAAK,CAAC,IAAI,CAAC,CAAC;IACvD,MAAM,KAAK,GAAG,IAAI,YAAY,CAAC,KAAK,CAAC,IAAI,CAAC,CAAC;IAE3C,MAAM,UAAU,GAAG,iBAAiB,CAAC,KAAK,EAAE,OAAO,CAAC,CAAC;IAErD,0BAA0B;IAC1B,KAAK,IAAI,MAAM,GAAG,CAAC,EAAE,MAAM,GAAG,KAAK,GAAG,CAAC,EAAE,MAAM,EAAE,EAAE,CAAC;QAClD,KAAK,IAAI,GAAG,GAAG,CAAC,EAAE,GAAG,GAAG,MAAM,GAAG,CAAC,EAAE,GAAG,EAAE,EAAE,CAAC;YAC1C,MAAM,gBAAgB,GAAG,SAAS,CAAC,QAAQ,CAAC,MAAM,EAAE,GAAG,EAAE,KAAK,EAAE,CAAC,CAAC,CAAC,CAAC;YACpE,MAAM,gBAAgB,GAAG,SAAS,CAAC,QAAQ,CAAC,MAAM,EAAE,GAAG,EAAE,KAAK,EAAE,CAAC,CAAC,CAAC,CAAC;YACpE,MAAM,SAAS,GAAG,YAAY,CAAC,gBAAgB,EAAE,gBAAgB,CAAC,CAAC;YACnE,MAAM,eAAe,GAAG,QAAQ,CAAC,QAAQ,CAAC,MAAM,EAAE,GAAG,EAAE,KAAK,EAAE,CAAC,CAAC,CAAC,CAAC;YAClE;YACE,aAAa;YACb,CAAC,SAAS,KAAK,CAAC;gBACd,eAAe,IAAI,QAAQ,CAAC,QAAQ,CAAC,MAAM,EAAE,GAAG,GAAG,CAAC,EAAE,KAAK,EAAE,CAAC,CAAC,CAAC;gBAChE,eAAe,IAAI,QAAQ,CAAC,QAAQ,CAAC,MAAM,EAAE,GAAG,GAAG,CAAC,EAAE,KAAK,EAAE,CAAC,CAAC,CAAC,CAAC;gBACnE,eAAe;gBACf,CAAC,SAAS,KAAK,CAAC;oBACd,eAAe;wBACb,QAAQ,CAAC,QAAQ,CAAC,MAAM,GAAG,CAAC,EAAE,GAAG,GAAG,CAAC,EAAE,KAAK,EAAE,CAAC,CAAC,CAAC;oBACnD,eAAe;wBACb,QAAQ,CAAC,QAAQ,CAAC,MAAM,GAAG,CAAC,EAAE,GAAG,GAAG,CAAC,EAAE,KAAK,EAAE,CAAC,CAAC,CAAC,CAAC;gBACtD,WAAW;gBACX,CAAC,SAAS,KAAK,CAAC;oBACd,eAAe,IAAI,QAAQ,CAAC,QAAQ,CAAC,MAAM,GAAG,CAAC,EAAE,GAAG,EAAE,KAAK,EAAE,CAAC,CAAC,CAAC;oBAChE,eAAe,IAAI,QAAQ,CAAC,QAAQ,CAAC,MAAM,GAAG,CAAC,EAAE,GAAG,EAAE,KAAK,EAAE,CAAC,CAAC,CAAC,CAAC;gBACnE,iBAAiB;gBACjB,CAAC,SAAS,KAAK,CAAC;oBACd,eAAe;wBACb,QAAQ,CAAC,QAAQ,CAAC,MAAM,GAAG,CAAC,EAAE,GAAG,GAAG,CAAC,EAAE,KAAK,EAAE,CAAC,CAAC,CAAC;oBACnD,eAAe,IAAI,QAAQ,CAAC,QAAQ,CAAC,MAAM,GAAG,CAAC,EAAE,GAAG,GAAG,CAAC,EAAE,KAAK,EAAE,CAAC,CAAC,CAAC,CAAC,EACvE,CAAC;gBACD,wCAAwC;gBACxC,iBAAiB,CAAC,QAAQ,CAAC,MAAM,EAAE,GAAG,EAAE,KAAK,EAAE,CAAC,CAAC,CAAC,GAAG,eAAe,CAAC;YACvE,CAAC;QACH,CAAC;IACH,CAAC;IAED,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,KAAK,GAAG,MAAM,EAAE,EAAE,CAAC,EAAE,CAAC;QACxC,MAAM,UAAU,GAAG,iBAAiB,CAAC,CAAC,CAAC,CAAC;QACxC,IAAI,WAAW,GAAG,CAAC,CAAC;QACpB,IAAI,UAAU,GAAG,QAAQ,EAAE,CAAC;YAC1B,WAAW,EAAE,CAAC;YACd,UAAU,CAAC,aAAa,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC;QACjC,CAAC;QACD,IAAI,UAAU,GAAG,QAAQ,EAAE,CAAC;YAC1B,WAAW,EAAE,CAAC;QAChB,CAAC;QAED,KAAK,CAAC,CAAC,CAAC,GAAG,WAAW,CAAC;IACzB,CAAC;IAED,yBAAyB;IACzB,IAAI,UAAU,EAAE,CAAC;QACf,IAAI,aAAa,GAAe,EAAE,CAAC;QACnC,KAAK,IAAI,MAAM,GAAG,CAAC,EAAE,MAAM,GAAG,KAAK,GAAG,CAAC,EAAE,EAAE,MAAM,EAAE,CAAC;YAClD,KAAK,IAAI,GAAG,GAAG,CAAC,EAAE,GAAG,GAAG,MAAM,GAAG,CAAC,EAAE,EAAE,GAAG,EAAE,CAAC;gBAC1C,IAAI,KAAK,CAAC,QAAQ,CAAC,MAAM,EAAE,GAAG,EAAE,KAAK,EAAE,CAAC,CAAC,CAAC,KAAK,CAAC,EAAE,CAAC;oBACjD,SAAS;gBACX,CAAC;gBAED,KAAK,EAAE,KACL,IAAI,UAAU,GAAG,MAAM,GAAG,CAAC,EAC3B,UAAU,GAAG,MAAM,GAAG,CAAC,EACvB,EAAE,UAAU,EACZ,CAAC;oBACD,KAAK,IAAI,OAAO,GAAG,GAAG,GAAG,CAAC,EAAE,OAAO,GAAG,GAAG,GAAG,CAAC,EAAE,EAAE,OAAO,EAAE,CAAC;wBACzD,IAAI,KAAK,CAAC,QAAQ,CAAC,UAAU,EAAE,OAAO,EAAE,KAAK,EAAE,CAAC,CAAC,CAAC,KAAK,CAAC,EAAE,CAAC;4BACzD,aAAa,CAAC,IAAI,CAAC,CAAC,MAAM,EAAE,GAAG,CAAC,CAAC,CAAC;4BAClC,UAAU,CAAC,QAAQ,CAAC,MAAM,EAAE,GAAG,EAAE,CAAC,EAAE,CAAC,CAAC,CAAC;4BACvC,MAAM,KAAK,CAAC;wBACd,CAAC;oBACH,CAAC;gBACH,CAAC;YACH,CAAC;QACH,CAAC;QAED,0BAA0B;QAC1B,OAAO,aAAa,CAAC,MAAM,GAAG,CAAC,EAAE,CAAC;YAChC,MAAM,SAAS,GAAG,EAAE,CAAC;YACrB,KAAK,MAAM,YAAY,IAAI,aAAa,EAAE,CAAC;gBACzC,KAAK,IAAI,CAAC,GAAG,CAAC,CAAC,EAAE,CAAC,GAAG,CAAC,EAAE,EAAE,CAAC,EAAE,CAAC;oBAC5B,KAAK,IAAI,CAAC,GAAG,CAAC,CAAC,EAAE,CAAC,GAAG,CAAC,EAAE,EAAE,CAAC,EAAE,CAAC;wBAC5B,IAAI,CAAC,KAAK,CAAC,IAAI,CAAC,KAAK,CAAC,EAAE,CAAC;4BACvB,SAAS;wBACX,CAAC;wBACD,MAAM,GAAG,GAAG,YAAY,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC;wBAChC,MAAM,MAAM,GAAG,YAAY,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC;wBACnC;wBACE,+BAA+B;wBAC/B,KAAK,CAAC,QAAQ,CAAC,MAAM,EAAE,GAAG,EAAE,KAAK,EAAE,CAAC,CAAC,CAAC,KAAK,CAAC;4BAC5C,UAAU,CAAC,QAAQ,CAAC,MAAM,EAAE,GAAG,EAAE,CAAC,CAAC,KAAK,CAAC,EACzC,CAAC;4BACD,SAAS,CAAC,IAAI,CAAC,CAAC,GAAG,EAAE,MAAM,CAAC,CAAC,CAAC;4BAC9B,UAAU,CAAC,QAAQ,CAAC,MAAM,EAAE,GAAG,EAAE,CAAC,EAAE,CAAC,CAAC,CAAC;wBACzC,CAAC;oBACH,CAAC;gBACH,CAAC;YACH,CAAC;YACD,aAAa,GAAG,SAAS,CAAC;QAC5B,CAAC;IACH,CAAC;IAED,OAAO,UAAU,CAAC;IAElB;;;;;;;;MAQE;AACJ,CAAC;AAED;;;;;GAKG;AACH,MAAM,UAAU,YAAY,CAAC,CAAS,EAAE,CAAS;IAC/C,OAAO,CAAC,IAAI,CAAC,KAAK,CAAC,IAAI,CAAC,KAAK,CAAC,CAAC,EAAE,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,IAAI,CAAC,EAAE,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC;AAChE,CAAC"}

package/lib/morphology/clearBorder.d.ts

@@ -0,0 +1,26 @@
+import type { Mask } from '../Mask.js';
+export interface ClearBorderOptions {
+    /**
+     * Consider pixels connected by corners?
+     * @default `false`
+     */
+    allowCorners?: boolean;
+    /**
+     * Image to which the resulting image has to be put.
+     */
+    out?: Mask;
+    /**
+     * Clear either white or black area touching the border.
+     * In practice it will invert the color.
+     * @default `'white'`
+     */
+    color?: 'white' | 'black';
+}
+/**
+ * Set the pixels connected to the border of the mask to zero. You can chose to allow corner connection of not with the `allowCorners` option.
+ * @param mask - The mask to process.
+ * @param options - Clear border options.
+ * @returns The image with cleared borders.
+ */
+export declare function clearBorder(mask: Mask, options?: ClearBorderOptions): Mask;
+//# sourceMappingURL=clearBorder.d.ts.map

package/lib/morphology/clearBorder.d.ts.map

@@ -0,0 +1 @@
+{"version":3,"file":"clearBorder.d.ts","sourceRoot":"","sources":["../../src/morphology/clearBorder.ts"],"names":[],"mappings":"AAAA,OAAO,KAAK,EAAE,IAAI,EAAE,MAAM,YAAY,CAAC;AAKvC,MAAM,WAAW,kBAAkB;IACjC;;;OAGG;IACH,YAAY,CAAC,EAAE,OAAO,CAAC;IACvB;;OAEG;IACH,GAAG,CAAC,EAAE,IAAI,CAAC;IACX;;;;OAIG;IACH,KAAK,CAAC,EAAE,OAAO,GAAG,OAAO,CAAC;CAC3B;AACD;;;;;GAKG;AACH,wBAAgB,WAAW,CACzB,IAAI,EAAE,IAAI,EACV,OAAO,GAAE,kBAAuB,GAC/B,IAAI,CASN"}

package/lib/morphology/clearBorder.js

@@ -0,0 +1,19 @@
+import { borderIterator } from '../utils/borderIterator.js';
+import { multipleFloodFill } from './multipleFloodFill.js';
+/**
+ * Set the pixels connected to the border of the mask to zero. You can chose to allow corner connection of not with the `allowCorners` option.
+ * @param mask - The mask to process.
+ * @param options - Clear border options.
+ * @returns The image with cleared borders.
+ */
+export function clearBorder(mask, options = {}) {
+    const { allowCorners = false, out, color = 'white' } = options;
+    return multipleFloodFill(mask, {
+        startPixels: borderIterator(mask),
+        startPixelValue: color === 'white' ? 1 : 0,
+        newPixelValue: color === 'white' ? 0 : 1,
+        allowCorners,
+        out,
+    });
+}
+//# sourceMappingURL=clearBorder.js.map

package/lib/morphology/clearBorder.js.map

@@ -0,0 +1 @@
+{"version":3,"file":"clearBorder.js","sourceRoot":"","sources":["../../src/morphology/clearBorder.ts"],"names":[],"mappings":"AACA,OAAO,EAAE,cAAc,EAAE,MAAM,4BAA4B,CAAC;AAE5D,OAAO,EAAE,iBAAiB,EAAE,MAAM,wBAAwB,CAAC;AAmB3D;;;;;GAKG;AACH,MAAM,UAAU,WAAW,CACzB,IAAU,EACV,UAA8B,EAAE;IAEhC,MAAM,EAAE,YAAY,GAAG,KAAK,EAAE,GAAG,EAAE,KAAK,GAAG,OAAO,EAAE,GAAG,OAAO,CAAC;IAC/D,OAAO,iBAAiB,CAAC,IAAI,EAAE;QAC7B,WAAW,EAAE,cAAc,CAAC,IAAI,CAAC;QACjC,eAAe,EAAE,KAAK,KAAK,OAAO,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;QAC1C,aAAa,EAAE,KAAK,KAAK,OAAO,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;QACxC,YAAY;QACZ,GAAG;KACJ,CAAC,CAAC;AACL,CAAC"}

package/lib/morphology/close.d.ts

@@ -0,0 +1,18 @@
+import { Image } from '../Image.js';
+import type { Mask } from '../Mask.js';
+export interface CloseOptions {
+    /**
+     * 3x3 matrix. The kernel can only have ones and zeros.
+     * Accessing a value: kernel[row][column].
+     * @default `[[1, 1, 1], [1, 1, 1], [1, 1, 1]]`
+     */
+    kernel?: number[][];
+    /**
+     * Number of iterations of the algorithm.
+     * @default `1`
+     */
+    iterations?: number;
+}
+export declare function close(image: Image, options?: CloseOptions): Image;
+export declare function close(image: Mask, options?: CloseOptions): Mask;
+//# sourceMappingURL=close.d.ts.map

package/lib/morphology/close.d.ts.map

@@ -0,0 +1 @@
+{"version":3,"file":"close.d.ts","sourceRoot":"","sources":["../../src/morphology/close.ts"],"names":[],"mappings":"AAAA,OAAO,EAAE,KAAK,EAAE,MAAM,aAAa,CAAC;AACpC,OAAO,KAAK,EAAE,IAAI,EAAE,MAAM,YAAY,CAAC;AAIvC,MAAM,WAAW,YAAY;IAC3B;;;;OAIG;IACH,MAAM,CAAC,EAAE,MAAM,EAAE,EAAE,CAAC;IACpB;;;OAGG;IACH,UAAU,CAAC,EAAE,MAAM,CAAC;CACrB;AAED,wBAAgB,KAAK,CAAC,KAAK,EAAE,KAAK,EAAE,OAAO,CAAC,EAAE,YAAY,GAAG,KAAK,CAAC;AACnE,wBAAgB,KAAK,CAAC,KAAK,EAAE,IAAI,EAAE,OAAO,CAAC,EAAE,YAAY,GAAG,IAAI,CAAC"}

package/lib/morphology/close.js

@@ -0,0 +1,33 @@
+import { Image } from '../Image.js';
+import { checkKernel } from '../utils/validators/checkKernel.js';
+import checkProcessable from '../utils/validators/checkProcessable.js';
+/**
+ * In mathematical morphology, the closing of a set A by a structuring element B is the erosion of the dilation of that set (Wikipedia).
+ * In image processing, closing is, together with opening, the basic workhorse of morphological noise removal.
+ * Opening removes small objects, while closing removes small holes.
+ * @see {@link http://docs.opencv.org/2.4/doc/tutorials/imgproc/opening_closing_hats/opening_closing_hats.html}
+ * @param image - Image to process.
+ * @param options - Close options.
+ * @returns Closed image.
+ */
+export function close(image, options = {}) {
+    const { kernel = [
+        [1, 1, 1],
+        [1, 1, 1],
+        [1, 1, 1],
+    ], iterations = 1, } = options;
+    if (image instanceof Image) {
+        checkProcessable(image, {
+            bitDepth: [1, 8, 16],
+            components: 1,
+            alpha: false,
+        });
+    }
+    checkKernel(kernel);
+    let newImage = image;
+    for (let i = 0; i < iterations; i++) {
+        newImage = newImage.dilate({ kernel }).erode({ kernel });
+    }
+    return newImage;
+}
+//# sourceMappingURL=close.js.map

package/lib/morphology/close.js.map

@@ -0,0 +1 @@
+{"version":3,"file":"close.js","sourceRoot":"","sources":["../../src/morphology/close.ts"],"names":[],"mappings":"AAAA,OAAO,EAAE,KAAK,EAAE,MAAM,aAAa,CAAC;AAEpC,OAAO,EAAE,WAAW,EAAE,MAAM,oCAAoC,CAAC;AACjE,OAAO,gBAAgB,MAAM,yCAAyC,CAAC;AAkBvE;;;;;;;;GAQG;AACH,MAAM,UAAU,KAAK,CACnB,KAAmB,EACnB,UAAwB,EAAE;IAE1B,MAAM,EACJ,MAAM,GAAG;QACP,CAAC,CAAC,EAAE,CAAC,EAAE,CAAC,CAAC;QACT,CAAC,CAAC,EAAE,CAAC,EAAE,CAAC,CAAC;QACT,CAAC,CAAC,EAAE,CAAC,EAAE,CAAC,CAAC;KACV,EACD,UAAU,GAAG,CAAC,GACf,GAAG,OAAO,CAAC;IAEZ,IAAI,KAAK,YAAY,KAAK,EAAE,CAAC;QAC3B,gBAAgB,CAAC,KAAK,EAAE;YACtB,QAAQ,EAAE,CAAC,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC;YACpB,UAAU,EAAE,CAAC;YACb,KAAK,EAAE,KAAK;SACb,CAAC,CAAC;IACL,CAAC;IAED,WAAW,CAAC,MAAM,CAAC,CAAC;IAEpB,IAAI,QAAQ,GAAG,KAAK,CAAC;IACrB,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,UAAU,EAAE,CAAC,EAAE,EAAE,CAAC;QACpC,QAAQ,GAAG,QAAQ,CAAC,MAAM,CAAC,EAAE,MAAM,EAAE,CAAC,CAAC,KAAK,CAAC,EAAE,MAAM,EAAE,CAAC,CAAC;IAC3D,CAAC;IACD,OAAO,QAAQ,CAAC;AAClB,CAAC"}

package/lib/morphology/dilate.d.ts

@@ -0,0 +1,18 @@
+import { Image } from '../Image.js';
+import { Mask } from '../Mask.js';
+export interface DilateOptions {
+    /**
+     * Matrix with odd dimensions (e.g. 1 by 3). The kernel can only have ones and zeros.
+     * Accessing a value: kernel[row][column].
+     * @default `[[1, 1, 1], [1, 1, 1], [1, 1, 1]]`
+     */
+    kernel?: number[][];
+    /**
+     * Number of iterations of the algorithm.
+     *  @default `1`
+     */
+    iterations?: number;
+}
+export declare function dilate(image: Image, options?: DilateOptions): Image;
+export declare function dilate(image: Mask, options?: DilateOptions): Mask;
+//# sourceMappingURL=dilate.d.ts.map

package/lib/morphology/dilate.d.ts.map

@@ -0,0 +1 @@
+{"version":3,"file":"dilate.d.ts","sourceRoot":"","sources":["../../src/morphology/dilate.ts"],"names":[],"mappings":"AAAA,OAAO,EAAE,KAAK,EAAE,MAAM,aAAa,CAAC;AACpC,OAAO,EAAE,IAAI,EAAE,MAAM,YAAY,CAAC;AAIlC,MAAM,WAAW,aAAa;IAC5B;;;;OAIG;IACH,MAAM,CAAC,EAAE,MAAM,EAAE,EAAE,CAAC;IACpB;;;OAGG;IACH,UAAU,CAAC,EAAE,MAAM,CAAC;CACrB;AAED,wBAAgB,MAAM,CAAC,KAAK,EAAE,KAAK,EAAE,OAAO,CAAC,EAAE,aAAa,GAAG,KAAK,CAAC;AACrE,wBAAgB,MAAM,CAAC,KAAK,EAAE,IAAI,EAAE,OAAO,CAAC,EAAE,aAAa,GAAG,IAAI,CAAC"}