@cornerstonejs/tools 1.41.0 → 1.42.0

package/dist/cjs/utilities/contours/generateContourSetsFromLabelmap.js

@@ -0,0 +1,124 @@
+"use strict";
+var __importDefault = (this && this.__importDefault) || function (mod) {
+    return (mod && mod.__esModule) ? mod : { "default": mod };
+};
+Object.defineProperty(exports, "__esModule", { value: true });
+exports.generateContourSetsFromLabelmap = void 0;
+const core_1 = require("@cornerstonejs/core");
+const ImageMarchingSquares_1 = __importDefault(require("@kitware/vtk.js/Filters/General/ImageMarchingSquares"));
+const DataArray_1 = __importDefault(require("@kitware/vtk.js/Common/Core/DataArray"));
+const ImageData_1 = __importDefault(require("@kitware/vtk.js/Common/DataModel/ImageData"));
+const mergePoints_1 = require("./mergePoints");
+const contourFinder_1 = require("./contourFinder");
+const SegmentationRepresentations_1 = __importDefault(require("../../enums/SegmentationRepresentations"));
+const { Labelmap } = SegmentationRepresentations_1.default;
+function generateContourSetsFromLabelmap({ segmentations }) {
+    var _a;
+    const { representationData, segments = [0, 1] } = segmentations;
+    const { volumeId: segVolumeId } = representationData[Labelmap];
+    const vol = core_1.cache.getVolume(segVolumeId);
+    if (!vol) {
+        console.warn(`No volume found for ${segVolumeId}`);
+        return;
+    }
+    const numSlices = vol.dimensions[2];
+    const imageVol = core_1.cache.getVolume(vol.referencedVolumeId);
+    if (!imageVol) {
+        console.warn(`No volume found for ${vol.referencedVolumeId}`);
+        return;
+    }
+    const segData = vol.imageData.getPointData().getScalars().getData();
+    const pixelsPerSlice = vol.dimensions[0] * vol.dimensions[1];
+    for (let z = 0; z < numSlices; z++) {
+        for (let y = 0; y < vol.dimensions[1]; y++) {
+            const index = y * vol.dimensions[0] + z * pixelsPerSlice;
+            segData[index] = 0;
+            segData[index + vol.dimensions[0] - 1] = 0;
+        }
+    }
+    const ContourSets = [];
+    const { FrameOfReferenceUID } = imageVol.metadata;
+    const numSegments = segments.length;
+    for (let segIndex = 0; segIndex < numSegments; segIndex++) {
+        const segment = segments[segIndex];
+        if (!segment) {
+            continue;
+        }
+        const sliceContours = [];
+        const scalars = DataArray_1.default.newInstance({
+            name: 'Scalars',
+            numberOfComponents: 1,
+            size: pixelsPerSlice * numSlices,
+            dataType: 'Uint8Array',
+        });
+        const { containedSegmentIndices } = segment;
+        for (let sliceIndex = 0; sliceIndex < numSlices; sliceIndex++) {
+            if (isSliceEmptyForSegment(sliceIndex, segData, pixelsPerSlice, segIndex)) {
+                continue;
+            }
+            const frameStart = sliceIndex * pixelsPerSlice;
+            try {
+                for (let i = 0; i < pixelsPerSlice; i++) {
+                    const value = segData[i + frameStart];
+                    if (value === segIndex || (containedSegmentIndices === null || containedSegmentIndices === void 0 ? void 0 : containedSegmentIndices.has(value))) {
+                        scalars.setValue(i + frameStart, 1);
+                    }
+                    else {
+                        scalars.setValue(i, 0);
+                    }
+                }
+                const mSquares = ImageMarchingSquares_1.default.newInstance({
+                    slice: sliceIndex,
+                });
+                const imageDataCopy = ImageData_1.default.newInstance();
+                imageDataCopy.shallowCopy(vol.imageData);
+                imageDataCopy.getPointData().setScalars(scalars);
+                mSquares.setInputData(imageDataCopy);
+                const cValues = [1];
+                mSquares.setContourValues(cValues);
+                mSquares.setMergePoints(false);
+                const msOutput = mSquares.getOutputData();
+                const reducedSet = (0, mergePoints_1.removeDuplicatePoints)(msOutput);
+                if ((_a = reducedSet.points) === null || _a === void 0 ? void 0 : _a.length) {
+                    const contours = (0, contourFinder_1.findContoursFromReducedSet)(reducedSet.lines);
+                    sliceContours.push({
+                        referencedImageId: imageVol.imageIds[sliceIndex],
+                        contours,
+                        polyData: reducedSet,
+                        FrameNumber: sliceIndex + 1,
+                        sliceIndex,
+                        FrameOfReferenceUID,
+                    });
+                }
+            }
+            catch (e) {
+                console.warn(sliceIndex);
+                console.warn(e);
+            }
+        }
+        const metadata = {
+            referencedImageId: imageVol.imageIds[0],
+            FrameOfReferenceUID,
+        };
+        const ContourSet = {
+            label: segment.label,
+            color: segment.color,
+            metadata,
+            sliceContours,
+        };
+        ContourSets.push(ContourSet);
+    }
+    return ContourSets;
+}
+exports.generateContourSetsFromLabelmap = generateContourSetsFromLabelmap;
+function isSliceEmptyForSegment(sliceIndex, segData, pixelsPerSlice, segIndex) {
+    const startIdx = sliceIndex * pixelsPerSlice;
+    const endIdx = startIdx + pixelsPerSlice;
+    for (let i = startIdx; i < endIdx; i++) {
+        if (segData[i] === segIndex) {
+            return false;
+        }
+    }
+    return true;
+}
+//# sourceMappingURL=generateContourSetsFromLabelmap.js.map

package/dist/cjs/utilities/contours/generateContourSetsFromLabelmap.js.map

