- Workshop presentation
Efficient calculation of g-factors for CG-SENSE in high dimensions: noise amplification in random undersampling
Journal of Cardiovascular Magnetic Resonancevolume 16, Article number: W28 (2014)
SENSE [1, 2] is one of the most used parallel imaging techniques. In , uniform undersampling was employed to efficiently reconstruct an unalised image, whereas in , a conjugate gradient-based method (CG-SENSE) was used for reconstruction with arbitrary trajectories. SENSE framework allows the calculation of g-factors, characterizing the noise amplification for a given k-space trajectory and coil configuration . However, calculation of g-factors for arbitrary trajectories in high dimensions is time-consuming . Furthermore, noise characteristics of random undersampling, used in compressed sensing, is not well-understood. In this work, we use a Monte-Carlo (MC) method for fast calculation of g-factors for CG-SENSE similar to [4, 5] and apply it to random Cartesian undersampling trajectories. Theory: SENSE involves a pre-whitening step [1, 2], thus without loss of generality, we assume white noise. SENSE reconstruction solves minm ||y - Em||2, where E is the system matrix, and y are the undersampled measurements. The g-factor for the kth voxel is given by gk = √([E*E]-1k,k [E*E]k,k). Inverting E*E is not feasible in high dimensions. Instead we note the gk corresponds to the kth diagonal of the reconstruction noise covariance matrix (for normalized coil sensitivities), where nrecon = (E*E)-1E*nmeas, and nmeas is measurement noise with identity covariance matrix. We calculate the sample correlation matrix using a MC approach (since sample mean goes to 0), as 1/(p-1)∑p nprecon (nprecon)* for p instances of nrecon. Note we only calculate and store the diagonal elements of this matrix, significantly increasing efficiency.
The MC method was first verified in a numerical simulation, where the g-factor was explicitly calculated for a 2D coil configuration, to determine how many MC simulations suffice. Whole-heart imaging was performed with an isotropic resolution of 1.3 mm using a 32-channel coil array. Two 4-fold accelerated acquisitions were performed, one with uniform undersampling (2 × 2 in the ky-kz plane) and one with random undersampling. Coil sensitivity maps were exported. Images were reconstructed using SENSE (for uniform) and CG-SENSE (for both). g-factors were also calculated with the proposed approach.
Figure 1 shows the results of numerical simulations, indicating the method converges in ~50 MC simulations. Figure 2 shows the reconstructions associated with the two undersampling patterns and reconstructions, and the corresponding g-factors respectively. The results exhibit the semi-convergence property for random undersampling but not for uniform. Furthermore, the g-factor for random undersampling is smaller at its convergent point than for uniform.
g-factors for random undersampling is better than those for uniform at high k-space dimensions and high acceleration rates.
Pruessmann : MRM. 1999
Pruessmann : MRM. 2001
Liu : ISMRM. 2008
Thunberg : MRI. 2007
Robson : MRM. 2008