image-js 1.6.1 → 1.7.0

package/lib/maskAnalysis/getContourMoore.js.map

@@ -1 +1 @@
-{"version":3,"file":"getContourMoore.js","sourceRoot":"","sources":["../../src/maskAnalysis/getContourMoore.ts"],"names":[],"mappings":"AAAA,OAAO,MAAM,MAAM,aAAa,CAAC;AAKjC,OAAO,EACL,UAAU,EACV,kBAAkB,GACnB,MAAM,oCAAoC,CAAC;AAY5C;;GAEG;AACH,MAAM,UAAU,GAAgB;IAC9B,EAAE,EAAE,EAAE,CAAC,CAAC,EAAE,EAAE,EAAE,CAAC,EAAE,EAAE,OAAO;IAC1B,EAAE,EAAE,EAAE,CAAC,CAAC,EAAE,EAAE,EAAE,CAAC,CAAC,EAAE,EAAE,UAAU;IAC9B,EAAE,EAAE,EAAE,CAAC,EAAE,EAAE,EAAE,CAAC,CAAC,EAAE,EAAE,KAAK;IACxB,EAAE,EAAE,EAAE,CAAC,EAAE,EAAE,EAAE,CAAC,CAAC,EAAE,EAAE,WAAW;IAC9B,EAAE,EAAE,EAAE,CAAC,EAAE,EAAE,EAAE,CAAC,EAAE,EAAE,QAAQ;IAC1B,EAAE,EAAE,EAAE,CAAC,EAAE,EAAE,EAAE,CAAC,EAAE,EAAE,aAAa;IAC/B,EAAE,EAAE,EAAE,CAAC,EAAE,EAAE,EAAE,CAAC,EAAE,EAAE,OAAO;IACzB,EAAE,EAAE,EAAE,CAAC,CAAC,EAAE,EAAE,EAAE,CAAC,EAAE,EAAE,YAAY;CAChC,CAAC;AACF,sCAAsC;AACtC,MAAM,oBAAoB,GAAG,IAAI,SAAS,CAAC,CAAC,CAAC,CAAC,IAAI,CAAC,CAAC,CAAC,CAAC,CAAC;AACvD,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,UAAU,CAAC,MAAM,EAAE,CAAC,EAAE,EAAE,CAAC;IAC3C,MAAM,EAAE,EAAE,EAAE,EAAE,EAAE,GAAG,UAAU,CAAC,CAAC,CAAC,CAAC;IACjC,oBAAoB,CAAC,CAAC,EAAE,GAAG,CAAC,CAAC,GAAG,CAAC,GAAG,CAAC,EAAE,GAAG,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC;AACpD,CAAC;AACD;;;;;GAKG;AACH,MAAM,UAAU,eAAe,CAAC,IAAU;IACxC,IAAI,KAAK,GAAG,CAAC,CAAC;IAEd,MAAM,OAAO,GAAY,EAAE,CAAC;IAC5B,MAAM,OAAO,GAAG,IAAI,UAAU,CAAC,IAAI,CAAC,KAAK,GAAG,IAAI,CAAC,MAAM,CAAC,CAAC;IACzD,OAAO,IAAI,CAAC,aAAa,CAAC,KAAK,CAAC,KAAK,CAAC,IAAI,KAAK,KAAK,IAAI,CAAC,IAAI,EAAE,CAAC;QAC9D,KAAK,EAAE,CAAC;IACV,CAAC;IAED,IAAI,KAAK,KAAK,IAAI,CAAC,IAAI;QAAE,OAAO,OAAO,CAAC;IACxC,MAAM,CAAC,GAAG,EAAE,GAAG,CAAC,GAAG,kBAAkB,CAAC,KAAK,EAAE,IAAI,CAAC,KAAK,CAAC,CAAC;IACzD,MAAM,aAAa,GAAG,EAAE,MAAM,EAAE,GAAG,EAAE,GAAG,EAAE,CAAC;IAC3C,IAAI,YAAY,GAAG,aAAa,CAAC;IAEjC,OAAO,CAAC,KAAK,CAAC,GAAG,CAAC,CAAC;IACnB,OAAO,CAAC,IAAI,CAAC,YAAY,CAAC,CAAC;IAE3B,MAAM,SAAS,GAAG,aAAa,CAAC,IAAI,EAAE,YAAY,EAAE,EAAE,MAAM,EAAE,GAAG,GAAG,CAAC,EAAE,GAAG,EAAE,CAAC,CAAC;IAC9E,IAAI,CAAC,SAAS,EAAE,CAAC;QACf,OAAO,OAAO,CAAC;IACjB,CAAC;IAED,IAAI,sBAAsB,GAAU,SAAS,CAAC,SAAS,CAAC;IAExD,IAAI,cAAc,GAAG,SAAS,CAAC,SAAS,CAAC;IACzC,YAAY,GAAG,SAAS,CAAC,SAAS,CAAC;IAEnC,OAAO,IAAI,EAAE,CAAC;QACZ,KAAK,GAAG,YAAY,CAAC,GAAG,GAAG,IAAI,CAAC,KAAK,GAAG,YAAY,CAAC,MAAM,CAAC;QAC5D,IAAI,CAAC,OAAO,CAAC,KAAK,CAAC,EAAE,CAAC;YACpB,OAAO,CAAC,KAAK,CAAC,GAAG,CAAC,CAAC;YACnB,OAAO,CAAC,IAAI,CAAC,YAAY,CAAC,CAAC;QAC7B,CAAC;QAED,MAAM,IAAI,GAAG,aAAa,CAAC,IAAI,EAAE,YAAY,EAAE,cAAc,CAAC,CAAC;QAC/D,MAAM,CAAC,EAAE,CAAC,IAAI,EAAE,iCAAiC,CAAC,CAAC;QACnD,cAAc,GAAG,IAAI,CAAC,SAAS,CAAC;QAChC,YAAY,GAAG,IAAI,CAAC,SAAS,CAAC;QAC9B,IACE,YAAY,CAAC,MAAM,KAAK,aAAa,CAAC,MAAM;YAC5C,YAAY,CAAC,GAAG,KAAK,aAAa,CAAC,GAAG,EACtC,CAAC;YACD,IACE,cAAc,CAAC,MAAM,KAAK,sBAAsB,CAAC,MAAM;gBACvD,cAAc,CAAC,GAAG,KAAK,sBAAsB,CAAC,GAAG,EACjD,CAAC;gBACD,MAAM;YACR,CAAC;iBAAM,CAAC;gBACN,sBAAsB,GAAG,cAAc,CAAC;YAC1C,CAAC;QACH,CAAC;IACH,CAAC;IAED,OAAO,OAAO,CAAC;AACjB,CAAC;AACD;;;;;;GAMG;AACH,SAAS,aAAa,CACpB,IAAU,EACV,OAAc,EACd,QAAe;IAEf,IAAI,SAAS,GAAG,QAAQ,CAAC;IAEzB,uCAAuC;IACvC,MAAM,MAAM,GAAG,QAAQ,CAAC,MAAM,GAAG,OAAO,CAAC,MAAM,CAAC;IAChD,MAAM,MAAM,GAAG,QAAQ,CAAC,GAAG,GAAG,OAAO,CAAC,GAAG,CAAC;IAE1C,MAAM,UAAU,GAAG,oBAAoB,CAAC,CAAC,MAAM,GAAG,CAAC,CAAC,GAAG,CAAC,GAAG,CAAC,MAAM,GAAG,CAAC,CAAC,CAAC,CAAC;IAEzE,2CAA2C;IAC3C,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,UAAU,CAAC,MAAM,EAAE,CAAC,EAAE,EAAE,CAAC;QAC3C,MAAM,EAAE,EAAE,EAAE,EAAE,EAAE,GAAG,UAAU,CAAC,EAAE,CAC9B,CAAC,UAAU,GAAG,CAAC,CAAC,GAAG,UAAU,CAAC,MAAM,CACxB,CAAC;QACf,MAAM,MAAM,GAAG,OAAO,CAAC,MAAM,GAAG,EAAE,CAAC;QACnC,MAAM,MAAM,GAAG,OAAO,CAAC,GAAG,GAAG,EAAE,CAAC;QAChC,IAAI,UAAU,CAAC,IAAI,EAAE,MAAM,EAAE,MAAM,CAAC,KAAK,CAAC,EAAE,CAAC;YAC3C,2DAA2D;YAC3D,0DAA0D;YAC1D,aAAa;YACb,OAAO;gBACL,SAAS;gBACT,SAAS,EAAE,EAAE,MAAM,EAAE,MAAM,EAAE,GAAG,EAAE,MAAM,EAAE;aAC3C,CAAC;QACJ,CAAC;aAAM,CAAC;YACN,SAAS,GAAG,EAAE,MAAM,EAAE,MAAM,EAAE,GAAG,EAAE,MAAM,EAAE,CAAC;QAC9C,CAAC;IACH,CAAC;IACD,uEAAuE;IACvE,qCAAqC;IACrC,OAAO,SAAS,CAAC;AACnB,CAAC"}
+{"version":3,"file":"getContourMoore.js","sourceRoot":"","sources":["../../src/maskAnalysis/getContourMoore.ts"],"names":[],"mappings":"AAEA,OAAO,EAAE,MAAM,EAAE,MAAM,+BAA+B,CAAC;AAEvD,OAAO,EACL,UAAU,EACV,kBAAkB,GACnB,MAAM,oCAAoC,CAAC;AAY5C;;GAEG;AACH,MAAM,UAAU,GAAgB;IAC9B,EAAE,EAAE,EAAE,CAAC,CAAC,EAAE,EAAE,EAAE,CAAC,EAAE,EAAE,OAAO;IAC1B,EAAE,EAAE,EAAE,CAAC,CAAC,EAAE,EAAE,EAAE,CAAC,CAAC,EAAE,EAAE,UAAU;IAC9B,EAAE,EAAE,EAAE,CAAC,EAAE,EAAE,EAAE,CAAC,CAAC,EAAE,EAAE,KAAK;IACxB,EAAE,EAAE,EAAE,CAAC,EAAE,EAAE,EAAE,CAAC,CAAC,EAAE,EAAE,WAAW;IAC9B,EAAE,EAAE,EAAE,CAAC,EAAE,EAAE,EAAE,CAAC,EAAE,EAAE,QAAQ;IAC1B,EAAE,EAAE,EAAE,CAAC,EAAE,EAAE,EAAE,CAAC,EAAE,EAAE,aAAa;IAC/B,EAAE,EAAE,EAAE,CAAC,EAAE,EAAE,EAAE,CAAC,EAAE,EAAE,OAAO;IACzB,EAAE,EAAE,EAAE,CAAC,CAAC,EAAE,EAAE,EAAE,CAAC,EAAE,EAAE,YAAY;CAChC,CAAC;AACF,sCAAsC;AACtC,MAAM,oBAAoB,GAAG,IAAI,SAAS,CAAC,CAAC,CAAC,CAAC,IAAI,CAAC,CAAC,CAAC,CAAC,CAAC;AACvD,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,UAAU,CAAC,MAAM,EAAE,CAAC,EAAE,EAAE,CAAC;IAC3C,MAAM,EAAE,EAAE,EAAE,EAAE,EAAE,GAAG,UAAU,CAAC,CAAC,CAAC,CAAC;IACjC,oBAAoB,CAAC,CAAC,EAAE,GAAG,CAAC,CAAC,GAAG,CAAC,GAAG,CAAC,EAAE,GAAG,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC;AACpD,CAAC;AACD;;;;;GAKG;AACH,MAAM,UAAU,eAAe,CAAC,IAAU;IACxC,IAAI,KAAK,GAAG,CAAC,CAAC;IAEd,MAAM,OAAO,GAAY,EAAE,CAAC;IAC5B,MAAM,OAAO,GAAG,IAAI,UAAU,CAAC,IAAI,CAAC,KAAK,GAAG,IAAI,CAAC,MAAM,CAAC,CAAC;IACzD,OAAO,IAAI,CAAC,aAAa,CAAC,KAAK,CAAC,KAAK,CAAC,IAAI,KAAK,KAAK,IAAI,CAAC,IAAI,EAAE,CAAC;QAC9D,KAAK,EAAE,CAAC;IACV,CAAC;IAED,IAAI,KAAK,KAAK,IAAI,CAAC,IAAI;QAAE,OAAO,OAAO,CAAC;IACxC,MAAM,CAAC,GAAG,EAAE,GAAG,CAAC,GAAG,kBAAkB,CAAC,KAAK,EAAE,IAAI,CAAC,KAAK,CAAC,CAAC;IACzD,MAAM,aAAa,GAAG,EAAE,MAAM,EAAE,GAAG,EAAE,GAAG,EAAE,CAAC;IAC3C,IAAI,YAAY,GAAG,aAAa,CAAC;IACjC,IAAI,CAAC,OAAO,CAAC,KAAK,CAAC,EAAE,CAAC;QACpB,OAAO,CAAC,KAAK,CAAC,GAAG,CAAC,CAAC;QACnB,OAAO,CAAC,IAAI,CAAC,YAAY,CAAC,CAAC;IAC7B,CAAC;IAED,MAAM,SAAS,GAAG,aAAa,CAAC,IAAI,EAAE,YAAY,EAAE,EAAE,MAAM,EAAE,GAAG,GAAG,CAAC,EAAE,GAAG,EAAE,CAAC,CAAC;IAC9E,IAAI,CAAC,SAAS,EAAE,CAAC;QACf,OAAO,OAAO,CAAC;IACjB,CAAC;IAED,IAAI,sBAAsB,GAAU,SAAS,CAAC,SAAS,CAAC;IAExD,IAAI,cAAc,GAAG,SAAS,CAAC,SAAS,CAAC;IACzC,YAAY,GAAG,SAAS,CAAC,SAAS,CAAC;IAEnC,OAAO,IAAI,EAAE,CAAC;QACZ,KAAK,GAAG,YAAY,CAAC,GAAG,GAAG,IAAI,CAAC,KAAK,GAAG,YAAY,CAAC,MAAM,CAAC;QAC5D,IAAI,CAAC,OAAO,CAAC,KAAK,CAAC,EAAE,CAAC;YACpB,OAAO,CAAC,KAAK,CAAC,GAAG,CAAC,CAAC;YACnB,OAAO,CAAC,IAAI,CAAC,YAAY,CAAC,CAAC;QAC7B,CAAC;QAED,MAAM,IAAI,GAAG,aAAa,CAAC,IAAI,EAAE,YAAY,EAAE,cAAc,CAAC,CAAC;QAC/D,MAAM,CAAC,IAAI,EAAE,iCAAiC,CAAC,CAAC;QAChD,cAAc,GAAG,IAAI,CAAC,SAAS,CAAC;QAChC,YAAY,GAAG,IAAI,CAAC,SAAS,CAAC;QAC9B,IACE,YAAY,CAAC,MAAM,KAAK,aAAa,CAAC,MAAM;YAC5C,YAAY,CAAC,GAAG,KAAK,aAAa,CAAC,GAAG,EACtC,CAAC;YACD,IACE,cAAc,CAAC,MAAM,KAAK,sBAAsB,CAAC,MAAM;gBACvD,cAAc,CAAC,GAAG,KAAK,sBAAsB,CAAC,GAAG,EACjD,CAAC;gBACD,MAAM;YACR,CAAC;iBAAM,CAAC;gBACN,sBAAsB,GAAG,cAAc,CAAC;YAC1C,CAAC;QACH,CAAC;IACH,CAAC;IAED,OAAO,OAAO,CAAC;AACjB,CAAC;AACD;;;;;;GAMG;AACH,SAAS,aAAa,CACpB,IAAU,EACV,OAAc,EACd,QAAe;IAEf,IAAI,SAAS,GAAG,QAAQ,CAAC;IAEzB,uCAAuC;IACvC,MAAM,MAAM,GAAG,QAAQ,CAAC,MAAM,GAAG,OAAO,CAAC,MAAM,CAAC;IAChD,MAAM,MAAM,GAAG,QAAQ,CAAC,GAAG,GAAG,OAAO,CAAC,GAAG,CAAC;IAE1C,MAAM,UAAU,GAAG,oBAAoB,CAAC,CAAC,MAAM,GAAG,CAAC,CAAC,GAAG,CAAC,GAAG,CAAC,MAAM,GAAG,CAAC,CAAC,CAAC,CAAC;IAEzE,2CAA2C;IAC3C,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,UAAU,CAAC,MAAM,EAAE,CAAC,EAAE,EAAE,CAAC;QAC3C,MAAM,EAAE,EAAE,EAAE,EAAE,EAAE,GAAG,UAAU,CAAC,EAAE,CAC9B,CAAC,UAAU,GAAG,CAAC,CAAC,GAAG,UAAU,CAAC,MAAM,CACxB,CAAC;QACf,MAAM,MAAM,GAAG,OAAO,CAAC,MAAM,GAAG,EAAE,CAAC;QACnC,MAAM,MAAM,GAAG,OAAO,CAAC,GAAG,GAAG,EAAE,CAAC;QAChC,IAAI,UAAU,CAAC,IAAI,EAAE,MAAM,EAAE,MAAM,CAAC,KAAK,CAAC,EAAE,CAAC;YAC3C,2DAA2D;YAC3D,0DAA0D;YAC1D,aAAa;YACb,OAAO;gBACL,SAAS;gBACT,SAAS,EAAE,EAAE,MAAM,EAAE,MAAM,EAAE,GAAG,EAAE,MAAM,EAAE;aAC3C,CAAC;QACJ,CAAC;aAAM,CAAC;YACN,SAAS,GAAG,EAAE,MAAM,EAAE,MAAM,EAAE,GAAG,EAAE,MAAM,EAAE,CAAC;QAC9C,CAAC;IACH,CAAC;IACD,uEAAuE;IACvE,qCAAqC;IACrC,OAAO,SAAS,CAAC;AACnB,CAAC"}