@@ -0,0 +1 @@
+{"version":3,"file":"generateContourSetsFromLabelmap.js","sourceRoot":"","sources":["../../../../src/utilities/contours/generateContourSetsFromLabelmap.ts"],"names":[],"mappings":";;;;;;AAAA,8CAAgE;AAChE,gHAA2F;AAC3F,sFAAiE;AACjE,2FAAsE;AAEtE,+CAAsD;AACtD,mDAA6D;AAC7D,0GAAkF;AAElF,MAAM,EAAE,QAAQ,EAAE,GAAG,qCAA2B,CAAC;AAEjD,SAAS,+BAA+B,CAAC,EAAE,aAAa,EAAE;;IACxD,MAAM,EAAE,kBAAkB,EAAE,QAAQ,GAAG,CAAC,CAAC,EAAE,CAAC,CAAC,EAAE,GAAG,aAAa,CAAC;IAChE,MAAM,EAAE,QAAQ,EAAE,WAAW,EAAE,GAAG,kBAAkB,CAAC,QAAQ,CAAC,CAAC;IAG/D,MAAM,GAAG,GAAG,YAAgB,CAAC,SAAS,CAAC,WAAW,CAAC,CAAC;IACpD,IAAI,CAAC,GAAG,EAAE;QACR,OAAO,CAAC,IAAI,CAAC,uBAAuB,WAAW,EAAE,CAAC,CAAC;QACnD,OAAO;KACR;IAED,MAAM,SAAS,GAAG,GAAG,CAAC,UAAU,CAAC,CAAC,CAAC,CAAC;IAGpC,MAAM,QAAQ,GAAG,YAAgB,CAAC,SAAS,CAAC,GAAG,CAAC,kBAAkB,CAAC,CAAC;IACpE,IAAI,CAAC,QAAQ,EAAE;QACb,OAAO,CAAC,IAAI,CAAC,uBAAuB,GAAG,CAAC,kBAAkB,EAAE,CAAC,CAAC;QAC9D,OAAO;KACR;IAID,MAAM,OAAO,GAAG,GAAG,CAAC,SAAS,CAAC,YAAY,EAAE,CAAC,UAAU,EAAE,CAAC,OAAO,EAAE,CAAC;IACpE,MAAM,cAAc,GAAG,GAAG,CAAC,UAAU,CAAC,CAAC,CAAC,GAAG,GAAG,CAAC,UAAU,CAAC,CAAC,CAAC,CAAC;IAE7D,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,SAAS,EAAE,CAAC,EAAE,EAAE;QAClC,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,GAAG,CAAC,UAAU,CAAC,CAAC,CAAC,EAAE,CAAC,EAAE,EAAE;YAC1C,MAAM,KAAK,GAAG,CAAC,GAAG,GAAG,CAAC,UAAU,CAAC,CAAC,CAAC,GAAG,CAAC,GAAG,cAAc,CAAC;YACzD,OAAO,CAAC,KAAK,CAAC,GAAG,CAAC,CAAC;YACnB,OAAO,CAAC,KAAK,GAAG,GAAG,CAAC,UAAU,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC;SAC5C;KACF;IAKD,MAAM,WAAW,GAAG,EAAE,CAAC;IAEvB,MAAM,EAAE,mBAAmB,EAAE,GAAG,QAAQ,CAAC,QAAQ,CAAC;IAElD,MAAM,WAAW,GAAG,QAAQ,CAAC,MAAM,CAAC;IACpC,KAAK,IAAI,QAAQ,GAAG,CAAC,EAAE,QAAQ,GAAG,WAAW,EAAE,QAAQ,EAAE,EAAE;QACzD,MAAM,OAAO,GAAG,QAAQ,CAAC,QAAQ,CAAC,CAAC;QAGnC,IAAI,CAAC,OAAO,EAAE;YACZ,SAAS;SACV;QAED,MAAM,aAAa,GAAG,EAAE,CAAC;QACzB,MAAM,OAAO,GAAG,mBAAY,CAAC,WAAW,CAAC;YACvC,IAAI,EAAE,SAAS;YACf,kBAAkB,EAAE,CAAC;YACrB,IAAI,EAAE,cAAc,GAAG,SAAS;YAChC,QAAQ,EAAE,YAAY;SACvB,CAAC,CAAC;QACH,MAAM,EAAE,uBAAuB,EAAE,GAAG,OAAO,CAAC;QAC5C,KAAK,IAAI,UAAU,GAAG,CAAC,EAAE,UAAU,GAAG,SAAS,EAAE,UAAU,EAAE,EAAE;YAE7D,IACE,sBAAsB,CAAC,UAAU,EAAE,OAAO,EAAE,cAAc,EAAE,QAAQ,CAAC,EACrE;gBACA,SAAS;aACV;YACD,MAAM,UAAU,GAAG,UAAU,GAAG,cAAc,CAAC;YAE/C,IAAI;gBAEF,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,cAAc,EAAE,CAAC,EAAE,EAAE;oBACvC,MAAM,KAAK,GAAG,OAAO,CAAC,CAAC,GAAG,UAAU,CAAC,CAAC;oBACtC,IAAI,KAAK,KAAK,QAAQ,KAAI,uBAAuB,aAAvB,uBAAuB,uBAAvB,uBAAuB,CAAE,GAAG,CAAC,KAAK,CAAC,CAAA,EAAE;wBAC5D,OAAe,CAAC,QAAQ,CAAC,CAAC,GAAG,UAAU,EAAE,CAAC,CAAC,CAAC;qBAC9C;yBAAM;wBACJ,OAAe,CAAC,QAAQ,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC;qBACjC;iBACF;gBAED,MAAM,QAAQ,GAAG,8BAAuB,CAAC,WAAW,CAAC;oBACnD,KAAK,EAAE,UAAU;iBAClB,CAAC,CAAC;gBAIH,MAAM,aAAa,GAAG,mBAAY,CAAC,WAAW,EAAE,CAAC;gBAEjD,aAAa,CAAC,WAAW,CAAC,GAAG,CAAC,SAAS,CAAC,CAAC;gBACzC,aAAa,CAAC,YAAY,EAAE,CAAC,UAAU,CAAC,OAAO,CAAC,CAAC;gBAGjD,QAAQ,CAAC,YAAY,CAAC,aAAa,CAAC,CAAC;gBACrC,MAAM,OAAO,GAAG,CAAC,CAAC,CAAC,CAAC;gBACpB,QAAQ,CAAC,gBAAgB,CAAC,OAAO,CAAC,CAAC;gBACnC,QAAQ,CAAC,cAAc,CAAC,KAAK,CAAC,CAAC;gBAG/B,MAAM,QAAQ,GAAG,QAAQ,CAAC,aAAa,EAAE,CAAC;gBAG1C,MAAM,UAAU,GAAG,IAAA,mCAAqB,EAAC,QAAQ,CAAC,CAAC;gBACnD,IAAI,MAAA,UAAU,CAAC,MAAM,0CAAE,MAAM,EAAE;oBAC7B,MAAM,QAAQ,GAAG,IAAA,0CAA0B,EAAC,UAAU,CAAC,KAAK,CAAC,CAAC;oBAE9D,aAAa,CAAC,IAAI,CAAC;wBACjB,iBAAiB,EAAE,QAAQ,CAAC,QAAQ,CAAC,UAAU,CAAC;wBAChD,QAAQ;wBACR,QAAQ,EAAE,UAAU;wBACpB,WAAW,EAAE,UAAU,GAAG,CAAC;wBAC3B,UAAU;wBACV,mBAAmB;qBACpB,CAAC,CAAC;iBACJ;aACF;YAAC,OAAO,CAAC,EAAE;gBACV,OAAO,CAAC,IAAI,CAAC,UAAU,CAAC,CAAC;gBACzB,OAAO,CAAC,IAAI,CAAC,CAAC,CAAC,CAAC;aACjB;SACF;QAED,MAAM,QAAQ,GAAG;YACf,iBAAiB,EAAE,QAAQ,CAAC,QAAQ,CAAC,CAAC,CAAC;YACvC,mBAAmB;SACpB,CAAC;QAEF,MAAM,UAAU,GAAG;YACjB,KAAK,EAAE,OAAO,CAAC,KAAK;YACpB,KAAK,EAAE,OAAO,CAAC,KAAK;YACpB,QAAQ;YACR,aAAa;SACd,CAAC;QAEF,WAAW,CAAC,IAAI,CAAC,UAAU,CAAC,CAAC;KAC9B;IAED,OAAO,WAAW,CAAC;AACrB,CAAC;AAeQ,0EAA+B;AAbxC,SAAS,sBAAsB,CAAC,UAAU,EAAE,OAAO,EAAE,cAAc,EAAE,QAAQ;IAC3E,MAAM,QAAQ,GAAG,UAAU,GAAG,cAAc,CAAC;IAC7C,MAAM,MAAM,GAAG,QAAQ,GAAG,cAAc,CAAC;IAEzC,KAAK,IAAI,CAAC,GAAG,QAAQ,EAAE,CAAC,GAAG,MAAM,EAAE,CAAC,EAAE,EAAE;QACtC,IAAI,OAAO,CAAC,CAAC,CAAC,KAAK,QAAQ,EAAE;YAC3B,OAAO,KAAK,CAAC;SACd;KACF;IAED,OAAO,IAAI,CAAC;AACd,CAAC"}

package/dist/cjs/utilities/contours/index.d.ts

@@ -0,0 +1,6 @@
+import contourFinder from './contourFinder';
+import mergePoints from './mergePoints';
+import detectContourHoles from './detectContourHoles';
+import { generateContourSetsFromLabelmap } from './generateContourSetsFromLabelmap';
+import AnnotationToPointData from './AnnotationToPointData';
+export { contourFinder, mergePoints, detectContourHoles, generateContourSetsFromLabelmap, AnnotationToPointData, };

package/dist/cjs/utilities/contours/index.js

@@ -0,0 +1,17 @@
+"use strict";
+var __importDefault = (this && this.__importDefault) || function (mod) {
+    return (mod && mod.__esModule) ? mod : { "default": mod };
+};
+Object.defineProperty(exports, "__esModule", { value: true });
+exports.AnnotationToPointData = exports.generateContourSetsFromLabelmap = exports.detectContourHoles = exports.mergePoints = exports.contourFinder = void 0;
+const contourFinder_1 = __importDefault(require("./contourFinder"));
+exports.contourFinder = contourFinder_1.default;
+const mergePoints_1 = __importDefault(require("./mergePoints"));
+exports.mergePoints = mergePoints_1.default;
+const detectContourHoles_1 = __importDefault(require("./detectContourHoles"));
+exports.detectContourHoles = detectContourHoles_1.default;
+const generateContourSetsFromLabelmap_1 = require("./generateContourSetsFromLabelmap");
+Object.defineProperty(exports, "generateContourSetsFromLabelmap", { enumerable: true, get: function () { return generateContourSetsFromLabelmap_1.generateContourSetsFromLabelmap; } });
+const AnnotationToPointData_1 = __importDefault(require("./AnnotationToPointData"));
+exports.AnnotationToPointData = AnnotationToPointData_1.default;
+//# sourceMappingURL=index.js.map

package/dist/cjs/utilities/contours/index.js.map

@@ -0,0 +1 @@
+{"version":3,"file":"index.js","sourceRoot":"","sources":["../../../../src/utilities/contours/index.ts"],"names":[],"mappings":";;;;;;AAAA,oEAA4C;AAQ1C,wBARK,uBAAa,CAQL;AAPf,gEAAwC;AAQtC,sBARK,qBAAW,CAQL;AAPb,8EAAsD;AAQpD,6BARK,4BAAkB,CAQL;AAPpB,uFAAoF;AAQlF,gHARO,iEAA+B,OAQP;AANjC,oFAA4D;AAO1D,gCAPK,+BAAqB,CAOL"}

package/dist/cjs/utilities/contours/mergePoints.d.ts

@@ -0,0 +1,8 @@
+export declare function removeDuplicatePoints(polyData: any, bypass?: boolean): {
+    points: any[];
+    lines: any[];
+};
+declare const _default: {
+    removeDuplicatePoints: typeof removeDuplicatePoints;
+};
+export default _default;

