Magnetic resonance imaging of ventilation and perfusion in the lung

• US 9,750,427 B2
• Filed: 12/07/2011
• Issued: 09/05/2017
• Est. Priority Date: 12/07/2010
• Status: Active Grant

• Current Assignee: Regents of the University of California (University of California)
• Original Assignee: Regents of the University of California (University of California)
• Inventors: Prisk, Gordon Kim, Hopkins, Susan Roberta, Buxton, Richard Bruce, Pereira De S, Rui Carlos, Theilmann, Rebecca Jean, Cronin, Matthew Vincent
• Primary Examiner(s): Chen, Tse
• Assistant Examiner(s): HOFFMAN, JOANNE M
• Application Number:
  US13/992,231
• Publication Number:
  US 20130338489A1
• Time in Patent Office: 2,099 Days
• Field of Search: None
• CPC Class Codes: A61B 5/055 involving electronic [EMR] or nuclear [NMR] magnetic resonance, e.g. magnetic resonance imaging
• CPC Class Codes: A61B 5/08 Detecting, measuring or recording devices for evaluating the respiratory organs A61B5/0205 takes precedence
• CPC Class Codes: A61B 5/0813 Measurement of pulmonary parameters by tracers, e.g. radioactive tracers
• CPC Class Codes: G01R 33/5601 involving use of a contrast agent for contrast manipulation, e.g. a paramagnetic, super-paramagnetic, ferromagnetic or hyperpolarised contrast agent
• CPC Class Codes: G01R 33/56366 Perfusion imaging

Abstract

Methods, devices, and systems are disclosed for implementing a fully quantitative non-injectable contrast proton MRI technique to measure spatial ventilation-perfusion (VA/Q) matching and spatial distribution of ventilation and perfusion. In one aspect, a method using MRI to characterize ventilation and perfusion in a lung includes acquiring an MR image of the lung having MR data in a voxel and obtaining a breathing frequency parameter, determining a water density value, a specific ventilation value, and a perfusion value in at least one voxel of the MR image based on the MR data and using the water density value to determine an air content value, and determining a ventilation-perfusion ratio value that is the product of the specific ventilation value, the air content value, the inverse of the perfusion value, and the breathing frequency.

Claims

1. A computer-implemented method using magnetic resonance imaging to characterize ventilation and perfusion in a lung, comprising:

acquiring, using a magnetic resonance imaging (MM) machine under the control of a processing unit, a magnetic resonance (MR) image of the lung that includes MR data in a voxel of the MR image, wherein the acquiring causes the MM machine to measure a breathing frequency value;

determining, by the processing unit a water density value and a perfusion value in at least one voxel of the MR image based on the MR data;

determining, by the processing unit from the water density value, an air content value in the at least one voxel;

determining, by the processing unit, a specific ventilation value in the at least one voxel based on a time delay between an onset of a stimulus to the lung and a response detected in the at least one voxel;

determining, by the processing unit a ventilation-perfusion ratio value that is the product of the specific ventilation value, the air content value, the inverse of the perfusion value, and the breathing frequency value; and

evaluating, by the processing unit, a spatial distribution of the specific ventilation value, perfusion value, and the ventilation-perfusion ratio by computing a standard deviation of the specific ventilation value and perfusion value;

wherein the specific ventilation value is computed by SVI*(1−

density)*f_B wherein SVI represents quantitative data from the MR image, (1−

density) represents the air content value, and f_B represents the breathing frequency value, thereby determining when the lung is diseased.

2. The method of claim 1, wherein the MR image is a 2D image or a 3D image.

3. The method of claim 1, wherein the determining the ventilation-perfusion ratio value is determined voxel-by-voxel.

4. The method of claim 1, wherein the acquiring the MR image uses oxygen as an inhaled contrast agent.

5. The method of claim 1, further comprising determining an alveolar ventilation value that is the product of the specific ventilation value, the air content value, and the breathing frequency value.

6. The method of claim 5, further comprising:

computing a first variance in the specific ventilation value, a second variance in the perfusion value, and a third variance in the ventilation-perfusion ratio value; and

generating a correlation coefficient from one or more of the first variance, the second variance, and the third variance.

7. The method of claim 1, further comprising producing an image that shows the spatial distribution of ventilation, perfusion, and ventilation-perfusion for each voxel associated with the acquired MR image.

8. The method of claim 1, wherein the spatial distribution of ventilation, perfusion, and ventilation-perfusion is used as an indicator of a diseased lung.

9. The method of claim 1, further comprising:

determining ventilation-perfusion matching in the lung from at least the ventilation-perfusion ratio value.

10. The method of claim 9, wherein ventilation-perfusion matching is used as an indicator of a diseased lung.

11. A magnetic resonance imaging system to characterize ventilation and perfusion in a lung, comprising:

a magnetic resonance imaging machine that acquires an MR image of the lung; and

a processing unit having a processor and a non-transitory computer-readable storage medium having instructions stored thereon that when executed by the processor cause the processing unit to at least;

