@cornerstonejs/tools 1.43.6 → 1.43.7

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

@@ -0,0 +1,4 @@
+import type { Types } from '@cornerstonejs/core';
+import type { vtkImageData } from '@kitware/vtk.js/Common/DataModel/ImageData';
+import { PointInShapeCallback } from './pointInShapeCallback';
+export default function pointInSurroundingSphereCallback(imageData: vtkImageData, circlePoints: [Types.Point3, Types.Point3], callback: PointInShapeCallback, viewport?: Types.IVolumeViewport): void;

package/dist/cjs/utilities/pointInSurroundingSphereCallback.js

@@ -0,0 +1,70 @@
+"use strict";
+var __importDefault = (this && this.__importDefault) || function (mod) {
+    return (mod && mod.__esModule) ? mod : { "default": mod };
+};
+Object.defineProperty(exports, "__esModule", { value: true });
+const core_1 = require("@cornerstonejs/core");
+const gl_matrix_1 = require("gl-matrix");
+const sphere_1 = require("./math/sphere");
+const pointInShapeCallback_1 = __importDefault(require("./pointInShapeCallback"));
+const boundingBox_1 = require("./boundingBox");
+const { transformWorldToIndex } = core_1.utilities;
+function pointInSurroundingSphereCallback(imageData, circlePoints, callback, viewport) {
+    const { boundsIJK, centerWorld, radiusWorld } = _getBounds(circlePoints, imageData, viewport);
+    const sphereObj = {
+        center: centerWorld,
+        radius: radiusWorld,
+    };
+    (0, pointInShapeCallback_1.default)(imageData, (pointLPS) => (0, sphere_1.pointInSphere)(sphereObj, pointLPS), callback, boundsIJK);
+}
+exports.default = pointInSurroundingSphereCallback;
+function _getBounds(circlePoints, imageData, viewport) {
+    const [bottom, top] = circlePoints;
+    const centerWorld = gl_matrix_1.vec3.fromValues((bottom[0] + top[0]) / 2, (bottom[1] + top[1]) / 2, (bottom[2] + top[2]) / 2);
+    const radiusWorld = gl_matrix_1.vec3.distance(bottom, top) / 2;
+    let boundsIJK;
+    if (!viewport) {
+        const centerIJK = transformWorldToIndex(imageData, centerWorld);
+        const spacings = imageData.getSpacing();
+        const minSpacing = Math.min(...spacings);
+        const maxRadiusIJK = Math.ceil(radiusWorld / minSpacing);
+        boundsIJK = [
+            [centerIJK[0] - maxRadiusIJK, centerIJK[0] + maxRadiusIJK],
+            [centerIJK[1] - maxRadiusIJK, centerIJK[1] + maxRadiusIJK],
+            [centerIJK[2] - maxRadiusIJK, centerIJK[2] + maxRadiusIJK],
+        ];
+        return {
+            boundsIJK,
+            centerWorld: centerWorld,
+            radiusWorld,
+        };
+    }
+    boundsIJK = _computeBoundsIJKWithCamera(imageData, viewport, circlePoints, centerWorld, radiusWorld);
+    return {
+        boundsIJK,
+        centerWorld: centerWorld,
+        radiusWorld,
+    };
+}
+function _computeBoundsIJKWithCamera(imageData, viewport, circlePoints, centerWorld, radiusWorld) {
+    const [bottom, top] = circlePoints;
+    const dimensions = imageData.getDimensions();
+    const camera = viewport.getCamera();
+    const viewUp = gl_matrix_1.vec3.fromValues(camera.viewUp[0], camera.viewUp[1], camera.viewUp[2]);
+    const viewPlaneNormal = gl_matrix_1.vec3.fromValues(camera.viewPlaneNormal[0], camera.viewPlaneNormal[1], camera.viewPlaneNormal[2]);
+    const viewRight = gl_matrix_1.vec3.create();
+    gl_matrix_1.vec3.cross(viewRight, viewUp, viewPlaneNormal);
+    const topLeftWorld = gl_matrix_1.vec3.create();
+    const bottomRightWorld = gl_matrix_1.vec3.create();
+    gl_matrix_1.vec3.scaleAndAdd(topLeftWorld, top, viewPlaneNormal, radiusWorld);
+    gl_matrix_1.vec3.scaleAndAdd(bottomRightWorld, bottom, viewPlaneNormal, -radiusWorld);
+    gl_matrix_1.vec3.scaleAndAdd(topLeftWorld, topLeftWorld, viewRight, -radiusWorld);
+    gl_matrix_1.vec3.scaleAndAdd(bottomRightWorld, bottomRightWorld, viewRight, radiusWorld);
+    const sphereCornersIJK = [
+        transformWorldToIndex(imageData, topLeftWorld),
+        (transformWorldToIndex(imageData, bottomRightWorld)),
+    ];
+    const boundsIJK = (0, boundingBox_1.getBoundingBoxAroundShape)(sphereCornersIJK, dimensions);
+    return boundsIJK;
+}
+//# sourceMappingURL=pointInSurroundingSphereCallback.js.map

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

