Vessel encoded arterial spin labeling using fourier encoding suitable for vascular territory mapping

Patent metadata

Field	Value
Publication number	US-9468394-B2
Application number	US-201313780323-A
Country	US
Kind code	B2
Filing date	Feb 28, 2013
Priority date	Mar 5, 2012
Publication date	Oct 18, 2016
Grant date	Oct 18, 2016

Abstract

Methods, systems, computer programs, circuits and workstations are configured to generate at least one two-dimensional color blood flow map that combines X and Y encoded primary and secondary peak data of 1-dimensional Inverse Fourier Transform images to visually indicate (i) anterior/posterior and right/left directional components of respective feeding arteries using defined colors and (ii) brightness to indicate an amount of blood flow associated with each voxel, with increased brightness associated with increased blood flow.

First claim

That which is claimed:

1. A method of mapping brain vascular perfusion using MRI comprising:
acquiring a set of multiple pseudo-continuous arterial spin label (PCASL) images with varying phase offsets,
acquiring a set of pseudo-continuous arterial spin label (PCASL) images with X directional encoding;
acquiring a set of PCASL images with Y directional encoding; electronically generating complex signals per voxel; then
electronically applying a one-dimensional inverse Fourier transform (IFT) per voxel of the X and Y encoded images; then
electronically evaluating each voxel of the X and Y encoded images using the applied one-dimensional IFT to determine X, Y locations of primary and secondary voxel intensity peaks to localize feeding arteries, wherein the primary voxel intensity peak is a highest voxel intensity and the secondary voxel intensity peak is a second highest voxel intensity; and
generating at least one two-dimensional color blood flow map that combines X and Y encoded primary and secondary peak data to visually indicate anterior/posterior and right/left directional components of respective feeding arteries using defined colors, the map further configured to visually indicate an amount of blood flow associated with each voxel using brightness of intensity, with increased brightness associated with increased blood flow.

2. The method of claim 1 , wherein the acquiring the PCASL encoded images is carried out by first acquiring PCASL images using a sinusoidal modulation with four phase offsets to form real and imaginary parts of a complex modulation and using gradients that are the same in X or Y directions across an entire tagging period, and wherein the applied IFT can resolve multiple feeding vessel locations in the corresponding X or Y directions.

3. The method of claim 1 , wherein the at least one color blood flow map is overlaid onto a time of flight image on a display and the color map includes visual indicia comprising circles that identifies each location of four arteries at a tagging plane.

4. The method of claim 1 , wherein the color map indicates the associated feeding artery per voxel in unique colors.

5. The method of claim 4 , wherein a right internal carotid artery is shown in blue, a left internal carotid artery is shown in red, a right vertebral artery is shown in cyan, and a left vertebral artery is shown in yellow.

6. The method of claim 1 , wherein the generating the at least one color blood flow image is carried out to generate a plurality of image slices that are concurrently presented in a panel or window on a display, and wherein the acquired images are non-contrast enhanced images obtained with a scan time of about 5 minutes or less.

7. The method of claim 2 , wherein the sinusoidal modulation is carried out so that the sinusoidal modulation is complex with an imaginary component of the modulated signal that is 90 degrees phase offset, wherein the four phase offsets (Δθ) are provided as: a real pair at 0°/180° and an imaginary pair at 90°/270°, wherein the gradients are carried out with gradient steps: n/(γFOV det ) n: −(N−1)/2, . . . , 0, . . . , (N−1)/2,
where γ is a gyromagnetic ratio, FOV det is a predefined detection FOV set at a tagging plane, and N is the number of encoding steps in a direction which defines a resolution of vessel localization in the color map.

8. The method of claim 7 , wherein to reduce gradient steps while preserving detection resolution, positive or negative gradient steps are skipped and synthesized by taking a complex conjugate of mirrored steps based on Hermitian symmetry after removing a phase component of DC signal from all acquired data which is acquired without gradient encoding.

9. The method of claim 1 , further comprising displaying the at least one map on a clinician workstation.

10. The method of claim 1 , wherein the acquiring steps are carried out using a processor associated with an MR Scanner.

11. The method of claim 1 , further comprising generating a series of mixing ratio color blood flow maps using a ratio of encoded primary and secondary peak data.

12. The method of claim 1 , wherein the IFT is applied per voxel after acquiring the PCASL images with phase offsets with respective X or Y encoding and before generating the color blood flow maps, and wherein the color blood flow maps are generated using two-dimensional locations of primary voxel intensity peaks identified by the applied IFT to resolve feeding artery location.

13. The method of claim 12 , wherein the generating step can identify multiple source arteries that feed blood to a single voxel using only the acquired PCASL images and the applied 1-D IFT.

14. The method of claim 1 , wherein the applied one-dimensional IFT converts a modulated complex signal into projected ASL signals onto the X or Y encoding direction to localize the feeding arteries.

15. The method of claim 1 , wherein the method further comprises electronically identifying locations of the primary voxel intensity peak and the secondary voxel intensity peak and correlating to a physical location in 2D or 3D space to represent an amount of blood flow and location of a primary perfusion component with peak locations of an individual voxel in two orthogonal directions to estimate 2D location of a feeding artery for the color blood flow map, and wherein the second highest intensity peak can identify when mixing due to multiple arterial sources is present.

16. An image processing circuit configured to carry out claim 1 .

17. A data processing circuit comprising at least one processor configured to: (i) acquire a set of multiple pseudo-continuous arterial spin label (PCASL) images with varying phase offsets, (ii) acquire a set of non-contrast enhanced pseudo-continuous arterial spin label (PCASL) images with X directional encoding; (iii) acquire a set of PCASL images with Y directional encoding; (iv) apply a one-dimensional inverse Fourier transform (IFT) per voxel of the X and Y encoded images; then (v) evaluate each voxel of the X and Y encoded images to determine two-dimensional locations of primary and secondary voxel intensity peaks based on the IFT, wherein the primary voxel intensity peak is a highest voxel intensity and the secondary voxel intensity peak is a second highest voxel intensity; and then (vi) generate at least one two-dimensional color blood flow map that combines X and Y encoded primary and secondary peak data to visually indicate (a) anterior/posterior and right/left directional components of respective feeding arteries using defined colors and (b) an amount of blood flow associated with each voxel using brightness, with increased brightness associated with increased blood flow.

18. The data processing circuit of claim 17 , wherein the circuit identifies locations of the primary voxel intensity peak and the secondary voxel intensity peak and correlates a physical location in 2D or 3D space to represent an amount of blood flow and location of a primary perfusion component with peak locations of an individual voxel in two orthogonal directions to estimate a 2D location of a feeding artery for the color blood flow map, and wherein a second highest intensity peak can identify when mixing due to multiple arterial sources is present.

19. A computer program product comprising non-transitory computer readable storage medium having a non-transitory computer readable program code embodied in the medium, the computer-readable program code configured to perform by at least one processor the following steps:
acquiring a set of multiple pseudo-c

Assignees

• Univ Wake Forest Health Sciences

Inventors

• Jung Youngkyoo

Classifications

• A61B2576/026

  for the brain · CPC title

• G06T2207/10088

  Magnetic resonance imaging [MRI] · CPC title

• A61B5/7425

  Displaying combinations of multiple images regardless of image source, e.g. displaying a reference anatomical image with a live image · CPC title

• A61B5/0263

  using NMR · CPC title

• G06T2207/30016

  Brain · CPC title