2023

Swin golden angle: a radial profile order for golden ratio sampling with navigator compatibility and eddy-current suppression
Zhongsen Li1, Aiqi Sun2, Shuai Wang1, Chuyu Liu1, Sirui Wu1, Haozhong Sun1, Rui Guo3, Haiyan Ding1, and Rui Li1
1Center for Biomedical Imaging Research, Tsinghua University, Beijing, China, 2Institute of Science and Technology for Brain-Inspired Intelligence, Fudan University, Shanghai, China, 3School of Medical Technology, Beijing Institute of Technology, Beijing, China

Synopsis

Keywords: Data Acquisition, Cardiovascular, golden angle radial, eddy-current suppresion, profile order

Golden angle radial trajectory suffers from eddy-current when combined with trueFISP sequence. If additional navigators are required, how to keep golden ratio sampling while maintain stable steady-state is still a problem to be solved. In this work, we propose a novel radial profile order, named "swin golden angle"(swinGA), which is able to achieve golden ratio sampling, stable trueFISP signal, and navigator acquisition simultaneously. We validated the proposed profile in static phantom study and in-vivo study for free-running cardiac imaging. The results show that the proposed swinGA significantly outperforms other sampling profiles and achieves good image quality.

Introduction

Golden angle (GA) radial trajectory[1], which uses a golden ratio angle increment in sampling the k-space lines, enables high scan efficiency and improves temporal resolution for dynamic imaging. However, directly applying GA to trueFISP sequence will induce image artifacts since the large change in gradient waveform between two readouts leads to non-smooth eddy-current effects and unstable signal[2]. GA using a tiny angle increment, which is named as tiny GA (tGA)[3], could mitigate this impaction. However, when the self-navigation acquisition is integrated, as shown in Figure 1, the trajectory will become complicated since the fixed readout direction of self-navigation will disrupt the smooth angle series.

In this work, we propose a novel radial profile order, which is able to simultaneously achieve golden ratio sampling, stable signal steady-state, and navigator acquisition.

Methods

When navigator acquisition is combined with golden-angle radial trueFISP sequence, there are two major considerations. The first is that the navigator readout should be inserted into the angle series in a smooth way, in order to reduce eddy-current. The second is that the navigator should be acquired along the same direction, because the acquired data is usually not symmetric with inversed readout direction due to system imperfections(e.g., gradient delay).
We design the swin golden angle(swinGA) profile to facilitate the combination of golden-angle radial trueFISP sequence and self-navigation. We first considered the trajectory without navigator as shown in Figure 2A. First, spokes of golden angle increment are generated for each frame. Second, the angles are taken modulus to 180$$$^{\circ}$$$ and mapped into range [0$$$^{\circ}$$$,180$$$^{\circ}$$$). Third, the spokes are sorted according to the modulus angles. For odd frames, spokes are sorted in ascending order. For even frames, spokes are sorted in descending order. After these three steps, the swinGA trajectory can be generated. The spokes sweep from 0$$$^{\circ}$$$ to 180$$$^{\circ}$$$ in odd frames and sweep backward in even frames, just like a swin.
Figure 2B displays how to generate swinGA trajectory with inserted navigators. The first and second steps are the same as Figure 2A. Then, the navigator spokes are inserted into the frame. Finally, all spokes (including GA and navigator) are sorted in the same way as Figure 2A.
In swinGA profile, since the navigators are inserted into GA patterns, the maximum azimuthal gap will be restricted by the fibonacci series[1], indicating the angle increment is smooth. Besides, because the spoke readout always starts from the same half-plane, the navigator quality is guaranteed.
The pseudo code for swinGA is provided in Figure 3 for generating the profile angles.

Experiments

TrueFISP sequences with five trajectories, including Cartesian(Cart), uniform radial(Uni), GA, tGA and swinGA, were implemented on a 3T clinical scanner (Philips IngeniaCX R5.7.1; Best, Netherlands). The sequence parameters are kept the same for all experiments: FOV=300mmx300mmx8mm, Resolution=1mmx1mmx8mm, TR/TE/FA=4.2ms/2.1ms/45$$$^{\circ}$$$, Duration=15s. For the navigator-interleaved acuiqsition, 5 imaging spokes and 2 navigator spokes(45$$$^{\circ}$$$ and 135$$$^{\circ}$$$) are acquired in a frame, which corresponds to a temporal resolution of 29.4ms. Besides, we choose N=5 for tGA sequence. For cartesian sequence, each frame also contains 5 imaging lines(random phase encoding) and 2 navigator lines(fixed as ky=-1 and ky=0).

Static Phantom Study
Scanning was performed on a static water phantom. After acquisition, only the imaging spokes (or phase encoding lines) were used to fill a fully-sampled k-space and reconstruct the phantom image to depict the eddy-current appearance in static objects.