@@ -0,0 +1 @@
+{"version":3,"file":"pointInSurroundingSphereCallback.js","sourceRoot":"","sources":["../../../src/utilities/pointInSurroundingSphereCallback.ts"],"names":[],"mappings":";;;;;AAAA,8CAA2D;AAI3D,yCAAiC;AACjC,0CAA8C;AAC9C,kFAEgC;AAEhC,+CAA0D;AAE1D,MAAM,EAAE,qBAAqB,EAAE,GAAG,gBAAO,CAAC;AAiB1C,SAAwB,gCAAgC,CACtD,SAAuB,EACvB,YAA0C,EAC1C,QAA8B,EAC9B,QAAgC;IAMhC,MAAM,EAAE,SAAS,EAAE,WAAW,EAAE,WAAW,EAAE,GAAG,UAAU,CACxD,YAAY,EACZ,SAAS,EACT,QAAQ,CACT,CAAC;IAEF,MAAM,SAAS,GAAG;QAChB,MAAM,EAAE,WAAW;QACnB,MAAM,EAAE,WAAW;KACpB,CAAC;IAEF,IAAA,8BAAoB,EAClB,SAAS,EACT,CAAC,QAAQ,EAAE,EAAE,CAAC,IAAA,sBAAa,EAAC,SAAS,EAAE,QAAQ,CAAC,EAChD,QAAQ,EACR,SAAS,CACV,CAAC;AACJ,CAAC;AA3BD,mDA2BC;AAED,SAAS,UAAU,CACjB,YAA0C,EAC1C,SAAuB,EACvB,QAAQ;IAMR,MAAM,CAAC,MAAM,EAAE,GAAG,CAAC,GAAG,YAAY,CAAC;IAGnC,MAAM,WAAW,GAAG,gBAAI,CAAC,UAAU,CACjC,CAAC,MAAM,CAAC,CAAC,CAAC,GAAG,GAAG,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,EACxB,CAAC,MAAM,CAAC,CAAC,CAAC,GAAG,GAAG,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,EACxB,CAAC,MAAM,CAAC,CAAC,CAAC,GAAG,GAAG,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CACzB,CAAC;IAGF,MAAM,WAAW,GAAG,gBAAI,CAAC,QAAQ,CAAC,MAAM,EAAE,GAAG,CAAC,GAAG,CAAC,CAAC;IAEnD,IAAI,SAAS,CAAC;IAEd,IAAI,CAAC,QAAQ,EAAE;QAUb,MAAM,SAAS,GAAG,qBAAqB,CACrC,SAAS,EACT,WAA2B,CAC5B,CAAC;QAEF,MAAM,QAAQ,GAAG,SAAS,CAAC,UAAU,EAAE,CAAC;QACxC,MAAM,UAAU,GAAG,IAAI,CAAC,GAAG,CAAC,GAAG,QAAQ,CAAC,CAAC;QAEzC,MAAM,YAAY,GAAG,IAAI,CAAC,IAAI,CAAC,WAAW,GAAG,UAAU,CAAC,CAAC;QAEzD,SAAS,GAAG;YACV,CAAC,SAAS,CAAC,CAAC,CAAC,GAAG,YAAY,EAAE,SAAS,CAAC,CAAC,CAAC,GAAG,YAAY,CAAC;YAC1D,CAAC,SAAS,CAAC,CAAC,CAAC,GAAG,YAAY,EAAE,SAAS,CAAC,CAAC,CAAC,GAAG,YAAY,CAAC;YAC1D,CAAC,SAAS,CAAC,CAAC,CAAC,GAAG,YAAY,EAAE,SAAS,CAAC,CAAC,CAAC,GAAG,YAAY,CAAC;SAC3D,CAAC;QAEF,OAAO;YACL,SAAS;YACT,WAAW,EAAE,WAA2B;YACxC,WAAW;SACZ,CAAC;KACH;IAED,SAAS,GAAG,2BAA2B,CACrC,SAAS,EACT,QAAQ,EACR,YAAY,EACZ,WAAW,EACX,WAAW,CACZ,CAAC;IAEF,OAAO;QACL,SAAS;QACT,WAAW,EAAE,WAA2B;QACxC,WAAW;KACZ,CAAC;AACJ,CAAC;AAED,SAAS,2BAA2B,CAClC,SAAS,EACT,QAAQ,EACR,YAAY,EACZ,WAAW,EACX,WAAW;IAEX,MAAM,CAAC,MAAM,EAAE,GAAG,CAAC,GAAG,YAAY,CAAC;IAEnC,MAAM,UAAU,GAAG,SAAS,CAAC,aAAa,EAAkB,CAAC;IAC7D,MAAM,MAAM,GAAG,QAAQ,CAAC,SAAS,EAAE,CAAC;IAKpC,MAAM,MAAM,GAAG,gBAAI,CAAC,UAAU,CAC5B,MAAM,CAAC,MAAM,CAAC,CAAC,CAAC,EAChB,MAAM,CAAC,MAAM,CAAC,CAAC,CAAC,EAChB,MAAM,CAAC,MAAM,CAAC,CAAC,CAAC,CACjB,CAAC;IACF,MAAM,eAAe,GAAG,gBAAI,CAAC,UAAU,CACrC,MAAM,CAAC,eAAe,CAAC,CAAC,CAAC,EACzB,MAAM,CAAC,eAAe,CAAC,CAAC,CAAC,EACzB,MAAM,CAAC,eAAe,CAAC,CAAC,CAAC,CAC1B,CAAC;IACF,MAAM,SAAS,GAAG,gBAAI,CAAC,MAAM,EAAE,CAAC;IAEhC,gBAAI,CAAC,KAAK,CAAC,SAAS,EAAE,MAAM,EAAE,eAAe,CAAC,CAAC;IAQ/C,MAAM,YAAY,GAAG,gBAAI,CAAC,MAAM,EAAE,CAAC;IACnC,MAAM,gBAAgB,GAAG,gBAAI,CAAC,MAAM,EAAE,CAAC;IAEvC,gBAAI,CAAC,WAAW,CAAC,YAAY,EAAE,GAAG,EAAE,eAAe,EAAE,WAAW,CAAC,CAAC;IAClE,gBAAI,CAAC,WAAW,CAAC,gBAAgB,EAAE,MAAM,EAAE,eAAe,EAAE,CAAC,WAAW,CAAC,CAAC;IAG1E,gBAAI,CAAC,WAAW,CAAC,YAAY,EAAE,YAAY,EAAE,SAAS,EAAE,CAAC,WAAW,CAAC,CAAC;IACtE,gBAAI,CAAC,WAAW,CAAC,gBAAgB,EAAE,gBAAgB,EAAE,SAAS,EAAE,WAAW,CAAC,CAAC;IAI7E,MAAM,gBAAgB,GAAG;QACT,qBAAqB,CAAC,SAAS,EAAgB,YAAY,CAAC;QAC5D,CACZ,qBAAqB,CAAC,SAAS,EAAgB,gBAAgB,CAAC,CACjE;KACF,CAAC;IAGF,MAAM,SAAS,GAAG,IAAA,uCAAyB,EAAC,gBAAgB,EAAE,UAAU,CAAC,CAAC;IAE1E,OAAO,SAAS,CAAC;AACnB,CAAC"}

