System and method for image based physiological monitoring of cardiovascular function

• US 8,577,441 B2
• Filed: 12/07/2006
• Issued: 11/05/2013
• Est. Priority Date: 12/08/2005
• Status: Expired due to Fees

• Current Assignee: Siemens Medical Solutions USA Incorporated (Siemens AG)
• Original Assignee: Siemens Medical Solutions USA Incorporated (Siemens AG)
• Inventors: Lorenz, Christine H, Maier, Corinna
• Primary Examiner(s): Cheng, Jacqueline
• Application Number:
  US11/635,101
• Publication Number:
  US 20070135705A1
• Time in Patent Office: 2,525 Days
• Field of Search: 600/410
• US Class Current: 600/410
• CPC Class Codes: A61B 5/055 involving electronic [EMR] or nuclear [NMR] magnetic resonance, e.g. magnetic resonance imaging
• CPC Class Codes: G06T 2207/10088 Magnetic resonance imaging [MRI]
• CPC Class Codes: G06T 2207/30048 Heart; Cardiac
• CPC Class Codes: G06T 7/12 Edge-based segmentation
• CPC Class Codes: G06T 7/149 involving deformable models, e.g. active contour models
• CPC Class Codes: G06T 7/20 Analysis of motion motion estimation for coding, decoding, compressing or decompressing digital video signals H04N19/43, H04N19/51

Abstract

A method for real time monitoring of cardiac function includes defining a plurality of planes through a heart from which to acquire heart imaging data, acquiring a sequence of 2-dimensional images from said selected planes, wherein each said image comprises a plurality of intensities defined on a domain of points on a 2D grid, selecting one or more of said sequence of images to provide imaging data for monitoring heart function, wherein said remaining unselected images are adapted for diagnostic purposes, and repeating said steps of acquiring a sequence of 2-dimensional images from said selected planes and selecting one or more of said sequence of images for monitoring heart function, wherein a time sequence of images is obtained, alternating between images for monitoring heart function and images for diagnostic purposes.

Claims

1. A computer implemented method for real time monitoring of cardiac function, the method performed by the computer comprising the steps of:

selecting a plane through a heart from which to acquire heart function data;

locating a center of a chamber of said heart;

acquiring a sequence of 2D magnetic resonance (MR) images from said plane over one or more cardiac cycles with a time resolution at least one order of magnitude less than the cardiac rate, wherein each of said images comprises a plurality of intensities defined on a domain of points on a 2D grid;

projecting a plurality of beams from said chamber center through each MR image as each MR image is acquired, each beam comprising a substantially collinear subset of points in said image;

sampling image intensities along periodically spaced points of each beam in each MR image;

obtaining image intensity in real-time as a function of time at corresponding beam points across the acquired MR images, wherein a time resolution of the image intensity function is at least that of the sequence of the 2D MR images;

detecting endocardial contours from said image intensity functions sampled at periodically spaced points of each of the plurality of beams in each MR image;

calculating distances between said detected endocardial contours, wherein the maximum distances as a function of time in each projection beam are taken as the end-diastolic (ED) diameter, and the minimum distances are taken as the end-systolic (ES) diameter; and

using said image intensity time function along each said beam to detect changes in values of a physiological parameter over time.

2. The method of claim 1, further comprising bilinearly interpolating additional intensities between sampled points along each said beam.

3. The method of claim 1, wherein said endocardial contour is detected using a horizontal 1D Canny filter.

4. The method of claim 1, wherein said endocardial contour is detected using an active contour model that minimizes an energy that is a weighted sum of external constraint forces determined by said image and internal constraint forces that depend on the shape of the contour.

5. The method of claim 1, further comprising calculating a chamber area from the area enclosed by the points of the endocardial contour, and calculating an ejection fraction from EDV-ESV/EDV, wherein EDV is the end-diastolic chamber area and ESV is the end-systolic chamber area.

6. The method of claim 5, wherein said chamber is the left ventricle.

7. The method of claim 1, further comprising detecting epicardial contours from said plurality of beams, and calculating a wall thickness as a function of time from said endocardial contours and said epicardial contours.

8. The method of claim 1, wherein detecting value changes comprises selecting a clinically relevant change threshold, maintaining a running average of said parameter values over a predetermined number of previous heartbeats, and comparing a parameter value at each heartbeat to either the running average value or a baseline value for said parameter.

9. The method of claim 1, wherein said acquired 2D MR images are part of a sequence of images that include images acquired for diagnostic purposes.

10. A computer implemented method for real time monitoring of cardiac function during an image guided interventional procedure, the method performed by the computer comprising the steps of:

defining a plurality of mutually orthogonal planes through a heart from which to acquire heart imaging data for acquiring a time series of 2-dimensional images from said selected planes, wherein each said image comprises a plurality of intensities defined on a domain of points on a 2D grid;

selecting a subset of said time series of images to provide imaging data for monitoring heart function, wherein said selected subset of images are to be acquired from a same plane, wherein said remaining unselected images are acquired from a different plane and are adapted for diagnostic purposes;

