Three-dimensional spatially resolved phase graph framework

Abstract

Purpose: An open-source spatially resolved phase graph framework is proposed for simulating arbitrary pulse sequences in the presence of piece-wise constant gradients with arbitrary orientations in three dimensions. It generalizes the extended phase graph algorithm for analysis of nonperiodic sequences while preserving its efficiency, and is able to estimate the signal modulation in the 3D spatial domain. Methods: The framework extends the recursive magnetization-evolution algorithm to account for anisotropic diffusion and exploits a novel 3D k-space grid-merging method to balance the computational effort and memory requirements against acceptable simulation errors. A new postsimulation module is proposed to track and visualize the signal evolution both in the k-space and in the image domain, which can be used for simulating image artifacts or finding frequency-response profiles. To illustrate the developed technique, three examples are presented: (1) fast off-resonance calculation for dictionary building in MR fingerprinting, (2) validation of a steady-state sequence with quasi-isotropic diffusion weighting, and (3) investigation of the magnetization evolution in PRESS-based spectroscopic imaging. Results: The grid-merging algorithm of the proposed framework demonstrates high calculation efficiency exemplified by frequency-response simulation of pseudo steady-state or diffusion-weighted steady-state sequences. It further helps to visualize the signal evolution in PRESS-based sequences. Conclusions: The proposed simulation framework has been validated based on several different example applications for analyzing signal evolution in the frequency and spatial domain.