package/dist/esm/utilities/pointInSurroundingSphereCallback.js

@@ -0,0 +1,64 @@
+import { utilities as csUtils } from '@cornerstonejs/core';
+import { vec3 } from 'gl-matrix';
+import { pointInSphere } from './math/sphere';
+import pointInShapeCallback from './pointInShapeCallback';
+import { getBoundingBoxAroundShape } from './boundingBox';
+const { transformWorldToIndex } = csUtils;
+export default function pointInSurroundingSphereCallback(imageData, circlePoints, callback, viewport) {
+    const { boundsIJK, centerWorld, radiusWorld } = _getBounds(circlePoints, imageData, viewport);
+    const sphereObj = {
+        center: centerWorld,
+        radius: radiusWorld,
+    };
+    pointInShapeCallback(imageData, (pointLPS) => pointInSphere(sphereObj, pointLPS), callback, boundsIJK);
+}
+function _getBounds(circlePoints, imageData, viewport) {
+    const [bottom, top] = circlePoints;
+    const centerWorld = vec3.fromValues((bottom[0] + top[0]) / 2, (bottom[1] + top[1]) / 2, (bottom[2] + top[2]) / 2);
+    const radiusWorld = vec3.distance(bottom, top) / 2;
+    let boundsIJK;
+    if (!viewport) {
+        const centerIJK = transformWorldToIndex(imageData, centerWorld);
+        const spacings = imageData.getSpacing();
+        const minSpacing = Math.min(...spacings);
+        const maxRadiusIJK = Math.ceil(radiusWorld / minSpacing);
+        boundsIJK = [
+            [centerIJK[0] - maxRadiusIJK, centerIJK[0] + maxRadiusIJK],
+            [centerIJK[1] - maxRadiusIJK, centerIJK[1] + maxRadiusIJK],
+            [centerIJK[2] - maxRadiusIJK, centerIJK[2] + maxRadiusIJK],
+        ];
+        return {
+            boundsIJK,
+            centerWorld: centerWorld,
+            radiusWorld,
+        };
+    }
+    boundsIJK = _computeBoundsIJKWithCamera(imageData, viewport, circlePoints, centerWorld, radiusWorld);
+    return {
+        boundsIJK,
+        centerWorld: centerWorld,
+        radiusWorld,
+    };
+}
+function _computeBoundsIJKWithCamera(imageData, viewport, circlePoints, centerWorld, radiusWorld) {
+    const [bottom, top] = circlePoints;
+    const dimensions = imageData.getDimensions();
+    const camera = viewport.getCamera();
+    const viewUp = vec3.fromValues(camera.viewUp[0], camera.viewUp[1], camera.viewUp[2]);
+    const viewPlaneNormal = vec3.fromValues(camera.viewPlaneNormal[0], camera.viewPlaneNormal[1], camera.viewPlaneNormal[2]);
+    const viewRight = vec3.create();
+    vec3.cross(viewRight, viewUp, viewPlaneNormal);
+    const topLeftWorld = vec3.create();
+    const bottomRightWorld = vec3.create();
+    vec3.scaleAndAdd(topLeftWorld, top, viewPlaneNormal, radiusWorld);
+    vec3.scaleAndAdd(bottomRightWorld, bottom, viewPlaneNormal, -radiusWorld);
+    vec3.scaleAndAdd(topLeftWorld, topLeftWorld, viewRight, -radiusWorld);
+    vec3.scaleAndAdd(bottomRightWorld, bottomRightWorld, viewRight, radiusWorld);
+    const sphereCornersIJK = [
+        transformWorldToIndex(imageData, topLeftWorld),
+        (transformWorldToIndex(imageData, bottomRightWorld)),
+    ];
+    const boundsIJK = getBoundingBoxAroundShape(sphereCornersIJK, dimensions);
+    return boundsIJK;
+}
+//# sourceMappingURL=pointInSurroundingSphereCallback.js.map