acquiring said time series of images over one or more heart cycles, wherein an image is acquired from each of the plurality of planes at each acquisition time, alternating between acquiring images for monitoring heart function and acquiring images for diagnostic purposes, wherein every n^th image acquired during each heart cycle is acquired for monitoring, wherein n is greater than or equal to 2, wherein image intensities are extracted from corresponding points across the acquired monitoring images as each monitoring image is acquired to define said image intensities as a function of time; and

using said image intensity time functions along each said beam to detect changes in values of a physiological parameter over time.

11. The method of claim 10, further comprising, for each monitoring image,

selecting a center of a chamber of said heart,

projecting a plurality of beams into said image through said chamber center, each beam comprising a substantially collinear subset of points in said image, and

sampling image intensities along each point of said beam, to obtain a series of time measurements wherein an image intensity along each said beam is plotted as a function of time.

12. The method of claim 11, further comprising, analyzing said image-intensity time behavior over a pre-determined time period after intravenous administration of a contrast agent.

13. The method of claim 10, wherein said images are magnetic resonance (MR) phase images, and

further comprising plotting a MR phase encoded for velocity as a function of time, and deriving performance parameters from said phase-time plot, including one or more of time to peak flow, percentage regurgitant flow, and peak flow velocity.

14. The method of claim 13, where each said phase image is a 1D image created by using a pencil beam 2D radio-frequency excitation to derive measurements from a column of tissue.

15. The method of claim 13, where each said phase image is a 2D image created from a velocity encoded MR image.

16. A program storage device readable by a computer, tangibly embodying a non-transitory program of instructions executable by the computer to perform the method steps for real time monitoring of cardiac function, said method comprising the steps of:

selecting a plane through a heart from which to acquire heart function data;

locating a center of a chamber of said heart;

acquiring a sequence of 2D magnetic resonance (MR) images from said plane over one or more cardiac cycles with a time resolution at least one order of magnitude less than the cardiac rate, wherein each of said images comprises a plurality of intensities defined on a domain of points on a 2D grid;

projecting a plurality of beams from said chamber center through each MR image as each MR image is acquired, each beam comprising a substantially collinear subset of points in said image;

sampling image intensities along periodically spaced points of each beam in each MR image;

obtaining image intensity in real-time as a function of time at corresponding beam points across the acquired MR images, wherein a time resolution of the image intensity function is at least that of the sequence of the 2D MR images;

detecting endocardial contours from said image intensity functions sampled at periodically spaced points of each of the plurality of beams in each MR image;

calculating distances between said detected endocardial contours, wherein the maximum distances as a function of time in each projection beam are taken as the end-diastolic (ED) diameter, and the minimum distances are taken as the end-systolic (ES) diameter; and

using said image intensity time function along each said beam to detect changes in values of a physiological parameter over time.

17. The computer readable program storage device of claim 16, the method further comprising bilinearly interpolating additional intensities between sampled points along each said beam.

18. The computer readable program storage device of claim 16, wherein said endocardial contour is detected using a horizontal 1D Canny filter.

19. The computer readable program storage device of claim 16, wherein said endocardial contour is detected using an active contour model that minimizes an energy that is a weighted sum of external constraint forces determined by said image and internal constraint forces that depend on the shape of the contour.

20. The computer readable program storage device of claim 16, the method further comprising calculating a chamber area from the area enclosed by the points of the endocardial contour, and calculating an ejection fraction from EDV-ESV/EDV, wherein EDV is the end-diastolic chamber area and ESV is the end-systolic chamber area.

21. The computer readable program storage device of claim 20, wherein said chamber is the left ventricle.

22. The computer readable program storage device of claim 16, the method further comprising detecting epicardial contours from said plurality of beams, and calculating a wall thickness as a function of time from said endocardial contours and said epicardial contours.

23. The computer readable program storage device of claim 16, wherein detecting value changes comprises selecting a clinically relevant change threshold, maintaining a running average of said parameter values over a predetermined number of previous heartbeats, and comparing a parameter value at each heartbeat to either the running average value or a baseline value for said parameter.

24. The computer readable program storage device of claim 16, wherein said acquired 2D MR images are part of a sequence of images that include images acquired for diagnostic purposes.

25. The computer readable program storage device of claim 16, the method further comprising analyzing said image-intensity time behavior over a pre-determined time period after intravenous administration of a contrast agent.

26. The computer readable program storage device of claim 16, wherein said images are magnetic resonance (MR) phase images, and the method further comprising plotting a MR phase encoded for velocity as a function of time, and deriving performance parameters from said phase-time plot, including one or more of time to peak flow, percentage regurgitant flow, and peak flow velocity.

27. The computer readable program storage device of claim 26, where each said phase image is a 1D image created by using a pencil beam 2D radio-frequency excitation to derive measurements from a column of tissue.

28. The computer readable program storage device of claim 26, where each said phase image is a 2D image created from a velocity encoded MR image.