Press n or j to go to the next uncovered block, b, p or k for the previous block.
| 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60 61 62 63 64 65 66 67 68 69 70 71 72 73 74 75 76 77 78 79 80 81 82 83 84 85 86 87 88 89 90 91 92 93 94 95 96 97 98 99 100 101 102 103 104 105 106 107 108 109 110 111 112 113 114 115 116 117 118 119 120 121 122 123 124 125 126 127 128 129 130 131 132 133 134 135 136 137 138 139 140 141 142 143 144 145 146 147 148 149 150 151 152 153 154 155 156 157 158 159 160 161 162 163 164 165 166 167 168 169 170 171 | 30x 30x 9x 21x 30x 30x 30x 30x 30x 30x 30x 30x 30x 30x 30x 30x 30x 30x 30x 30x 30x 30x 30x 30x 30x 39x 39x 948x 948x 24116x 24116x 21225x 21225x 21225x 21225x 21225x 21225x 24116x 948x 948x 39x 39x 30x | import { vec3 } from 'gl-matrix';
import type { Types } from '@cornerstonejs/core';
import type { vtkImageData } from '@kitware/vtk.js/Common/DataModel/ImageData';
import BoundsIJK from '../types/BoundsIJK';
export type PointInShape = {
value: number;
index: number;
pointIJK: vec3;
pointLPS: vec3 | number[];
};
export type PointInShapeCallback = ({
value,
index,
pointIJK,
pointLPS,
}: {
value: number;
index: number;
pointIJK: vec3;
pointLPS: vec3;
}) => void;
export type ShapeFnCriteria = (pointLPS: vec3, pointIJK: vec3) => boolean;
/**
* For each point in the image (If boundsIJK is not provided, otherwise, for each
* point in the provided bounding box), It runs the provided callback IF the point
* passes the provided criteria to be inside the shape (which is defined by the
* provided pointInShapeFn)
*
* @param imageData - The image data object.
* @param dimensions - The dimensions of the image.
* @param pointInShapeFn - A function that takes a point in LPS space and returns
* true if the point is in the shape and false if it is not.
* @param callback - A function that will be called for
* every point in the shape.
* @param boundsIJK - The bounds of the volume in IJK coordinates.
*/
export default function pointInShapeCallback(
imageData: vtkImageData | Types.CPUImageData,
pointInShapeFn: ShapeFnCriteria,
callback?: PointInShapeCallback,
boundsIJK?: BoundsIJK
): Array<PointInShape> {
let iMin, iMax, jMin, jMax, kMin, kMax;
let scalarData;
const { numComps } = imageData as any;
// if getScalarData is a method on imageData
if ((imageData as Types.CPUImageData).getScalarData) {
scalarData = (imageData as Types.CPUImageData).getScalarData();
} else {
scalarData = (imageData as vtkImageData)
.getPointData()
.getScalars()
.getData();
}
Iif (!scalarData) {
console.warn('No scalar data found for imageData', imageData);
return;
}
const dimensions = imageData.getDimensions();
Iif (!boundsIJK) {
iMin = 0;
iMax = dimensions[0];
jMin = 0;
jMax = dimensions[1];
kMin = 0;
kMax = dimensions[2];
} else {
[[iMin, iMax], [jMin, jMax], [kMin, kMax]] = boundsIJK;
}
const start = vec3.fromValues(iMin, jMin, kMin);
const direction = imageData.getDirection();
const rowCosines = direction.slice(0, 3);
const columnCosines = direction.slice(3, 6);
const scanAxisNormal = direction.slice(6, 9);
const spacing = imageData.getSpacing();
const [rowSpacing, columnSpacing, scanAxisSpacing] = spacing;
// @ts-ignore will be fixed in vtk-master
const worldPosStart = imageData.indexToWorld(start);
const rowStep = vec3.fromValues(
rowCosines[0] * rowSpacing,
rowCosines[1] * rowSpacing,
rowCosines[2] * rowSpacing
);
const columnStep = vec3.fromValues(
columnCosines[0] * columnSpacing,
columnCosines[1] * columnSpacing,
columnCosines[2] * columnSpacing
);
const scanAxisStep = vec3.fromValues(
scanAxisNormal[0] * scanAxisSpacing,
scanAxisNormal[1] * scanAxisSpacing,
scanAxisNormal[2] * scanAxisSpacing
);
const xMultiple =
numComps ||
scalarData.length / dimensions[2] / dimensions[1] / dimensions[0];
const yMultiple = dimensions[0] * xMultiple;
const zMultiple = dimensions[1] * yMultiple;
const pointsInShape: Array<PointInShape> = [];
const currentPos = vec3.clone(worldPosStart);
for (let k = kMin; k <= kMax; k++) {
const startPosJ = vec3.clone(currentPos);
for (let j = jMin; j <= jMax; j++) {
const startPosI = vec3.clone(currentPos);
for (let i = iMin; i <= iMax; i++) {
const pointIJK: Types.Point3 = [i, j, k];
// The current world position (pointLPS) is now in currentPos
if (pointInShapeFn(currentPos as Types.Point3, pointIJK)) {
const index = k * zMultiple + j * yMultiple + i * xMultiple;
let value;
Iif (xMultiple > 2) {
value = [
scalarData[index],
scalarData[index + 1],
scalarData[index + 2],
];
} else {
value = scalarData[index];
}
pointsInShape.push({
value,
index,
pointIJK,
pointLPS: currentPos.slice(),
});
Eif (callback) {
callback({ value, index, pointIJK, pointLPS: currentPos });
}
}
// Increment currentPos by rowStep for the next iteration
vec3.add(currentPos, currentPos, rowStep);
}
// Reset currentPos to the start of the next J line and increment by columnStep
vec3.copy(currentPos, startPosI);
vec3.add(currentPos, currentPos, columnStep);
}
// Reset currentPos to the start of the next K slice and increment by scanAxisStep
vec3.copy(currentPos, startPosJ);
vec3.add(currentPos, currentPos, scanAxisStep);
}
return pointsInShape;
}
|