DBP3:Utah:RegSegPipeline

From NAMIC Wiki
Jump to: navigation, search
Home < DBP3:Utah:RegSegPipeline
  back to DBP3 home

The CARMA DBP: MRI-based study and treatment of atrial fibrillation

Alex Zaitsev, Dominik Meier, Ron Kikinis

Pilot Studies on a Registration & Segmentation Pipeline & Workflow

Overall processing steps are (order tentative)

  1. N4 bias field correction for the MRI (surface coils):
    1. run on entire image gives some benefit that may be improved with masking: again the dominant intensity dropoff from the surface coil occurs along the chest wall and ribcage. Even if that is not the structure of interest, it is the low-freq. variation the bias correction algorithm is searching for, and masking that out can be counter-productive: via masking we may end up with a smoother image, but the intensity variations removed were not caused by the coil but are actually true signal.
    2. Module used: N4 ITK; Parameters: convergence: 1e-5, iterations: 50,40,30,20, shrink factor: 3
  2. registration MRA>cMRI
    1. the MRA contains the same FOV and has surrounding structures (liver, chest, spine etc) visible also, despite lower intensities. A global affine is thus not necessarily going to benefit from masking the heart, unless the relative motion of the heart becomes the dominant reason for misalignment.
    2. Module used: BRAINSfit
      1. tried masking with both BrainsFit and RobustMultires modules. Both failed to provide better alignment with masking.
      2. images have good alignment to begin with, need reliable fiducial to evaluate registration. E..g. descending aorta? non-cardiac structures?
  3. registration follow-up -> baseline
    1. most reliably done on the post contrast MRI.
    2. DOF at least 12, because image is captured at different phases in the breathing/cardiac cycle, some low-res BSpline should be beneficial
    3. Module used: BRAINSfit
    4. the IS FOV can differ, e.g. how much of the liver is included. If the two exams differ significantly (>30%) in that content, a prior cropping step is suggested to better match image content before registration. Resampling to isotropic voxel size at this stage is also advantageous.
  4. ROI definition (manual box ROI or automated via atlas)
  5. segmentation of LA from MRA -> inner wall
    1. as a dynamic image the MRA contains significant spread and likely requires interactive segmentation/thresholding to yield a satisfactory LA volume
    2. Module used: Editor: thresholding or thresholding within Volumes thresholding option within Display tab, use iron colormap & low alpha setting to check for ventricular wall borders.
    1. cropping and island removal
  1. LA wall segmentation
    1. very small structure, most reliably done manually direct. Starting with automation may yield more effort on post-edits
    2. Module used: Editor: manual outline
  2. segmentation of enhancement within LA wall: intensity statistics. An atlas-based set of intensity distributions may be more meaningful here than a simple Otsu, because both amount and location of enhancement is unknown and can in theory be 0.