PDE-based spatial smoothing: a practical demonstration of impacts on MRI brain extraction, tissue segmentation and registration

Spatial smoothing is typically used to denoise magnetic resonance imaging (MRI) data. Gaussian smoothing kernels, associated with heat equations or isotropic diffusion (ISD), are widely adopted for this purpose because of their easy implementation and efficient computation, but despite these advantages, Gaussian smoothing kernels blur the edges, curvature and texture of images. To overcome these issues, researchers have proposed anisotropic diffusion (ASD) and non-local means [i.e., diffusion (NLD)] kernels. However, these new filtering paradigms are rarely applied to MRI analyses. In the current study, using real degraded MRI data, we demonstrated the effect of denoising using ISD, ASD and NLD kernels. Furthermore, we evaluated their impact on three common preprocessing steps of MRI data analysis: brain extraction, segmentation and registration. Results suggest that NLD-based spatial smoothing is most effective at improving the quality of MRI data preprocessing and thus should become the new standard method of smoothing in MRI data processing.