package/dist/cjs/utilities/contours/mergePoints.js

@@ -0,0 +1,77 @@
+"use strict";
+Object.defineProperty(exports, "__esModule", { value: true });
+exports.removeDuplicatePoints = void 0;
+function ptInArray(array, pt) {
+    let index = -1;
+    for (let i = 0; i < array.length; i++) {
+        if (isSamePoint(pt, array[i])) {
+            index = i;
+        }
+    }
+    return index;
+}
+function isSamePoint(ptA, ptB) {
+    if (ptA[0] == ptB[0] && ptA[1] == ptB[1] && ptA[2] == ptB[2]) {
+        return true;
+    }
+    else {
+        return false;
+    }
+}
+function replacePointIndexReferences(linesArray, oldIndex, newIndex) {
+    for (let i = 0; i < linesArray.length; i++) {
+        const line = linesArray[i];
+        if (line.a == oldIndex) {
+            line.a = newIndex;
+        }
+        else if (line.b == oldIndex) {
+            line.b = newIndex;
+        }
+    }
+}
+function removeDuplicatePoints(polyData, bypass = false) {
+    const points = polyData.getPoints();
+    const lines = polyData.getLines();
+    const pointsArray = [];
+    for (let i = 0; i < points.getNumberOfPoints(); i++) {
+        const pt = points.getPoint(i).slice();
+        pointsArray.push(pt);
+    }
+    const linesArray = [];
+    for (let i = 0; i < lines.getNumberOfCells(); i++) {
+        const cell = lines.getCell(i * 3).slice();
+        const a = cell[0];
+        const b = cell[1];
+        const line = {
+            a,
+            b,
+        };
+        linesArray.push(line);
+    }
+    if (bypass) {
+        return { points: pointsArray, lines: linesArray };
+    }
+    const newPoints = [];
+    for (let i = 0; i < pointsArray.length; i++) {
+        const pt = pointsArray[i];
+        let index = ptInArray(newPoints, pt);
+        if (index >= 0) {
+            replacePointIndexReferences(linesArray, i, index);
+        }
+        else {
+            index = newPoints.length;
+            newPoints.push(pt);
+            replacePointIndexReferences(linesArray, i, index);
+        }
+    }
+    const newLines = [];
+    linesArray.forEach((line) => {
+        if (line.a != line.b) {
+            newLines.push(line);
+        }
+    });
+    return { points: newPoints, lines: newLines };
+}
+exports.removeDuplicatePoints = removeDuplicatePoints;
+exports.default = { removeDuplicatePoints };
+//# sourceMappingURL=mergePoints.js.map

package/dist/cjs/utilities/contours/mergePoints.js.map

@@ -0,0 +1 @@
+{"version":3,"file":"mergePoints.js","sourceRoot":"","sources":["../../../../src/utilities/contours/mergePoints.ts"],"names":[],"mappings":";;;AAMA,SAAS,SAAS,CAAC,KAAK,EAAE,EAAE;IAC1B,IAAI,KAAK,GAAG,CAAC,CAAC,CAAC;IACf,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,KAAK,CAAC,MAAM,EAAE,CAAC,EAAE,EAAE;QACrC,IAAI,WAAW,CAAC,EAAE,EAAE,KAAK,CAAC,CAAC,CAAC,CAAC,EAAE;YAC7B,KAAK,GAAG,CAAC,CAAC;SACX;KACF;IACD,OAAO,KAAK,CAAC;AACf,CAAC;AAQD,SAAS,WAAW,CAAC,GAAG,EAAE,GAAG;IAC3B,IAAI,GAAG,CAAC,CAAC,CAAC,IAAI,GAAG,CAAC,CAAC,CAAC,IAAI,GAAG,CAAC,CAAC,CAAC,IAAI,GAAG,CAAC,CAAC,CAAC,IAAI,GAAG,CAAC,CAAC,CAAC,IAAI,GAAG,CAAC,CAAC,CAAC,EAAE;QAC5D,OAAO,IAAI,CAAC;KACb;SAAM;QACL,OAAO,KAAK,CAAC;KACd;AACH,CAAC;AAQD,SAAS,2BAA2B,CAAC,UAAU,EAAE,QAAQ,EAAE,QAAQ;IACjE,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,UAAU,CAAC,MAAM,EAAE,CAAC,EAAE,EAAE;QAC1C,MAAM,IAAI,GAAG,UAAU,CAAC,CAAC,CAAC,CAAC;QAC3B,IAAI,IAAI,CAAC,CAAC,IAAI,QAAQ,EAAE;YACtB,IAAI,CAAC,CAAC,GAAG,QAAQ,CAAC;SACnB;aAAM,IAAI,IAAI,CAAC,CAAC,IAAI,QAAQ,EAAE;YAC7B,IAAI,CAAC,CAAC,GAAG,QAAQ,CAAC;SACnB;KACF;AACH,CAAC;AASD,SAAgB,qBAAqB,CAAC,QAAQ,EAAE,MAAM,GAAG,KAAK;IAC5D,MAAM,MAAM,GAAG,QAAQ,CAAC,SAAS,EAAE,CAAC;IACpC,MAAM,KAAK,GAAG,QAAQ,CAAC,QAAQ,EAAE,CAAC;IAElC,MAAM,WAAW,GAAG,EAAE,CAAC;IACvB,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,MAAM,CAAC,iBAAiB,EAAE,EAAE,CAAC,EAAE,EAAE;QACnD,MAAM,EAAE,GAAG,MAAM,CAAC,QAAQ,CAAC,CAAC,CAAC,CAAC,KAAK,EAAE,CAAC;QACtC,WAAW,CAAC,IAAI,CAAC,EAAE,CAAC,CAAC;KACtB;IACD,MAAM,UAAU,GAAG,EAAE,CAAC;IACtB,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,KAAK,CAAC,gBAAgB,EAAE,EAAE,CAAC,EAAE,EAAE;QACjD,MAAM,IAAI,GAAG,KAAK,CAAC,OAAO,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,KAAK,EAAE,CAAC;QAE1C,MAAM,CAAC,GAAG,IAAI,CAAC,CAAC,CAAC,CAAC;QAClB,MAAM,CAAC,GAAG,IAAI,CAAC,CAAC,CAAC,CAAC;QAClB,MAAM,IAAI,GAAG;YACX,CAAC;YACD,CAAC;SACF,CAAC;QACF,UAAU,CAAC,IAAI,CAAC,IAAI,CAAC,CAAC;KACvB;IAED,IAAI,MAAM,EAAE;QACV,OAAO,EAAE,MAAM,EAAE,WAAW,EAAE,KAAK,EAAE,UAAU,EAAE,CAAC;KACnD;IAGD,MAAM,SAAS,GAAG,EAAE,CAAC;IACrB,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,WAAW,CAAC,MAAM,EAAE,CAAC,EAAE,EAAE;QAC3C,MAAM,EAAE,GAAG,WAAW,CAAC,CAAC,CAAC,CAAC;QAC1B,IAAI,KAAK,GAAG,SAAS,CAAC,SAAS,EAAE,EAAE,CAAC,CAAC;QAErC,IAAI,KAAK,IAAI,CAAC,EAAE;YAEd,2BAA2B,CAAC,UAAU,EAAE,CAAC,EAAE,KAAK,CAAC,CAAC;SACnD;aAAM;YACL,KAAK,GAAG,SAAS,CAAC,MAAM,CAAC;YACzB,SAAS,CAAC,IAAI,CAAC,EAAE,CAAC,CAAC;YACnB,2BAA2B,CAAC,UAAU,EAAE,CAAC,EAAE,KAAK,CAAC,CAAC;SACnD;KACF;IAGD,MAAM,QAAQ,GAAG,EAAE,CAAC;IACpB,UAAU,CAAC,OAAO,CAAC,CAAC,IAAI,EAAE,EAAE;QAC1B,IAAI,IAAI,CAAC,CAAC,IAAI,IAAI,CAAC,CAAC,EAAE;YACpB,QAAQ,CAAC,IAAI,CAAC,IAAI,CAAC,CAAC;SACrB;IACH,CAAC,CAAC,CAAC;IAEH,OAAO,EAAE,MAAM,EAAE,SAAS,EAAE,KAAK,EAAE,QAAQ,EAAE,CAAC;AAChD,CAAC;AAnDD,sDAmDC;AAED,kBAAe,EAAE,qBAAqB,EAAE,CAAC"}

package/dist/cjs/utilities/index.d.ts

@@ -14,6 +14,7 @@ import scroll from './scroll';
 import roundNumber from './roundNumber';
 import { pointToString } from './pointToString';
 import annotationFrameRange from './annotationFrameRange';
+import * as contours from './contours';
 import * as segmentation from './segmentation';
 import * as drawing from './drawing';
 import * as math from './math';
@@ -31,4 +32,4 @@ import * as dynamicVolume from './dynamicVolume';
 import * as polyDataUtils from './polyData/utils';
 import * as voi from './voi';
 import { triggerEvent } from '@cornerstonejs/core';
