View article

In diffusion MRI (dMRI) domain, many High Angular Resolution Diffusion Imaging (HARDI) methods were proposed to estimate Ensemble Average Propagator (EAP) and Orientation Distribution Function (ODF). They normally need many samples, which limits their applications. Some Compressive Sensing (CS) based methods were proposed to estimate ODF in Q-Ball Imaging (QBI) from limited samples. However EAP estimation is much more difficult than ODF in QBI. Recently Spherical Polar Fourier Imaging (SPFI) was proposed to represent diffusion signal using Spherical Polar Fourier (SPF) basis without specific assumption on diffusion signals and analytically obtain EAP and ODF via the Fourier dual SPF (dSPF) basis from arbitrarily sampled signal. Normally the coefficients of SPF basis are estimated via Least Square with weighted L2 norm regularization (L2-SPFI). However, L2-SPFI needs a truncated basis to avoid overfitting, which brings some estimation errors. By considering the Fourier relationship between EAP and signal and the Fourier basis pair provided in SPFI, we propose a novel EAP estimation method, named L1-SPFI, to estimate EAP from limited samples using CS technique, and favorably compare it to the classical L2-SPFI method. L1-SPFI estimates the coefficients in SPFI using least square with weighted L1 norm regularization. The weights are designed to enhance the sparsity. L1-SPFI significantly accelerates the ordinary CS based Fourier reconstruction method. This is performed by using SPF basis pair in CS estimation process which avoids the numerical Fourier transform in each iteration step. By considering high …