Compressed sensing high resolution functional magnetic resonance imaging

Patent metadata

Field	Value
Publication number	US-11357402-B2
Application number	US-202016859267-A
Country	US
Kind code	B2
Filing date	Apr 27, 2020
Priority date	Aug 31, 2015
Publication date	Jun 14, 2022
Grant date	Jun 14, 2022

Abstract

The present disclosure provides methods and systems for high-resolution functional magnetic resonance imaging (fMRI), including real-time high-resolution functional MRI methods and systems.

First claim

What is claimed is:

1. A method comprising:
acquiring image data of a target area in a subject using a variable density spiral (VDS) trajectory, wherein the VDS trajectory includes a stack of VDS trajectories, wherein a change in a radial direction versus a change in a rate in an angular direction of the VDS trajectories is determined at least in part by a number of interleaves and an effective field of view, wherein the field of view follows an exponential function and an angle between each interleaf of the number of interleaves is varied; and
producing an image of the target area in the subject based on the image data.

2. The method of claim 1 , wherein producing the image comprises reconstructing the image data using a wavelet spatial regularization.

3. The method of claim 2 , wherein the wavelet spatial regularization includes a three dimensional Dubechies-4 discrete wavelet transform.

4. The method of claim 3 , wherein the image is produced by taking an inverse of the discrete wavelet transform.

5. The method of claim 1 , wherein the number of interleaves decreases as k z increases in a k-space of the image data.

6. The method of claim 1 , wherein at least some of the interleaves on an outer portion of a k-space of the image data are randomly skipped.

7. The method of claim 6 , wherein the interleaves on the outer portion of the k-space are randomly skipped following a Gaussian distribution.

8. A method comprising:
acquiring image data of a target area in a subject using a stack of multi-interleaf variable density spiral (VDS) trajectories, wherein a total number of interleaves follows a Laplacian distribution and a center of a k-space is more densely sampled than an outer portion of the k-space; and
producing an image of the target area in the subject based on the image data.

9. The method of claim 8 further comprising:
acquiring two phase-cycled images of the target area; and
combining the two phased-cycled images by maximum intensity projection.

10. The method of claim 8 , wherein producing the image comprises iteratively reconstructing the image from the image data using an L1 regularized cost function.

11. The method of claim 10 , wherein reconstructing comprises regularizing a temporal domain and a spatial domain by a discrete cosine transform.

12. The method of claim 10 , wherein the L1 regularized cost function is solved by a gradient descent method.

13. The method of claim 8 , wherein producing the image comprises reconstructing the image from the image data using one or more regularization parameters.

14. The method of claim 13 , wherein the one or more regularization parameters includes at least one of a contrast to noise ratio, an active volume within a designed active region of the image data, mean F statistic value, normalized root mean squared error, or peak hemodynamic response function amplitude.

15. A functional magnetic resonance imaging (fMRI) system, the system comprising:
a receiver configured to acquire image data of a target area in a subject using a variable density spiral (VDS) trajectory, wherein the VDS trajectory includes a stack of VDS trajectories, wherein a change in a radial direction versus a change in a rate in an angular direction of the VDS trajectories is determined by a number of interleaves and an effective field of view, wherein the field of view follows an exponential function and an angle between each interleaf of the number of interleaves is varied; and
a processor configured to produce an image of the target area in the subject based on the acquired image data.

16. The system of claim 15 , further comprising a coil configured to apply a balanced steady state free precession (b-SSFP) sequence to the target area in the subject.

17. The system of claim 15 , wherein the processor is further configured to analyze the image data using a fast iterative shrinkage thresholding algorithm (FISTA).

18. The system of claim 15 , wherein the processor is configured to produce the image during acquisition of the image data by the receiver or immediately after acquisition of the image data by the receiver.

19. The system of claim 15 , further comprising a coil configured to apply an excitation waveform to the target area in the subject.

20. The system of claim 15 , wherein the processor is further configured to analyze the image data using a sparsifying transform.

Assignees

• Univ Leland Stanford Junior
• Univ California

Inventors

• Lee Jin Hyung
• Fang Zhongnan

Classifications

• A61B5/055 Primary

  involving electronic [EMR] or nuclear [NMR] magnetic resonance, e.g. magnetic resonance imaging · CPC title

• G01R33/5614

  using a fully balanced steady-state free precession [bSSFP] pulse sequence, e.g. trueFISP · CPC title

• G01R33/56308

  Characterization of motion or flow; Dynamic imaging · CPC title

• A61B5/0263

  using NMR · CPC title

• A61B5/0042 Primary

  for the brain · CPC title