command the magnetic resonance imaging machine to acquire the MR image that includes MR data in each voxel of the MR image,

process the MR data to determine a water density value and a perfusion value in at least one voxel of the MR image, wherein the water density value is used to determine an air content value,

measure a breathing frequency of the lung;

determine a specific ventilation value in the at least one voxel based on a time delay between an onset of a stimulus to the lung and a response detected in the at least one voxel;

determine a ventilation-perfusion ratio value that is the product of the specific ventilation value, the air content value, the inverse of the perfusion value, and the breathing frequency; and

evaluate a spatial distribution of the specific ventilation value, perfusion value, and the ventilation-perfusion ratio by computing a standard deviation of the specific ventilation value and perfusion value;

wherein the specific ventilation value is computed by SVI*(1−

density)*f_B wherein SVI represents quantitative data from the MR image, (1−

density) represents the air content value, and f_B represents the breathing frequency value, thereby determining when the lung is diseased.

12. The magnetic resonance imaging system of claim 11, wherein the MRI system uses oxygen as an inhaled contrast agent to acquire the MR image of the lung.

13. The magnetic resonance imaging system of claim 11, wherein the instructions, when executed by the processor, further cause the processing unit to at least:

determine an alveolar ventilation value that is the product of the specific ventilation value, the air content value, and the breathing frequency.

14. The magnetic resonance imaging system of claim 13, wherein the instructions, when executed by the processor, further cause the processing unit to at least:

compute a first variance in the specific ventilation value, a second variance in the perfusion value, and a third variance in the ventilation-perfusion ratio; and

generating a correlation coefficient from one or more of the first variance, the second variance, and the third variance.

15. The magnetic resonance imaging system of claim 11, wherein the instructions, when executed by the processor, further cause the processing unit to at least:

produce an image that shows the spatial distribution of ventilation, perfusion, and ventilation-perfusion for each voxel associated with the acquired MR image.

16. The magnetic resonance imaging system of claim 11, wherein the processing unit uses the spatial distribution of ventilation, perfusion, and ventilation-perfusion as an indicator of a diseased lung.

17. The magnetic resonance imaging system of claim 11, wherein the instructions, when executed by the processor, further cause the processing unit to at least:

determine ventilation-perfusion matching in the lung from at least the ventilation-perfusion ratio value.

18. The magnetic resonance imaging system of claim 17, wherein the processing unit uses ventilation-perfusion matching as an indicator of a diseased lung.

19. A computer program product comprising a non-transitory computer-readable storage medium having instructions stored thereon, the instructions comprising:

code for acquiring, using a magnetic resonance imaging (MRI) machine, an MR image of a lung that includes MR data in a voxel of the MR image;

code for determining a water density value and a perfusion value in at least one voxel of the MR image based on the MR data, wherein the water density value is used to determine an air content value;

code for measuring a breathing frequency value;

code for determining a specific ventilation value in the at least one voxel based on a time delay between an onset of a stimulus to the lung and a response detected in the at least one voxel;

code for determining a ventilation-perfusion ratio value that is the product of the specific ventilation value, the air content value, the inverse of the perfusion value, and the breathing frequency value; and

code for evaluating a spatial distribution of the specific ventilation value, perfusion value, and the ventilation-perfusion ratio by computing a standard deviation of the specific ventilation value and perfusion value,

wherein the specific ventilation value is computed by SVI*(1−

density)*f_B wherein SVI represents quantitative data from the MR image, (1−

density) represents the air content value, and f_B represents the breathing frequency value, thereby determining when the lung is diseased.

20. The non-transitory computer-readable storage medium of claim 19, wherein the code for acquiring the MRI image of the lung includes code that uses oxygen as an inhaled contrast agent.

21. The non-transitory computer-readable storage medium of claim 19, further comprising code for determining an alveolar ventilation value that is the product of the specific ventilation value, the air content value, and the breathing frequency.

22. The non-transitory computer-readable storage medium of claim 21, wherein code for determining the alveolar ventilation value includes code for computing a first variance in the specific ventilation value, a second variance in the perfusion value, and a third variance in the ventilation-perfusion ratio;

and

generating a correlation coefficient from one or more of the first variance, the second variance, and the third variance.

23. The non-transitory computer-readable storage medium of claim 19, wherein code for determining the alveolar ventilation value further comprises code for producing an image that shows the spatial distribution of ventilation, perfusion, and ventilation-perfusion for each voxel associated with the acquired MR image.

24. The non-transitory computer-readable storage medium of claim 19, wherein code for determining the alveolar ventilation value further comprises code for using the spatial distribution of ventilation, perfusion, and ventilation-perfusion to indicate a diseased lung.

25. The non-transitory computer-readable storage medium of claim 19, further comprising code for determining ventilation-perfusion matching in the lung from at least the ventilation-perfusion ratio value.

26. The non-transitory computer-readable storage medium of claim 25, wherein code for determining ventilation-perfusion matching further comprises code for using ventilation-perfusion matching to indicate a diseased lung.