Dynamic In-vivo Study
Two healthy volunteers were recruited. To exclude irrelevant factors, breath-holding is required. A PS-model based low-rank reconstruction algorithm is used for reconstructing the image for both cartesian and radial trajectory[4]. A retrospective ECG-synchronized CINE scan is additionally performed for comparison, which shares quite similar sequence parameters.

Results

The phantom study results are shown in Figure 4. It can be seen that GA always suffers from severe eddy-current artifacts. Cart, Uni, tGA keep good image quality without navigator, but display residual eddy-current artifacts when navigators are inserted. SwinGA presents the minimum artifacts with navigator acquisition.
The in-vivo study results are shown in Figure 5. It can be seen that reconstruction fails for Uni, which is due to its low sampling efficiency (spokes are not evenly distributed), leading to highly ill-posed inverse problem. GA suffers from severe eddy-current artifacts. The tGA performs slightly better than GA, but the artifacts are still very obvious. Cart also suffers from eddy current, which obscures the image details at the heart region. The image produced by swinGA displays the best quality, and is highly consistent with retrospectively ECG-synchronized CINE image.

Discussion and Conclusion

In this work, we proposed "swin golden angle", which is able to simultaneously achieve golden ratio sampling, stable trueFISP signal, and robust navigator acquisition. Our validations indicated that the proposed SwinGA could effectively suppress the artifacts in cardiac imaging compared to other GA trajectories when self-navigation signals was sampled. Adding some fixed navigator readouts into GA sampling may be beneficial, because the dynamic information could be directly characterized by the navigators. Besides, since navigators do not overlap with GA spokes, all spokes can be utilized as imaging data without sacrificing the scan efficiency. SwinGA may contribute to the development of more advanced cardiac imaging techniques.

Acknowledgements

No acknowledgement found.

References

[1]. Winkelmann, S., Schaeffter, T., Koehler, T., Eggers, H., & Doessel, O. (2006). An optimal radial profile order based on the Golden Ratio for time-resolved MRI. IEEE transactions on medical imaging, 26(1), 68-76.

[2]. Wundrak, S., Paul, J., Ulrici, J., Hell, E., & Rasche, V. (2014). A small surrogate for the golden angle in time-resolved radial MRI based on generalized fibonacci sequences. IEEE transactions on medical imaging, 34(6), 1262-1269.

[3]. Scheffler, K., & Hennig, J. (2003). Eddy current optimized phase encoding schemes to reduce artifacts in balanced SSFP imaging.

[4]. Li, Z., Sun, A., Liu, C., Wei, H., Wang, S., Fu, M., & Li, R. (2022). Accelerated partial separable model using dimension-reduced optimization technique for ultra-fast cardiac MRI. arXiv preprint arXiv:2210.00493.

Figures

Figure 1. Illustration of three common radial profile orders: uniform(Uni), golden angle(GA) and tiny golden angle(tGA). For the default radial profile order, tGA enjoys both high sampling efficiency (spokes are more evenly distributed) and eddy-current suppresion (the angle increment between spokes are smooth). However, when navigator readouts are inserted, large angle jump can occur in all three radial profiles, which will incur eddy-current artifacts.

Figure 2. Illustration of the proposed swin golden angle(swinGA) radial profile order. First, spokes of golden angle increment are generated for each frame. Second, the angles are taken modulus and mapped into range [0, 180). Third, the spokes are sorted according to the angles. For odd frames, the spokes are sorted in ascending order. For even frames, the spokes are sorted in descending order. These three steps generate a swin-like profile order in (A). When navigator is required, the only modification is to insert the navigators into the spokes to be sorted, as shown in (B).

Figure 3. The pseudo code for the algorithm to produce the swinGA radial profile order. The bold font indicates that the variable is an array.

Figure 4. Phantom experiment results. It can be seen that GA always suffers from severe eddy-current artifacts. Cart, Uni, tGA keep good image quality without navigator, but display residual eddy-current artifacts when navigators are inserted. SwinGA presents the minimum artifacts with navigator acquisition.

Figure 5. In-vivo experiment results on healthy subjects. It can be seen that reconstruction fails for Uni, which is due to the low sampling efficiency (spokes are not evenly distributed), leading to highly ill-posed inverse problem. GA suffers from severe eddy-current artifacts. The tGA performs slightly better than GA, but the artifacts are still very obvious. Cart also suffers from eddy-current, which obscures the image details at the heart region. The image produced by swinGA displays the best quality, and is highly consistent with retrospectively ECG-synchronized CINE image.

Proc. Intl. Soc. Mag. Reson. Med. 31 (2023)
2023
DOI: https://doi.org/10.58530/2023/2023