1
|
Consagra W, Ning L, Rathi Y. Neural orientation distribution fields for estimation and uncertainty quantification in diffusion MRI. Med Image Anal 2024; 93:103105. [PMID: 38377728 DOI: 10.1016/j.media.2024.103105] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/17/2023] [Revised: 12/13/2023] [Accepted: 02/05/2024] [Indexed: 02/22/2024]
Abstract
Inferring brain connectivity and structure in-vivo requires accurate estimation of the orientation distribution function (ODF), which encodes key local tissue properties. However, estimating the ODF from diffusion MRI (dMRI) signals is a challenging inverse problem due to obstacles such as significant noise, high-dimensional parameter spaces, and sparse angular measurements. In this paper, we address these challenges by proposing a novel deep-learning based methodology for continuous estimation and uncertainty quantification of the spatially varying ODF field. We use a neural field (NF) to parameterize a random series representation of the latent ODFs, implicitly modeling the often ignored but valuable spatial correlation structures in the data, and thereby improving efficiency in sparse and noisy regimes. An analytic approximation to the posterior predictive distribution is derived which can be used to quantify the uncertainty in the ODF estimate at any spatial location, avoiding the need for expensive resampling-based approaches that are typically employed for this purpose. We present empirical evaluations on both synthetic and real in-vivo diffusion data, demonstrating the advantages of our method over existing approaches.
Collapse
Affiliation(s)
- William Consagra
- Psychiatry Neuroimaging Laboratory, Brigham and Women's Hospital, Harvard Medical School, 399 Revolution Drive, Boston, 02215, MA, United States.
| | - Lipeng Ning
- Psychiatry Neuroimaging Laboratory, Brigham and Women's Hospital, Harvard Medical School, 399 Revolution Drive, Boston, 02215, MA, United States
| | - Yogesh Rathi
- Psychiatry Neuroimaging Laboratory, Brigham and Women's Hospital, Harvard Medical School, 399 Revolution Drive, Boston, 02215, MA, United States
| |
Collapse
|
2
|
Guo L, Lyu J, Zhang Z, Shi J, Feng Q, Feng Y, Gao M, Zhang X. A Joint Framework for Denoising and Estimating Diffusion Kurtosis Tensors Using Multiple Prior Information. IEEE TRANSACTIONS ON MEDICAL IMAGING 2022; 41:308-319. [PMID: 34520348 DOI: 10.1109/tmi.2021.3112515] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 06/13/2023]
Abstract
Diffusion kurtosis imaging (DKI) has been shown to be valuable in a wide range of neuroscientific and clinical applications. However, reliable estimation of DKI tensors is often compromised by noise, especially for the kurtosis tensor (KT). Here, we propose a joint denoising and estimating framework that integrates multiple sources of prior information, including nonlocal structural self-similarity (NSS), local spatial smoothness (LSS), physical relevance (PR) of the DKI model, and noise characteristics of magnitude diffusion MRI (dMRI) images for improved estimation of DKI tensors. The local and nonlocal spatial smoothing constraints are complementary to each other, making the proposed framework highly effective in reducing the noise fluctuations on DKI tensors, especially KT. As an additional refinement, we propose to impose a physically relevant constraint within our joint denoising and estimation framework. We further adopt the first-moment noise-corrected fitting model (M1NCM) to remove the noncentral χ -distribution noise bias. The effectiveness of integrating multiple sources of priors into the joint framework is verified by comparing the proposed M1NCM-NSS-LSS-PR method with various versions of M1NCM-based estimators and two state-of-the-art methods. Results show that the proposed method outperformed the compared methods in simulations and in-vivo dMRI datasets of both spatially stationary and nonstationary noise distributions. The in-vivo experiments also show that the proposed M1NCM-NSS-LSS-PR method was robust to the number of diffusion directions.
Collapse
|
3
|
Chu CY, Sun CY, Kuai ZX, Yang F, Zhu YM. Structure Prior Constrained Estimation of Human Cardiac Diffusion Tensors. IEEE Trans Biomed Eng 2019; 66:3220-3230. [PMID: 30843792 DOI: 10.1109/tbme.2019.2902381] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/10/2022]
Abstract
OBJECTIVE The purpose of this paper is to increase the accuracy of human cardiac diffusion tensor (DT) estimation in diffusion magnetic resonance imaging (dMRI) with a few diffusion gradient directions. METHODS A structure prior constrained (SPC) method is proposed. The method consists in introducing two regularizers in the conventional nonlinear least squares estimator. The two regularizers penalize the dissimilarity between neighboring DTs and the difference between estimated and prior fiber orientations, respectively. A novel numerical solution is presented to ensure the positive definite estimation. RESULTS Experiments on ex vivo human cardiac data show that the SPC method is able to well estimate DTs at most voxels, and is superior to state-of-the-art methods in terms of the mean errors of principal eigenvector, second eigenvector, helix angle, transverse angle, fractional anisotropy, and mean diffusivity. CONCLUSION The SPC method is a practical and reliable alternative to current denoising- or regularization-based methods for the estimation of human cardiac DT. SIGNIFICANCE The SPC method is able to accurately estimate human cardiac DTs in dMRI with a few diffusion gradient directions.
Collapse
|
4
|
|
5
|
Snow WM, Dale R, O'Brien-Moran Z, Buist R, Peirson D, Martin M, Albensi BC. In Vivo Detection of Gray Matter Neuropathology in the 3xTg Mouse Model of Alzheimer's Disease with Diffusion Tensor Imaging. J Alzheimers Dis 2018; 58:841-853. [PMID: 28505976 DOI: 10.3233/jad-170136] [Citation(s) in RCA: 18] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/15/2023]
Abstract
A diagnosis of Alzheimer's disease (AD), a neurodegenerative disorder accompanied by severe functional and cognitive decline, is based on clinical findings, with final confirmation of the disease at autopsy by the presence of amyloid-β (Aβ) plaques and neurofibrillary tangles. Given that microstructural brain alterations occur years prior to clinical symptoms, efforts to detect brain changes early could significantly enhance our ability to diagnose AD sooner. Diffusion tensor imaging (DTI), a type of MRI that characterizes the magnitude, orientation, and anisotropy of the diffusion of water in tissues, has been used to infer neuropathological changes in vivo. Its utility in AD, however, is still under investigation. The current study used DTI to examine brain regions susceptible to AD-related pathology; the cerebral cortex, entorhinal cortex, and hippocampus, in 12-14-month-old 3xTg AD mice that possess both Aβ plaques and neurofibrillary tangles. Mean diffusivity did not differ between 3xTg and control mice in any region. Decreased fractional anisotropy (p < 0.01) and axial diffusivity (p < 0.05) were detected only in the hippocampus, in which both congophilic Aβ plaques and hyperphosphorylated tau accumulation, consistent with neurofibrillary tangle formation, were detected. Pathological tau accumulation was seen in the cortex. The entorhinal cortex was largely spared from AD-related neuropathology. This is the first study to demonstrate DTI abnormalities in gray matter in a mouse model of AD in which both pathological hallmarks are present, suggesting the feasibility of DTI as a non-invasive means of detecting brain pathology in vivo in early-stage AD.
Collapse
Affiliation(s)
- Wanda M Snow
- Division of Neurodegenerative Disorders, St. Boniface Hospital Albrechtsen Research Centre, Winnipeg, MB, Canada.,Department of Pharmacology & Therapeutics, University of Manitoba, Winnipeg, MB, Canada
| | - Ryan Dale
- Division of Neurodegenerative Disorders, St. Boniface Hospital Albrechtsen Research Centre, Winnipeg, MB, Canada.,Department of Pharmacology & Therapeutics, University of Manitoba, Winnipeg, MB, Canada
| | | | - Richard Buist
- Department of Radiology, University of Manitoba, Winnipeg, MB, Canada
| | - Danial Peirson
- Division of Neurodegenerative Disorders, St. Boniface Hospital Albrechtsen Research Centre, Winnipeg, MB, Canada.,Department of Pharmacology & Therapeutics, University of Manitoba, Winnipeg, MB, Canada
| | - Melanie Martin
- Department of Pharmacology & Therapeutics, University of Manitoba, Winnipeg, MB, Canada.,Department of Physics, University of Winnipeg, Winnipeg, MB, Canada.,Department of Radiology, University of Manitoba, Winnipeg, MB, Canada
| | - Benedict C Albensi
- Division of Neurodegenerative Disorders, St. Boniface Hospital Albrechtsen Research Centre, Winnipeg, MB, Canada.,Department of Pharmacology & Therapeutics, University of Manitoba, Winnipeg, MB, Canada
| |
Collapse
|
6
|
Recent Computational Advances in Denoising for Magnetic Resonance Diffusional Kurtosis Imaging (DKI). J Indian Inst Sci 2017. [DOI: 10.1007/s41745-017-0036-2] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
|
7
|
Matsuda KM, Lopes-Calcas A, Honke ML, O'Brien-Moran Z, Buist R, West M, Martin M. Ex vivo tissue imaging for radiology-pathology correlation: a pilot study with a small bore 7-T MRI in a rare pigmented ganglioglioma exhibiting complex MR signal characteristics associated with melanin and hemosiderin. J Med Imaging (Bellingham) 2017; 4:036001. [PMID: 28924575 PMCID: PMC5596201 DOI: 10.1117/1.jmi.4.3.036001] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/31/2017] [Accepted: 08/09/2017] [Indexed: 12/02/2022] Open
Abstract
To advance magnetic resonance imaging (MRI) technologies further for in vivo tissue characterization with histopathologic validation, we investigated the feasibility of ex vivo tissue imaging of a surgically removed human brain tumor as a comprehensive approach for radiology–pathology correlation in histoanatomically identical fashion in a rare case of pigmented ganglioglioma with complex paramagnetic properties. Pieces of surgically removed ganglioglioma, containing melanin and hemosiderin pigments, were imaged with a small bore 7-T MRI scanner to obtain T1-, T2-, and T2*-weighted image and diffusion tensor imaging (DTI). Corresponding histopathological slides were prepared for routine hematoxylin and eosin stain and special stains for melanin and iron/hemosiderin to correlate with MRI signal characteristics. Furthermore, mean diffusivity (MD) maps were generated from DTI data and correlated with cellularity using image analysis. While the presence of melanin was difficult to interpret in in vivo MRI with certainty due to concomitant hemosiderin pigments and calcium depositions, ex vivo tissue imaging clearly demonstrated pieces of tissue exhibiting the characteristic MR signal pattern for melanin with pathologic confirmation in a histoanatomically identical location. There was also concordant correlation between MD and cellularity. Although it is still in an initial phase of development, ex vivo tissue imaging is a promising approach, which offers radiology–pathology correlation in a straightforward and comprehensive manner.
Collapse
Affiliation(s)
- Kant M Matsuda
- University of Manitoba, Max Rady College of Medicine, Department of Pathology, Rady Faculty of Health Sciences, Winnipeg, Manitoba, Canada.,Health Sciences Centre Winnipeg, Department of Pathology, Diagnostic Services of Manitoba, Winnipeg, Manitoba, Canada.,Memorial Sloan-Kettering Cancer Center, Department of Pathology, New York, New York, United States
| | - Ana Lopes-Calcas
- University of Manitoba, Max Rady College of Medicine, Department of Pathology, Rady Faculty of Health Sciences, Winnipeg, Manitoba, Canada
| | - Michael L Honke
- University of Winnipeg, Department of Physics, Winnipeg, Manitoba, Canada
| | - Zoe O'Brien-Moran
- University of Winnipeg, Department of Physics, Winnipeg, Manitoba, Canada
| | - Richard Buist
- University of Manitoba, Max Rady College of Medicine, Department of Radiology, Rady Faculty of Health Sciences, Winnipeg, Manitoba, Canada
| | - Michael West
- University of Manitoba, Max Rady College of Medicine, Department of Neurosurgery, Rady Faculty of Health Sciences, Winnipeg, Manitoba, Canada
| | - Melanie Martin
- University of Winnipeg, Department of Physics, Winnipeg, Manitoba, Canada.,University of Manitoba, Max Rady College of Medicine, Department of Radiology, Rady Faculty of Health Sciences, Winnipeg, Manitoba, Canada
| |
Collapse
|
8
|
Denoise diffusion-weighted images using higher-order singular value decomposition. Neuroimage 2017; 156:128-145. [PMID: 28416450 DOI: 10.1016/j.neuroimage.2017.04.017] [Citation(s) in RCA: 21] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/31/2016] [Revised: 02/22/2017] [Accepted: 04/06/2017] [Indexed: 11/21/2022] Open
Abstract
Noise usually affects the reliability of quantitative analysis in diffusion-weighted (DW) magnetic resonance imaging (MRI), especially at high b-values and/or high spatial resolution. Higher-order singular value decomposition (HOSVD) has recently emerged as a simple, effective, and adaptive transform to exploit sparseness within multidimensional data. In particular, the patch-based HOSVD denoising has demonstrated superb performance when applied to T1-, T2-, and proton density-weighted MRI data. In this study, we aim to investigate the feasibility of denoising DW data using the HOSVD transform. With the low signal-to-noise ratio in typical DW data, the patch-based HOSVD denoising suffers from stripe artifacts in homogeneous regions because of the HOSVD bases learned from the noisy patches. To address this problem, we propose a novel denoising method. It first introduces a global HOSVD-based denoising as a prefiltering stage to guide the subsequent patch-based HOSVD denoising stage. The HOSVD bases from the patch groups in prefiltered images are then used to transform the noisy patch groups in original DW data. Experiments were performed using simulated and in vivo DW data. Results show that the proposed method significantly reduces stripe artifacts compared with conventional patch-based HOSVD denoising methods, and outperforms two state-of-the-art denoising methods in terms of denoising quality and diffusion parameters estimation.
Collapse
|
9
|
Rao S, Ibrahim JG, Cheng J, Yap PT, Zhu H. SR-HARDI: Spatially Regularizing High Angular Resolution Diffusion Imaging. J Comput Graph Stat 2016; 25:1195-1211. [PMID: 27974868 DOI: 10.1080/10618600.2015.1105750] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/22/2022]
Abstract
High angular resolution diffusion imaging (HARDI) has recently been of great interest in mapping the orientation of intra-voxel crossing fibers, and such orientation information allows one to infer the connectivity patterns prevalent among different brain regions and possible changes in such connectivity over time for various neurodegenerative and neuropsychiatric diseases. The aim of this paper is to propose a penalized multi-scale adaptive regression model (PMARM) framework to spatially and adaptively infer the orientation distribution function (ODF) of water diffusion in regions with complex fiber configurations. In PMARM, we reformulate the HARDI imaging reconstruction as a weighted regularized least-squares regression (WRLSR) problem. Similarity and distance weights are introduced to account for spatial smoothness of HARDI, while preserving the unknown discontinuities (e.g., edges between white matter and grey matter) of HARDI. The L1 penalty function is introduced to ensure the sparse solutions of ODFs, while a scaled L1 weighted estimator is calculated to correct the bias introduced by the L1 penalty at each voxel. In PMARM, we integrate the multiscale adaptive regression models (Li et al., 2011), the propagation-separation method (Polzehl and Spokoiny, 2000), and Lasso (least absolute shrinkage and selection operator) (Tibshirani, 1996) to adaptively estimate ODFs across voxels. Experimental results indicate that PMARM can reduce the angle detection errors on fiber crossing area and provide more accurate reconstruction than standard voxel-wise methods.
Collapse
Affiliation(s)
- Shangbang Rao
- Department of Biostatistics and Biomedical Research Imaging Center University of North Carolina at Chapel Hill Chapel Hill, NC 27599, USA
| | - Joseph G Ibrahim
- Department of Biostatistics and Biomedical Research Imaging Center University of North Carolina at Chapel Hill Chapel Hill, NC 27599, USA
| | - Jian Cheng
- Department of Biostatistics and Biomedical Research Imaging Center University of North Carolina at Chapel Hill Chapel Hill, NC 27599, USA
| | - Pew-Thian Yap
- Department of Biostatistics and Biomedical Research Imaging Center University of North Carolina at Chapel Hill Chapel Hill, NC 27599, USA
| | - Hongtu Zhu
- Department of Biostatistics and Biomedical Research Imaging Center University of North Carolina at Chapel Hill Chapel Hill, NC 27599, USA
| |
Collapse
|
10
|
Mazumder R, Clymer BD, Mo X, White RD, Kolipaka A. Adaptive anisotropic gaussian filtering to reduce acquisition time in cardiac diffusion tensor imaging. Int J Cardiovasc Imaging 2016; 32:921-34. [PMID: 26843150 DOI: 10.1007/s10554-016-0848-6] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Submit a Manuscript] [Subscribe] [Scholar Register] [Received: 10/05/2015] [Accepted: 01/28/2016] [Indexed: 10/22/2022]
Abstract
Diffusion tensor imaging (DTI) is used to quantify myocardial fiber orientation based on helical angles (HA). Accurate HA measurements require multiple excitations (NEX) and/or several diffusion encoding directions (DED). However, increasing NEX and/or DED increases acquisition time (TA). Therefore, in this study, we propose to reduce TA by implementing a 3D adaptive anisotropic Gaussian filter (AAGF) on the DTI data acquired from ex-vivo healthy and infarcted porcine hearts. DTI was performed on ex-vivo hearts [9-healthy, 3-myocardial infarction (MI)] with several combinations of DED and NEX. AAGF, mean (AVF) and median filters (MF) were applied on the primary eigenvectors of the diffusion tensor prior to HA estimation. The performance of AAGF was compared against AVF and MF. Root mean square error (RMSE), concordance correlation-coefficients and Bland-Altman's technique was used to determine optimal combination of DED and NEX that generated the best HA maps in the least possible TA. Lastly, the effect of implementing AAGF on the infarcted porcine hearts was also investigated. RMSE in HA estimation for AAGF was lower compared to AVF or MF. Post-filtering (AAGF) fewer DED and NEX were required to achieve HA maps with similar integrity as those obtained from higher NEX and/or DED. Pathological alterations caused in HA orientation in the MI model were preserved post-filtering (AAGF). Our results demonstrate that AAGF reduces TA without affecting the integrity of the myocardial microstructure.
Collapse
Affiliation(s)
- Ria Mazumder
- Department of Electrical and Computer Engineering, The Ohio State University, 205 Dreese Laboratories, 2015 Neil Avenue, Columbus, OH, 43210, USA.,Department of Radiology, The Ohio State University, Room 460, 395 West 12th Avenue, 4th Floor, Columbus, OH, 43210, USA
| | - Bradley D Clymer
- Department of Electrical and Computer Engineering, The Ohio State University, 205 Dreese Laboratories, 2015 Neil Avenue, Columbus, OH, 43210, USA
| | - Xiaokui Mo
- Department of Biomedical Informatics, Center for Biostatistics, Room 320D, Lincoln Tower, 1800 Cannon Drive, Columbus, OH, 43210, USA
| | - Richard D White
- Department of Radiology, The Ohio State University, Room 460, 395 West 12th Avenue, 4th Floor, Columbus, OH, 43210, USA.,Division of Cardiovascular Medicine, Department of Internal Medicine, The Ohio State University, 244 Davis Heart and Lung Research Institute, 473 W. 12th Avenue, Columbus, OH, 43210, USA
| | - Arunark Kolipaka
- Department of Radiology, The Ohio State University, Room 460, 395 West 12th Avenue, 4th Floor, Columbus, OH, 43210, USA. .,Division of Cardiovascular Medicine, Department of Internal Medicine, The Ohio State University, 244 Davis Heart and Lung Research Institute, 473 W. 12th Avenue, Columbus, OH, 43210, USA.
| |
Collapse
|