feat(seg): Improve Labelmap Statistics, Interpolation, and Threshold Configuration

common/reviews/api/core.api.md

@@ -275,6 +275,9 @@ type BoundsLPS = [Point3, Point3, Point3];
 // @public (undocumented)
 export const cache: Cache_2;
 
+// @public (undocumented)
+function calculateSpacingBetweenImageIds(imageIds: string[]): number;
+
 // @public (undocumented)
 function calculateViewportsSpatialRegistration(viewport1: StackViewport | IVolumeViewport, viewport2: StackViewport | IVolumeViewport): void;
 
@@ -1191,6 +1194,9 @@ interface FlipDirection {
     flipVertical?: boolean;
 }
 
+// @public (undocumented)
+function fnv1aHash(str: string): string;
+
 // @public (undocumented)
 class FrameRange {
     // (undocumented)
@@ -3963,6 +3969,7 @@ declare namespace utilities {
         scaleRGBTransferFunction as scaleRgbTransferFunction,
         triggerEvent,
         imageIdToURI,
+        fnv1aHash,
         metadataProvider as calibratedPixelSpacingMetadataProvider,
         clamp,
         uuidv4,
@@ -4047,7 +4054,8 @@ declare namespace utilities {
         clip,
         transformWorldToIndexContinuous,
         createSubVolume,
-        getVolumeDirectionVectors
+        getVolumeDirectionVectors,
+        calculateSpacingBetweenImageIds
     }
 }
 export { utilities }

common/reviews/api/labelmap-interpolation.api.md

@@ -0,0 +1,18 @@
+## API Report File for "@cornerstonejs/labelmap-interpolation"
+
+> Do not edit this file. It is a report generated by [API Extractor](https://api-extractor.com/).
+
+```ts
+
+// @public (undocumented)
+export function interpolate({ segmentationId, segmentIndex, configuration, }: {
+    segmentationId: string;
+    segmentIndex: number;
+    configuration?: MorphologicalContourInterpolationOptions & {
+        preview: boolean;
+    };
+}): Promise<void>;
+
+// (No @packageDocumentation comment for this package)
+
+```

common/reviews/api/tools.api.md

@@ -757,9 +757,10 @@ class BasicStatsCalculator_2 extends Calculator {
         unit: string;
     }) => NamedStatistics;
     // (undocumented)
-    static statsCallback: ({ value: newValue, pointLPS }: {
+    static statsCallback: ({ value: newValue, pointLPS, pointIJK, }: {
         value: any;
         pointLPS?: any;
+        pointIJK?: any;
     }) => void;
     // (undocumented)
     static statsInit(options: {
@@ -1004,6 +1005,10 @@ enum ChangeTypes {
 
 // @public (undocumented)
 enum ChangeTypes_2 {
+    // (undocumented)
+    COMPUTE_STATISTICS = "Computing Statistics",
+    // (undocumented)
+    INTERPOLATE_LABELMAP = "Interpolating Labelmap",
     // (undocumented)
     POLYSEG_CONTOUR_TO_LABELMAP = "Converting Contour to Labelmap",
     // (undocumented)
@@ -2819,12 +2824,18 @@ function getNormal2(polyline: Types_2.Point2[]): Types_2.Point3;
 // @public (undocumented)
 function getNormal3(polyline: Types_2.Point3[]): Types_2.Point3;
 
+// @public (undocumented)
+function getOrCreateSegmentationVolume(segmentationId: any): any;
+
 // @public (undocumented)
 function getOrientationStringLPS(vector: Types_2.Point3): string;
 
 // @public (undocumented)
 function getPixelValueUnits(modality: string, imageId: string, options: pixelUnitsOptions): string;
 
+// @public (undocumented)
+function getPixelValueUnitsImageId(imageId: string, options: pixelUnitsOptions): string;
+
 // @public (undocumented)
 function getPoint(points: any, idx: any): Types_2.Point3;
 
@@ -2903,6 +2914,13 @@ function getStackSegmentationImageIdsForViewport(viewportId: string, segmentatio
 // @public (undocumented)
 function getState(annotation?: Annotation): AnnotationStyleStates;
 
+// @public (undocumented)
+function getStatistics({ segmentationId, segmentIndices, viewportId, }: {
+    segmentationId: string;
+    segmentIndices: number[] | number;
+    viewportId: string;
+}): Promise<any>;
+
 // @public (undocumented)
 function getStyle<T extends SegmentationRepresentations>(specifier: SpecifierWithType<T>): StyleForType<T>;
 
@@ -3512,12 +3530,11 @@ export class LabelmapBaseTool extends BaseTool {
         segmentColor: Types_2.Color;
     };
     // (undocumented)
-    protected getEditData({ viewport, representationData, segmentsLocked, segmentationId, volumeOperation, }: {
+    protected getEditData({ viewport, representationData, segmentsLocked, segmentationId, }: {
         viewport: any;
         representationData: any;
         segmentsLocked: any;
         segmentationId: any;
-        volumeOperation?: boolean;
     }): EditDataReturnType;
     // (undocumented)
     protected getOperationData(element?: any): ModifiedLabelmapToolOperationData;
@@ -3574,6 +3591,7 @@ type LabelmapToolOperationData = {
     viewPlaneNormal: number[];
     viewUp: number[];
     strategySpecificConfiguration: any;
+    activeStrategy: string;
     points: Types_2.Point3[];
     voxelManager: any;
     override: {
@@ -5281,7 +5299,9 @@ declare namespace segmentation_2 {
         getBrushToolInstances,
         growCut,
         LabelmapMemo,
-        IslandRemoval
+        IslandRemoval,
+        getOrCreateSegmentationVolume,
+        getStatistics
     }
 }
 
@@ -5889,6 +5909,8 @@ type Statistics = {
     label?: string;
     value: number | number[];
     unit: null | string;
+    pointIJK?: Types_2.Point3;
+    pointLPS?: Types_2.Point3;
 };
 
 // @public (undocumented)
@@ -5931,6 +5953,10 @@ enum StrategyCallbacks {
     // (undocumented)
     CreateIsInThreshold = "createIsInThreshold",
     // (undocumented)
+    EnsureImageVolumeFor3DManipulation = "ensureImageVolumeFor3DManipulation",
+    // (undocumented)
+    EnsureSegmentationVolumeFor3DManipulation = "ensureSegmentationVolumeFor3DManipulation",
+    // (undocumented)
     Fill = "fill",
     // (undocumented)
     GetStatistics = "getStatistics",
@@ -6711,6 +6737,7 @@ declare namespace utilities {
         getCalibratedProbeUnitsAndValue,
         getCalibratedAspect,
         getPixelValueUnits,
+        getPixelValueUnitsImageId,
         segmentation_2 as segmentation,
         contours,
         triggerAnnotationRenderForViewportIds,

lerna.json

@@ -6,7 +6,8 @@
     "packages/adapters",
     "packages/nifti-volume-loader",
     "packages/dicomImageLoader",
-    "packages/ai"
+    "packages/ai",
+    "packages/labelmap-interpolation"
   ],
   "npmClient": "yarn"
 }

package.json

@@ -20,11 +20,12 @@
     "build:esm": "npx lerna run build:esm --stream",
     "watch": "npx lerna watch -- lerna run build --scope=$LERNA_PACKAGE_NAME --include-dependents",
     "build:update-api:ai": "cd packages/ai && npm run build:update-api",
+    "build:update-api:labelmap-interpolation": "cd packages/labelmap-interpolation && npm run build:update-api",
     "build:update-api:core": "cd packages/core && npm run build:update-api",
     "build:update-api:tools": "cd packages/tools && npm run build:update-api",
     "build:update-api:nifti": "cd packages/nifti-volume-loader && npm run build:update-api",
     "build:update-api:dicomImageLoader": "cd packages/dicomImageLoader && npm run build:update-api",
-    "build:update-api": "npm run build && npm run build:update-api:ai && npm run build:update-api:core && npm run build:update-api:tools && npm run build:update-api:nifti && npm run build:update-api:dicomImageLoader",
+    "build:update-api": "npm run build && npm run build:update-api:ai && npm run build:update-api:core && npm run build:update-api:tools && npm run build:update-api:nifti && npm run build:update-api:dicomImageLoader && npm run build:update-api:labelmap-interpolation",
     "clean": "npx lerna run clean --stream",
     "clean:deep": "npx lerna run clean:deep --stream",
     "example": "node ./utils/ExampleRunner/example-runner-cli.js",
@@ -184,5 +185,6 @@
     "not ie < 11",
     "not op_mini all"
   ],
-  "packageManager": "[email protected]+sha512.a6b2f7906b721bba3d67d4aff083df04dad64c399707841b7acf00f6b133b7ac24255f2652fa22ae3534329dc6180534e98d17432037ff6fd140556e2bb3137e"
+  "packageManager": "[email protected]+sha512.a6b2f7906b721bba3d67d4aff083df04dad64c399707841b7acf00f6b133b7ac24255f2652fa22ae3534329dc6180534e98d17432037ff6fd140556e2bb3137e",
+  "dependencies": {}
 }

packages/core/src/cache/cache.ts

@@ -16,6 +16,7 @@ import imageIdToURI from '../utilities/imageIdToURI';
 import eventTarget from '../eventTarget';
 import Events from '../enums/Events';
 import { ImageQualityStatus } from '../enums';
+import fnv1aHash from '../utilities/fnv1aHash';
 
 const ONE_GB = 1073741824;
 
@@ -33,13 +34,45 @@ class Cache {
   private readonly _imageCache = new Map<string, ICachedImage>();
   // used to store volume data (3d)
   private readonly _volumeCache = new Map<string, ICachedVolume>();
+  // used to store the reverse lookup from imageIds to volumeId
+  private readonly _imageIdsToVolumeIdCache = new Map<string, string>();
   // Todo: contour for now, but will be used for surface, etc.
   private readonly _geometryCache = new Map<string, ICachedGeometry>();
 
   private _imageCacheSize = 0;
   private _maxCacheSize = 3 * ONE_GB;
   private _geometryCacheSize = 0;
 
+  /**
+   * Generates a deterministic volume ID from a list of image IDs
+   * @param imageIds - Array of image IDs
+   * @returns A deterministic volume ID
+   */
+  public generateVolumeId(imageIds: string[]): string {
+    const imageURIs = imageIds.map(imageIdToURI).sort();
+
+    let combinedHash = 0x811c9dc5;
+    for (const id of imageURIs) {
+      const idHash = fnv1aHash(id);
+      for (let i = 0; i < idHash.length; i++) {
+        combinedHash ^= idHash.charCodeAt(i);
+        combinedHash +=
+          (combinedHash << 1) +
+          (combinedHash << 4) +
+          (combinedHash << 7) +
+          (combinedHash << 8) +
+          (combinedHash << 24);
+      }
+    }
+    return `volume-${(combinedHash >>> 0).toString(36)}`;
+  }
+
+  public getImageIdsForVolumeId(volumeId: string): string[] {
+    return Array.from(this._imageIdsToVolumeIdCache.entries())
+      .filter(([_, id]) => id === volumeId)
+      .map(([key]) => key);
+  }
+
   /**
    * Set the maximum cache Size
    *

packages/core/src/utilities/VoxelManager.ts

@@ -893,7 +893,20 @@ export default class VoxelManager<T> {
           const sliceData = new SliceDataConstructor(sliceSize);
           // @ts-ignore
           sliceData.set(scalarData.subarray(sliceStart, sliceEnd));
-          imageVoxelManager.scalarData = sliceData;
+
+          // Instead of directly assigning scalarData, use TypedArray's set method
+          // previously here we were using imageVoxelManager.scalarData = sliceData
+          // which had some weird side effects
+          if (imageVoxelManager.scalarData) {
+            imageVoxelManager.scalarData.set(sliceData);
+            // Ensure the voxel manager knows about the changes
+            imageVoxelManager.modifiedSlices.add(sliceIndex);
+          } else {
+            // Fallback to individual updates if scalarData is not directly accessible
+            for (let i = 0; i < sliceSize; i++) {
+              imageVoxelManager.setAtIndex(i, sliceData[i]);
+            }
+          }
 
           // Update min/max values for this slice
           for (let i = 0; i < sliceData.length; i++) {