Effects of Non-Local Diffusion on Structural MRI Preprocessing and Default Network Mapping: Statistical Comparisons with Isotropic/Anisotropic Diffusion

Neuroimaging community usually employs spatial smoothing to denoise magnetic resonance imaging (MRI) data, e.g., Gaussian smoothing kernels. Such an isotropic diffusion (ISD) based smoothing is widely adopted for denoising purpose due to its easy implementation and efficient computation. Beyond these advantages, Gaussian smoothing kernels tend to blur the edges, curvature and texture of images. Researchers have proposed anisotropic diffusion (ASD) and non-local diffusion (NLD) kernels. We recently demonstrated the effect of these new filtering paradigms on preprocessing real degraded MRI images from three individual subjects. Here, to further systematically investigate the effects at a group level, we collected both structural and functional MRI data from 23 participants. We first evaluated the three smoothing strategies' impact on brain extraction, segmentation and registration. Finally, we investigated how they affect subsequent mapping of default network based on resting-state functional MRI (R-fMRI) data. Our findings suggest that NLD-based spatial smoothing maybe more effective and reliable at improving the quality of both MRI data preprocessing and default network mapping. We thus recommend NLD may become a promising method of smoothing structural MRI images of R-fMRI pipeline.

[1]  P. Coupé,et al.  Impact of Non-local Means filtering on Brain Tissue Segmentation , 2010 .

[2]  B. Biswal,et al.  The resting brain: unconstrained yet reliable. , 2009, Cerebral cortex.

[3]  Hae Yong Kim,et al.  Robust anisotropic diffusion to produce enhanced statistical parametric map from noisy fMRI , 2005, Comput. Vis. Image Underst..

[4]  R. Buckner,et al.  Functional-Anatomic Fractionation of the Brain's Default Network , 2010, Neuron.

[5]  Christian Windischberger,et al.  Toward discovery science of human brain function , 2010, Proceedings of the National Academy of Sciences.

[6]  D. Schacter,et al.  The Brain's Default Network , 2008, Annals of the New York Academy of Sciences.

[7]  Guy Gilboa,et al.  Nonlocal Operators with Applications to Image Processing , 2008, Multiscale Model. Simul..

[8]  Judith D. Singer,et al.  Using SAS PROC MIXED to Fit Multilevel Models, Hierarchical Models, and Individual Growth Models , 1998 .

[9]  P. Lions,et al.  Axioms and fundamental equations of image processing , 1993 .

[10]  Gregory Camilli,et al.  Application of a Method of Estimating DIF for Polytomous Test Items , 1999 .

[11]  Alan C. Evans,et al.  Growing Together and Growing Apart: Regional and Sex Differences in the Lifespan Developmental Trajectories of Functional Homotopy , 2010, The Journal of Neuroscience.

[12]  M. Jenkinson Non-linear registration aka Spatial normalisation , 2007 .