package/lib/maskAnalysis/getConvexityDefects.d.ts

@@ -0,0 +1,17 @@
+import type { Point } from '../index_full.ts';
+export interface ConvexityDefectsOptions {
+    /**
+     * Depth threshold to remove smaller convexity defaults.
+     * @default 0
+     */
+    depthThreshold?: number;
+}
+/**
+ * Find border points that are farthest from its convex hull lines. Inspired by openCV's convexityDefects.
+ * @param borderPoints - Mask's border points.
+ * @param hullPoints - Mask's convex hull points.
+ * @param options - Convexity defects options.
+ * @returns Array of convexity defects.
+ */
+export declare function getConvexityDefects(borderPoints: Point[], hullPoints: Point[], options?: ConvexityDefectsOptions): Point[];
+//# sourceMappingURL=getConvexityDefects.d.ts.map

package/lib/maskAnalysis/getConvexityDefects.d.ts.map

@@ -0,0 +1 @@
+{"version":3,"file":"getConvexityDefects.d.ts","sourceRoot":"","sources":["../../src/maskAnalysis/getConvexityDefects.ts"],"names":[],"mappings":"AAAA,OAAO,KAAK,EAAE,KAAK,EAAE,MAAM,kBAAkB,CAAC;AAE9C,MAAM,WAAW,uBAAuB;IACtC;;;OAGG;IACH,cAAc,CAAC,EAAE,MAAM,CAAC;CACzB;AAED;;;;;;GAMG;AACH,wBAAgB,mBAAmB,CACjC,YAAY,EAAE,KAAK,EAAE,EACrB,UAAU,EAAE,KAAK,EAAE,EACnB,OAAO,GAAE,uBAA4B,GACpC,KAAK,EAAE,CAsET"}

