Hi everyone! First time posting here, so apologies if this is unclear, I can provide extra details if needed, but I am currently working on my first big project in slicer segmenting cephalic glands in snake heads. The anatomical planes of the scans were not aligned with the correct axes to start off with, so I’ve done a quick transform to rotate the heads into the proper position, then resampled the scalar volume using the lanczos interpolation method. Pretty much just following the steps of this video: Tutorial: Re-orientation of 3D volumes in 3DSlicer by Dr Jeremy Shaw - YouTube
Having the anatomical planes be aligned with the axes of the program is pretty important, because many of these cephalic head glands and other bits used for orientation are quite difficult to identify without easily visible symmetries and familiar planes through which to orientate oneself in the digital 3D space (if that makes sense).
The issue is, after completely segmenting two snake heads, I’ve had the issue of the third head ending up with stripe artifacts during segmentation that make the process much more difficult (at least, I think that’s what these are called. Because I’ve rotated the original volume in all three axes, the result is less the stripes shown here and more triangles and other non-quadrilateral shapes, thereby making incomplete segmentation on an individual layer and bleeding through into other nearby layers). I’ve tried several solutions to this thus far, to little or no success.
- Cropping the volume. Didn’t solve the problem.
- Clicking the small button to the right of ‘Source Volume’ in segment editor, changing the source volume to the cropped version, turning on isotropic spacing, and the oversampling factor. The oversampling factor needs to be turned off to 2.5 or 3 in order to completely eliminate the artifacts, but this seems to quite quickly result in the program running out of memory and crashing, and each mark takes about a minute to process. Not ideal.
- I am aware of the Rotate To Volume Plane button in the top bar of the slice windows, but that overall is counterproductive, because when I do it, it simply just seems to undo all the transformation progress I’ve made, turning the slices back into their original orientations.
For additional help, I’ve been creating a WIP workflow to share with others working in the same lab, so these are the exact steps I’ve been following.
Files were loaded into 3DSlicer using the SlicerMorph Extension. From the modules menu, SlicerMorph → Input and Output → ImageStacks. Then, the desired reconstructed file was selected for Filename Pattern, and images other than the first in the file seemed to work better. Note that the reconstructed file tended to be nested inside the original scan file. Adjust the spacing to the pixel size (convert µm to mm), then load a preview. In the Volume Rendering module, select the eye symbol next to Volume and Display ROI, and select a preset to display in the 3D window, adjusting the shift to get a clear image. Then, adjust the ROI to be as close around the intended area to analyze as possible, reducing the file size, then return to ImageStacks to load the model at half or even full size.
Orientation of the Model
- In the Volume Rendering module, make the volume and the Region of Interest (ROI) visible by clicking the ‘eye’ buttons next to each of them
- Go the Data module. In the second tab, ‘Transform Hierarchy’, right click and select ‘Add Transform’. This should create a transformation with a name like ‘LinearTransform’. Then, back in the first tab, ‘Subject Hierarchy’, right click the grid icon at the far right of the module window for the volume being worked with (created during the previous ImageStacks step) and change the transform select the transformation created (for example, ‘LinearTransform’, in my case).
- Go to the Transforms module, and move the file you are wanding to orient to the right window by selecting it and clicking the green arrow. Then, use the Rotation (and Translation, if necessary) tools to move the volume so that the coronal, sagittal, and transverse planes of the specimen line up with the coronal, sagittal, and transverse planes of the program. When this is complete, click the file in the Transformed window, and click the ‘Solidify Transforms’ button below the two green arrows to permanently apply the transformations to the figure.
- Go to the Resample Scalar Volume module (you will have to look it up, select the lanczos option, the input volume to be the transformed one, and the output volume to be ‘Create New Volume As’. Name this with something you’ll remember, then hit apply. This might take a few minutes.
- After that, you may note that the ROI in the slice windows will likely no longer be a rectangle, and the slices will not be the squared off segments you rotated the volume to get. To fix this, hover over the top left arrow of each window, and change the plane from ‘Reformatted’ to ‘Axial,’ ‘Sagittal’, or ‘Coronal’. This should fix the issue.
- After this, cropping is important to prevent oblique slices from forming and complicating the segmentation process. Go to the “Crop Volume” Module, set the input volume to the output of your resampled scalar volume, select ‘Create New AnnotationROI As’ and name it, and select ‘Create New Volume As’ and name it. The bounding box will likely be diagonal, so return to the Data module, right click the grid icon which should be tilted, and set the transformation to ‘None’. This should bring the box in line. Then, set the white bounding box ROI around the figure, be sure the interpolator is set to linear, and click ‘Apply’