Magnetic resonance angiography with vessel segmentation

• US 6,381,486 B1
• Filed: 05/18/1999
• Issued: 04/30/2002
• Est. Priority Date: 01/08/1999
• Status: Expired due to Term

• Current Assignee: Wisconsin Alumni Research Foundation (University of Wisconsin System)
• Original Assignee: Wisconsin Alumni Research Foundation (University of Wisconsin System)
• Inventors: Mazaheri, Yousef, Mistretta, Charles A., Grist, Thomas M., Carroll, Timothy J.
• Primary Examiner(s): Lateef, Marvin M.
• Assistant Examiner(s): Lin, Jeoyuh
• Application Number:
  US09/314,227
• Time in Patent Office: 1,078 Days
• Field of Search: 324/309, 324/318, 324/307, 358/111, 128/653.3, 600/420, 600/410, 600/413, 600/431, 600/419, 600/415, 334/309, 334/306
• US Class Current: 600/420
• CPC Class Codes: G01R 33/5619 by temporal sharing of data, e.g. keyhole, block regional interpolation scheme for k-Space [BRISK]
• CPC Class Codes: G01R 33/5635 Angiography, e.g. contrast-enhanced angiography [CE-MRA] or time-of-flight angiography [TOF-MRA]

Abstract

A magnetic resonance angiogram (MRA) is acquired using a contrast enhancement method in which a series of low resolution NMR images are rapidly acquired during a time resolved phase of the examination in which the contrast bolus makes a first pass through the arteries and veins. Additional, high resolution NMR image data is acquired in a subsequent steady-state phase of the examination from which a high resolution NMR image is reconstructed. Segmentation of arteries, veins and background in the high resolution image is performed using information in the low resolution NMR images.

Claims

1. A method for producing a magnetic resonance angiogram of a patient using an MRI system, the steps comprising:

a) injecting the patient with a contrast agent which flows into a region of interest in the patient;

b) acquiring NMR image data for a series of NMR time course images during a time resolved phase of an examination, during which the contrast agent makes a first passage through the patient'"'"'s arteries in the region of interest;

c) producing from the NMR time course images an arterial enhancement reference curve that indicates the signal level enhancement of image voxels depicting arteries during the time resolved phase;

d) acquiring NMR image data during a steady state phase of the examination which follows the time resolved phase;

e) reconstructing an image after the examination using NMR data acquired during the examination;

f) producing the magnetic resonance angiogram by segmenting the image reconstructed in step e) using the arterial enhancement reference curve.

2. The method as recited in claim 1 in which the segmentation includes reducing the brightness of voxels in the image which do not depict arteries.

3. The method as recited in claim 1 which includes producing from the NMR time course images a venous enhancement reference curve that indicates the signal enhancement of images voxels depicting veins during the time resolved phase of the examination, and the segmenting in step f employs the venous enhancement reference curve.

4. The method as recited in claim 1 in which the segmenting in step e) includes comparing the brightness levels of each image voxel during the time resolved phase of the examination with corresponding points in the arterial enhancement reference curve and modifying the brightness of the voxel in the magnetic resonance angiogram as a function of this comparison.

5. The method as recited in claim 1 in which the segmenting in step e) includes:

comparing each voxel signal level enhancement in the NMR time course images with the arterial enhancement reference curve and producing a correlation value SSD_A for each such voxel;

producing from the NMR time course images a venous enhancement reference curve that indicates the signal level enhancement of image voxels depicting veins during the time resolved phase of the examination;

comparing each voxel signal level enhancement in the NMR time course images with the venous enhancement reference curve and producing a correlation value SSD_V for each such voxel; and

each corresponding voxel in the reconstructed image is characterized by evaluating the correlation values SSD_A and SSD_V.

6. The method as recited in claim 1 in which the signal levels of voxels in the series of NMR time course images are compared with the arterial enhancement reference curve to characterize image voxels depicting arteries as arterial voxel, and the segmenting of the image reconstructed in step e) includes selecting for inclusion in the magnetic resonance angiograms images voxels therein which correspond to arterial voxels.

7. The method as recited in claim 6 in which the segmenting of the image reconstructed in step e) includes suppressing the signal level of voxels therein which are not characterized as arterial voxels.