-export { math, planar, viewportFilters, drawing, debounce, dynamicVolume, throttle, orientation, isObject, touch, triggerEvent, calibrateImageSpacing, getCalibratedLengthUnits, getCalibratedAreaUnits, getCalibratedScale, segmentation, triggerAnnotationRenderForViewportIds, triggerAnnotationRender, pointInShapeCallback, pointInSurroundingSphereCallback, getAnnotationNearPoint, getAnnotationNearPointOnEnabledElement, jumpToSlice, viewport, cine, clip, boundingBox, rectangleROITool, planarFreehandROITool, stackPrefetch, stackContextPrefetch, scroll, roundNumber, pointToString, polyDataUtils, voi, annotationFrameRange, };
+export { math, planar, viewportFilters, drawing, debounce, dynamicVolume, throttle, orientation, isObject, touch, triggerEvent, calibrateImageSpacing, getCalibratedLengthUnits, getCalibratedAreaUnits, getCalibratedScale, segmentation, contours, triggerAnnotationRenderForViewportIds, triggerAnnotationRender, pointInShapeCallback, pointInSurroundingSphereCallback, getAnnotationNearPoint, getAnnotationNearPointOnEnabledElement, jumpToSlice, viewport, cine, clip, boundingBox, rectangleROITool, planarFreehandROITool, stackPrefetch, stackContextPrefetch, scroll, roundNumber, pointToString, polyDataUtils, voi, annotationFrameRange, };

package/dist/cjs/utilities/index.js

@@ -26,7 +26,7 @@ var __importDefault = (this && this.__importDefault) || function (mod) {
     return (mod && mod.__esModule) ? mod : { "default": mod };
 };
 Object.defineProperty(exports, "__esModule", { value: true });
-exports.annotationFrameRange = exports.voi = exports.polyDataUtils = exports.pointToString = exports.roundNumber = exports.scroll = exports.stackContextPrefetch = exports.stackPrefetch = exports.planarFreehandROITool = exports.rectangleROITool = exports.boundingBox = exports.clip = exports.cine = exports.viewport = exports.jumpToSlice = exports.getAnnotationNearPointOnEnabledElement = exports.getAnnotationNearPoint = exports.pointInSurroundingSphereCallback = exports.pointInShapeCallback = exports.triggerAnnotationRender = exports.triggerAnnotationRenderForViewportIds = exports.segmentation = exports.getCalibratedScale = exports.getCalibratedAreaUnits = exports.getCalibratedLengthUnits = exports.calibrateImageSpacing = exports.triggerEvent = exports.touch = exports.isObject = exports.orientation = exports.throttle = exports.dynamicVolume = exports.debounce = exports.drawing = exports.viewportFilters = exports.planar = exports.math = void 0;
+exports.annotationFrameRange = exports.voi = exports.polyDataUtils = exports.pointToString = exports.roundNumber = exports.scroll = exports.stackContextPrefetch = exports.stackPrefetch = exports.planarFreehandROITool = exports.rectangleROITool = exports.boundingBox = exports.clip = exports.cine = exports.viewport = exports.jumpToSlice = exports.getAnnotationNearPointOnEnabledElement = exports.getAnnotationNearPoint = exports.pointInSurroundingSphereCallback = exports.pointInShapeCallback = exports.triggerAnnotationRender = exports.triggerAnnotationRenderForViewportIds = exports.contours = exports.segmentation = exports.getCalibratedScale = exports.getCalibratedAreaUnits = exports.getCalibratedLengthUnits = exports.calibrateImageSpacing = exports.triggerEvent = exports.touch = exports.isObject = exports.orientation = exports.throttle = exports.dynamicVolume = exports.debounce = exports.drawing = exports.viewportFilters = exports.planar = exports.math = void 0;
 const getAnnotationNearPoint_1 = require("./getAnnotationNearPoint");
 Object.defineProperty(exports, "getAnnotationNearPoint", { enumerable: true, get: function () { return getAnnotationNearPoint_1.getAnnotationNearPoint; } });
 Object.defineProperty(exports, "getAnnotationNearPointOnEnabledElement", { enumerable: true, get: function () { return getAnnotationNearPoint_1.getAnnotationNearPointOnEnabledElement; } });
@@ -62,6 +62,8 @@ const pointToString_1 = require("./pointToString");
 Object.defineProperty(exports, "pointToString", { enumerable: true, get: function () { return pointToString_1.pointToString; } });
 const annotationFrameRange_1 = __importDefault(require("./annotationFrameRange"));
 exports.annotationFrameRange = annotationFrameRange_1.default;
+const contours = __importStar(require("./contours"));
+exports.contours = contours;
 const segmentation = __importStar(require("./segmentation"));
 exports.segmentation = segmentation;
 const drawing = __importStar(require("./drawing"));

package/dist/cjs/utilities/index.js.map

@@ -1 +1 @@
-{"version":3,"file":"index.js","sourceRoot":"","sources":["../../../src/utilities/index.ts"],"names":[],"mappings":";;;;;;;;;;;;;;;;;;;;;;;;;;;;;AAAA,qEAGkC;AAkEhC,uGApEA,+CAAsB,OAoEA;AACtB,uHApEA,+DAAsC,OAoEA;AAhExC,0DAAkC;AA+ChC,mBA/CK,kBAAQ,CA+CL;AA9CV,0DAAkC;AAgDhC,mBAhDK,kBAAQ,CAgDL;AA/CV,0DAAkC;AAiDhC,mBAjDK,kBAAQ,CAiDL;AAhDV,kDAA0B;AAiExB,eAjEK,cAAI,CAiEL;AAhEN,oFAA4D;AAkD1D,gCAlDK,+BAAqB,CAkDL;AAjDvB,6DAI8B;AA8C5B,yGAjDA,6CAAwB,OAiDA;AACxB,uGAjDA,2CAAsB,OAiDA;AACtB,mGAjDA,uCAAkB,OAiDA;AA/CpB,oHAA4F;AAiD1F,gDAjDK,+CAAqC,CAiDL;AAhDvC,wFAAgE;AAiD9D,kCAjDK,iCAAuB,CAiDL;AAhDzB,yEAAiD;AAqD/C,sBArDK,qBAAW,CAqDL;AAnDb,kFAA0D;AA+CxD,+BA/CK,8BAAoB,CA+CL;AA9CtB,0GAAkF;AA+ChF,2CA/CK,0CAAgC,CA+CL;AA9ClC,sDAA8B;AA0D5B,iBA1DK,gBAAM,CA0DL;AAzDR,gEAAwC;AA0DtC,sBA1DK,qBAAW,CA0DL;AAzDb,mDAAgD;AA0D9C,8FA1DO,6BAAa,OA0DP;AAzDf,kFAA0D;AA4DxD,+BA5DK,8BAAoB,CA4DL;AAzDtB,6DAA+C;AAoC7C,oCAAY;AAnCd,mDAAqC;AAuBnC,0BAAO;AAtBT,6CAA+B;AAmB7B,oBAAI;AAlBN,iDAAmC;AAmBjC,wBAAM;AAlBR,mEAAqD;AAmBnD,0CAAe;AAlBjB,2DAA6C;AAuB3C,kCAAW;AAtBb,6CAA+B;AAuC7B,oBAAI;AAtCN,2DAA6C;AAwC3C,kCAAW;AAvCb,+EAAiE;AAyC/D,sDAAqB;AAxCvB,qEAAuD;AAuCrD,4CAAgB;AAtClB,mDAAsE;AAwCpE,8FAxCO,6BAAa,OAwCP;AACb,qGAzCsB,oCAAoB,OAyCtB;AAxCtB,qDAAuC;AAiCrC,4BAAQ;AAhCV,+CAAiC;AAkB/B,sBAAK;AAjBP,+DAAiD;AAa/C,sCAAa;AAZf,gEAAkD;AAyChD,sCAAa;AAxCf,2CAA6B;AAyC3B,kBAAG;AAtCL,8CAAmD;AAajD,6FAbO,mBAAY,OAaP"}
+{"version":3,"file":"index.js","sourceRoot":"","sources":["../../../src/utilities/index.ts"],"names":[],"mappings":";;;;;;;;;;;;;;;;;;;;;;;;;;;;;AAAA,qEAGkC;AAoEhC,uGAtEA,+CAAsB,OAsEA;AACtB,uHAtEA,+DAAsC,OAsEA;AAlExC,0DAAkC;AAgDhC,mBAhDK,kBAAQ,CAgDL;AA/CV,0DAAkC;AAiDhC,mBAjDK,kBAAQ,CAiDL;AAhDV,0DAAkC;AAkDhC,mBAlDK,kBAAQ,CAkDL;AAjDV,kDAA0B;AAmExB,eAnEK,cAAI,CAmEL;AAlEN,oFAA4D;AAmD1D,gCAnDK,+BAAqB,CAmDL;AAlDvB,6DAI8B;AA+C5B,yGAlDA,6CAAwB,OAkDA;AACxB,uGAlDA,2CAAsB,OAkDA;AACtB,mGAlDA,uCAAkB,OAkDA;AAhDpB,oHAA4F;AAmD1F,gDAnDK,+CAAqC,CAmDL;AAlDvC,wFAAgE;AAmD9D,kCAnDK,iCAAuB,CAmDL;AAlDzB,yEAAiD;AAuD/C,sBAvDK,qBAAW,CAuDL;AArDb,kFAA0D;AAiDxD,+BAjDK,8BAAoB,CAiDL;AAhDtB,0GAAkF;AAiDhF,2CAjDK,0CAAgC,CAiDL;AAhDlC,sDAA8B;AA4D5B,iBA5DK,gBAAM,CA4DL;AA3DR,gEAAwC;AA4DtC,sBA5DK,qBAAW,CA4DL;AA3Db,mDAAgD;AA4D9C,8FA5DO,6BAAa,OA4DP;AA3Df,kFAA0D;AA8DxD,+BA9DK,8BAAoB,CA8DL;AA3DtB,qDAAuC;AAsCrC,4BAAQ;AArCV,6DAA+C;AAoC7C,oCAAY;AAnCd,mDAAqC;AAuBnC,0BAAO;AAtBT,6CAA+B;AAmB7B,oBAAI;AAlBN,iDAAmC;AAmBjC,wBAAM;AAlBR,mEAAqD;AAmBnD,0CAAe;AAlBjB,2DAA6C;AAuB3C,kCAAW;AAtBb,6CAA+B;AAwC7B,oBAAI;AAvCN,2DAA6C;AAyC3C,kCAAW;AAxCb,+EAAiE;AA0C/D,sDAAqB;AAzCvB,qEAAuD;AAwCrD,4CAAgB;AAvClB,mDAAsE;AAyCpE,8FAzCO,6BAAa,OAyCP;AACb,qGA1CsB,oCAAoB,OA0CtB;AAzCtB,qDAAuC;AAkCrC,4BAAQ;AAjCV,+CAAiC;AAkB/B,sBAAK;AAjBP,+DAAiD;AAa/C,sCAAa;AAZf,gEAAkD;AA0ChD,sCAAa;AAzCf,2CAA6B;AA0C3B,kBAAG;AAvCL,8CAAmD;AAajD,6FAbO,mBAAY,OAaP"}