package/lib/maskAnalysis/getConvexityDefects.js

@@ -0,0 +1,75 @@
+/**
+ * Find border points that are farthest from its convex hull lines. Inspired by openCV's convexityDefects.
+ * @param borderPoints - Mask's border points.
+ * @param hullPoints - Mask's convex hull points.
+ * @param options - Convexity defects options.
+ * @returns Array of convexity defects.
+ */
+export function getConvexityDefects(borderPoints, hullPoints, options = {}) {
+    const { depthThreshold = 0 } = options;
+    const defects = [];
+    if (hullPoints.length === 0) {
+        throw new RangeError('No hull points were defined for convexity defects detection.');
+    }
+    if (borderPoints.length === 0) {
+        throw new RangeError('No border points were defined for convexity defects detection.');
+    }
+    const lastHullPoint = hullPoints.at(-1);
+    let currBorderIndex = borderPoints.findIndex((bp) => limitReached(bp, lastHullPoint));
+    if (currBorderIndex === -1) {
+        throw new RangeError('Could not find a border point matching the convex hull endpoint.');
+    }
+    for (let i = hullPoints.length - 1; i >= 0; i--) {
+        const currHullPoint = hullPoints.at(i);
+        const nextHullPoint = hullPoints.at(i - 1);
+        if (!currHullPoint || !nextHullPoint)
+            continue;
+        const v0x = nextHullPoint.column - currHullPoint.column;
+        const v0y = nextHullPoint.row - currHullPoint.row;
+        const edgeLenSq = v0x * v0x + v0y * v0y;
+        let maxDepth = 0;
+        let maxDepthIndex = -1;
+        const segmentStart = currBorderIndex;
+        let reached = false;
+        for (let j = 0; j < borderPoints.length; j++) {
+            const checkIndex = (segmentStart + j) % borderPoints.length;
+            const bp = borderPoints[checkIndex];
+            if (limitReached(bp, nextHullPoint)) {
+                currBorderIndex = (checkIndex + 1) % borderPoints.length;
+                reached = true;
+                break;
+            }
+            const v1x = bp.column - currHullPoint.column;
+            const v1y = bp.row - currHullPoint.row;
+            const t = Math.max(0, Math.min(1, (v1x * v0x + v1y * v0y) / edgeLenSq));
+            const closestX = currHullPoint.column + t * v0x;
+            const closestY = currHullPoint.row + t * v0y;
+            const depth = Math.hypot(bp.column - closestX, bp.row - closestY);
+            if (depth > maxDepth) {
+                maxDepth = depth;
+                maxDepthIndex = checkIndex;
+            }
+        }
+        if (!reached) {
+            throw new RangeError('Could not reach the next hull point while scanning border points; hull and border may be inconsistent.');
+        }
+        if (maxDepth > depthThreshold) {
+            const defectPoint = borderPoints[maxDepthIndex];
+            if (defectPoint)
+                defects.push(defectPoint);
+        }
+    }
+    return defects;
+}
+/**
+ * Checks if borderPoint reached convex hull's point.
+ * @param borderPoint - Border point.
+ * @param hullPoint - Endpoint of convex hull segment.
+ * @returns whether endpoint is reached.
+ */
+function limitReached(borderPoint, hullPoint) {
+    const tolerance = 1;
+    return (Math.abs(borderPoint.column - hullPoint.column) <= tolerance &&
+        Math.abs(borderPoint.row - hullPoint.row) <= tolerance);
+}
+//# sourceMappingURL=getConvexityDefects.js.map

