Difference between revisions of "Initial atlas construction workflow"

From NAMIC Wiki
Jump to: navigation, search
m (moved Current atlas construction workflow to Initial atlas construction workflow: The atlas construction workflow has been revised. New workflow has been adopted.)
 
(3 intermediate revisions by 2 users not shown)
Line 32: Line 32:
 
1.3. Slicer Histogram Matching module
 
1.3. Slicer Histogram Matching module
  
1.4. FSL/Flirt
+
1.4. FSL-Flirt
  
 
1.5. Slicer Mask, Editor modules (verify the mask encloses actual ICC!)
 
1.5. Slicer Mask, Editor modules (verify the mask encloses actual ICC!)
Line 40: Line 40:
 
|  
 
|  
  
2.1. Use the skullstriped image from Step 1.2 as the template to repeat Steps 1.3-1.4.
+
2.1. Use the skullstriped image from Step 1.2 as the template to repeat Step 1.4.
  
 
2.2. Perform non-rigid alignment of the subjects affinely registered to the template.
 
2.2. Perform non-rigid alignment of the subjects affinely registered to the template.
Line 51: Line 51:
 
2.2. Diffeomorphic demons
 
2.2. Diffeomorphic demons
  
2.3. Some tool in Slicer under develop
+
2.3. Some tool in Slicer under development
  
 
|-
 
|-
Line 60: Line 60:
  
 
3.2. Slightly smooth the segmentations Styner et al. recommend kernel of 0.4 mm variance.
 
3.2. Slightly smooth the segmentations Styner et al. recommend kernel of 0.4 mm variance.
 +
 +
3.3. Probabilistic atlas is the result of smoothing of the hard segmentation of the averaged template with Gaussian kernel.
  
 
|  
 
|  
Line 65: Line 67:
 
3.1. Kmeans + manual
 
3.1. Kmeans + manual
  
3.2. Gaussian smoothing module in Slicer, 0.6 and 0.65 kernel size used for Valentino and Tommy, respectively.
+
3.2. Gaussian smoothing module in Slicer, 0.6 and 0.65 kernel size used for Valentino and Tommy, respectively. '''AF: I do not understand what is the relation of Valentino and Tommy to this, if the operation was performed on hard segmentation from the average template'''
 +
 
 +
3.3. Slicer Gaussian smoothing module
  
*We did not perform back propagation due to the large structural differences in the dataset.
+
*We did not perform back propagation as in Styner07 due to the large structural differences in the dataset.
  
 
|}
 
|}
 
== Standing questions ==
 
 
* when to perform bias correction? After Step 1.4?
 
* deep brain structures -- similar approach for identification in subjects as for GM/WM/CSF in Step 3.1?
 
* Styner et al. used T2 to segment CSF. We only have T1
 
  
 
== References ==
 
== References ==

Latest revision as of 17:51, 3 December 2009

Home < Initial atlas construction workflow

Content

This page describes the current working pipeline of atlas construction for vervet data.

Data

10 vervet subjects (2 subjects possibly have unacceptable image quality), T1 sequence

Workflow / Tools

Workflow step Description Tools
1. Preparation

1.1. Select the subject that is best oriented as the template image

1.2. Manually identify the ICC in the template subject.

1.3. Perform intensity calibration for all other subjects to the template.

1.4. 12 DOF registration of each of the subjects to the template.

1.5. Skull-strip each subject using dilated ICC of the template.

1.1. n/a

1.2. manual in SNAP.

1.3. Slicer Histogram Matching module

1.4. FSL-Flirt

1.5. Slicer Mask, Editor modules (verify the mask encloses actual ICC!)

2. Atlas construction

2.1. Use the skullstriped image from Step 1.2 as the template to repeat Step 1.4.

2.2. Perform non-rigid alignment of the subjects affinely registered to the template.

2.3. Compute the atlas as the average.

2.1 see above

2.2. Diffeomorphic demons

2.3. Some tool in Slicer under development

3. Probabilistic atlas construction

3.1. Segment WM/GM/CSF from the averaged template, manually edit to ensure accuracy

3.2. Slightly smooth the segmentations Styner et al. recommend kernel of 0.4 mm variance.

3.3. Probabilistic atlas is the result of smoothing of the hard segmentation of the averaged template with Gaussian kernel.

3.1. Kmeans + manual

3.2. Gaussian smoothing module in Slicer, 0.6 and 0.65 kernel size used for Valentino and Tommy, respectively. AF: I do not understand what is the relation of Valentino and Tommy to this, if the operation was performed on hard segmentation from the average template

3.3. Slicer Gaussian smoothing module

  • We did not perform back propagation as in Styner07 due to the large structural differences in the dataset.

References

  1. M. Styner, R. Knickmeyer, S. Joshi, C. Coe, S. J. Short, and J. Gilmore. Automatic brain segmentation in rhesus monkeys. Proc SPIE Medical Imaging Conference, Proc SPIE Vol 6512 Medical Imaging 2007, pp 65122L-1 - 65122L-8 pdf
  2. Balci S.K., Golland P., Wells W.M. Non-rigid Groupwise Registration using B-Spline Deformation Model. Insight Journal - 2007 MICCAI Open Science Workshop. link
  3. Previous descriptions of the atlas construction workflow: Summary by Ginger Li, Description off BSL atlas page