package/dist/cjs/utilities/segmentation/contourAndFindLargestBidirectional.d.ts

@@ -0,0 +1 @@
+export default function contourAndFindLargestBidirectional(segmentation: any): any;

package/dist/cjs/utilities/segmentation/contourAndFindLargestBidirectional.js

@@ -0,0 +1,31 @@
+"use strict";
+var __importDefault = (this && this.__importDefault) || function (mod) {
+    return (mod && mod.__esModule) ? mod : { "default": mod };
+};
+Object.defineProperty(exports, "__esModule", { value: true });
+const contours_1 = require("../contours");
+const SegmentationRepresentations_1 = __importDefault(require("../../enums/SegmentationRepresentations"));
+const findLargestBidirectional_1 = __importDefault(require("./findLargestBidirectional"));
+const EPSILON = 1e-2;
+const { Labelmap } = SegmentationRepresentations_1.default;
+function contourAndFindLargestBidirectional(segmentation) {
+    const contours = (0, contours_1.generateContourSetsFromLabelmap)({
+        segmentations: segmentation,
+    });
+    if (!(contours === null || contours === void 0 ? void 0 : contours.length) || !contours[0].sliceContours.length) {
+        return;
+    }
+    const { representationData, segments = [
+        null,
+        { label: 'Unspecified', color: null, containedSegmentIndices: null },
+    ], } = segmentation;
+    const { volumeId: segVolumeId } = representationData[Labelmap];
+    const segmentIndex = segments.findIndex((it) => !!it);
+    if (segmentIndex === -1) {
+        return;
+    }
+    segments[segmentIndex].segmentIndex = segmentIndex;
+    return (0, findLargestBidirectional_1.default)(contours[0], segVolumeId, segments[segmentIndex]);
+}
+exports.default = contourAndFindLargestBidirectional;
+//# sourceMappingURL=contourAndFindLargestBidirectional.js.map

package/dist/cjs/utilities/segmentation/contourAndFindLargestBidirectional.js.map

@@ -0,0 +1 @@
+{"version":3,"file":"contourAndFindLargestBidirectional.js","sourceRoot":"","sources":["../../../../src/utilities/segmentation/contourAndFindLargestBidirectional.ts"],"names":[],"mappings":";;;;;AAAA,0CAA8D;AAC9D,0GAAkF;AAClF,0FAAkE;AAElE,MAAM,OAAO,GAAG,IAAI,CAAC;AACrB,MAAM,EAAE,QAAQ,EAAE,GAAG,qCAA2B,CAAC;AAYjD,SAAwB,kCAAkC,CAAC,YAAY;IACrE,MAAM,QAAQ,GAAG,IAAA,0CAA+B,EAAC;QAC/C,aAAa,EAAE,YAAY;KAC5B,CAAC,CAAC;IAEH,IAAI,CAAC,CAAA,QAAQ,aAAR,QAAQ,uBAAR,QAAQ,CAAE,MAAM,CAAA,IAAI,CAAC,QAAQ,CAAC,CAAC,CAAC,CAAC,aAAa,CAAC,MAAM,EAAE;QAC1D,OAAO;KACR;IAED,MAAM,EACJ,kBAAkB,EAClB,QAAQ,GAAG;QACT,IAAI;QACJ,EAAE,KAAK,EAAE,aAAa,EAAE,KAAK,EAAE,IAAI,EAAE,uBAAuB,EAAE,IAAI,EAAE;KACrE,GACF,GAAG,YAAY,CAAC;IACjB,MAAM,EAAE,QAAQ,EAAE,WAAW,EAAE,GAAG,kBAAkB,CAAC,QAAQ,CAAC,CAAC;IAE/D,MAAM,YAAY,GAAG,QAAQ,CAAC,SAAS,CAAC,CAAC,EAAE,EAAE,EAAE,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC;IACtD,IAAI,YAAY,KAAK,CAAC,CAAC,EAAE;QACvB,OAAO;KACR;IACD,QAAQ,CAAC,YAAY,CAAC,CAAC,YAAY,GAAG,YAAY,CAAC;IACnD,OAAO,IAAA,kCAAwB,EAC7B,QAAQ,CAAC,CAAC,CAAC,EACX,WAAW,EACX,QAAQ,CAAC,YAAY,CAAC,CACvB,CAAC;AACJ,CAAC;AA5BD,qDA4BC"}

package/dist/cjs/utilities/segmentation/createBidirectionalToolData.d.ts

@@ -0,0 +1,14 @@
+import type { Types } from '@cornerstonejs/core';
+import type { Annotation } from '../../types/AnnotationTypes';
+export declare type BidirectionalData = {
+    majorAxis: [Types.Point3, Types.Point3];
+    minorAxis: [Types.Point3, Types.Point3];
+    maxMajor: number;
+    maxMinor: number;
+    segmentIndex: number;
+    label?: string;
+    color?: string | number[];
+    referencedImageId: string;
+    FrameOfReferenceUID: string;
+};
+export default function createBidirectionalToolData(bidirectionalData: BidirectionalData, viewport: any): Annotation;

package/dist/cjs/utilities/segmentation/createBidirectionalToolData.js

@@ -0,0 +1,43 @@
+"use strict";
+Object.defineProperty(exports, "__esModule", { value: true });
+function createBidirectionalToolData(bidirectionalData, viewport) {
+    const { majorAxis, minorAxis, label = '', FrameOfReferenceUID, referencedImageId, } = bidirectionalData;
+    const [major0, major1] = majorAxis;
+    const [minor0, minor1] = minorAxis;
+    const { viewUp, viewPlaneNormal } = viewport.getCamera();
+    const points = [major0, major1, minor0, minor1];
+    const bidirectionalToolData = {
+        highlighted: true,
+        invalidated: true,
+        metadata: {
+            toolName: 'Bidirectional',
+            viewPlaneNormal,
+            viewUp,
+            FrameOfReferenceUID,
+            referencedImageId,
+        },
+        data: {
+            handles: {
+                points,
+                textBox: {
+                    hasMoved: false,
+                    worldPosition: [0, 0, 0],
+                    worldBoundingBox: {
+                        topLeft: [0, 0, 0],
+                        topRight: [0, 0, 0],
+                        bottomLeft: [0, 0, 0],
+                        bottomRight: [0, 0, 0],
+                    },
+                },
+                activeHandleIndex: null,
+            },
+            label,
+            cachedStats: {},
+        },
+        isLocked: false,
+        isVisible: true,
+    };
+    return bidirectionalToolData;
+}
+exports.default = createBidirectionalToolData;
+//# sourceMappingURL=createBidirectionalToolData.js.map

package/dist/cjs/utilities/segmentation/createBidirectionalToolData.js.map