package/lib/maskAnalysis/getConvexityDefects.js.map

@@ -0,0 +1 @@
+{"version":3,"file":"getConvexityDefects.js","sourceRoot":"","sources":["../../src/maskAnalysis/getConvexityDefects.ts"],"names":[],"mappings":"AAUA;;;;;;GAMG;AACH,MAAM,UAAU,mBAAmB,CACjC,YAAqB,EACrB,UAAmB,EACnB,UAAmC,EAAE;IAErC,MAAM,EAAE,cAAc,GAAG,CAAC,EAAE,GAAG,OAAO,CAAC;IACvC,MAAM,OAAO,GAAG,EAAE,CAAC;IAEnB,IAAI,UAAU,CAAC,MAAM,KAAK,CAAC,EAAE,CAAC;QAC5B,MAAM,IAAI,UAAU,CAClB,8DAA8D,CAC/D,CAAC;IACJ,CAAC;IACD,IAAI,YAAY,CAAC,MAAM,KAAK,CAAC,EAAE,CAAC;QAC9B,MAAM,IAAI,UAAU,CAClB,gEAAgE,CACjE,CAAC;IACJ,CAAC;IAED,MAAM,aAAa,GAAG,UAAU,CAAC,EAAE,CAAC,CAAC,CAAC,CAAU,CAAC;IAEjD,IAAI,eAAe,GAAG,YAAY,CAAC,SAAS,CAAC,CAAC,EAAE,EAAE,EAAE,CAClD,YAAY,CAAC,EAAE,EAAE,aAAa,CAAC,CAChC,CAAC;IAEF,IAAI,eAAe,KAAK,CAAC,CAAC,EAAE,CAAC;QAC3B,MAAM,IAAI,UAAU,CAClB,kEAAkE,CACnE,CAAC;IACJ,CAAC;IACD,KAAK,IAAI,CAAC,GAAG,UAAU,CAAC,MAAM,GAAG,CAAC,EAAE,CAAC,IAAI,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC;QAChD,MAAM,aAAa,GAAG,UAAU,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC;QACvC,MAAM,aAAa,GAAG,UAAU,CAAC,EAAE,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC;QAC3C,IAAI,CAAC,aAAa,IAAI,CAAC,aAAa;YAAE,SAAS;QAC/C,MAAM,GAAG,GAAG,aAAa,CAAC,MAAM,GAAG,aAAa,CAAC,MAAM,CAAC;QACxD,MAAM,GAAG,GAAG,aAAa,CAAC,GAAG,GAAG,aAAa,CAAC,GAAG,CAAC;QAClD,MAAM,SAAS,GAAG,GAAG,GAAG,GAAG,GAAG,GAAG,GAAG,GAAG,CAAC;QACxC,IAAI,QAAQ,GAAG,CAAC,CAAC;QACjB,IAAI,aAAa,GAAG,CAAC,CAAC,CAAC;QACvB,MAAM,YAAY,GAAG,eAAe,CAAC;QACrC,IAAI,OAAO,GAAG,KAAK,CAAC;QACpB,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,YAAY,CAAC,MAAM,EAAE,CAAC,EAAE,EAAE,CAAC;YAC7C,MAAM,UAAU,GAAG,CAAC,YAAY,GAAG,CAAC,CAAC,GAAG,YAAY,CAAC,MAAM,CAAC;YAC5D,MAAM,EAAE,GAAG,YAAY,CAAC,UAAU,CAAC,CAAC;YACpC,IAAI,YAAY,CAAC,EAAE,EAAE,aAAa,CAAC,EAAE,CAAC;gBACpC,eAAe,GAAG,CAAC,UAAU,GAAG,CAAC,CAAC,GAAG,YAAY,CAAC,MAAM,CAAC;gBACzD,OAAO,GAAG,IAAI,CAAC;gBACf,MAAM;YACR,CAAC;YAED,MAAM,GAAG,GAAG,EAAE,CAAC,MAAM,GAAG,aAAa,CAAC,MAAM,CAAC;YAC7C,MAAM,GAAG,GAAG,EAAE,CAAC,GAAG,GAAG,aAAa,CAAC,GAAG,CAAC;YAEvC,MAAM,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,CAAC,EAAE,IAAI,CAAC,GAAG,CAAC,CAAC,EAAE,CAAC,GAAG,GAAG,GAAG,GAAG,GAAG,GAAG,GAAG,CAAC,GAAG,SAAS,CAAC,CAAC,CAAC;YACxE,MAAM,QAAQ,GAAG,aAAa,CAAC,MAAM,GAAG,CAAC,GAAG,GAAG,CAAC;YAChD,MAAM,QAAQ,GAAG,aAAa,CAAC,GAAG,GAAG,CAAC,GAAG,GAAG,CAAC;YAC7C,MAAM,KAAK,GAAG,IAAI,CAAC,KAAK,CAAC,EAAE,CAAC,MAAM,GAAG,QAAQ,EAAE,EAAE,CAAC,GAAG,GAAG,QAAQ,CAAC,CAAC;YAClE,IAAI,KAAK,GAAG,QAAQ,EAAE,CAAC;gBACrB,QAAQ,GAAG,KAAK,CAAC;gBACjB,aAAa,GAAG,UAAU,CAAC;YAC7B,CAAC;QACH,CAAC;QACD,IAAI,CAAC,OAAO,EAAE,CAAC;YACb,MAAM,IAAI,UAAU,CAClB,wGAAwG,CACzG,CAAC;QACJ,CAAC;QACD,IAAI,QAAQ,GAAG,cAAc,EAAE,CAAC;YAC9B,MAAM,WAAW,GAAG,YAAY,CAAC,aAAa,CAAC,CAAC;YAChD,IAAI,WAAW;gBAAE,OAAO,CAAC,IAAI,CAAC,WAAW,CAAC,CAAC;QAC7C,CAAC;IACH,CAAC;IAED,OAAO,OAAO,CAAC;AACjB,CAAC;AACD;;;;;GAKG;AACH,SAAS,YAAY,CAAC,WAAkB,EAAE,SAAgB;IACxD,MAAM,SAAS,GAAG,CAAC,CAAC;IACpB,OAAO,CACL,IAAI,CAAC,GAAG,CAAC,WAAW,CAAC,MAAM,GAAG,SAAS,CAAC,MAAM,CAAC,IAAI,SAAS;QAC5D,IAAI,CAAC,GAAG,CAAC,WAAW,CAAC,GAAG,GAAG,SAAS,CAAC,GAAG,CAAC,IAAI,SAAS,CACvD,CAAC;AACJ,CAAC"}

