Method and system for identifying optimal image within a series of images that depict a moving organ

• US 7,587,074 B2
• Filed: 07/14/2004
• Issued: 09/08/2009
• Est. Priority Date: 07/21/2003
• Status: Expired due to Fees

• Current Assignee: Medtronic Incorporated (Medtronic Plc)
• Original Assignee: Paieon Inc.
• Inventors: Zarkh, Michael, Klaiman, Moshe
• Primary Examiner(s): Bali; Vikkram
• Assistant Examiner(s): Allison; Andrae S
• Application Number:
  US10/565,224
• Publication Number:
  US 20060188135A1
• Time in Patent Office: 1,882 Days
• Field of Search: 382128-132, 600420-427, 600/431, 378/69, 128/922
• US Class Current: 382/128
• CPC Class Codes: A61B 5/7285 for synchronising or triggering a physiological measurement or image acquisition with a physiological event or waveform, e.g. an ECG signal A61B6/5205, A61B6/541, A61B8/5207, A61B8/543 take precedence
• CPC Class Codes: A61B 6/503 for diagnosis of the heart
• CPC Class Codes: A61B 6/504 for diagnosis of blood vessels, e.g. by angiography
• CPC Class Codes: A61B 6/541 involving acquisition triggered by a physiological signal
• CPC Class Codes: G06T 11/005 Specific pre-processing for tomographic reconstruction, e.g. calibration, source positioning, rebinning, scatter correction, retrospective gating
• CPC Class Codes: G06T 2211/404 Angiography
• CPC Class Codes: G06T 2211/412 Dynamic
• 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 and system for quantifying a cyclic motion within a series of images depicting a moving object subject to composite motion containing a cyclic component and a non-cyclic component of lower frequency than the cyclic component. Composite motion is computed as well as the non-cyclic component as the integral of motion over a motion cycle. The non-cyclic component is subtracted from the composite motion so as to obtain the cyclic component

Claims

1. A method for obtaining a cyclic motion within a series of images depicting a moving object subject to composite motion containing a cyclic motion component having a cyclic period and a non-cyclic consistent component of a lower frequency than the cyclic motion component, the method comprising:

a computer performing the steps of;

(a) computing the composite motion between at least one pair of successive images, the composite motion represented by at least one vector;

(b) computing the non-cyclic consistent component as an integral of the composite motion over the cyclic period;

(c) computing a proportional part of the non cyclic consistent component for each of the at least one pair of successive images; and

(d) for each of the at least one pair of successive images, subtracting the proportional part of the non-cyclic consistent component from the composite motion so as to obtain the cyclic motion component,

wherein the series of images comprises an at least one series of N images acquired during a cyclic period, each frame having an index i within the cyclic period, i=1 . . . N, and wherein the proportional part of the non cyclic component for each of the at least one pair of successive images i-1 and i is determined by dividing the non cyclic component by N and multiplying by i-1.

2. The method according to claim 1, wherein the cyclic period of the cyclic motion component is computed using spectral analysis.

3. The method according to claim 1, wherein the composite motion is determined by optical flow.

4. The method according to claim 1, wherein the composite motion is determined using phase correlation of said images.

5. The method according to claim 1, where cyclic motion values are used for evaluating performance of a body organ.

6. The method according to claim 4, when used in a cardiac application to evaluate heart performance.

7. The method according to claim 6, when used for Ejection Fraction analysis.

8. The method according to claim 6, when used for Left Ventricular analysis.

9. The method according to claim 6, when used for Wall Motion analysis.

10. A method for identifying an image depicting an event associated with cyclic motion, the method comprising:

(a) computing the cyclic motion according to the method of claim 1;

(b) using a graphical representation of the cyclic motion to identify all images matching said event; and

(c) selecting one of said images.

11. The method according to claim 10, wherein the selected image is closest to a predetermined approximation.

12. The method according to claim 10, wherein the event is least motion.

13. The method according to claim 12, for selecting angiographic images to participate in three-dimensional reconstruction of coronary vessels.

14. The method according to claim 13, including deriving cycle period and approximation for least-motion image from an analysis of an electro cardiogram (ECG) signal.

15. The method according to claim 13, including distinguishing the end diastole instance from the end-systole instance by the state of coronary vessel—

maximal spreading versus minimal spreading, respectively.

16. The method according to any one of claim 1 when used for selecting optimal image or images for Quantitative Coronary (QCA) analysis.

17. The method according to any one of claim 1 when used for selecting optimal image or images for Intra Vascular Ultra Sound (IVUS) analysis.

18. The method according to any one of claim 1 when used for selecting optimal image or images for Left Ventricular LVA) analysis.

19. The method according to claim 1 when used for selecting optimal image or images for Wall Motion analysis.

20. The method according to any one of claim 1 when used for Computerized Tomography (CT)reconstruction.

21. The method according to any one of claim 1 when used for Magnetic Resonance Imaging (MRI) reconstruction.

22. The method according to any one of claim 1 when used for Positron Emission Tomography (PET) reconstruction.

23. A system for obtaining a cyclic motion within a series of images depicting a moving object subject to composite motion containing a cyclic motion component having a cyclic period and a non-cyclic consistent component of a lower frequency than the cyclic motion component, the system comprising:

a composite motion unit for computing the composite motion between at least one pair of successive images, the composite motion represented by at least one vector;

a non-cyclic motion unit for computing the non-cyclic consistent component as an integral of the composite motion over the cyclic period;

a proportional part unit for computing a proportional part of the non cyclic motion component for each of the at least one pair of successive images; and

a subtraction unit for subtracting the proportional part of the non-cyclic consistent component from the composite motion occurring between each of the at least one pair of successive images, so as to obtain the cyclic motion component,

wherein the series of images comprises an at least one series of N images acquired during a cyclic period, each frame having an index i within the cyclic period, i=1 . . . N, and wherein the proportional part of the non cyclic component for each of the at least one pair of successive images i-1 and i is determined by dividing the non cyclic component by N and multiplying by i-1 , wherein at least one of the above units is implemented in a hardware or hardware software combination.

24. A system for identifying an image depicting an event associated with cyclic motion, the system comprising:

the cyclic motion system of claim 23 for computing the cyclic motion and deriving data representative of a graphical representation thereof,

an image identification unit responsive to said data representative of a graphical representation of the cyclic motion for identifying all images matching said event, and

an image selection unit for selecting one of said images.

25. The system according to claim 24, wherein the image identification unit is adapted to identify minimal cyclic motion.

26. The system according to claim 25, wherein the image selection unit is adapted to select angiographic images to participate in three-dimensional reconstruction of coronary vessels.

27. The system according to claim 26, including an Electro Cardiogram (ECG) analyzer for deriving cycle period and approximation for least-motion image from an analysis of an ECG signal.

28. The system according to claim 26, including an image processing unit coupled to the image selection unit for distinguishing the end-diastole instance from the end-systole instance by the state of coronary vessel—

maximal spreading versus minimal spreading, respectively.