@@ -0,0 +1 @@
+{"version":3,"file":"createBidirectionalToolData.js","sourceRoot":"","sources":["../../../../src/utilities/segmentation/createBidirectionalToolData.ts"],"names":[],"mappings":";;AAmBA,SAAwB,2BAA2B,CACjD,iBAAoC,EACpC,QAAQ;IAER,MAAM,EACJ,SAAS,EACT,SAAS,EACT,KAAK,GAAG,EAAE,EACV,mBAAmB,EACnB,iBAAiB,GAClB,GAAG,iBAAiB,CAAC;IACtB,MAAM,CAAC,MAAM,EAAE,MAAM,CAAC,GAAG,SAAS,CAAC;IACnC,MAAM,CAAC,MAAM,EAAE,MAAM,CAAC,GAAG,SAAS,CAAC;IAEnC,MAAM,EAAE,MAAM,EAAE,eAAe,EAAE,GAAG,QAAQ,CAAC,SAAS,EAAE,CAAC;IACzD,MAAM,MAAM,GAAG,CAAC,MAAM,EAAE,MAAM,EAAE,MAAM,EAAE,MAAM,CAAC,CAAC;IAChD,MAAM,qBAAqB,GAAG;QAC5B,WAAW,EAAE,IAAI;QACjB,WAAW,EAAE,IAAI;QACjB,QAAQ,EAAE;YACR,QAAQ,EAAE,eAAe;YACzB,eAAe;YACf,MAAM;YACN,mBAAmB;YACnB,iBAAiB;SAClB;QACD,IAAI,EAAE;YACJ,OAAO,EAAE;gBACP,MAAM;gBACN,OAAO,EAAE;oBACP,QAAQ,EAAE,KAAK;oBACf,aAAa,EAAE,CAAC,CAAC,EAAE,CAAC,EAAE,CAAC,CAAiB;oBACxC,gBAAgB,EAAE;wBAChB,OAAO,EAAE,CAAC,CAAC,EAAE,CAAC,EAAE,CAAC,CAAiB;wBAClC,QAAQ,EAAE,CAAC,CAAC,EAAE,CAAC,EAAE,CAAC,CAAiB;wBACnC,UAAU,EAAE,CAAC,CAAC,EAAE,CAAC,EAAE,CAAC,CAAiB;wBACrC,WAAW,EAAE,CAAC,CAAC,EAAE,CAAC,EAAE,CAAC,CAAiB;qBACvC;iBACF;gBACD,iBAAiB,EAAE,IAAI;aACxB;YACD,KAAK;YACL,WAAW,EAAE,EAAE;SAChB;QACD,QAAQ,EAAE,KAAK;QACf,SAAS,EAAE,IAAI;KAChB,CAAC;IACF,OAAO,qBAAqB,CAAC;AAC/B,CAAC;AAhDD,8CAgDC"}

package/dist/cjs/utilities/segmentation/findLargestBidirectional.d.ts

@@ -0,0 +1 @@
+export default function findLargestBidirectional(contours: any, segVolumeId: string, segment: any): any;

package/dist/cjs/utilities/segmentation/findLargestBidirectional.js

@@ -0,0 +1,94 @@
+"use strict";
+Object.defineProperty(exports, "__esModule", { value: true });
+const gl_matrix_1 = require("gl-matrix");
+const isLineInSegment_1 = require("./isLineInSegment");
+const EPSILON = 1e-2;
+function findLargestBidirectional(contours, segVolumeId, segment) {
+    const { sliceContours } = contours;
+    const { segmentIndex, containedSegmentIndices } = segment;
+    let maxBidirectional;
+    const isInSegment = (0, isLineInSegment_1.createIsInSegment)(segVolumeId, segmentIndex, containedSegmentIndices);
+    for (const sliceContour of sliceContours) {
+        const bidirectional = createBidirectionalForSlice(sliceContour, isInSegment, maxBidirectional);
+        if (!bidirectional) {
+            continue;
+        }
+        maxBidirectional = bidirectional;
+    }
+    if (maxBidirectional) {
+        Object.assign(maxBidirectional, segment);
+    }
+    return maxBidirectional;
+}
+exports.default = findLargestBidirectional;
+function createBidirectionalForSlice(sliceContour, isInSegment, currentMax = { maxMajor: 0, maxMinor: 0 }) {
+    const { points } = sliceContour.polyData;
+    const { maxMinor: currentMaxMinor, maxMajor: currentMaxMajor } = currentMax;
+    let maxMajor = currentMaxMajor * currentMaxMajor;
+    let maxMinor = currentMaxMinor * currentMaxMinor;
+    let maxMajorPoints;
+    for (let index1 = 0; index1 < points.length; index1++) {
+        for (let index2 = index1 + 1; index2 < points.length; index2++) {
+            const point1 = points[index1];
+            const point2 = points[index2];
+            const distance2 = gl_matrix_1.vec3.sqrDist(point1, point2);
+            if (distance2 < maxMajor) {
+                continue;
+            }
+            if (distance2 - EPSILON < maxMajor + EPSILON && maxMajorPoints) {
+                continue;
+            }
+            if (!isInSegment.testCenter(point1, point2)) {
+                continue;
+            }
+            if (!(0, isLineInSegment_1.isLineInSegment)(point1, point2, isInSegment)) {
+                continue;
+            }
+            maxMajor = distance2 - EPSILON;
+            maxMajorPoints = [index1, index2];
+            maxMinor = 0;
+        }
+    }
+    if (!maxMajorPoints) {
+        return;
+    }
+    maxMajor = Math.sqrt(maxMajor + EPSILON);
+    const handle0 = points[maxMajorPoints[0]];
+    const handle1 = points[maxMajorPoints[1]];
+    const unitMajor = gl_matrix_1.vec3.sub(gl_matrix_1.vec3.create(), handle0, handle1);
+    gl_matrix_1.vec3.scale(unitMajor, unitMajor, 1 / maxMajor);
+    let maxMinorPoints;
+    for (let index1 = 0; index1 < points.length; index1++) {
+        for (let index2 = index1 + 1; index2 < points.length; index2++) {
+            const point1 = points[index1];
+            const point2 = points[index2];
+            const distance2 = gl_matrix_1.vec3.sqrDist(point1, point2);
+            if (distance2 <= maxMinor) {
+                continue;
+            }
+            const delta = gl_matrix_1.vec3.sub(gl_matrix_1.vec3.create(), point1, point2);
+            const dot = Math.abs(gl_matrix_1.vec3.dot(delta, unitMajor)) / Math.sqrt(distance2);
+            if (dot > EPSILON) {
+                continue;
+            }
+            if (!isInSegment.testCenter(point1, point2)) {
+                continue;
+            }
+            if (!(0, isLineInSegment_1.isLineInSegment)(point1, point2, isInSegment)) {
+                continue;
+            }
+            maxMinor = distance2;
+            maxMinorPoints = [index1, index2];
+        }
+    }
+    if (!maxMinorPoints) {
+        return;
+    }
+    maxMinor = Math.sqrt(maxMinor);
+    const handle2 = points[maxMinorPoints[0]];
+    const handle3 = points[maxMinorPoints[1]];
+    const bidirectional = Object.assign({ majorAxis: [handle0, handle1], minorAxis: [handle2, handle3], maxMajor,
+        maxMinor }, sliceContour);
+    return bidirectional;
+}
+//# sourceMappingURL=findLargestBidirectional.js.map

package/dist/cjs/utilities/segmentation/findLargestBidirectional.js.map

