Extreme MRI: Large‐scale volumetric dynamic imaging from continuous non‐gated acquisitions

• Published in: Magnetic resonance in medicine Vol. 84; no. 4; pp. 1763 - 1780
• Main Authors: Ong, Frank, Zhu, Xucheng, Cheng, Joseph Y., Johnson, Kevin M., Larson, Peder E. Z., Vasanawala, Shreyas S., Lustig, Michael
• Format: Journal Article
• Language: English
• Published: United States Wiley Subscription Services, Inc 01.10.2020
• Subjects: Algorithms; Computation; Cones; Contrast Media; Data acquisition; DCE‐MRI, multiscale low rank; Gating; Image acquisition; Image Enhancement; Image Interpretation, Computer-Assisted; Image Processing, Computer-Assisted; Image reconstruction; Imaging, Three-Dimensional; Lung - diagnostic imaging; Magnetic Resonance Imaging; Memory management; Optimization; Pulmonary MRI; Resource management; stochastic optimization; Stochasticity; volumetric dynamic MRI; Pulmonary MRI; DCE-MRI, multiscale low rank; volumetric dynamic MRI; stochastic optimization
• ISSN: 0740-3194; 1522-2594
• DOI: 10.1002/mrm.28235

Purpose To develop a framework to reconstruct large‐scale volumetric dynamic MRI from rapid continuous and non‐gated acquisitions, with applications to pulmonary and dynamic contrast‐enhanced (DCE) imaging. Theory and Methods The problem considered here requires recovering 100 gigabytes of dynamic volumetric image data from a few gigabytes of k‐space data, acquired continuously over several minutes. This reconstruction is vastly under‐determined, heavily stressing computing resources as well as memory management and storage. To overcome these challenges, we leverage intrinsic three‐dimensional (3D) trajectories, such as 3D radial and 3D cones, with ordering that incoherently cover time and k‐space over the entire acquisition. We then propose two innovations: (a) A compressed representation using multiscale low‐rank matrix factorization that constrains the reconstruction problem, and reduces its memory footprint. (b) Stochastic optimization to reduce computation, improve memory locality, and minimize communications between threads and processors. We demonstrate the feasibility of the proposed method on DCE imaging acquired with a golden‐angle ordered 3D cones trajectory and pulmonary imaging acquired with a bit‐reversed ordered 3D radial trajectory. We compare it with “soft‐gated" dynamic reconstruction for DCE and respiratory‐resolved reconstruction for pulmonary imaging. Results The proposed technique shows transient dynamics that are not seen in gating‐based methods. When applied to datasets with irregular, or non‐repetitive motions, the proposed method displays sharper image features. Conclusions We demonstrated a method that can reconstruct massive 3D dynamic image series in the extreme undersampling and extreme computation setting.