Robust ODF smoothing for accurate estimation of fiber orientation.

Q-ball imaging was presented as a model free, linear and multimodal diffusion sensitive approach to reconstruct diffusion orientation distribution function (ODF) using diffusion weighted MRI data. The ODFs are widely used to estimate the fiber orientations. However, the smoothness constraint was proposed to achieve a balance between the angular resolution and noise stability for ODF constructs. Different regularization methods were proposed for this purpose. However, these methods are not robust and quite sensitive to the global regularization parameter. Although, numerical methods such as L-curve test are used to define a globally appropriate regularization parameter, it cannot serve as a universal value suitable for all regions of interest. This may result in over smoothing and potentially end up in neglecting an existing fiber population. In this paper, we propose to include an interpolation step prior to the spherical harmonic decomposition. This interpolation based approach is based on Delaunay triangulation provides a reliable, robust and accurate smoothing approach. This method is easy to implement and does not require other numerical methods to define the required parameters. Also, the fiber orientations estimated using this approach are more accurate compared to other common approaches.

Maximum likelihood estimation for Rician distributed data in analytical q-ball imaging.

Analytical q-ball imaging is widely used for reconstruction of orientation distribution function (ODF) using diffusion weighted MRI data. Estimating the spherical harmonic coefficients is a critical step in this method. Least squares (LS) is widely used for this purpose assuming the noise to be additive Gaussian. However, Rician noise is considered as a more appropriate model to describe noise in MR signal. Therefore, the current estimation techniques are valid only for high SNRs with Gaussian distribution approximating the Rician distribution. The aim of this study is to present an estimation approach considering the actual distribution of the data to provide reliable results particularly for the case of low SNR values. Maximum likelihood (ML) is investigated as a more effective estimation method. However, no closed form estimator is presented as the estimator becomes nonlinear for the noise assumption of the Rician distribution. Consequently, the results of LS estimator is used as an initial guess and the more refined answer is achieved using iterative numerical methods. According to the results, the ODFs reconstructed from low SNR data are in close agreement with ODFs reconstructed from high SNRs when Rician distribution is considered. Also, the error between the estimated and actual fiber orientations was compared using ML and LS estimator. In low SNRs, ML estimator achieves less error compared to the LS estimator.