@@ -0,0 +1 @@
+{"version":3,"file":"findLargestBidirectional.js","sourceRoot":"","sources":["../../../../src/utilities/segmentation/findLargestBidirectional.ts"],"names":[],"mappings":";;AAAA,yCAAiC;AAEjC,uDAAuE;AAGvE,MAAM,OAAO,GAAG,IAAI,CAAC;AAqBrB,SAAwB,wBAAwB,CAC9C,QAAQ,EACR,WAAmB,EACnB,OAAO;IAEP,MAAM,EAAE,aAAa,EAAE,GAAG,QAAQ,CAAC;IACnC,MAAM,EAAE,YAAY,EAAE,uBAAuB,EAAE,GAAG,OAAO,CAAC;IAC1D,IAAI,gBAAgB,CAAC;IACrB,MAAM,WAAW,GAAG,IAAA,mCAAiB,EACnC,WAAW,EACX,YAAY,EACZ,uBAAuB,CACxB,CAAC;IACF,KAAK,MAAM,YAAY,IAAI,aAAa,EAAE;QACxC,MAAM,aAAa,GAAG,2BAA2B,CAC/C,YAAY,EACZ,WAAW,EACX,gBAAgB,CACjB,CAAC;QACF,IAAI,CAAC,aAAa,EAAE;YAClB,SAAS;SACV;QACD,gBAAgB,GAAG,aAAa,CAAC;KAClC;IACD,IAAI,gBAAgB,EAAE;QACpB,MAAM,CAAC,MAAM,CAAC,gBAAgB,EAAE,OAAO,CAAC,CAAC;KAC1C;IACD,OAAO,gBAAgB,CAAC;AAC1B,CAAC;AA5BD,2CA4BC;AAWD,SAAS,2BAA2B,CAClC,YAAY,EACZ,WAAW,EACX,UAAU,GAAG,EAAE,QAAQ,EAAE,CAAC,EAAE,QAAQ,EAAE,CAAC,EAAE;IAEzC,MAAM,EAAE,MAAM,EAAE,GAAG,YAAY,CAAC,QAAQ,CAAC;IACzC,MAAM,EAAE,QAAQ,EAAE,eAAe,EAAE,QAAQ,EAAE,eAAe,EAAE,GAAG,UAAU,CAAC;IAC5E,IAAI,QAAQ,GAAG,eAAe,GAAG,eAAe,CAAC;IACjD,IAAI,QAAQ,GAAG,eAAe,GAAG,eAAe,CAAC;IACjD,IAAI,cAAc,CAAC;IACnB,KAAK,IAAI,MAAM,GAAG,CAAC,EAAE,MAAM,GAAG,MAAM,CAAC,MAAM,EAAE,MAAM,EAAE,EAAE;QACrD,KAAK,IAAI,MAAM,GAAG,MAAM,GAAG,CAAC,EAAE,MAAM,GAAG,MAAM,CAAC,MAAM,EAAE,MAAM,EAAE,EAAE;YAC9D,MAAM,MAAM,GAAG,MAAM,CAAC,MAAM,CAAC,CAAC;YAC9B,MAAM,MAAM,GAAG,MAAM,CAAC,MAAM,CAAC,CAAC;YAC9B,MAAM,SAAS,GAAG,gBAAI,CAAC,OAAO,CAAC,MAAM,EAAE,MAAM,CAAC,CAAC;YAC/C,IAAI,SAAS,GAAG,QAAQ,EAAE;gBACxB,SAAS;aACV;YACD,IAAI,SAAS,GAAG,OAAO,GAAG,QAAQ,GAAG,OAAO,IAAI,cAAc,EAAE;gBAG9D,SAAS;aACV;YACD,IAAI,CAAC,WAAW,CAAC,UAAU,CAAC,MAAM,EAAE,MAAM,CAAC,EAAE;gBAG3C,SAAS;aACV;YACD,IAAI,CAAC,IAAA,iCAAe,EAAC,MAAM,EAAE,MAAM,EAAE,WAAW,CAAC,EAAE;gBAEjD,SAAS;aACV;YACD,QAAQ,GAAG,SAAS,GAAG,OAAO,CAAC;YAC/B,cAAc,GAAG,CAAC,MAAM,EAAE,MAAM,CAAC,CAAC;YAClC,QAAQ,GAAG,CAAC,CAAC;SACd;KACF;IACD,IAAI,CAAC,cAAc,EAAE;QACnB,OAAO;KACR;IAED,QAAQ,GAAG,IAAI,CAAC,IAAI,CAAC,QAAQ,GAAG,OAAO,CAAC,CAAC;IACzC,MAAM,OAAO,GAAG,MAAM,CAAC,cAAc,CAAC,CAAC,CAAC,CAAC,CAAC;IAC1C,MAAM,OAAO,GAAG,MAAM,CAAC,cAAc,CAAC,CAAC,CAAC,CAAC,CAAC;IAC1C,MAAM,SAAS,GAAG,gBAAI,CAAC,GAAG,CAAC,gBAAI,CAAC,MAAM,EAAE,EAAE,OAAO,EAAE,OAAO,CAAC,CAAC;IAC5D,gBAAI,CAAC,KAAK,CAAC,SAAS,EAAE,SAAS,EAAE,CAAC,GAAG,QAAQ,CAAC,CAAC;IAE/C,IAAI,cAAc,CAAC;IAEnB,KAAK,IAAI,MAAM,GAAG,CAAC,EAAE,MAAM,GAAG,MAAM,CAAC,MAAM,EAAE,MAAM,EAAE,EAAE;QACrD,KAAK,IAAI,MAAM,GAAG,MAAM,GAAG,CAAC,EAAE,MAAM,GAAG,MAAM,CAAC,MAAM,EAAE,MAAM,EAAE,EAAE;YAC9D,MAAM,MAAM,GAAG,MAAM,CAAC,MAAM,CAAC,CAAC;YAC9B,MAAM,MAAM,GAAG,MAAM,CAAC,MAAM,CAAC,CAAC;YAC9B,MAAM,SAAS,GAAG,gBAAI,CAAC,OAAO,CAAC,MAAM,EAAE,MAAM,CAAC,CAAC;YAC/C,IAAI,SAAS,IAAI,QAAQ,EAAE;gBACzB,SAAS;aACV;YACD,MAAM,KAAK,GAAG,gBAAI,CAAC,GAAG,CAAC,gBAAI,CAAC,MAAM,EAAE,EAAE,MAAM,EAAE,MAAM,CAAC,CAAC;YAEtD,MAAM,GAAG,GAAG,IAAI,CAAC,GAAG,CAAC,gBAAI,CAAC,GAAG,CAAC,KAAK,EAAE,SAAS,CAAC,CAAC,GAAG,IAAI,CAAC,IAAI,CAAC,SAAS,CAAC,CAAC;YACxE,IAAI,GAAG,GAAG,OAAO,EAAE;gBACjB,SAAS;aACV;YAED,IAAI,CAAC,WAAW,CAAC,UAAU,CAAC,MAAM,EAAE,MAAM,CAAC,EAAE;gBAG3C,SAAS;aACV;YACD,IAAI,CAAC,IAAA,iCAAe,EAAC,MAAM,EAAE,MAAM,EAAE,WAAW,CAAC,EAAE;gBACjD,SAAS;aACV;YACD,QAAQ,GAAG,SAAS,CAAC;YACrB,cAAc,GAAG,CAAC,MAAM,EAAE,MAAM,CAAC,CAAC;SACnC;KACF;IAED,IAAI,CAAC,cAAc,EAAE;QAEnB,OAAO;KACR;IACD,QAAQ,GAAG,IAAI,CAAC,IAAI,CAAC,QAAQ,CAAC,CAAC;IAC/B,MAAM,OAAO,GAAG,MAAM,CAAC,cAAc,CAAC,CAAC,CAAC,CAAC,CAAC;IAC1C,MAAM,OAAO,GAAG,MAAM,CAAC,cAAc,CAAC,CAAC,CAAC,CAAC,CAAC;IAE1C,MAAM,aAAa,GAAG,gBACpB,SAAS,EAAE,CAAC,OAAO,EAAE,OAAO,CAAC,EAC7B,SAAS,EAAE,CAAC,OAAO,EAAE,OAAO,CAAC,EAC7B,QAAQ;QACR,QAAQ,IACL,YAAY,CACK,CAAC;IACvB,OAAO,aAAa,CAAC;AACvB,CAAC"}

package/dist/cjs/utilities/segmentation/index.d.ts

@@ -9,4 +9,7 @@ import floodFill from './floodFill';
 import { getBrushSizeForToolGroup, setBrushSizeForToolGroup } from './brushSizeForToolGroup';
 import { getBrushThresholdForToolGroup, setBrushThresholdForToolGroup } from './brushThresholdForToolGroup';
 import thresholdSegmentationByRange from './thresholdSegmentationByRange';
-export { thresholdVolumeByRange, createMergedLabelmapForIndex, isValidRepresentationConfig, getDefaultRepresentationConfig, createLabelmapVolumeForViewport, rectangleROIThresholdVolumeByRange, triggerSegmentationRender, floodFill, getBrushSizeForToolGroup, setBrushSizeForToolGroup, getBrushThresholdForToolGroup, setBrushThresholdForToolGroup, thresholdSegmentationByRange, };
+import contourAndFindLargestBidirectional from './contourAndFindLargestBidirectional';
+import createBidirectionalToolData from './createBidirectionalToolData';
+import segmentContourAction from './segmentContourAction';
+export { thresholdVolumeByRange, createMergedLabelmapForIndex, isValidRepresentationConfig, getDefaultRepresentationConfig, createLabelmapVolumeForViewport, rectangleROIThresholdVolumeByRange, triggerSegmentationRender, floodFill, getBrushSizeForToolGroup, setBrushSizeForToolGroup, getBrushThresholdForToolGroup, setBrushThresholdForToolGroup, thresholdSegmentationByRange, contourAndFindLargestBidirectional, createBidirectionalToolData, segmentContourAction, };

package/dist/cjs/utilities/segmentation/index.js

@@ -3,7 +3,7 @@ var __importDefault = (this && this.__importDefault) || function (mod) {
     return (mod && mod.__esModule) ? mod : { "default": mod };
 };
 Object.defineProperty(exports, "__esModule", { value: true });
