1
|
Cheng J, Shen D, Yap PT, Basser PJ. Single- and Multiple-Shell Uniform Sampling Schemes for Diffusion MRI Using Spherical Codes. IEEE TRANSACTIONS ON MEDICAL IMAGING 2018; 37:185-199. [PMID: 28952937 PMCID: PMC5867228 DOI: 10.1109/tmi.2017.2756072] [Citation(s) in RCA: 9] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 06/07/2023]
Abstract
In diffusion MRI (dMRI), a good sampling scheme is important for efficient acquisition and robust reconstruction. Diffusion weighted signal is normally acquired on single or multiple shells in q-space. Signal samples are typically distributed uniformly on different shells to make them invariant to the orientation of structures within tissue, or the laboratory coordinate frame. The Electrostatic Energy Minimization (EEM) method, originally proposed for single shell sampling scheme in dMRI, was recently generalized to multi-shell schemes, called Generalized EEM (GEEM). GEEM has been successfully used in the Human Connectome Project (HCP). However, EEM does not directly address the goal of optimal sampling, i.e., achieving large angular separation between sampling points. In this paper, we propose a more natural formulation, called Spherical Code (SC), to directly maximize the minimal angle between different samples in single or multiple shells. We consider not only continuous problems to design single or multiple shell sampling schemes, but also discrete problems to uniformly extract sub-sampled schemes from an existing single or multiple shell scheme, and to order samples in an existing scheme. We propose five algorithms to solve the above problems, including an incremental SC (ISC), a sophisticated greedy algorithm called Iterative Maximum Overlap Construction (IMOC), an 1-Opt greedy method, a Mixed Integer Linear Programming (MILP) method, and a Constrained Non-Linear Optimization (CNLO) method. To our knowledge, this is the first work to use the SC formulation for single or multiple shell sampling schemes in dMRI. Experimental results indicate that SC methods obtain larger angular separation and better rotational invariance than the state-of-the-art EEM and GEEM. The related codes and a tutorial have been released in DMRITool.
Collapse
|
2
|
Koay CG, Yeh PH, Ollinger JM, İrfanoğlu MO, Pierpaoli C, Basser PJ, Oakes TR, Riedy G. Tract Orientation and Angular Dispersion Deviation Indicator (TOADDI): A framework for single-subject analysis in diffusion tensor imaging. Neuroimage 2015; 126:151-63. [PMID: 26638985 DOI: 10.1016/j.neuroimage.2015.11.046] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/22/2015] [Revised: 11/05/2015] [Accepted: 11/18/2015] [Indexed: 11/19/2022] Open
Abstract
The purpose of this work is to develop a framework for single-subject analysis of diffusion tensor imaging (DTI) data. This framework is termed Tract Orientation and Angular Dispersion Deviation Indicator (TOADDI) because it is capable of testing whether an individual tract as represented by the major eigenvector of the diffusion tensor and its corresponding angular dispersion are significantly different from a group of tracts on a voxel-by-voxel basis. This work develops two complementary statistical tests based on the elliptical cone of uncertainty, which is a model of uncertainty or dispersion of the major eigenvector of the diffusion tensor. The orientation deviation test examines whether the major eigenvector from a single subject is within the average elliptical cone of uncertainty formed by a collection of elliptical cones of uncertainty. The shape deviation test is based on the two-tailed Wilcoxon-Mann-Whitney two-sample test between the normalized shape measures (area and circumference) of the elliptical cones of uncertainty of the single subject against a group of controls. The False Discovery Rate (FDR) and False Non-discovery Rate (FNR) were incorporated in the orientation deviation test. The shape deviation test uses FDR only. TOADDI was found to be numerically accurate and statistically effective. Clinical data from two Traumatic Brain Injury (TBI) patients and one non-TBI subject were tested against the data obtained from a group of 45 non-TBI controls to illustrate the application of the proposed framework in single-subject analysis. The frontal portion of the superior longitudinal fasciculus seemed to be implicated in both tests (orientation and shape) as significantly different from that of the control group. The TBI patients and the single non-TBI subject were well separated under the shape deviation test at the chosen FDR level of 0.0005. TOADDI is a simple but novel geometrically based statistical framework for analyzing DTI data. TOADDI may be found useful in single-subject, graph-theoretic and group analyses of DTI data or DTI-based tractography techniques.
Collapse
Affiliation(s)
- Cheng Guan Koay
- National Intrepid Center of Excellence (NICoE), Bethesda, MD, USA; Section on Tissue Biophysics and Biomimetics, NICHD, National Institutes of Health, Bethesda, MD, USA; NorthTide Group, LLC, USA.
| | - Ping-Hong Yeh
- National Intrepid Center of Excellence (NICoE), Bethesda, MD, USA; The Henry M. Jackson Foundation for the Advancement of Military Medicine, Bethesda, MD, USA
| | - John M Ollinger
- National Intrepid Center of Excellence (NICoE), Bethesda, MD, USA
| | - M Okan İrfanoğlu
- The Henry M. Jackson Foundation for the Advancement of Military Medicine, Bethesda, MD, USA; Section on Tissue Biophysics and Biomimetics, NICHD, National Institutes of Health, Bethesda, MD, USA
| | - Carlo Pierpaoli
- Section on Tissue Biophysics and Biomimetics, NICHD, National Institutes of Health, Bethesda, MD, USA
| | - Peter J Basser
- Section on Tissue Biophysics and Biomimetics, NICHD, National Institutes of Health, Bethesda, MD, USA
| | - Terrence R Oakes
- National Intrepid Center of Excellence (NICoE), Bethesda, MD, USA
| | - Gerard Riedy
- National Intrepid Center of Excellence (NICoE), Bethesda, MD, USA; National Capital Neuroimaging Consortium, Bethesda, MD, USA
| |
Collapse
|
3
|
Pasha Hosseinbor A, Chung MK, Koay CG, Schaefer SM, van Reekum CM, Schmitz LP, Sutterer M, Alexander AL, Davidson RJ. 4D hyperspherical harmonic (HyperSPHARM) representation of surface anatomy: a holistic treatment of multiple disconnected anatomical structures. Med Image Anal 2015; 22:89-101. [PMID: 25828650 DOI: 10.1016/j.media.2015.02.004] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/11/2014] [Revised: 01/19/2015] [Accepted: 02/20/2015] [Indexed: 11/19/2022]
Abstract
Image-based parcellation of the brain often leads to multiple disconnected anatomical structures, which pose significant challenges for analyses of morphological shapes. Existing shape models, such as the widely used spherical harmonic (SPHARM) representation, assume topological invariance, so are unable to simultaneously parameterize multiple disjoint structures. In such a situation, SPHARM has to be applied separately to each individual structure. We present a novel surface parameterization technique using 4D hyperspherical harmonics in representing multiple disjoint objects as a single analytic function, terming it HyperSPHARM. The underlying idea behind HyperSPHARM is to stereographically project an entire collection of disjoint 3D objects onto the 4D hypersphere and subsequently simultaneously parameterize them with the 4D hyperspherical harmonics. Hence, HyperSPHARM allows for a holistic treatment of multiple disjoint objects, unlike SPHARM. In an imaging dataset of healthy adult human brains, we apply HyperSPHARM to the hippocampi and amygdalae. The HyperSPHARM representations are employed as a data smoothing technique, while the HyperSPHARM coefficients are utilized in a support vector machine setting for object classification. HyperSPHARM yields nearly identical results as SPHARM, as will be shown in the paper. Its key advantage over SPHARM lies computationally; HyperSPHARM possess greater computational efficiency than SPHARM because it can parameterize multiple disjoint structures using much fewer basis functions and stereographic projection obviates SPHARM's burdensome surface flattening. In addition, HyperSPHARM can handle any type of topology, unlike SPHARM, whose analysis is confined to topologically invariant structures.
Collapse
Affiliation(s)
- A Pasha Hosseinbor
- Waisman Laboratory for Brain Imaging and Behavior, University of Wisconsin-Madison, Madison, WI, USA.
| | - Moo K Chung
- Waisman Laboratory for Brain Imaging and Behavior, University of Wisconsin-Madison, Madison, WI, USA; Department of Biostatistics and Medical Informatics, University of Wisconsin-Madison, Madison, WI, USA
| | - Cheng Guan Koay
- Department of Medical Physics, University of Wisconsin-Madison, Madison, WI, USA
| | - Stacey M Schaefer
- Waisman Laboratory for Brain Imaging and Behavior, University of Wisconsin-Madison, Madison, WI, USA
| | - Carien M van Reekum
- Department of Psychology, University of Reading, Reading, Berkshire RG6 6UR, UK
| | - Lara Peschke Schmitz
- Waisman Laboratory for Brain Imaging and Behavior, University of Wisconsin-Madison, Madison, WI, USA
| | - Matt Sutterer
- Waisman Laboratory for Brain Imaging and Behavior, University of Wisconsin-Madison, Madison, WI, USA
| | - Andrew L Alexander
- Waisman Laboratory for Brain Imaging and Behavior, University of Wisconsin-Madison, Madison, WI, USA; Department of Medical Physics, University of Wisconsin-Madison, Madison, WI, USA; Department of Psychiatry, University of Wisconsin-Madison, Madison, WI, USA
| | - Richard J Davidson
- Waisman Laboratory for Brain Imaging and Behavior, University of Wisconsin-Madison, Madison, WI, USA; Department of Psychology, University of Wisconsin-Madison, Madison, WI, USA
| |
Collapse
|