THOMAS: Thalamus Optimized Multi-Atlas Segmentation at 3T
Jason Su1, Thomas Tourdias2, Manojkumar Saranathan3, Pejman Ghanouni4, and Brian Rutt4

1Electrical Engineering, Stanford University, Stanford, CA, United States, 2Neuroradiology, Bordeaux University Hospital, Bordeaux, France, 3Radiology, University of Arizona, Tucson, AZ, United States, 4Radiology, Stanford University, Stanford, CA, United States

Synopsis

The efficacy of the Thalamus Optimized Multi-Atlas Segmentation (THOMAS) algorithm for segmentation of thalamic nuclei with white-matter-nulled MP-RAGE images is studied in 3T and 7T variants of the image contrast. 5 subjects are evaluated at both field strengths and ground truth manual delineations of nuclei are performed on the 7T images. We demonstrate that the algorithm performs as well on 3T images as on 7T within a dice coefficient of ±0.1 as evaluated against the ground truth. This indicates that THOMAS can now reach a much wider audience of interested groups.

Target Audience

Scientists interested in automatic segmentation of the thalamus and its nuclei, esp. for the study of atrophy or treatment targeting.

Purpose

The white matter nulled MP-RAGE (WMnMPRAGE) contrast at 7T has enabled detailed manual delineation of thalamic nuclei1 by a trained neuroradiologist guided by the Morel atlas2 as well as the creation of a fully automatic segmentation pipeline called THOMAS (Thalamus Optimized Multi-Atlas Segmentation), the first of its kind to segment the 12 thalamic nuclei at 1mm3 resolution with validation against expert-defined labels.3 However, the availability of 7T scanners is limited. In this work we investigate the application of THOMAS to WMnMPRAGE images acquired at 3T. The performance of the algorithm is evaluated over 5 subjects scanned at 3T and 7T against ground truth manual delineations from 7T images.

Methods

Following informed consent, 6 essential tremor patients were scanned at 7T (Discovery MR950, GE Healthcare) using a 32-channel head coil (Nova Medical) with WMnMPRAGE scan parameters: TS 6s, TI 680ms, TR 10ms, BW 12kHz, flip angle 4°, FOV 18cm, 180x180x200 matrix, slice thickness 1mm, ARC parallel imaging 1.5x1.5 (2D radial fanbeam), scan time 5.5 minutes.4 Using these images, the whole thalami and 12 thalamic nuclei were manually delineated with a high degree of reproducibility to form the ground truth.1 WMnMPRAGE images were also acquired at 3T (Discovery MR750, GE Healthcare) using an 8-channel head coil: TS 4.5s, TI 500ms, TR 10ms, BW 9kHz, flip angle 7°, FOV 18cm, 180x180x200 matrix, slice thickness 1mm, no parallel imaging, scan time ~10 min with prospective motion correction. THOMAS was applied without any modification to all WMn images to automatically segment the thalamus and its nuclei. The pipeline involves non-linearly registering a library of prior thalamic atlases to the target subject and treating these as potential candidate segmentations. These candidates are combined together and their differences resolved with a label fusion by Wang and Yushkevich.5 The 7T images were registered to the 3T images using FSL’s FLIRT with 7 degrees-of-freedom. Special care was taken to ignore signal artifacts around the sinuses which were present at both field strengths but more severe at 7T in the registration. The results were evaluated against manually delineated truth labels using Dice’s coefficient. In the 3T case, the ground truth labels were transferred from 7T using the registration. One subject was removed from the analysis due to poor image quality from coil sensitivity fall off.

Results

Segmentation time for a single subject was 30-35 minutes on a 12-core 2.66 GHz Intel Xeon independent of field strength. Fig. 1 shows the performance of THOMAS for the whole thalamus and each of the 12 nuclei at 3T and 7T. Fig. 2 shows different in the dice between 3T and 7T. Notably, it appears that THOMAS works as well with 3T WMnMPRAGE as with 7T images. The difference in dice is generally within ±0.1 between the two on an individual subject basis (Fig. 2). Fig. 3 aggregates the Dice across all 12 nuclei and computes a kernel density estimate of the probability of segmentation performance at each field strength. A Wilcoxon signed-rank test has a p-value of 0.531, which indicates that the null hypothesis that the differences between the two are symmetric about 0 cannot be rejected. Example segmentations are depicted in Fig. 4.

Discussion/Conclusion

The THOMAS algorithm can successfully segment thalamic nuclei in both 3T and 7T WMnMPRAGE images with essentially equivalent quality, when evaluated against manual delineations. In one subject, THOMAS performed exceptionally well on the 3T images compared to 7T, we speculate that this may be due to improved B1+ uniformity at 3T. The initial step in the pipeline does bias-correct the images to remove any gross image non-uniformities but this may not be sufficient for 7T images and does not resolve non-uniformities in image contrast caused by B1+ non-uniformity. For these reasons, 3T WMnMPRAGE may eventually prove to be preferred over 7T WMnMPRAGE despite the improved signal-to-noise ratio at 7T. Reliable automatic segmentation of the whole thalamus and its nuclei is critically required for any large-scale study of thalamic changes with disease. It may also play an important role for treatment planning, e.g. in thalamic ablation for essential tremor or deep-brain stimulation for a variety of disorders. Demonstration that our method successfully operates at 3T means that it can find applicability in many more studies. Ongoing work will evaluate segmentation performance using a variety of 3T WMnMPRAGE volumes collected over a several year time period across different sites and scanner manufacturers.

Acknowledgements

Research support from NIH P41 EB015891, NIH 1 S10 RR026351-01Al, GE Healthcare, BioX-Sanofi BIOSTAR program.

References

[1]Tourdias et al. Neuroimage. 2013 Sep 7;84C:534-545 [2]Niemann et al. Neuroimage. 2000 Dec;12(6):601-16. [3]Su et al. Proceedings of ISMRM. 2015 May 30;6231. [4]Saranathan et al. Magn Reson Med. 2014 May 29. [5]Wang and Yushkevich. Front Neuroinform. 2013 Nov 22;7:27.

Figures

A box plot showing the performance of THOMAS with 3T (blue) and 7T (red) WMnMPRAGE images. The white line indicates the median of the group of 5 subjects, the box indicates the interquartile range, and the whiskers the total range.

The difference in dice between the automatic segmentation computed on a per-subject basis and displayed as a box plot. The red line indicates the median of the group of 5 subjects, the box indicates the interquartile range, and the whiskers up to the total range.

An example segmentation from 3T (left) and 7T (right) with the predicted ROIs as solid colors and the ground truth in yellow outline.



Proc. Intl. Soc. Mag. Reson. Med. 24 (2016)
4328