An augmented Lagrangian based compressed sensing reconstruction for non-Cartesian magnetic resonance imaging without gridding and regridding at every iteration

• Published in: PloS one 2014-09, Vol.9 (9), p.e107107-e107107
• Main Authors: Akçakaya, Mehmet, Nam, Seunghoon, Basha, Tamer A, Kawaji, Keigo, Tarokh, Vahid, Nezafat, Reza
• Format: Article
• Language: English
• Subjects: Adult; Algorithms; Complexity; Computation; Computer applications; Convergence; Data Compression; Female; Fourier transforms; Humans; Image acquisition; Image Processing, Computer-Assisted; Image quality; Image reconstruction; Imaging, Three-Dimensional; Iterative methods; Lagrange multiplier; Magnetic resonance; Magnetic resonance imaging; Magnetic Resonance Imaging - methods; Male; Mathematical programming; Medical imaging; Medical schools; Medicine; Medicine and Health Sciences; Middle Aged; NMR; Nuclear magnetic resonance; Phantoms, Imaging; Sharpness; Time Factors; Trajectories; Young Adult
• ISSN: 1932-6203; 1932-6203
• DOI: 10.1371/journal.pone.0107107

Non-Cartesian trajectories are used in a variety of fast imaging applications, due to the incoherent image domain artifacts they create when undersampled. While the gridding technique is commonly utilized for reconstruction, the incoherent artifacts may be further removed using compressed sensing (CS). CS reconstruction is typically done using conjugate-gradient (CG) type algorithms, which require gridding and regridding to be performed at every iteration. This leads to a large computational overhead that hinders its applicability. We sought to develop an alternative method for CS reconstruction that only requires two gridding and one regridding operation in total, irrespective of the number of iterations. This proposed technique is evaluated on phantom images and whole-heart coronary MRI acquired using 3D radial trajectories, and compared to conventional CS reconstruction using CG algorithms in terms of quantitative vessel sharpness, vessel length, computation time, and convergence rate. Both CS reconstructions result in similar vessel length (P = 0.30) and vessel sharpness (P = 0.62). The per-iteration complexity of the proposed technique is approximately 3-fold lower than the conventional CS reconstruction (17.55 vs. 52.48 seconds in C++). Furthermore, for in-vivo datasets, the convergence rate of the proposed technique is faster (60±13 vs. 455±320 iterations) leading to a ∼23-fold reduction in reconstruction time. The proposed reconstruction provides images of similar quality to the conventional CS technique in terms of removing artifacts, but at a much lower computational complexity.