Systems and methods for three-dimensional spiral perfusion imaging

1. A method for perfusion imaging of a subject, comprising:

acquiring perfusion imaging data associated with the heart of a subject, the acquiring comprising applying an imaging pulse sequence with a three-dimensional stack-of-spirals trajectory;

reconstructing perfusion images from the acquired perfusion imaging data, wherein the reconstructing comprises parallel imaging and motion-guided compressed sensing;

determining absolute perfusion values from the reconstructed perfusion images wherein the determining comprises transforming a volumetric imaging dataset associated with the reconstructed perfusion images into a set of time-intensity curves representative of time-intensity relationships; and

generating a quantitative volumetric perfusion flow map based on the determined absolute perfusion values;

wherein the absolute perfusion values are determined on a pixel-wise or segmental basis, from the time-intensity relationships, to quantify myocardial blood flow of the heart of the subject.

2. The method of claim 1 , wherein the acquiring of the perfusion imaging data is performed during first pass of a contrast agent.

3. The method of claim 1 , wherein the imaging pulse sequence uses a plurality of variable density spiral waveforms.

4. The method of claim 3 , wherein the plurality of variable density spiral waveforms are used for each of a plurality of partitions configured such that volumetric image data for a whole ventricle or the whole heart of the subject is obtained.

5. The method of claim 1 , wherein the imaging pulse sequence uses spiral waveforms arranged in a linear stack.

6. The method of claim 1 , wherein the imaging pulse sequence uses spiral waveforms arranged in a spiral-radial configuration in which spiral planes are organized into a spherical pattern.

7. The method of claim 1 , wherein the spiral planes are arranged with a random radial and azimuthal angle covering a sphere in k-space.

8. The method of claim 1 , wherein the imaging pulse sequence uses spiral interleaves that are uniformly spaced in the k-z direction with an angularly uniform pattern in time.

9. The method of claim 1 , wherein acquiring the perfusion imaging data by the imaging pulse sequence comprises acquiring magnetic resonance data associated with the heart of the subject.

10. The method of claim 1 , wherein the parallel imaging and motion-guided compressed sensing uses Block Low-rank Sparsity with Motion guidance (BLOSM) combined with sensitivity encoding (SENSE).

11. The method of claim 1 , wherein determining the absolute perfusion values on the pixel-wise or segmental basis comprises using at least one of Fermi-function deconvolution, multi-compartment Kety modeling, distributed-parameter models, and model-independent deconvolution.

12. The method of claim 1 , further comprising acquiring one or more two-dimensional spiral images to quantify an arterial input function (AIF) for quantification of myocardial perfusion.

13. The method of claim 1 , further comprising acquiring two-dimensional proton density (PD) images for normalization of an arterial input function (AIF) to convert signal intensity of acquired AIF images into gadolinium concentration to perform quantification of perfusion.

14. A system for perfusion imaging of a subject, comprising:

a data acquisition device configured to acquire perfusion imaging data associated with the heart of a subject, the acquiring comprising applying an imaging pulse sequence with a three-dimensional stack-of-spirals trajectory; and

one or more processors configured to:

reconstruct perfusion images from the acquired perfusion imaging data, wherein the reconstructing comprises parallel imaging and motion-guided compressed sensing,

determine absolute perfusion values from the reconstructed perfusion images, wherein the absolute perfusion values are determined by transforming a volumetric imaging dataset associated with the reconstructed perfusion images into a set of time-intensity curves representative of time-intensity relationships, and

generate a quantitative volumetric perfusion flow map based on the determined absolute perfusion values;

wherein the absolute perfusion values are determined on a pixel-wise or segmental basis, from the time-intensity relationships, to quantify the myocardial blood flow.

15. The system of claim 14 , wherein determining the absolute perfusion values on the pixel-wise or segmental basis comprises using at least one of Fermi-function deconvolution, multi-compartment Kety modeling, distributed-parameter models, and model-independent deconvolution.

16. The system of claim 14 , wherein the data acquisition device comprises a magnetic resonance imaging device.

17. The system of claim 14 , wherein the acquiring of the perfusion imaging data is performed during first pass of a contrast agent.

18. The system of claim 14 , wherein the imaging pulse sequence uses a plurality of variable density spiral waveforms.

19. The system of claim 18 , wherein the plurality of variable density spiral waveforms are used for each of a plurality of partitions configured such that volumetric image data for a whole ventricle or the whole heart of the subject is obtained.

20. The system of claim 14 , wherein the imaging pulse sequence uses spiral waveforms arranged in a linear stack.

21. The system of claim 14 , wherein the imaging pulse sequence uses spiral waveforms arranged in a spiral-radial configuration in which spiral planes are organized into a spherical pattern.

22. The system of claim 14 , wherein the spiral planes are arranged with a random radial and azimuthal angle covering a sphere in k-space.

