Differences between revisions 19 and 20
Revision 19 as of 2015-11-30 19:54:08
Size: 8456
Comment:
Revision 20 as of 2015-11-30 21:48:45
Size: 10224
Comment:
Deletions are marked like this. Additions are marked like this.
Line 27: Line 27:
== Segmentation using Slicer == == Slicer Procedures ==
Line 41: Line 41:
Load MR image data (NifTI): Load MRI data:
Line 61: Line 61:
 * Select Segmented label map image (NifTI, .nii)  * Select Segmented label map image(s) (NifTI, .nii)
Line 63: Line 63:
 * Click/check 'Show Options'
 * Click/check '!LabelMap'
Line 81: Line 83:
 * Rename the triangulated surface as desired  * Rename the triangulated surface as desired, see naming convention: ["Specifications/DataManagement"]
Line 84: Line 86:
=== Viewing Options ===

In any of the 2D slice viewer windows, using the default mouse cursor pointer tool, click and drag up/down and left/right to change contrast and brightness of MRI
Line 92: Line 98:

 * Editor Module
  * Applying Growth Cut Effect (region growing algorithm using cellular automata):
   * Use Draw Effect or Paint Effect to highlight anatomy of interest (foreground) using specific label (label1).
   * Use Draw Effect or Paint Effect to highlight area around anatomy of interest (background) with different label (label2).
   * Use Paint Effect to "cap" boundaries of anatomy of interest at two slices (slice before anatomy of interest appears on MRI and slice after anatomy of interest appears on MRI).
    * See figures below for example of segmentation using GrowthCut Effect with foreground (green) and background (yellow).
 * From Modules drop-down menu, select Editor
 * If first time entering Editor Module, click Apply to select the default !GenericAnatomyColors to create a new Merge Volume
  NOTE: if not first time, make sure label map displayed (i.e. being edited) in 2D slice viewers is the same as the Merge Volume.
 * Perform the following steps on roughly equally-spaced slices in each image plane orientation as appropriate (coronal, sagittal, axial).
  * Use Paint Effect (under 'Edit Selected Label Map' tab, top row of icons, 3rd icon from left) [with an appropriate/largest radius] to highlight anatomy of interest (foreground) using a specific label (e.g. label1).
  * Use Paint Effect to highlight area around anatomy of interest (background) with a different label (e.g. label2).
  * Use Paint Effect to "cap" boundaries of anatomy of interest (using background label) to prevent foreground leakage (before/after anatomy of interest appears).
 * Select !GrowCutEffect icon under 'Edit Selected Label Map' tab, bottom row of icons, 3rd icon from left.
 * Click Apply button. [window graphics will freeze until grow cut algorithm is complete.

 * See figures below for example of segmentation using GrowthCut Effect with foreground (green) and background (yellow).
Line 102: Line 113:
attachment:gc4.png
Line 104: Line 114:
   * Repeat for five slices at varying depths for each view (coronal, sagittal, axial).
Line 107: Line 116:
   * Manually segmentation to clean up results of Growth Cut Effect.
   *** For more information on Growth Cut Effect refer to the 3DSlicerWiki [http://www.slicer.org/slicerWiki/index.php/Modules:GrowCutSegmentation-Documentation-3.6].

attachment:gc4.png

   * For more information on Growth Cut Effect refer to the 3DSlicerWiki [http://www.slicer.org/slicerWiki/index.php/Modules:GrowCutSegmentation-Documentation-3.6].
Line 110: Line 121:
  * Generating raw volume from label for anatomy of interest.
   * In Make Model Effect, choose model name, deselect smooth model and select apply.
   * Creates vtk file.
   * To save as surface representation, save vtk file as STL file.
Line 134: Line 140:
 * From the Modules drop-down menu, select 'Surface Models' -> 'Label Map Smoothing'
 * Set Gaussian Smoothing Parameters, Sigma based on image resolution (e.g. 0.7 for cartilage MRI, resolution = 0.35 x 0.35 x 0.7 mm)
 * From 'Input Volume' drop-down menu, select desired label map
 * From 'Output Volume' drop-down menu, select 'Create and rename new Volume'
 * Enter new name as desired, click 'OK'
 * Click Apply

==== Generate Triangulated Surface from Label Map ====

 * From the Modules drop-down menu, select 'Surface Models' -> 'Model Maker'
  * Under 'IO' tab, from 'Input Volume' drop-down menu, select desired Label Map
  * In Model drop-down menu, either 'Create and rename new Hierarchy' (e.g. default 'Models'), or select an existing Model Hierarchy.
   NOTE: all models can be generated under the same or separate model hierarchies
  * Under 'Create Model' tab, specify a desired 'Model Name' for the select Label Map.
  * If no smoothing is desired, under 'Model Maker Parameters set 'Smooth' to zero.
  * Click Apply button.
   
 * From Editor:
  * Select the !MakeModelEffect
  * Set desired Label Map as the Merge Volume
  * If no smoothing is desired, deselect 'smooth model'
  * Click Apply button.
Line 150: Line 178:
Input: cartilage MRI INPUT: cartilage MRI
Line 161: Line 189:
Output labels: FMB (femur), TBB (tibia), FBB (fibula), PTB (patella). OUTPUT LABELS: FMB (femur), TBB (tibia), FBB (fibula), PTB (patella).
Line 165: Line 193:
INPUT: cartilage MRI
Line 167: Line 197:
Tissue labels include FMC (femoral), TBC-M (medial tibial), TBC-L (lateral tibial), PTC (patellar). OUTPUT LABELS: FMC (femoral), TBC-M (medial tibial), TBC-L (lateral tibial), PTC (patellar).

TableOfContents()

Target Outcome

This specification targets at volumetric reconstruction of a tissue of interest, specifically the definition of the boundaries of the tissue,

  • as an image volume, and
  • as a surface representation.

Prerequisites

Include(Infrastructure/AuxiliarySoftware, "Infrastructure", 2, from="= Image Analysis & Segmentation =", to="= Geometry & Mesh Generation & Manipulation =")

Previous Protocols

For more details, see ["Specifications/ExperimentationAnatomicalImaging"].

Protocols

Input

Set(s) of MRI in NIFTI format

Slicer Procedures

Overview

  • Download desired image volume for the desired specimen (NifTI, .nii).
  • Load the desired image volume into Slicer.
  • Segment the tissue of interest using the following detailed segmentation procedures can be segmented by manually painting over its boundaries or using automated or semi-automated algorithms available in the image segmentation & analysis software.

    • NOTES:
    • Save the segmented label map regularly to avoid losing tedious work through program crashes.ation in process should be saved in the relevant repository area specific to the specimen and segmentation results. The file format should be NIFTI and file naming convention should rely on tissue labeling described in ["Specifications/DataManagement"], e.g., oks001-fmb-01.nii, indicating the segmentation for femur bone of the Open Knee(s) - Specimen 1, the appended number indicates the segmentation label (to accommodate multiple segmentations of the same volume by others).
    • Segmented volume can be overlayed on other image sets to evaluate boundaries, however it only be further edited on image volumes with the same resolution as the image from which it was originally defined.
  • After completion of the segmentation, the volume should be exported as an image volume (NIFTI) and a triangulated surface (without smoothing) in STL format, e.g. oks001-fmb-01.stl

Image/Data Input/Output Procedures

Load MRI data:

  • NifTI:
    • File -> Add Data -> [Select MRI file in NifTI format]

    • Open drop-down menu from toolbar for list of modules.
    • Select Editor to edit the labelmap (a volume representation of highlighted regions).
    • One can manually segment bones for every slice, or use Growth Cut Effect and Threshold Effect.
  • DICOM:
    • File -> DICOM, or click DCM icon

    • Click Import
    • Select directory containing DICOM image slices (.IMA)
    • Click 'Add link'
    • CLick 'OK' after images load
    • Select desired Series in DICOM Browser window
    • Click 'Load'

Load segmented label map image data (NifTI):

  • File -> Add Data

  • Click 'Choose File(s) to Add'
  • Select Segmented label map image(s) (NifTI, .nii)
  • Click 'Open'
  • Click/check 'Show Options'
  • Click/check 'LabelMap'

  • NOTE: Link 2D slice views in drop-down tack menu to ensure same label map is displayed in all views.
  • NOTE: In Editor, make sure 'Merge Volume' is set to label map if applying grow-cut or other editor methods.

Save segmented label map (as NifTI):

  • File -> Save, or Save icon

  • Deselect all modified/selected files in the 'File Name' field
  • Select/Check the desired segmented label map volume (NRRD, .nrrd)
  • Set the 'File Format' to NifTI (.nii)
  • Rename the label map volume as desired, see naming convention: ["Specifications/DataManagement"]
  • Click 'Save'

Save triangulated surface (a.k.a. Slicer Model):

  • File -> Save, or Save icon

  • Deselect all modified/selected files in the 'File Name' field
  • Select/Check the desired triangulated surface (VTK, .vtk)
  • Set the 'File Format' to STL (.stl)
  • Rename the triangulated surface as desired, see naming convention: ["Specifications/DataManagement"]
  • Click 'Save'

Viewing Options

In any of the 2D slice viewer windows, using the default mouse cursor pointer tool, click and drag up/down and left/right to change contrast and brightness of MRI

Segmentation Procedures

Manual Segmentation

Insert Katie's instructions here!

Grow-Cut Segmentation

  • From Modules drop-down menu, select Editor
  • If first time entering Editor Module, click Apply to select the default GenericAnatomyColors to create a new Merge Volume

    • NOTE: if not first time, make sure label map displayed (i.e. being edited) in 2D slice viewers is the same as the Merge Volume.
  • Perform the following steps on roughly equally-spaced slices in each image plane orientation as appropriate (coronal, sagittal, axial).
    • Use Paint Effect (under 'Edit Selected Label Map' tab, top row of icons, 3rd icon from left) [with an appropriate/largest radius] to highlight anatomy of interest (foreground) using a specific label (e.g. label1).
    • Use Paint Effect to highlight area around anatomy of interest (background) with a different label (e.g. label2).
    • Use Paint Effect to "cap" boundaries of anatomy of interest (using background label) to prevent foreground leakage (before/after anatomy of interest appears).
  • Select GrowCutEffect icon under 'Edit Selected Label Map' tab, bottom row of icons, 3rd icon from left.

  • Click Apply button. [window graphics will freeze until grow cut algorithm is complete.
  • See figures below for example of segmentation using GrowthCut Effect with foreground (green) and background (yellow).

attachment:gc1.png attachment:gc2.png attachment:gc3.png

  • Use Growth Cut Effect to complete segmentation of bone. (Apply)
  • Remove background by using Change Label Effect; select label2 as input and background [0] as output. (Apply)

attachment:gc4.png

  • Tips and Shortcuts using 3DSlicer:
    • In Editor Mode:
      • 'c' to open list of labels, type # of label id and enter to switch labels quickly
      • 'z' to undo
      • 'y' to redo
      • 'd' for Draw Effect
      • 'p' for Paint Effect
      • 'e' to erase from label
      • arrow keys to move forward/background a slice
      • To change transparency of label in views or switch label to outline view.
        • Hover over thumb tack in top left of view.
        • Select '>>' button.

          • Change percent opacity by changing value near eye icon: 0.0 (transparent) to 1.0 (opaque).
          • Click on button next to eye icon to switch to outline view.

Image Processing Procedures

Label Map Smoothing

  • From the Modules drop-down menu, select 'Surface Models' -> 'Label Map Smoothing'

  • Set Gaussian Smoothing Parameters, Sigma based on image resolution (e.g. 0.7 for cartilage MRI, resolution = 0.35 x 0.35 x 0.7 mm)
  • From 'Input Volume' drop-down menu, select desired label map
  • From 'Output Volume' drop-down menu, select 'Create and rename new Volume'
  • Enter new name as desired, click 'OK'
  • Click Apply

Generate Triangulated Surface from Label Map

  • From the Modules drop-down menu, select 'Surface Models' -> 'Model Maker'

    • Under 'IO' tab, from 'Input Volume' drop-down menu, select desired Label Map
    • In Model drop-down menu, either 'Create and rename new Hierarchy' (e.g. default 'Models'), or select an existing Model Hierarchy.
      • NOTE: all models can be generated under the same or separate model hierarchies
    • Under 'Create Model' tab, specify a desired 'Model Name' for the select Label Map.
    • If no smoothing is desired, under 'Model Maker Parameters set 'Smooth' to zero.
    • Click Apply button.
  • From Editor:
    • Select the MakeModelEffect

    • Set desired Label Map as the Merge Volume
    • If no smoothing is desired, deselect 'smooth model'
    • Click Apply button.

Tissue-Specific Proceduress

Registration Markers

For Registration Marker segmentation, the following technique was used:

  • Use the level tracing effect to fill in the registration markers in the beginning, middle, and end of it's respective appearances within the slices, for each view.
  • Next fill in any gaps with the paint tool.
  • Finally, on those slices, outline the area surrounding the markers with a different color and use the grow cut effect to fill in the segmentation (as described above).

Volume labels include bone type and marker location, e.g. REG-FM, where REG stands for registration, F for femur, M for medial. In result, potential labels are REG-FM, REG-FL, REG-FA, REG-FP, REG-TM, REG-TL, REG-TA, REG-TP, REG-PS, REG-PL, REG-PM.

Bones

INPUT: cartilage MRI

Bone segmenatation procedure:

  1. Perform Grow Cut segmentation procedure on desired bone
  2. Perform Label Map Smoothing to remove Grow Cut boundary noise
  3. Perform manual segmentation to more accurately define bone boundary
  4. Iteratively repeat last two steps until boundary is as desired NOTE: Cortical bone will appear black in MR images, so outer edge of black bone region defines the bone surface.

OUTPUT LABELS: FMB (femur), TBB (tibia), FBB (fibula), PTB (patella).

Cartilage

INPUT: cartilage MRI

Cartilage segmentation should be informed by bone segmentation, i.e. subtraction of bone volume to trim back surface of cartilage.

OUTPUT LABELS: FMC (femoral), TBC-M (medial tibial), TBC-L (lateral tibial), PTC (patellar).

Menisci

Menisci segmentation should be informed by cartilage segmentation, i.e. subtraction of menisci volume to match cartilage boundaries.

Tissue labels include MNS-M (medial), MNS-L (lateral).

Ligaments & Tendons

Ligament segmentation should be informed by bone segmentation, i.e. subtraction of bone volume to trim insertion areas.

Tissue labels include ACL, PCL, LCL, PTL, QAT, etc.

Output

  • Volume of tissue of interest as a binary image aligned with original MRI coordinate system (raw, without filtering and smoothing)
  • Surface representation of tissue of interest in STL format in MRI coordinate system (raw, without filtering and smoothing)

Specifications/ImageSegmentation (last edited 2020-09-04 17:07:52 by snehalkc)