[13]  Jean-Michel Morel,et al.  A non-local algorithm for image denoising , 2005, 2005 IEEE Computer Society Conference on Computer Vision and Pattern Recognition (CVPR'05).

[14]  Xiu-Xia Xing,et al.  PDE-based spatial smoothing: a practical demonstration of impacts on MRI brain extraction, tissue segmentation and registration. , 2011, Magnetic resonance imaging.

[15]  Pierrick Coupé,et al.  Author manuscript, published in "Journal of Magnetic Resonance Imaging 2010;31(1):192-203" DOI: 10.1002/jmri.22003 Adaptive Non-Local Means Denoising of MR Images with Spatially Varying Noise Levels , 2010 .

[16]  Paul M. Thompson,et al.  A Parameterization-Based Numerical Method for Isotropic and Anisotropic Diffusion Smoothing on Non-Flat Surfaces , 2009, IEEE Transactions on Image Processing.

[17]  Bharat B. Biswal,et al.  The oscillating brain: Complex and reliable , 2010, NeuroImage.

[18]  Bin Dong,et al.  Level Set Based Nonlocal Surface Restoration , 2008, Multiscale Model. Simul..

[19]  Joachim Weickert,et al.  Anisotropic diffusion in image processing , 1996 .

[20]  Moo K. Chung,et al.  Deformation-based surface morphometry applied to gray matter deformation , 2003, NeuroImage.

[21]  Xi-Nian Zuo,et al.  Reliable intrinsic connectivity networks: Test–retest evaluation using ICA and dual regression approach , 2010, NeuroImage.

[22]  Isabelle Bloch,et al.  A primal sketch of the cortex mean curvature: a morphogenesis based approach to study the variability of the folding patterns , 2003, IEEE Transactions on Medical Imaging.

[23]  Ayse Pinar Saygin,et al.  Smoothing and cluster thresholding for cortical surface-based group analysis of fMRI data , 2006, NeuroImage.

[24]  Yong He,et al.  Graph Theoretical Analysis of Functional Brain Networks: Test-Retest Evaluation on Short- and Long-Term Resting-State Functional MRI Data , 2011, PloS one.

[25]  Guido Gerig,et al.  Nonlinear anisotropic filtering of MRI data , 1992, IEEE Trans. Medical Imaging.

[26]  Pierre Kornprobst,et al.  Mathematical problems in image processing - partial differential equations and the calculus of variations , 2010, Applied mathematical sciences.

[27]  G. Aubert,et al.  Mathematical Problems in Image Processing: Partial Differential Equations and the Calculus of Variations (Applied Mathematical Sciences) , 2006 .

[28]  Stephen M Smith,et al.  Fast robust automated brain extraction , 2002, Human brain mapping.

[29]  Pierrick Coupé,et al.  An Optimized Blockwise Nonlocal Means Denoising Filter for 3-D Magnetic Resonance Images , 2008, IEEE Transactions on Medical Imaging.

[30]  Alan C. Evans,et al.  Cortical thickness analysis examined through power analysis and a population simulation , 2005, NeuroImage.

[31]  Jean-Michel Morel,et al.  Image Denoising Methods. A New Nonlocal Principle , 2010, SIAM Rev..

[32]  O. Sporns,et al.  Network centrality in the human functional connectome. , 2012, Cerebral cortex.

[33]  D Le Bihan,et al.  Detection of fMRI activation using Cortical Surface Mapping , 2001, Human brain mapping.

[34]  José V. Manjón,et al.  Improved estimates of partial volume coefficients from noisy brain MRI using spatial context , 2010, NeuroImage.

[35]  Enmin Song,et al.  Denoising 3D MR images by the enhanced non-local means filter for Rician noise. , 2010, Magnetic resonance imaging.

[36]  Paul H. C. Eilers,et al.  Enhancing scatterplots with smoothed densities , 2004, Bioinform..

[37]  Pierrick Coupé,et al.  Rician Noise Removal by Non-Local Means Filtering for Low Signal-to-Noise Ratio MRI: Applications to DT-MRI , 2008, MICCAI.

[38]  Stephen M. Smith,et al.  Segmentation of brain MR images through a hidden Markov random field model and the expectation-maximization algorithm , 2001, IEEE Transactions on Medical Imaging.

[39]  José V. Manjón,et al.  MRI denoising using Non-Local Means , 2008, Medical Image Anal..

[40]  Jean-Michel Morel,et al.  A Review of Image Denoising Algorithms, with a New One , 2005, Multiscale Model. Simul..

[41]  A. Anderson,et al.  Reduction of noise in diffusion tensor images using anisotropic smoothing , 2005, Magnetic resonance in medicine.

[42]  Mark W. Woolrich,et al.  Advances in functional and structural MR image analysis and implementation as FSL , 2004, NeuroImage.

[43]  Anders M. Dale,et al.  Cortical Surface-Based Analysis I. Segmentation and Surface Reconstruction , 1999, NeuroImage.