23. The system of claim 14 , wherein the imaging pulse sequence uses spiral interleaves that are uniformly spaced in the k-z direction with an angularly uniform pattern in time.

24. The system of claim 14 , wherein acquiring the perfusion imaging data by the imaging pulse sequence comprises acquiring magnetic resonance data associated with the heart of the subject.

25. The system of claim 14 , wherein the parallel imaging and motion-guided compressed sensing uses Block Low-rank Sparsity with Motion guidance (BLOSM) combined with sensitivity encoding (SENSE).

26. The system of claim 14 , wherein the data acquisition device and one or more processors are further configured to acquire one or more two-dimensional spiral images to quantify an arterial input function (AIF) for quantification of myocardial perfusion.

27. The system of claim 14 , wherein the data acquisition device and one or more processors are further configured to acquire two-dimensional proton density (PD) images for normalization of an arterial input function (AIF) to convert signal intensity of acquired AIF images into gadolinium concentration to perform quantification of perfusion.

28. A method for magnetic resonance imaging of a subject, comprising:

acquiring magnetic resonance data associated with an area of interest of a subject, the acquiring comprising applying an imaging pulse sequence with a three-dimensional stack-of-spirals trajectory;

reconstructing images corresponding to the area of interest from the acquired magnetic resonance data, wherein the reconstructing comprises parallel imaging and motion-guided compressed sensing using Block Low-rank Sparsity with Motion guidance (BLOSM);

quantifying, from the reconstructed images, values associated with a physiological activity in the area of interest, wherein the quantifying is based on time-intensity relationships; and

generating, from the quantified values, a quantitative volumetric map representing the physiological activity;

wherein the values are determined on a pixel-wise or segmental basis, from the time-intensity relationships, to quantify myocardial blood flow.

29. The method of claim 28 , wherein the physiological activity comprises a fluid flow in the area of interest.

30. The method of claim 28 , wherein the area of interest of the subject comprises at least a part of the heart of the subject.

31. The method of claim 30 , wherein the quantified values associated with the physiological activity area associated with absolute perfusion values.

32. The method of claim 28 , wherein quantifying the values using the time-intensity relationships comprises transforming a volumetric imaging dataset associated with the reconstructed images into a set of time-intensity curves.

33. The method of claim 28 , wherein generating the quantitative volumetric map comprises generating a volumetric perfusion flow map based on the quantified values.

34. The method of claim 28 , wherein the acquiring of the magnetic resonance data is performed during first pass of a contrast agent introduced in the area of interest.

35. The method of claim 28 , wherein the imaging pulse sequence uses a plurality of variable density spiral waveforms.

36. The method of claim 35 , wherein the plurality of variable density spiral waveforms are used for each of a plurality of partitions configured such that volumetric image data for the entire area of interest is obtained.

37. The method of claim 28 , wherein the imaging pulse sequence uses spiral waveforms arranged in a linear stack.

38. The method of claim 28 , wherein the imaging pulse sequence uses spiral waveforms arranged in a spiral-radial configuration in which spiral planes are organized into a spherical pattern.

39. The method of claim 28 , wherein the spiral planes are arranged with a random radial and azimuthal angle covering a sphere in k-space.

40. The method of claim 28 , wherein the imaging pulse sequence uses spiral interleaves that are uniformly spaced in the k-z direction with an angularly uniform pattern in time.

41. The method of claim 28 , wherein determining the absolute perfusion values on the pixel-wise or segmental basis comprises using at least one of Fermi-function deconvolution, multi-compartment Kety modeling, distributed-parameter models, and model-independent deconvolution.

42. The method of claim 28 , further comprising acquiring one or more two-dimensional spiral images to quantify an arterial input function (AIF) for quantification of myocardial perfusion.

43. The method of claim 28 , further comprising acquiring two-dimensional proton density (PD) images for normalization of an arterial input function (AIF) to convert signal intensity of acquired AIF images into gadolinium concentration to perform quantification of perfusion.

44. A non-transitory computer-readable medium having stored instructions that, when executed by one or more processors, cause a computing device to perform functions that comprise:

acquiring perfusion imaging data associated with the heart of a subject, the acquiring comprising applying an imaging pulse sequence with a three-dimensional stack-of-spirals trajectory;

reconstructing perfusion images from the acquired perfusion imaging data, wherein the reconstructing comprises parallel imaging and motion-guided compressed sensing;

determining absolute perfusion values from the reconstructed perfusion images, wherein the determining comprises transforming a volumetric imaging dataset associated with the reconstructed perfusion images into a set of time-intensity curves representative of time-intensity relationships; and

generating a quantitative volumetric perfusion flow map based on the determined absolute perfusion values;

wherein the absolute perfusion values are determined on a pixel-wise or segmental basis, from the time-intensity relationships, to quantify myocardial blood flow of the heart of the subject.