package/lib/maskAnalysis/getExternalContour.d.ts

@@ -7,5 +7,5 @@ export type GetExternalContourOptions = Omit<GetBorderPointsOptions, 'innerBorer
  * @param options - GetExternalContourOptions.
  * @returns Array of contour points.
  */
-export declare function getExternalContour(mask: Mask, options?: GetExternalContourOptions): import("../index_full.ts").Point[];
+export declare function getExternalContour(mask: Mask, options?: GetExternalContourOptions): import("../index.ts").Point[];
 //# sourceMappingURL=getExternalContour.d.ts.map

package/lib/maskAnalysis/getExternalContour.d.ts.map

@@ -1 +1 @@
-{"version":3,"file":"getExternalContour.d.ts","sourceRoot":"","sources":["../../src/maskAnalysis/getExternalContour.ts"],"names":[],"mappings":"AAAA,OAAO,KAAK,EAAE,IAAI,EAAE,MAAM,YAAY,CAAC;AAIvC,OAAO,KAAK,EAAE,sBAAsB,EAAE,MAAM,yBAAyB,CAAC;AAEtE,MAAM,MAAM,yBAAyB,GAAG,IAAI,CAC1C,sBAAsB,EACtB,aAAa,CACd,CAAC;AACF;;;;;GAKG;AACH,wBAAgB,kBAAkB,CAChC,IAAI,EAAE,IAAI,EACV,OAAO,GAAE,yBAA8B,sCAQxC"}
+{"version":3,"file":"getExternalContour.d.ts","sourceRoot":"","sources":["../../src/maskAnalysis/getExternalContour.ts"],"names":[],"mappings":"AAAA,OAAO,KAAK,EAAE,IAAI,EAAE,MAAM,YAAY,CAAC;AAIvC,OAAO,KAAK,EAAE,sBAAsB,EAAE,MAAM,yBAAyB,CAAC;AAEtE,MAAM,MAAM,yBAAyB,GAAG,IAAI,CAC1C,sBAAsB,EACtB,aAAa,CACd,CAAC;AACF;;;;;GAKG;AACH,wBAAgB,kBAAkB,CAChC,IAAI,EAAE,IAAI,EACV,OAAO,GAAE,yBAA8B,iCAQxC"}

package/lib/maskAnalysis/index.d.ts

@@ -1,2 +1,3 @@
 export * from './maskAnalysis.types.js';
+export * from './getConvexityDefects.ts';
 //# sourceMappingURL=index.d.ts.map

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

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

package/lib/maskAnalysis/index.js

@@ -1,2 +1,3 @@
 export * from './maskAnalysis.types.js';
+export * from "./getConvexityDefects.js";
 //# sourceMappingURL=index.js.map

package/lib/maskAnalysis/index.js.map

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

package/lib/roi/getExternalContour.d.ts

@@ -6,5 +6,5 @@ import type { Roi } from './Roi.ts';
  * @param options - GetExternalContourOptions.
  * @returns Array of contour points.
  */
