Quantitative perfusion analysis

• US 9,406,146 B2
• Filed: 06/23/2010
• Issued: 08/02/2016
• Est. Priority Date: 06/30/2009
• Status: Active Grant

• Current Assignee: Koninklijke Philips Electronics N.V. (Koninklijke Philips N.V.)
• Original Assignee: Koninklijke Philips N.V.
• Inventors: Wiemker, Rafael, Buelow, Thomas
• Primary Examiner(s): Koziol, Stephen R
• Assistant Examiner(s): CHOI, TIMOTHY WING HO
• Application Number:
  US13/380,610
• Publication Number:
  US 20120093390A1
• Time in Patent Office: 2,232 Days
• US Class Current: 1/1
• CPC Class Codes: G06T 2207/10096 Dynamic contrast-enhanced magnetic resonance imaging [DCE-MRI]
• CPC Class Codes: G06T 2207/20104 Interactive definition of region of interest [ROI]
• CPC Class Codes: G06T 2207/30104 Vascular flow; Blood flow; Perfusion
• CPC Class Codes: G06T 7/0016 involving temporal comparison
• CPC Class Codes: G06T 7/20 Analysis of motion motion estimation for coding, decoding, compressing or decompressing digital video signals H04N19/43, H04N19/51
• CPC Class Codes: G06T 7/41 based on statistical description of texture

Abstract

A system is disclosed for quantitative analysis of perfusion images comprising image elements having intensity values associated therewith. The system comprises a frequency distribution computing subsystem (1) for computing a plurality of frequency distributions of the intensity values of at least part of the images. The system comprises a perfusion information extractor (2) for extracting information relating to perfusion from the plurality of frequency distributions. The perfusion information extractor (2) comprises a shift detector (3) for detecting a shift of the intensity values of the frequency distribution. The perfusion information extractor (2) is arranged for extracting the information relating to perfusion, based on the detected shift. A user interface element (8) enables a user to indicate a boundary between the core region and the rim region by a single degree of freedom. A vesselness subsystem (9) associates a vesselness value with an image element.

Claims

1. A system for quantitative analysis of perfusion images comprising image elements having intensity values associated therewith, comprising

a frequency distribution computing subsystem including a processor configured to compute a plurality of frequency distributions of the intensity values of at least part of the perfusion images, wherein the perfusion images include a plurality of images of a same region of interest at different time points over a predetermined period of time and include voxels or pixels with intensity values corresponding to a concentration of contrast agent uptake in tissue of the region of interest over the predetermined period of time;

and

a perfusion information extractor including the processor further configured to compare the plurality of frequency distributions with a single reference frequency distribution to detect respective shifts of respective frequency distributions relative to the reference frequency and extract information relating to perfusion from the plurality of frequency distributions based on the shifts, wherein the information relating to perfusion includes quantitative perfusion information about the contrast agent uptake.

2. The system according to claim 1, wherein the processor compares at least two frequency distributions computed for different images for detecting a shift of the intensity values of the frequency distributions extract the information relating to perfusion, based on the detected shift.

3. The system according to claim 2, wherein the processor comprises a cross-correlator for cross-correlating the at least two frequency distributions to obtain a cross correlation curve, and the processor detects the shift, based on the cross-correlation curve.

4. The system according to claim 3, the wherein the processor further comprises a peak detector for detecting a peak in the cross-correlation curve.

5. The system according to claim 2, the perfusion information extractor comprising a shift curve generator for combining the shifts detected for respective frequency distributions to form a shift curve.

6. The system according to claim 1, further comprising a region-identifying subsystem including the processor further configured to identify a region of interest in the images, the frequency distribution computing subsystem being arranged for computing the frequency distributions for the region of interest in the images.

7. The system according to claim 6, the region-identifying subsystem being arranged for identifying at least part of a rim region surrounding a core region of a lesion as the region of interest.

8. The system according to claim 7, the region-identifying subsystem comprising a user interface element for enabling a user to indicate a boundary between the core region and the rim region by a single degree of freedom, the single degree of freedom being related to a size of the core region.

9. The system according to claim 1, further comprising a vesselness subsystem including the processor further configured to associate a vesselness value with an image element, the vesselness value representing a likelihood that the image element belongs to a vessel, the frequency distribution computing subsystem being arranged for weighting a contribution of an intensity value of the image element to the frequency distribution, based on the vesselness value.

10. The system according to claim 1, the information relating to perfusion comprising at least one of:

an arrival time of a contrast agent, an amount of contrast agent flowing through a region, a maximum enhancement, a time of maximum enhancement, a maximum upslope or intensity-value time derivative, a time of maximum upslope, a fitting of an analytical perfusion curve with free parameters to the observed time curve.

11. The system according to claim 1, wherein the perfusion images are a result of a dynamic scan performed before and after administering a contrast agent.

12. The system according to claim 1, further comprising:

acquiring data through a time series of scans before and after administering a contrast agent; and

generating the perfusion images by reconstructing the acquired data.

13. The system according to claim 11, wherein the perfusion images are successive contrast images of the dynamic scan.

14. The system according to claim 2, where the perfusion information extractor detects contrast agent inflow by detecting a shift in the intensity values of the frequency distributions of two compared frequency distributions of the frequency distributions computed for the different images.

15. The system according to claim 2, where the perfusion information extractor detects a volume of contrast agent based on an amount of the shift in the intensity values and further based on one of a time duration between a shift corresponding to a decrease of the intensity values and a corresponding subsequent increase of the intensity values or an increase of the intensity values followed by a decrease of the intensity values.

16. The system according to claim 1, further comprising:

computing a spatial core-to-rim-curve of the tissue that indicates contrast agent uptake with respect to a native scan as a function of a distance from a core region.

17. The system according to claim 16, further comprising:

partitioning the tissue into the core region and a rim region, and computing the spatial core-to-rim-curve based on the partitioning.

18. A system for quantitative analysis of perfusion images comprising image elements having intensity values associated therewith, comprising

a frequency distribution computing subsystem including a processor configured to compute a plurality of frequency distributions of the intensity values of at least part of the perfusion images, wherein the perfusion images include a plurality of images of a same region of interest at different time points over a predetermined period of time and include voxels or pixels with intensity values corresponding to a concentration of contrast agent uptake in tissue of the region of interest over the predetermined period of time;

and

a perfusion information extractor including the processor further configured to compare at least two frequency distributions computed for different images for detecting a shift of the intensity values of the frequency distribution, extract information relating to perfusion from the plurality of frequency distributions based on the detected shift, wherein the information relating to perfusion includes quantitative perfusion information about the contrast agent uptake, and combine shifts detected for respective frequency distributions to form a shift curve.

19. A system for quantitative analysis of perfusion images comprising image elements having intensity values associated therewith, comprising

a frequency distribution computing subsystem including a processor configured to compute a plurality of frequency distributions of the intensity values of at least part of the perfusion images, wherein the perfusion images include a plurality of images of a same region of interest at different time points over a predetermined period of time and include voxels or pixels with intensity values corresponding to a concentration of contrast agent uptake in tissue of the region of interest over the predetermined period of time;

a perfusion information extractor including the processor further configured to extract information relating to perfusion from the plurality of frequency distributions, wherein the information relating to perfusion includes quantitative perfusion information about the contrast agent uptake; and

a vesselness subsystem including the processor further configured to associate a vesselness value with an image element, the vesselness value representing a likelihood that the image element belongs to a vessel, the frequency distribution computing subsystem being arranged for weighting a contribution of an intensity value of the image element to the frequency distribution, based on the vesselness value.