Multi-echo pseudo-golden angle stack of stars thermometry with high spatial and temporal resolution using k-space weighted image contrast
A method for producing an image of a subject using a magnetic resonance imaging (MRI) system includes acquiring a series of echo signals by sampling k-space along radial lines that each pass through the center of k-space. Each projection of the radial lines is divided into multiple echoes and successive projections are spaced by a predetermined angular distance. The series of echo signals are reconstructed into a plurality of images, wherein each image corresponds to a distinct echo signal.
1. A method for producing an image of a subject using a magnetic resonance imaging (MRI) system, the method comprising:
acquiring a series of echo signals by sampling k-space along radial lines that each pass through the center of k-space, wherein each projection of the radial lines is divided into multiple echoes and successive projections are spaced by a predetermined angular distance;
reconstructing the series of echo signals into a plurality of images, wherein each image corresponds to a distinct echo signal;
for each echo signal, calculating a proton resonance frequency temperature change in the subject using phase information from a first echo signal acquired in the series of echo signals as a reference phase.
2. The method of claim 1 , wherein the predetermined angular distance is a pseudo golden angle.
3. The method of claim 1 , wherein the predetermined angular distance is a factor that depends on the golden angle.
4. The method of claim 1 , wherein the predetermined angular distance is a multiple of a Fibonacci number.
5. The method of claim 1 , further comprising:
calculating an initial signal magnitude value based on the series of echo signals; and
calculating a transverse relaxation time value based on the series of echo signals.
6. The method of claim 1 , wherein the series of projections comprises (i) a series of inner projections spaced according to a first angle corresponding to a ratio of a first Fibonacci number to a maximum Fibonacci number and (ii) a series of outer projections spaced according to a second angle corresponding to a ratio of a second Fibonacci number to the maximum Fibonacci number.
7. The method of claim 1 , wherein the plurality of echo signals is reconstructed into a plurality of images using a k-space weighted image contrast (KWIC) reconstruction process.
8. The method of claim 7 , wherein the KWIC reconstruction process uses a sliding window that advances over multiple projections for each echo signal.
9. The method of claim 8 , wherein the sliding window is asymmetric in time and k-space center is reconstructed at the end of the sliding window.
10. The method of claim 1 , further comprising:
measuring a plurality of respiration phase offsets by linearly fitting the slope of the phase at k-space center from the series of echo signals as a function of TE using linear regression; and
performing respiration correction on the series of echo signals using the plurality of respiration phase offsets.
11. The method of claim 10 , wherein the MRI system comprises a plurality of coils that independently measure the slope and the method further comprises:
averaging the slope of the phase across all of the coils over multiple respiratory cycles to obtain a baseline slope; and
subtracting the baseline slope from the phase slope measured by each coil.
12. A method for performing thermometry using a magnetic resonance imaging (MRI) system, the method comprising:
acquiring a series of echo signals by sampling k-space along each projection in a series of projections using a stack of stars sequence, wherein successive projections in the series of projections are spaced apart by an angular distance equal to a pseudo golden angle; and
calculating proton resonance frequency temperature change in the subject using phase information from the series of echo signals.
13. The method of claim 12 , further comprising:
reconstructing the plurality of echo signals into a plurality of images using a KWIC reconstruction process, wherein each image corresponds to a distinct echo signal.
14. The method of claim 13 , wherein the KWIC reconstruction process uses a sliding window that advances over multiple projections for each echo signal and the sliding window is asymmetric in time and k-space center is reconstructed at the end of the sliding window.
15. The method of claim 12 , wherein the proton resonance frequency temperature change is calculated independently for each echo signal using phase information from a first echo signal acquired in the series of echo signals as a reference phase.
16. The method of claim 12 , wherein the proton resonance frequency temperature change is calculated independently for each echo signal using phase information corresponding to a reference signal acquired at the same rotation from k-space center as the echo signal as a reference phase.
17. The method of claim 12 , wherein the proton resonance frequency temperature change is calculated across the series of echo signals using phase information from a first echo signal acquired in the series of echo signals.
18. The method of claim 12 , wherein the proton resonance frequency temperature change is calculated across the series of echo signals and the proton resonance frequency temperature change for each individual echo signal is calculated using a reference signal acquired at the same rotation from k-space center as the individual echo signal as a reference phase.
19. The method of claim 12 , further comprising:
calculating an initial signal magnitude value based on the series of echo signals; and
calculating a transverse relaxation time value based on the series of echo signals.