-export declare function getExternalContour(roi: Roi, options?: GetExternalContourOptions): import("../index_full.ts").Point[];
+export declare function getExternalContour(roi: Roi, options?: GetExternalContourOptions): import("../index.ts").Point[];
 //# sourceMappingURL=getExternalContour.d.ts.map

package/lib/roi/getExternalContour.d.ts.map

@@ -1 +1 @@
-{"version":3,"file":"getExternalContour.d.ts","sourceRoot":"","sources":["../../src/roi/getExternalContour.ts"],"names":[],"mappings":"AAAA,OAAO,KAAK,EAAE,yBAAyB,EAAE,MAAM,uCAAuC,CAAC;AAEvF,OAAO,KAAK,EAAE,GAAG,EAAE,MAAM,UAAU,CAAC;AAEpC;;;;;GAKG;AACH,wBAAgB,kBAAkB,CAChC,GAAG,EAAE,GAAG,EACR,OAAO,CAAC,EAAE,yBAAyB,sCAIpC"}
+{"version":3,"file":"getExternalContour.d.ts","sourceRoot":"","sources":["../../src/roi/getExternalContour.ts"],"names":[],"mappings":"AAAA,OAAO,KAAK,EAAE,yBAAyB,EAAE,MAAM,uCAAuC,CAAC;AAEvF,OAAO,KAAK,EAAE,GAAG,EAAE,MAAM,UAAU,CAAC;AAEpC;;;;;GAKG;AACH,wBAAgB,kBAAkB,CAChC,GAAG,EAAE,GAAG,EACR,OAAO,CAAC,EAAE,yBAAyB,iCAIpC"}

package/package.json

