title: Non-Negative Spherical Deconvolution (NNSD) for estimation of fiber Orientation Distribution Function in single-/multi-shell diffusion MRI.
abstract: Spherical Deconvolution (SD) is commonly used for estimating fiber Orientation Distribution Functions (fODFs) from diffusion-weighted signals. Existing SD methods can be classified into two categories: 1) Continuous Representation based SD (CR-SD), where typically Spherical Harmonic (SH) representation is used for convenient analytical solutions, and 2) Discrete Representation based SD (DR-SD), where the signal profile is represented by a discrete set of basis functions uniformly oriented on the unit sphere. A feasible fODF should be non-negative and should integrate to unity throughout the unit sphere S(2). However, to our knowledge, most existing SH-based SD methods enforce non-negativity only on discretized points and not the whole continuum of S(2). Maximum Entropy SD (MESD) and Cartesian Tensor Fiber Orientation Distributions (CT-FOD) are the only SD methods that ensure non-negativity throughout the unit sphere. They are however computational intensive and are susceptible to errors caused by numerical spherical integration. Existing SD methods are also known to overestimate the number of fiber directions, especially in regions with low anisotropy. DR-SD introduces additional error in peak detection owing to the angular discretization of the unit sphere. This paper proposes a SD framework, called Non-Negative SD (NNSD), to overcome all the limitations above. NNSD is significantly less susceptible to the false-positive peaks, uses SH representation for efficient analytical spherical deconvolution, and allows accurate peak detection throughout the whole unit sphere. We further show that NNSD and most existing SD methods can be extended to work on multi-shell data by introducing a three-dimensional fiber response function. We evaluated NNSD in comparison with Constrained SD (CSD), a quadratic programming variant of CSD, MESD, and an L1-norm regularized non-negative least-squares DR-SD. Experiments on synthetic and real single-/multi-shell data indicate that NNSD improves estimation performance in terms of mean difference of angles, peak detection consistency, and anisotropy contrast between isotropic and anisotropic regions.
authors.affiliation: Department of Radiology and BRIC, The University of North Carolina at Chapel Hill, USA. Electronic address: jian_cheng@med.unc.edu.
authors.forename: Jian
authors.initials: J
authors.lastname: Cheng
authors.affiliation: Athena Project-Team, INRIA Sophia Antipolis-Méditerranée, France.
authors.forename: Rachid
authors.initials: R
authors.lastname: Deriche
authors.affiliation: Center for Computational Medicine, LIAMA, Institute of Automation, Chinese Academy of Sciences, China.
authors.forename: Tianzi
authors.initials: T
authors.lastname: Jiang
authors.affiliation: Department of Radiology and BRIC, The University of North Carolina at Chapel Hill, USA.
authors.forename: Dinggang
authors.initials: D
authors.lastname: Shen
authors.affiliation: Department of Radiology and BRIC, The University of North Carolina at Chapel Hill, USA. Electronic address: ptyap@med.unc.edu.
authors.forename: Pew-Thian
authors.initials: PT
authors.lastname: Yap
mesh_terms: D016880:Anisotropy; D001921:Brain / Q000033:anatomy & histology*; D038524:Diffusion Magnetic Resonance Imaging / Q000379:methods*; D006801:Humans; D007091:Image Processing, Computer-Assisted / Q000379:methods*; D008962:Models, Theoretical; D009413:Nerve Fibers, Myelinated
keywords: Diffusion MRI; Fiber Orientation Distribution Function; Non-negativity constraint; Spherical deconvolution; Spherical harmonics
references.citation: Neuroimage. 2004 Nov;23(3):1176-85
references.pmid: 15528117
references.citation: IEEE Trans Med Imaging. 1993;12(2):328-33
references.pmid: 18218422
references.citation: Magn Reson Med. 2005 Nov;54(5):1194-206
references.pmid: 16161109
references.citation: Magn Reson Med. 2002 Oct;48(4):577-82
references.pmid: 12353272
references.citation: Neuroimage. 2012 Sep;62(3):1622-33
references.pmid: 22613620
references.citation: Magn Reson Med. 2007 Sep;58(3):497-510
references.pmid: 17763358
references.citation: Med Image Anal. 2009 Oct;13(5):715-29
references.pmid: 19665917
references.citation: Med Image Comput Comput Assist Interv. 2010;13(Pt 1):582-9
references.pmid: 20879278
references.citation: Magn Reson Med. 2005 Dec;54(6):1377-86
references.pmid: 16247738
references.citation: Neuroimage. 2011 Feb 1;54(3):1812-22
references.pmid: 20970508
references.citation: Neuroimage. 2013 Oct 15;80:125-43
references.pmid: 23702418
references.citation: Med Image Comput Comput Assist Interv. 2014;17(Pt 3):281-8
references.pmid: 25320810
references.citation: Med Image Anal. 2011 Aug;15(4):603-21
references.pmid: 20685153
references.citation: Med Image Comput Comput Assist Interv. 2011;14(Pt 2):66-73
references.pmid: 21995014
references.citation: Med Image Anal. 2011 Aug;15(4):369-96
references.pmid: 21397549
references.citation: Magn Reson Med. 2012 Dec;68(6):1846-55
references.pmid: 22334356
references.citation: Neuroimage. 2014 Nov 1;101:750-64
references.pmid: 25108182
references.citation: Biophys J. 1994 Jan;66(1):259-67
references.pmid: 8130344
references.citation: IEEE Trans Med Imaging. 2014 Aug;33(8):1627-40
references.pmid: 24801775
references.citation: Magn Reson Med. 2010 Aug;64(2):554-66
references.pmid: 20535807
references.citation: Magn Reson Med. 2004 Dec;52(6):1358-72
references.pmid: 15562495
references.citation: Med Image Comput Comput Assist Interv. 2012;15(Pt 2):322-30
references.pmid: 23286064
references.citation: Med Image Comput Comput Assist Interv. 2009;12(Pt 1):911-8
references.pmid: 20426075
references.citation: Neuroimage. 2012 Feb 1;59(3):2045-56
references.pmid: 22019883
references.citation: Med Image Comput Comput Assist Interv. 2009;12(Pt 1):877-85
references.pmid: 20426071
references.citation: IEEE Trans Med Imaging. 2007 Nov;26(11):1464-71
references.pmid: 18041262
references.citation: Inf Process Med Imaging. 2005;19:76-87
references.pmid: 17354686
references.citation: Med Image Anal. 2012 Aug;16(6):1121-9
references.pmid: 22890050
references.citation: IEEE Trans Med Imaging. 2012 Nov;31(11):2035-49
references.pmid: 22711770
references.citation: IEEE Trans Med Imaging. 2011 Oct;30(10):1829-40
references.pmid: 21571607
references.citation: Magn Reson Med. 2012 Mar;67(3):844-55
references.pmid: 22183751
references.citation: Med Image Comput Comput Assist Interv. 2010;13(Pt 1):648-56
references.pmid: 20879286
references.citation: Med Image Comput Comput Assist Interv. 2008;11(Pt 2):122-30
references.pmid: 18982597
references.citation: PLoS One. 2011;6(9):e24678
references.pmid: 21966364
references.citation: Med Image Comput Comput Assist Interv. 2012;15(Pt 2):313-21
references.pmid: 23286063
references.citation: Neuroimage. 2012 Feb 1;59(3):2175-86
references.pmid: 22019877
references.citation: Neuroimage. 2007 May 1;35(4):1459-72
references.pmid: 17379540
references.citation: Neuroimage. 2010 Jan 15;49(2):1446-58
references.pmid: 19781650
references.citation: Hum Brain Mapp. 2013 Sep;34(9):2089-102
references.pmid: 22461221
references.citation: Magn Reson Med. 2006 Jul;56(1):104-17
references.pmid: 16755539
references.citation: Med Image Comput Comput Assist Interv. 2011;14(Pt 2):98-106
references.pmid: 21995018
references.citation: IEEE Trans Biomed Eng. 2007 Mar;54(3):462-72
references.pmid: 17355058
references.citation: Med Image Comput Comput Assist Interv. 2010;13(Pt 1):590-7
references.pmid: 20879279