-exports.thresholdSegmentationByRange = exports.setBrushThresholdForToolGroup = exports.getBrushThresholdForToolGroup = exports.setBrushSizeForToolGroup = exports.getBrushSizeForToolGroup = exports.floodFill = exports.triggerSegmentationRender = exports.rectangleROIThresholdVolumeByRange = exports.createLabelmapVolumeForViewport = exports.getDefaultRepresentationConfig = exports.isValidRepresentationConfig = exports.createMergedLabelmapForIndex = exports.thresholdVolumeByRange = void 0;
+exports.segmentContourAction = exports.createBidirectionalToolData = exports.contourAndFindLargestBidirectional = exports.thresholdSegmentationByRange = exports.setBrushThresholdForToolGroup = exports.getBrushThresholdForToolGroup = exports.setBrushSizeForToolGroup = exports.getBrushSizeForToolGroup = exports.floodFill = exports.triggerSegmentationRender = exports.rectangleROIThresholdVolumeByRange = exports.createLabelmapVolumeForViewport = exports.getDefaultRepresentationConfig = exports.isValidRepresentationConfig = exports.createMergedLabelmapForIndex = exports.thresholdVolumeByRange = void 0;
 const thresholdVolumeByRange_1 = __importDefault(require("./thresholdVolumeByRange"));
 exports.thresholdVolumeByRange = thresholdVolumeByRange_1.default;
 const rectangleROIThresholdVolumeByRange_1 = __importDefault(require("./rectangleROIThresholdVolumeByRange"));
@@ -28,4 +28,10 @@ Object.defineProperty(exports, "getBrushThresholdForToolGroup", { enumerable: tr
 Object.defineProperty(exports, "setBrushThresholdForToolGroup", { enumerable: true, get: function () { return brushThresholdForToolGroup_1.setBrushThresholdForToolGroup; } });
 const thresholdSegmentationByRange_1 = __importDefault(require("./thresholdSegmentationByRange"));
 exports.thresholdSegmentationByRange = thresholdSegmentationByRange_1.default;
+const contourAndFindLargestBidirectional_1 = __importDefault(require("./contourAndFindLargestBidirectional"));
+exports.contourAndFindLargestBidirectional = contourAndFindLargestBidirectional_1.default;
+const createBidirectionalToolData_1 = __importDefault(require("./createBidirectionalToolData"));
+exports.createBidirectionalToolData = createBidirectionalToolData_1.default;
+const segmentContourAction_1 = __importDefault(require("./segmentContourAction"));
+exports.segmentContourAction = segmentContourAction_1.default;
 //# sourceMappingURL=index.js.map

package/dist/cjs/utilities/segmentation/index.js.map

@@ -1 +1 @@
-{"version":3,"file":"index.js","sourceRoot":"","sources":["../../../../src/utilities/segmentation/index.ts"],"names":[],"mappings":";;;;;;AAAA,sFAA8D;AAmB5D,iCAnBK,gCAAsB,CAmBL;AAlBxB,8GAAsF;AAuBpF,6CAvBK,4CAAkC,CAuBL;AAtBpC,kGAA0E;AAkBxE,uCAlBK,sCAA4B,CAkBL;AAjB9B,gGAAwE;AAkBtE,sCAlBK,qCAA2B,CAkBL;AAjB7B,sGAA8E;AAkB5E,yCAlBK,wCAA8B,CAkBL;AAjBhC,wGAAgF;AAkB9E,0CAlBK,yCAA+B,CAkBL;AAjBjC,2EAAwE;AAmBtE,0GAnBO,qDAAyB,OAmBP;AAlB3B,4DAAoC;AAmBlC,oBAnBK,mBAAS,CAmBL;AAlBX,mEAGiC;AAgB/B,yGAlBA,gDAAwB,OAkBA;AACxB,yGAlBA,gDAAwB,OAkBA;AAhB1B,6EAGsC;AAcpC,8GAhBA,0DAA6B,OAgBA;AAC7B,8GAhBA,0DAA6B,OAgBA;AAd/B,kGAA0E;AAexE,uCAfK,sCAA4B,CAeL"}
+{"version":3,"file":"index.js","sourceRoot":"","sources":["../../../../src/utilities/segmentation/index.ts"],"names":[],"mappings":";;;;;;AAAA,sFAA8D;AAsB5D,iCAtBK,gCAAsB,CAsBL;AArBxB,8GAAsF;AA0BpF,6CA1BK,4CAAkC,CA0BL;AAzBpC,kGAA0E;AAqBxE,uCArBK,sCAA4B,CAqBL;AApB9B,gGAAwE;AAqBtE,sCArBK,qCAA2B,CAqBL;AApB7B,sGAA8E;AAqB5E,yCArBK,wCAA8B,CAqBL;AApBhC,wGAAgF;AAqB9E,0CArBK,yCAA+B,CAqBL;AApBjC,2EAAwE;AAsBtE,0GAtBO,qDAAyB,OAsBP;AArB3B,4DAAoC;AAsBlC,oBAtBK,mBAAS,CAsBL;AArBX,mEAGiC;AAmB/B,yGArBA,gDAAwB,OAqBA;AACxB,yGArBA,gDAAwB,OAqBA;AAnB1B,6EAGsC;AAiBpC,8GAnBA,0DAA6B,OAmBA;AAC7B,8GAnBA,0DAA6B,OAmBA;AAjB/B,kGAA0E;AAkBxE,uCAlBK,sCAA4B,CAkBL;AAjB9B,8GAAsF;AAkBpF,6CAlBK,4CAAkC,CAkBL;AAjBpC,gGAAwE;AAkBtE,sCAlBK,qCAA2B,CAkBL;AAjB7B,kFAA0D;AAkBxD,+BAlBK,8BAAoB,CAkBL"}

package/dist/cjs/utilities/segmentation/isLineInSegment.d.ts

@@ -0,0 +1,9 @@
+import type { Types } from '@cornerstonejs/core';
+import { vec3 } from 'gl-matrix';
+export default function isLineInSegment(point1: Types.Point3, point2: Types.Point3, isInSegment: any): boolean;
+declare function createIsInSegment(segVolumeId: string, segmentIndex: number, containedSegmentIndices?: Set<number>): {
+    testCenter: (point1: any, point2: any) => boolean;
+    toIJK: (point: any) => vec3;
+    testIJK: (ijk: any) => boolean;
+};
+export { createIsInSegment, isLineInSegment };

package/dist/cjs/utilities/segmentation/isLineInSegment.js

@@ -0,0 +1,55 @@
+"use strict";
+Object.defineProperty(exports, "__esModule", { value: true });
+exports.isLineInSegment = exports.createIsInSegment = void 0;
+const core_1 = require("@cornerstonejs/core");
+const gl_matrix_1 = require("gl-matrix");
+function isLineInSegment(point1, point2, isInSegment) {
+    const ijk1 = isInSegment.toIJK(point1);
+    const ijk2 = isInSegment.toIJK(point2);
+    const testPoint = gl_matrix_1.vec3.create();
+    const { testIJK } = isInSegment;
+    const delta = gl_matrix_1.vec3.sub(gl_matrix_1.vec3.create(), ijk1, ijk2);
+    const testSize = Math.round(Math.max(...delta.map(Math.abs)));
+    if (testSize < 2) {
+        return true;
+    }
+    const unitDelta = gl_matrix_1.vec3.scale(gl_matrix_1.vec3.create(), delta, 1 / testSize);
+    for (let i = 1; i < testSize; i++) {
+        gl_matrix_1.vec3.scaleAndAdd(testPoint, ijk2, unitDelta, i);
+        if (!testIJK(testPoint)) {
+            return false;
+        }
+    }
+    return true;
+}
+exports.default = isLineInSegment;
+exports.isLineInSegment = isLineInSegment;
+function createIsInSegment(segVolumeId, segmentIndex, containedSegmentIndices) {
+    const vol = core_1.cache.getVolume(segVolumeId);
+    if (!vol) {
+        console.warn(`No volume found for ${segVolumeId}`);
+        return;
+    }
+    const segData = vol.imageData.getPointData().getScalars().getData();
+    const width = vol.dimensions[0];
+    const pixelsPerSlice = width * vol.dimensions[1];
+    return {
+        testCenter: (point1, point2) => {
+            const point = gl_matrix_1.vec3.add(gl_matrix_1.vec3.create(), point1, point2).map((it) => it / 2);
+            const ijk = vol.imageData.worldToIndex(point).map(Math.round);
+            const [i, j, k] = ijk;
+            const index = i + j * width + k * pixelsPerSlice;
+            const value = segData[index];
+            return value === segmentIndex || (containedSegmentIndices === null || containedSegmentIndices === void 0 ? void 0 : containedSegmentIndices.has(value));
+        },
+        toIJK: (point) => vol.imageData.worldToIndex(point),
+        testIJK: (ijk) => {
+            const [i, j, k] = ijk;
+            const index = Math.round(i) + Math.round(j) * width + Math.round(k) * pixelsPerSlice;
+            const value = segData[index];
+            return value === segmentIndex || (containedSegmentIndices === null || containedSegmentIndices === void 0 ? void 0 : containedSegmentIndices.has(value));
+        },
+    };
+}
+exports.createIsInSegment = createIsInSegment;
+//# sourceMappingURL=isLineInSegment.js.map