package/dist/esm/utilities/pointInSurroundingSphereCallback.js.map

@@ -0,0 +1 @@
+{"version":3,"file":"pointInSurroundingSphereCallback.js","sourceRoot":"","sources":["../../../src/utilities/pointInSurroundingSphereCallback.ts"],"names":[],"mappings":"AAAA,OAAO,EAAE,SAAS,IAAI,OAAO,EAAE,MAAM,qBAAqB,CAAC;AAI3D,OAAO,EAAE,IAAI,EAAE,MAAM,WAAW,CAAC;AACjC,OAAO,EAAE,aAAa,EAAE,MAAM,eAAe,CAAC;AAC9C,OAAO,oBAEN,MAAM,wBAAwB,CAAC;AAEhC,OAAO,EAAE,yBAAyB,EAAE,MAAM,eAAe,CAAC;AAE1D,MAAM,EAAE,qBAAqB,EAAE,GAAG,OAAO,CAAC;AAiB1C,MAAM,CAAC,OAAO,UAAU,gCAAgC,CACtD,SAAuB,EACvB,YAA0C,EAC1C,QAA8B,EAC9B,QAAgC;IAMhC,MAAM,EAAE,SAAS,EAAE,WAAW,EAAE,WAAW,EAAE,GAAG,UAAU,CACxD,YAAY,EACZ,SAAS,EACT,QAAQ,CACT,CAAC;IAEF,MAAM,SAAS,GAAG;QAChB,MAAM,EAAE,WAAW;QACnB,MAAM,EAAE,WAAW;KACpB,CAAC;IAEF,oBAAoB,CAClB,SAAS,EACT,CAAC,QAAQ,EAAE,EAAE,CAAC,aAAa,CAAC,SAAS,EAAE,QAAQ,CAAC,EAChD,QAAQ,EACR,SAAS,CACV,CAAC;AACJ,CAAC;AAED,SAAS,UAAU,CACjB,YAA0C,EAC1C,SAAuB,EACvB,QAAQ;IAMR,MAAM,CAAC,MAAM,EAAE,GAAG,CAAC,GAAG,YAAY,CAAC;IAGnC,MAAM,WAAW,GAAG,IAAI,CAAC,UAAU,CACjC,CAAC,MAAM,CAAC,CAAC,CAAC,GAAG,GAAG,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,EACxB,CAAC,MAAM,CAAC,CAAC,CAAC,GAAG,GAAG,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,EACxB,CAAC,MAAM,CAAC,CAAC,CAAC,GAAG,GAAG,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CACzB,CAAC;IAGF,MAAM,WAAW,GAAG,IAAI,CAAC,QAAQ,CAAC,MAAM,EAAE,GAAG,CAAC,GAAG,CAAC,CAAC;IAEnD,IAAI,SAAS,CAAC;IAEd,IAAI,CAAC,QAAQ,EAAE;QAUb,MAAM,SAAS,GAAG,qBAAqB,CACrC,SAAS,EACT,WAA2B,CAC5B,CAAC;QAEF,MAAM,QAAQ,GAAG,SAAS,CAAC,UAAU,EAAE,CAAC;QACxC,MAAM,UAAU,GAAG,IAAI,CAAC,GAAG,CAAC,GAAG,QAAQ,CAAC,CAAC;QAEzC,MAAM,YAAY,GAAG,IAAI,CAAC,IAAI,CAAC,WAAW,GAAG,UAAU,CAAC,CAAC;QAEzD,SAAS,GAAG;YACV,CAAC,SAAS,CAAC,CAAC,CAAC,GAAG,YAAY,EAAE,SAAS,CAAC,CAAC,CAAC,GAAG,YAAY,CAAC;YAC1D,CAAC,SAAS,CAAC,CAAC,CAAC,GAAG,YAAY,EAAE,SAAS,CAAC,CAAC,CAAC,GAAG,YAAY,CAAC;YAC1D,CAAC,SAAS,CAAC,CAAC,CAAC,GAAG,YAAY,EAAE,SAAS,CAAC,CAAC,CAAC,GAAG,YAAY,CAAC;SAC3D,CAAC;QAEF,OAAO;YACL,SAAS;YACT,WAAW,EAAE,WAA2B;YACxC,WAAW;SACZ,CAAC;KACH;IAED,SAAS,GAAG,2BAA2B,CACrC,SAAS,EACT,QAAQ,EACR,YAAY,EACZ,WAAW,EACX,WAAW,CACZ,CAAC;IAEF,OAAO;QACL,SAAS;QACT,WAAW,EAAE,WAA2B;QACxC,WAAW;KACZ,CAAC;AACJ,CAAC;AAED,SAAS,2BAA2B,CAClC,SAAS,EACT,QAAQ,EACR,YAAY,EACZ,WAAW,EACX,WAAW;IAEX,MAAM,CAAC,MAAM,EAAE,GAAG,CAAC,GAAG,YAAY,CAAC;IAEnC,MAAM,UAAU,GAAG,SAAS,CAAC,aAAa,EAAkB,CAAC;IAC7D,MAAM,MAAM,GAAG,QAAQ,CAAC,SAAS,EAAE,CAAC;IAKpC,MAAM,MAAM,GAAG,IAAI,CAAC,UAAU,CAC5B,MAAM,CAAC,MAAM,CAAC,CAAC,CAAC,EAChB,MAAM,CAAC,MAAM,CAAC,CAAC,CAAC,EAChB,MAAM,CAAC,MAAM,CAAC,CAAC,CAAC,CACjB,CAAC;IACF,MAAM,eAAe,GAAG,IAAI,CAAC,UAAU,CACrC,MAAM,CAAC,eAAe,CAAC,CAAC,CAAC,EACzB,MAAM,CAAC,eAAe,CAAC,CAAC,CAAC,EACzB,MAAM,CAAC,eAAe,CAAC,CAAC,CAAC,CAC1B,CAAC;IACF,MAAM,SAAS,GAAG,IAAI,CAAC,MAAM,EAAE,CAAC;IAEhC,IAAI,CAAC,KAAK,CAAC,SAAS,EAAE,MAAM,EAAE,eAAe,CAAC,CAAC;IAQ/C,MAAM,YAAY,GAAG,IAAI,CAAC,MAAM,EAAE,CAAC;IACnC,MAAM,gBAAgB,GAAG,IAAI,CAAC,MAAM,EAAE,CAAC;IAEvC,IAAI,CAAC,WAAW,CAAC,YAAY,EAAE,GAAG,EAAE,eAAe,EAAE,WAAW,CAAC,CAAC;IAClE,IAAI,CAAC,WAAW,CAAC,gBAAgB,EAAE,MAAM,EAAE,eAAe,EAAE,CAAC,WAAW,CAAC,CAAC;IAG1E,IAAI,CAAC,WAAW,CAAC,YAAY,EAAE,YAAY,EAAE,SAAS,EAAE,CAAC,WAAW,CAAC,CAAC;IACtE,IAAI,CAAC,WAAW,CAAC,gBAAgB,EAAE,gBAAgB,EAAE,SAAS,EAAE,WAAW,CAAC,CAAC;IAI7E,MAAM,gBAAgB,GAAG;QACT,qBAAqB,CAAC,SAAS,EAAgB,YAAY,CAAC;QAC5D,CACZ,qBAAqB,CAAC,SAAS,EAAgB,gBAAgB,CAAC,CACjE;KACF,CAAC;IAGF,MAAM,SAAS,GAAG,yBAAyB,CAAC,gBAAgB,EAAE,UAAU,CAAC,CAAC;IAE1E,OAAO,SAAS,CAAC;AACnB,CAAC"}

