An improved EM algorithm for solving MLE in constrained diffusion kurtosis imaging of human brain

TL;DR

This paper introduces an improved EM algorithm for maximum likelihood estimation in diffusion kurtosis imaging, effectively handling physical constraints and noise, with demonstrated success on synthetic and real human brain data.

Contribution

It develops a novel EM-based estimation scheme for DKI that incorporates physical constraints and noise modeling, enhancing accuracy and computational efficiency.

Findings

The method outperforms existing algorithms in accuracy.

It effectively enforces positivity and bounds on kurtosis.

Demonstrated on synthetic and real brain data.

Abstract

The displacement distribution of a water molecular is characterized mathematically as Gaussianity without considering potential diffusion barriers and compartments. However, this is not true in real scenario: most biological tissues are comprised of cell membranes, various intracellular and extracellular spaces, and of other compartments, where the water diffusion is referred to have a non-Gaussian distribution. Diffusion kurtosis imaging (DKI), recently considered to be one sensitive biomarker, is an extension of diffusion tensor imaging, which quantifies the degree of non-Gaussianity of the diffusion. This work proposes an efficient scheme of maximum likelihood estimation (MLE) in DKI: we start from the Rician noise model of the signal intensities. By augmenting a Von-Mises distributed latent phase variable, the Rician likelihood is transformed to a tractable joint density without loss of generality. A fast computational method, an expectation-maximization (EM) algorithm for MLE is proposed in DKI. To guarantee the physical relevance of the diffusion kurtosis we apply the ternary quartic (TQ) parametrization to utilize its positivity, which imposes the upper bound to the kurtosis. A Fisher-scoring method is used for achieving fast convergence of the individual diffusion compartments. In addition, we use the barrier method to constrain the lower bound to the kurtosis. The proposed estimation scheme is conducted on both synthetic and real data with an objective of healthy human brain. We compared the method with the other popular ones with promising performance shown in the results.