@@ -1,6 +1,6 @@
 {
   "name": "image-js",
-  "version": "1.6.1",
+  "version": "1.7.0",
   "license": "MIT",
   "description": "Image processing and manipulation in JavaScript",
   "author": "Michaël Zasso",
@@ -108,7 +108,7 @@
     "rimraf": "^6.1.2",
     "tailwindcss": "^4.1.18",
     "typescript": "~5.9.3",
-    "vite": "^7.3.0",
+    "vite": "^8.1.0",
     "vitest": "^4.0.16"
   },
   "repository": {

package/src/Stack.ts

@@ -1,6 +1,4 @@
-import type { BitDepth } from 'fast-png';
-
-import type { Image } from './Image.js';
+import type { BitDepth, Image } from './Image.js';
 import type { HistogramOptions } from './compute/index.js';
 import {
   histogram,

package/src/align/affineTransfrom/getAffineTransform.ts

@@ -60,9 +60,8 @@ export interface GetAffineTransformOptions {
   debug?: boolean;
   /**
    * Path of the debug image.
-   * @default `${import.meta.dirname}/montage.png`
    */
-  debugImagePath?: string;
+  debugImagePath?: string | URL;
 }
 
 export interface AffineTransform {
@@ -131,7 +130,7 @@ export function getAffineTransform(
     destinationOrigin = { column: 0, row: 0 },
     maxRansacNbIterations,
     debug = false,
-    debugImagePath = `${import.meta.dirname}/montage.png`,
+    debugImagePath,
   } = options;
 
   if (source.colorModel !== ImageColorModel.GREY) {
@@ -242,6 +241,9 @@ export function getAffineTransform(
 
   // create debug image
   if (debug) {
+    if (!debugImagePath) {
+      throw new RangeError('Debug image file path is not specified.');
+    }
     const montage = new Montage(source, destination, {
       disposition: MontageDisposition.VERTICAL,
     });

package/src/featureMatching/keypoints/getEigenvaluesForScore.ts

@@ -1,51 +1,53 @@
-import { EigenvalueDecomposition, Matrix, WrapperMatrix1D } from 'ml-matrix';
+import { EigenvalueDecomposition, Matrix } from 'ml-matrix';
 
 import type { Image } from '../../Image.ts';
+import { rawDirectConvolution } from '../../filters/convolution.ts';
 import type { Point } from '../../index_full.ts';
 import { SOBEL_X, SOBEL_Y } from '../../utils/constants/kernels.js';
-
 /**
  * A function that calculates eigenvalues to calculate feature score for Harris and Shi-Tomasi algorithms.
  * @param image - Image take data from.
  * @param origin - Center of the window, where the corner should be.
- * @param windowSize - Size of the window, where data should be scanned.
+ * @param cropSize - Size of the window, where data should be scanned.
  * @returns Array of two eigenvalues.
  */
 export function getEigenvaluesForScore(
   image: Image,
   origin: Point,
-  windowSize = 7,
+  cropSize = 5,
 ) {
-  if (!(windowSize % 2)) {
+  if (!(cropSize % 2)) {
     throw new TypeError('windowSize must be an odd integer');
   }
-
+  const kernelRadius = (SOBEL_X.length - 1) / 2;
+  const windowRadius = (cropSize - 1) / 2;
+  const padded = cropSize + 2 * kernelRadius;
   const cropOrigin = {
-    row: origin.row - (windowSize - 1) / 2,
-    column: origin.column - (windowSize - 1) / 2,
+    row: origin.row - windowRadius - kernelRadius,
+    column: origin.column - windowRadius - kernelRadius,
   };
   const window = image.crop({
     origin: cropOrigin,
-    width: windowSize,
-    height: windowSize,
-  });
-  const xDerivative = window.gradientFilter({ kernelX: SOBEL_X });
-  const yDerivative = window.gradientFilter({ kernelY: SOBEL_Y });
-
-  const xMatrix = new WrapperMatrix1D(xDerivative.getRawImage().data, {
-    rows: xDerivative.height,
+    width: padded,
+    height: padded,
   });
-  const yMatrix = new WrapperMatrix1D(yDerivative.getRawImage().data, {
-    rows: yDerivative.height,
-  });
-
-  const xx = xMatrix.mmul(xMatrix);
-  const xy = yMatrix.mmul(xMatrix);
-  const yy = yMatrix.mmul(yMatrix);
-
-  const xxSum = xx.sum();
-  const xySum = xy.sum();
-  const yySum = yy.sum();
+  const xDerivative = rawDirectConvolution(window, SOBEL_X);
+  const yDerivative = rawDirectConvolution(window, SOBEL_Y);
+
+  let xxSum = 0;
+  let xySum = 0;
+  let yySum = 0;
+
+  for (let i = kernelRadius; i < window.height - kernelRadius; i++) {
+    for (let j = kernelRadius; j < window.width - kernelRadius; j++) {
+      const idx = i * window.width + j;
+      const gx = xDerivative[idx];
+      const gy = yDerivative[idx];
+      xxSum += gx * gx;
+      xySum += gx * gy;
+      yySum += gy * gy;
+    }
+  }
 
   const structureTensor = new Matrix([
     [xxSum, xySum],

package/src/featureMatching/keypoints/getFastKeypoints.ts

@@ -41,6 +41,10 @@ export interface GetFastKeypointsOptions extends IsFastKeypointOptions {
    * @default `'FAST'`
    */
   scoreAlgorithm?: 'HARRIS' | 'FAST' | 'TOMASI';
+  /**
+   * Value between 0 and 1 that signifies a percentage from maximum score found in the image. Points can be filtered based on this relative maximum score.
+   */
+  qualityThreshold?: number;
   /**
    * Options for the Harris score computation.
    */
@@ -72,7 +76,7 @@ export function getFastKeypoints(
   options: GetFastKeypointsOptions = {},
 ): FastKeypoint[] {
   const { fastRadius = 3, scoreAlgorithm = 'FAST', scoreOptions } = options;
-
+  const windowSize = scoreOptions?.windowSize ?? 7;
   const circlePoints = getCirclePoints(fastRadius);
   const compassPoints = getCompassPoints(fastRadius);
 
@@ -81,6 +85,7 @@ export function getFastKeypoints(
     nbContiguousPixels = (3 / 4) * circlePoints.length,
     threshold = 20,
     nonMaxSuppression = true,
+    qualityThreshold = 0,
   } = options;
 
   checkProcessable(image, {
@@ -108,13 +113,14 @@ export function getFastKeypoints(
     .with('TOMASI', () => tomasiScore)
     .exhaustive();
 
+  const padding = scoreAlgorithm === 'FAST' ? 0 : (windowSize + 1) / 2;
   const allKeypoints: FastKeypoint[] = [];
 
   const scoreArray = new Float64Array(image.size).fill(
     Number.NEGATIVE_INFINITY,
   );
-  for (let row = 0; row < image.height; row++) {
-    for (let column = 0; column < image.width; column++) {
+  for (let row = padding; row < image.height - padding; row++) {
+    for (let column = padding; column < image.width - padding; column++) {
       const corner = { row, column };
       if (
         isFastKeypoint(corner, image, circlePoints, compassPoints, {
@@ -157,8 +163,20 @@ export function getFastKeypoints(
       }
     }
   }
+  let maxScore = -Infinity;
+  for (const k of keypoints) {
+    if (k.score > maxScore) maxScore = k.score;
+  }
+  const filtered: FastKeypoint[] = [];
+  // threshold by default negative infinity
+  const filterThreshold = qualityThreshold * maxScore;
 
-  keypoints.sort((a, b) => b.score - a.score);
+  for (const k of keypoints) {
+    if (k.score >= filterThreshold) {
+      k.score /= maxScore;
+      filtered.push(k);
+    }
+  }
 
-  return keypoints.slice(0, maxNbFeatures);
+  return filtered.toSorted((a, b) => b.score - a.score).slice(0, maxNbFeatures);
 }

package/src/featureMatching/keypoints/getHarrisScore.ts

@@ -24,7 +24,7 @@ export function getHarrisScore(
   origin: Point,
   options: GetHarrisScoreOptions = {},
 ): number {
-  const { harrisConstant = 0.04, windowSize } = options;
+  const { harrisConstant = 0.04, windowSize = 5 } = options;
   const eigenValues = getEigenvaluesForScore(image, origin, windowSize);
 
   return (

package/src/featureMatching/keypoints/getShiTomasiScore.ts

@@ -4,11 +4,10 @@ import type { Point } from '../../index_full.ts';
 import { getEigenvaluesForScore } from './getEigenvaluesForScore.js';
 
 export interface GetShiTomasiScoreOptions {
-  qualityLevel?: number;
   /**
    * Size of the window to compute the Harris score.
    * Should be an odd number so that the window can be centered on the corner.
-   * @default `7`
+   * @default `5`
    */
   windowSize?: number;
 }
@@ -29,7 +28,7 @@ export function getShiTomasiScore(
   origin: Point,
   options: GetShiTomasiScoreOptions = {},
 ): number {
-  const { windowSize = 7 } = options;
+  const { windowSize = 5 } = options;
 
   const eigenValues = getEigenvaluesForScore(image, origin, windowSize);
   return Math.min(eigenValues[0], eigenValues[1]);

package/src/featureMatching/visualize/drawKeypoints.ts

@@ -1,5 +1,7 @@
 import type { Image } from '../../Image.js';
 import type { Point } from '../../geometry/index.js';
+import type { DrawTextOptions } from '../../index_full.ts';
+import { drawText } from '../../index_full.ts';
 import { sum } from '../../utils/geometry/points.js';
 import { getOutputImage } from '../../utils/getOutputImage.js';
 import type { GetColorsOptions } from '../featureMatching.types.js';
@@ -47,7 +49,7 @@ export interface DrawKeypointsOptions {
   /**
    * Options for the coloring of the keypoints depending on their score (useful if showScore = true).
    */
-  showScoreOptions?: GetColorsOptions;
+  showScoreOptions?: GetColorsOptions & DrawTextOptions;
 }
 
 export interface DrawOrientedKeypointsOptions extends DrawKeypointsOptions {
@@ -89,7 +91,8 @@ export function drawKeypoints(
   } = options;
   let { maxNbKeypoints = keypoints.length } = options;
   const { strokeColor = [255, 0, 0] } = options;
-
+  const fontColor = showScoreOptions?.fontColor ?? [0, 0, 255];
+  const font = showScoreOptions?.font ?? '20px Helvetica';
   if (maxNbKeypoints > keypoints.length) {
     maxNbKeypoints = keypoints.length;
   }
@@ -118,6 +121,7 @@ export function drawKeypoints(
       strokeColor: keypointColor,
       out: newImage,
     });
+
     if (
       isOrientedFastKeypoint(keypoint) &&
       (options as DrawOrientedKeypointsOptions).showOrientation
@@ -135,7 +139,20 @@ export function drawKeypoints(
       });
     }
   }
-
+  if (showScore) {
+    drawText(
+      newImage,
+      keypoints.map((k) => ({
+        content: (k.score * 100).toFixed(2),
+        position: sum(k.origin, origin),
+        font,
+        fontColor: [...fontColor, 255],
+      })),
+      {
+        out: newImage,
+      },
+    );
+  }
   return newImage;
 }
 

package/src/load/decodeBmp.ts

@@ -4,6 +4,10 @@ import { Image } from '../Image.ts';
 import { Mask } from '../Mask.ts';
 import type { ImageColorModel } from '../utils/constants/colorModels.ts';
 
+import type { Resolution } from './load.types.ts';
+
+type DecodedBmp = ReturnType<typeof decode>;
+
 /**
  * Decode a BMP. See the fast-bmp npm module.
  * @param data - The data to decode.
@@ -32,6 +36,24 @@ export function decodeBmp(data: Uint8Array): Image {
     return new Image(decodedData.width, decodedData.height, {
       colorModel,
       data: decodedData.data as Uint8Array,
+      resolution: getBmpResolution(decodedData),
     });
   }
 }
+
+/**
+ * Gets image's resolution from its parsed data. BMP stores resolution in
+ * pixels per meter.
+ * @param bmp - Parsed BMP image.
+ * @returns Object with resolution data if it is defined.
+ */
+function getBmpResolution(bmp: DecodedBmp): Resolution | undefined {
+  if (!bmp.xPixelsPerMeter || !bmp.yPixelsPerMeter) {
+    return undefined;
+  }
+  return {
+    x: bmp.xPixelsPerMeter,
+    y: bmp.yPixelsPerMeter,
+    unit: 'meter',
+  };
+}

package/src/maskAnalysis/getContourMoore.ts

@@ -1,7 +1,6 @@
-import assert from 'node:assert';
-
 import type { Mask } from '../Mask.ts';
 import type { Point } from '../utils/geometry/points.ts';
+import { assert } from '../utils/validators/assert.ts';
 
 import {
   getBitSafe,
@@ -56,9 +55,10 @@ export function getContourMoore(mask: Mask): Point[] {
   const [row, col] = getCoordsFromIndex(index, mask.width);
   const startingPoint = { column: col, row };
   let currentPoint = startingPoint;
-
-  visited[index] = 1;
-  contour.push(currentPoint);
+  if (!visited[index]) {
+    visited[index] = 1;
+    contour.push(currentPoint);
+  }
 
   const firstNext = findNextPoint(mask, currentPoint, { column: col - 1, row });
   if (!firstNext) {
@@ -78,7 +78,7 @@ export function getContourMoore(mask: Mask): Point[] {
     }
 
     const next = findNextPoint(mask, currentPoint, backtrackPoint);
-    assert.ok(next, 'Next border point is undefined.');
+    assert(next, 'Next border point is undefined.');
     backtrackPoint = next.prevPoint;
     currentPoint = next.currPoint;
     if (

package/src/maskAnalysis/getConvexityDefects.ts

@@ -0,0 +1,105 @@
+import type { Point } from '../index_full.ts';
+
+export interface ConvexityDefectsOptions {
+  /**
+   * Depth threshold to remove smaller convexity defaults.
+   * @default 0
+   */
+  depthThreshold?: number;
+}
+
+/**
+ * Find border points that are farthest from its convex hull lines. Inspired by openCV's convexityDefects.
+ * @param borderPoints - Mask's border points.
+ * @param hullPoints - Mask's convex hull points.
+ * @param options - Convexity defects options.
+ * @returns Array of convexity defects.
+ */
+export function getConvexityDefects(
+  borderPoints: Point[],
+  hullPoints: Point[],
+  options: ConvexityDefectsOptions = {},
+): Point[] {
+  const { depthThreshold = 0 } = options;
+  const defects = [];
+
+  if (hullPoints.length === 0) {
+    throw new RangeError(
+      'No hull points were defined for convexity defects detection.',
+    );
+  }
+  if (borderPoints.length === 0) {
+    throw new RangeError(
+      'No border points were defined for convexity defects detection.',
+    );
+  }
+
+  const lastHullPoint = hullPoints.at(-1) as Point;
+
+  let currBorderIndex = borderPoints.findIndex((bp) =>
+    limitReached(bp, lastHullPoint),
+  );
+
+  if (currBorderIndex === -1) {
+    throw new RangeError(
+      'Could not find a border point matching the convex hull endpoint.',
+    );
+  }
+  for (let i = hullPoints.length - 1; i >= 0; i--) {
+    const currHullPoint = hullPoints.at(i);
+    const nextHullPoint = hullPoints.at(i - 1);
+    if (!currHullPoint || !nextHullPoint) continue;
+    const v0x = nextHullPoint.column - currHullPoint.column;
+    const v0y = nextHullPoint.row - currHullPoint.row;
+    const edgeLenSq = v0x * v0x + v0y * v0y;
+    let maxDepth = 0;
+    let maxDepthIndex = -1;
+    const segmentStart = currBorderIndex;
+    let reached = false;
+    for (let j = 0; j < borderPoints.length; j++) {
+      const checkIndex = (segmentStart + j) % borderPoints.length;
+      const bp = borderPoints[checkIndex];
+      if (limitReached(bp, nextHullPoint)) {
+        currBorderIndex = (checkIndex + 1) % borderPoints.length;
+        reached = true;
+        break;
+      }
+
+      const v1x = bp.column - currHullPoint.column;
+      const v1y = bp.row - currHullPoint.row;
+
+      const t = Math.max(0, Math.min(1, (v1x * v0x + v1y * v0y) / edgeLenSq));
+      const closestX = currHullPoint.column + t * v0x;
+      const closestY = currHullPoint.row + t * v0y;
+      const depth = Math.hypot(bp.column - closestX, bp.row - closestY);
+      if (depth > maxDepth) {
+        maxDepth = depth;
+        maxDepthIndex = checkIndex;
+      }
+    }
+    if (!reached) {
+      throw new RangeError(
+        'Could not reach the next hull point while scanning border points; hull and border may be inconsistent.',
+      );
+    }
+    if (maxDepth > depthThreshold) {
+      const defectPoint = borderPoints[maxDepthIndex];
+      if (defectPoint) defects.push(defectPoint);
+    }
+  }
+
+  return defects;
+}
+/**
+ * Checks if borderPoint reached convex hull's point.
+ * @param borderPoint - Border point.
+ * @param hullPoint - Endpoint of convex hull segment.
+ * @returns whether endpoint is reached.
+ */
+function limitReached(borderPoint: Point, hullPoint: Point): boolean {
+  const tolerance = 1;
+  return (
+    Math.abs(borderPoint.column - hullPoint.column) <= tolerance &&
+    Math.abs(borderPoint.row - hullPoint.row) <= tolerance
+  );
+}

package/src/maskAnalysis/index.ts

@@ -1 +1,2 @@
 export * from './maskAnalysis.types.js';
+export * from './getConvexityDefects.ts';