package/dist/types/utilities/pointInSurroundingSphereCallback.d.ts

@@ -0,0 +1,5 @@
+import type { Types } from '@cornerstonejs/core';
+import type { vtkImageData } from '@kitware/vtk.js/Common/DataModel/ImageData';
+import { PointInShapeCallback } from './pointInShapeCallback';
+export default function pointInSurroundingSphereCallback(imageData: vtkImageData, circlePoints: [Types.Point3, Types.Point3], callback: PointInShapeCallback, viewport?: Types.IVolumeViewport): void;
+//# sourceMappingURL=pointInSurroundingSphereCallback.d.ts.map

package/dist/types/utilities/pointInSurroundingSphereCallback.d.ts.map

@@ -0,0 +1 @@
+{"version":3,"file":"pointInSurroundingSphereCallback.d.ts","sourceRoot":"","sources":["../../../src/utilities/pointInSurroundingSphereCallback.ts"],"names":[],"mappings":"AACA,OAAO,KAAK,EAAE,KAAK,EAAE,MAAM,qBAAqB,CAAC;AAEjD,OAAO,KAAK,EAAE,YAAY,EAAE,MAAM,4CAA4C,CAAC;AAG/E,OAA6B,EAC3B,oBAAoB,EACrB,MAAM,wBAAwB,CAAC;AAqBhC,MAAM,CAAC,OAAO,UAAU,gCAAgC,CACtD,SAAS,EAAE,YAAY,EACvB,YAAY,EAAE,CAAC,KAAK,CAAC,MAAM,EAAE,KAAK,CAAC,MAAM,CAAC,EAC1C,QAAQ,EAAE,oBAAoB,EAC9B,QAAQ,CAAC,EAAE,KAAK,CAAC,eAAe,GAC/B,IAAI,CAsBN"}