8. The method as recited in claim 1 in which the image reconstructed in step e) is formed using NMR image data acquired during the time resolved phase and during the steady-state phase of the examination.

9. A method for producing a magnetic resonance angiogram of a patient using an MRI system, the steps comprising:

a) injecting the patient with a contrast agent which flows into a region of interest in the patient;

b) acquiring NMR central k-space data for a series of NMR time course images during a time resolved phase of an examination, during which the contrast agent makes a first passage through the patient'"'"'s arteries in the region of interest;

c) acquiring NMR k-space data during a steady state phase of the examination which follows the time resolved phase;

d) characterizing voxels in the NMR time course images that depict arteries by comparing the signal level enhancement of each voxel caused by the contrast agent with the signal level of an arterial enhancement reference curve produced from the NMR time course images and indicating signal level enhancement of image voxels depicting arteries during the time resolved phase;

e) reconstructing an image using NMR k-space data acquired during the steady-state phase of the examination; and

f) segmenting the reconstructed image by suppressing signals from voxels that have not been characterized as arterial.

10. The method as recited in claim 9 in which NMR k-space data is repeatedly acquired during the steady-state phase of the examination to produce multiple k-space samples and the multiple k-space samples are combined to form average k-space samples.

11. The method as recited in claim 10 in which the multiple k-space samples are combined to form weighted average k-space samples.

12. The method as recited in claim 9 in which step b) is performed by:

repeatedly performing a first NMR pulse sequence which acquires NMR central k-space data;

and step c) is performed by;

repeatedly performing a second NMR pulse sequence which acquires NMR k-space data from which a high resolution image may be reconstructed.

13. The method as recited in claim 12 in which the first NMR pulse sequence includes sampling k-space in a spiral trajectory.

14. The method as recited in claim 12 in which the first NMR pulse sequence includes sampling k-space from a projection angle and the projection angle is changed for each repetition of the first NMR pulse sequence.

15. The method as recited in claim 12 in which the first NMR pulse sequence includes applying phase encoding gradients which sample k-space in a centric order as the first NMR pulse sequence is repeated.

16. The method as recited in claim 9 which includes:

g) combining data from the segmented reconstructed image with each of a plurality of the NMR time course images to produce enhanced NMR time course images.

17. The method as recited in claim 16 in which step g) includes:

i) Fourier transforming the segmented reconstructed image to produce a segmented k-space data set;

ii) combining peripheral k-space data acquired in step b); and

iii) reconstructing the plurality of NMR time course images from the k-space data combined in step ii).

18. The method as recited in claim 9 in which the contrast agent is an intra-vascular material which enhances the NMR signals in the patient'"'"'s vasculature.

19. The method as recited in claim 9 which includes:

h) characterizing voxels in the NMR time course images that depict veins by analyzing the signal level enhancement of each voxel caused by the contrast agent.

20. The method as recited in claim 19 in which the comparing in step d) includes:

correlating each voxel signal level enhancement with the arterial enhancement reference curve and producing a correlation value SSD_A for each voxel;

the analyzing in step h) includes;

producing from the NMR time course images a venous enhancement reference curve that indicates the signal level enhancement of image voxels depicting veins during the time resolved phase of the examination;

correlating each voxel signal level enhancement with the venous enhancement reference curve and producing a correlation value SSD_V for each voxel; and

each voxel is characterized by evaluating the correlation values SSD_A and SSD_V.

21. The method as recited in claim 20 which includes:

evaluating the spatial location of the voxel relative to other voxels characterized as arterial.

22. The method as recited in claim 20 which includes:

evaluating the level of the voxel signal relative to the signal level of surrounding voxels.

23. The method as recited in claim 10 in which step f) includes:

i) producing a binary mask which indicates the voxels characterized as depicting arteries; and

ii) combining the binary mask with the reconstructed image to suppress the signals from non-arterial voxel.

24. The method as recited in claim 23 in which the non-arterial voxel signals are set to zero.

25. The method as recited in claim 9 in which step d) includes calculating a weighting factor for each voxel in the NMR time course images;

and step f) includes multiplying each voxel by its corresponding weighting factor.