package/package.json

@@ -1,6 +1,6 @@
 {
   "name": "@cornerstonejs/tools",
-  "version": "1.43.6",
+  "version": "1.43.7",
   "description": "Cornerstone3D Tools",
   "main": "src/index.ts",
   "types": "dist/types/index.d.ts",
@@ -29,7 +29,7 @@
     "webpack:watch": "webpack --mode development --progress --watch  --config ./.webpack/webpack.dev.js"
   },
   "dependencies": {
-    "@cornerstonejs/core": "^1.43.6",
+    "@cornerstonejs/core": "^1.43.7",
     "comlink": "^4.4.1",
     "lodash.clonedeep": "4.5.0",
     "lodash.get": "^4.4.2"
@@ -53,5 +53,5 @@
     "type": "individual",
     "url": "https://ohif.org/donate"
   },
-  "gitHead": "66da8823e7582d11c671c88902b2b2be8d6f9617"
+  "gitHead": "2d0c559929aa0524bf303922230ded50ac0ef1df"
 }

package/src/utilities/pointInSurroundingSphereCallback.ts

@@ -0,0 +1,188 @@
+import { utilities as csUtils } from '@cornerstonejs/core';
+import type { Types } from '@cornerstonejs/core';
+
+import type { vtkImageData } from '@kitware/vtk.js/Common/DataModel/ImageData';
+import { vec3 } from 'gl-matrix';
+import { pointInSphere } from './math/sphere';
+import pointInShapeCallback, {
+  PointInShapeCallback,
+} from './pointInShapeCallback';
+import { BoundsIJK } from '../types';
+import { getBoundingBoxAroundShape } from './boundingBox';
+
+const { transformWorldToIndex } = csUtils;
+
+/**
+ * Given an imageData, and the great circle top and bottom points of a sphere,
+ * this function will run the callback for each point of the imageData that is
+ * within the sphere defined by the great circle points. If the viewport
+ * is provided, region of interest will be an accurate approximation of the
+ * sphere (using viewport camera), and the resulting performance will be
+ * better.
+ *
+ * @privateRemarks great circle also known as orthodrome is the intersection of
+ * the sphere and the plane that passes through the center of the sphere
+ *
+ * @param imageData - The volume imageData
+ * @param circlePoints - bottom and top points of the great circle in world coordinates
+ * @param callback - A callback function that will be called for each point in the shape.
+ */
+export default function pointInSurroundingSphereCallback(
+  imageData: vtkImageData,
+  circlePoints: [Types.Point3, Types.Point3],
+  callback: PointInShapeCallback,
+  viewport?: Types.IVolumeViewport
+): void {
+  // We can run the sphere equation to determine if a point is inside
+  // the sphere; however, since the imageData dimensions can be quite large, we
+  // can narrow down the search by estimating the bounds of the sphere in index
+  // space.
+  const { boundsIJK, centerWorld, radiusWorld } = _getBounds(
+    circlePoints,
+    imageData,
+    viewport
+  );
+
+  const sphereObj = {
+    center: centerWorld,
+    radius: radiusWorld,
+  };
+
+  pointInShapeCallback(
+    imageData,
+    (pointLPS) => pointInSphere(sphereObj, pointLPS),
+    callback,
+    boundsIJK
+  );
+}
+
+function _getBounds(
+  circlePoints: [Types.Point3, Types.Point3],
+  imageData: vtkImageData,
+  viewport
+): {
+  boundsIJK: BoundsIJK;
+  centerWorld: Types.Point3;
+  radiusWorld: number;
+} {
+  const [bottom, top] = circlePoints;
+
+  // Sphere center in world
+  const centerWorld = vec3.fromValues(
+    (bottom[0] + top[0]) / 2,
+    (bottom[1] + top[1]) / 2,
+    (bottom[2] + top[2]) / 2
+  );
+
+  // sphere radius in world
+  const radiusWorld = vec3.distance(bottom, top) / 2;
+
+  let boundsIJK;
+
+  if (!viewport) {
+    // If no viewport is provide (no camera), we can estimate the bounding box
+    // of the sphere in index space.
+    // This is done by calculating the maximum value for radius in the index
+    // space (since the radius is in world space, we need to convert it to index, and
+    // each dimensions can have a different scale factor). Therefore, by finding
+    // the minimum spacing value in the imageData, we can calculate the maximum
+    // radius in index space and use that to calculate the bounds of the sphere
+    // This will not be accurate, but it is a good first approximation.
+    // sphere center in index
+    const centerIJK = transformWorldToIndex(
+      imageData,
+      centerWorld as Types.Point3
+    );
+
+    const spacings = imageData.getSpacing();
+    const minSpacing = Math.min(...spacings);
+
+    const maxRadiusIJK = Math.ceil(radiusWorld / minSpacing);
+
+    boundsIJK = [
+      [centerIJK[0] - maxRadiusIJK, centerIJK[0] + maxRadiusIJK],
+      [centerIJK[1] - maxRadiusIJK, centerIJK[1] + maxRadiusIJK],
+      [centerIJK[2] - maxRadiusIJK, centerIJK[2] + maxRadiusIJK],
+    ];
+
+    return {
+      boundsIJK,
+      centerWorld: centerWorld as Types.Point3,
+      radiusWorld,
+    };
+  }
+
+  boundsIJK = _computeBoundsIJKWithCamera(
+    imageData,
+    viewport,
+    circlePoints,
+    centerWorld,
+    radiusWorld
+  );
+
+  return {
+    boundsIJK,
+    centerWorld: centerWorld as Types.Point3,
+    radiusWorld,
+  };
+}
+
+function _computeBoundsIJKWithCamera(
+  imageData,
+  viewport,
+  circlePoints,
+  centerWorld,
+  radiusWorld
+) {
+  const [bottom, top] = circlePoints;
+
+  const dimensions = imageData.getDimensions() as Types.Point3;
+  const camera = viewport.getCamera();
+
+  // Calculate viewRight from the camera, this will get used in order to
+  // calculate circles topLeft and bottomRight on different planes of intersection
+  // between sphere and viewPlane
+  const viewUp = vec3.fromValues(
+    camera.viewUp[0],
+    camera.viewUp[1],
+    camera.viewUp[2]
+  );
+  const viewPlaneNormal = vec3.fromValues(
+    camera.viewPlaneNormal[0],
+    camera.viewPlaneNormal[1],
+    camera.viewPlaneNormal[2]
+  );
+  const viewRight = vec3.create();
+
+  vec3.cross(viewRight, viewUp, viewPlaneNormal);
+
+  // we need to find the bounding box of the sphere in the image, e.g., the
+  // topLeftWorld and bottomRightWorld points of the bounding box.
+  // We go from the sphereCenter in the normal direction of amount radius, and
+  // we go left to find the topLeftWorld point of the bounding box. Next we go
+  // in the opposite direction and go right to find the bottomRightWorld point
+  // of the bounding box.
+  const topLeftWorld = vec3.create();
+  const bottomRightWorld = vec3.create();
+
+  vec3.scaleAndAdd(topLeftWorld, top, viewPlaneNormal, radiusWorld);
+  vec3.scaleAndAdd(bottomRightWorld, bottom, viewPlaneNormal, -radiusWorld);
+
+  // go in the direction of viewRight with the value of radius
+  vec3.scaleAndAdd(topLeftWorld, topLeftWorld, viewRight, -radiusWorld);
+  vec3.scaleAndAdd(bottomRightWorld, bottomRightWorld, viewRight, radiusWorld);
+
+  // convert the world coordinates to index coordinates
+
+  const sphereCornersIJK = [
+    <Types.Point3>transformWorldToIndex(imageData, <Types.Point3>topLeftWorld),
+    <Types.Point3>(
+      transformWorldToIndex(imageData, <Types.Point3>bottomRightWorld)
+    ),
+  ];
+
+  // get the bounding box of the sphere in the image
+  const boundsIJK = getBoundingBoxAroundShape(sphereCornersIJK, dimensions);
+
+  return boundsIJK;
+}