Method and apparatus for quantification of cardiac wall motion asynchrony

• US 7,974,694 B2
• Filed: 02/08/2006
• Issued: 07/05/2011
• Est. Priority Date: 07/19/2002
• Status: Active Grant

• Current Assignee: Cardiac Pacemakers Incorporated (Boston Scientific Corporation)
• Original Assignee: Cardiac Pacemakers Incorporated (Boston Scientific Corporation)
• Inventors: Huvelle, Etienne, Kramer, Andrew P.
• Primary Examiner(s): Patel; Niketa I
• Assistant Examiner(s): ALTER, ALYSSA MARGO
• Application Number:
  US11/275,978
• Publication Number:
  US 20060116731A1
• Time in Patent Office: 1,973 Days
• Field of Search: 607/24, 607/17, 600/527, 600/508, 600/439, 128/922, 128/923, 128/928
• US Class Current: 607/17
• CPC Class Codes: A61B 8/08 Detecting organic movements or changes, e.g. tumours, cysts, swellings A61B8/02 - A61B8/06 take precedence
• CPC Class Codes: A61N 1/36514 controlled by a physiological quantity other than heart potential, e.g. blood pressure controlled by two or more physical parameters A61N1/36585
• CPC Class Codes: A61N 1/36528 the parameter being measured by means of ultrasound

Abstract

A phase analysis technique provides for quantification of regional wall motion asynchrony from endocardial border contours generated from two-dimensional echocardiographic ventricular images. The technique produces results including a degree of radial ventricular asynchrony in heart failure patients with ventricular conduction delay to predict a magnitude of contractile function improvement with pacing therapy. Quantification of change in ventricular regional wall motion asynchrony in response to a therapy provides for a means to identify candidates to receive the therapy and quantitatively predict the benefit of the therapy. Quantification of changes in ventricular regional wall motion asynchrony in response to a sequence of therapies provides for a means to determine an approximately optimal therapy for an intended patient response.

Claims

1. A system for use with an echocardiographic sensor and an echocardiographic machine adapted to record an echocardiogram, the system comprising:

a controller including;

a therapy protocol synthesizer adapted to generate a sequence of therapy parameters; and

an automatic therapy protocol execution module adapted to time delivery of a sequence of therapies each using at least one parameter of the sequence of therapy parameters;

a therapy circuit, coupled to the controller, to deliver the sequence of therapies; and

a cardiac wall motion evaluation module adapted to receive the echocardiogram recorded by the echocardiographic machine during the delivery of the sequence of therapies and calculate parameters representative of degrees of cardiac wall motion asynchrony each associated with one therapy of the sequence of therapies based on a portion of the echocardiogram recorded during the delivery of the one therapy of the sequence of therapies;

wherein the cardiac wall motion evaluation module further comprises a therapy decision module coupled to the phase computer, the therapy decision module adapted to determine one or more approximately optimal pacing parameters based on the relative phase.

2. The system of claim 1, wherein the cardiac wall motion evaluation module comprises:

a curve generator adapted to generate first and second regional displacement curves based on magnitudes of cardiac regional wall displacement calculated over a cardiac cycle length interval for first and second cardiac regions, respectively, based on the portion of the echocardiogram recorded during the delivery of the one therapy of the sequence of therapies;

a curve offset module, coupled to the curve generator, to provide a common magnitude reference for the first and second wall displacement curves;

a curve averager, coupled to the curve offset module, to average each of the first and second regional displacement curves over a predetermined number of cardiac cycles;

a curve smoothing module, coupled to the curve averager, to smooth the averaged first and second regional displacement curves; and

a phase computer coupled to the curve smoothing module, the phase computer adapted to computer a relative phase representing a phase relationship between the first and second regional displacement curves based on a frequency analysis, the relative phase indicative of the degree of cardiac wall motion asynchrony associated with the one therapy of the sequence of pacing therapies.

3. The system of claim 2, wherein the cardiac wall motion evaluation module further comprises a therapy decision module coupled to the phase computer, the therapy decision module adapted to determine one or more approximately optimal pacing parameters based on the relative phase.

4. The system of claim 2, wherein the phase computer comprises:

a regional displacement phase calculator coupled to the curve smoothing module, the regional displacement phase calculator adapted to calculate first and second regional displacement phase angles based on Fourier transforms of the first and second regional displacement curves, respectively; and

a relative phase calculator, coupled to the regional displacement phase calculator, to calculate the relative phase being a difference between the first and second phase angles.

5. The system of claim 4, wherein the regional displacement phase calculator includes a septal wall displacement phase angle calculator to calculate a septal phase angle Φ

_S and a lateral wall displacement phase angle calculator to calculator a lateral phase angle Φ

_L.

6. The system of claim 5, wherein the relative phase calculator includes a difference calculator to calculate a phase difference:

Φ

_LS=Φ

_L−

Φ

_S.

7. The system of claim 1, wherein the therapy circuit comprises a cardiac pacing circuit, and the therapy protocol synthesizer comprises an atrioventricular delay (AVD) calculator to calculate a sequence of AVDs each being used in one therapy of the sequence of therapies.

8. The system of claim 7, further comprising an AVD selector to select an approximately optimal AVD from the sequence of AVDs based on the degrees of cardiac wall motion asynchrony.

9. The system of claim 8, further comprising an AVD selector to select an AVD associated with an approximately maximum ventricular resynchronization from the sequence of AVDs.

10. The system of claim 9, wherein the AVD selector is adapted to select an AVD associated with a minimum |Φ

_LS|.

11. The system of claim 8, further comprising an AVD selector to select an AVD associated with an approximately optimal unloading of stress on the lateral wall from the sequence of AVDs.

12. The system of claim 11, wherein the AVD selector is adapted to select an AVD associated with a negative Φ

_LS.

13. The system of claim 1, wherein the therapy circuit comprises a cardiac pacing circuit, and the therapy protocol synthesizer is adapted to generate a sequence of pacing sites each being used in one therapy of the sequence of therapies, the pacing sites each defined by parameters specifying one or more sites and timing of pacing pulse delivery with respect to each of the one or more sites.

14. The system of claim 13, further comprising a pacing site selector to select an approximately optimal pacing site or pacing site combination from the sequence of pacing sites based on the degrees of cardiac wall motion asynchrony.

15. The system of claim 14, further comprising a pacing site selector to select a pacing site or pacing site combination associated with an approximately maximum ventricular resynchronization from the sequence of pacing sites.

16. The system of claim 15, wherein the pacing site selector is adapted to select a pacing site or pacing site combination associated with a minimum |Φ

_LS|.

17. The system of claim 14, further comprising a pacing site selector to select a pacing site or pacing site combination associated with an approximately optimal unloading of stress on a ventricular wall from the sequence of pacing sites.

18. The system of claim 17, wherein the pacing site selector is adapted to select a pacing site or pacing site combination associated with a negative Φ

_LS.

19. A method, comprising:

delivering a sequence of therapies;

recording an echocardiogram during the delivery of the sequence of therapies;

computing degrees of ventricular asynchrony each associated with one therapy of the sequence of therapies based on a portion of the echocardiogram recorded during the delivery of the one therapy of the sequence of therapies; and

determining an approximately optimal therapy based on the computed degrees of ventricular asynchrony.

20. The method of claim 19, wherein computing degrees of ventricular asynchrony each associated with one therapy of the sequence of therapies comprises:

calculating regional displacement magnitudes for septal and lateral regions over one cardiac cycle based on cardiac wall motion contours detected from the portion of the echocardiogram recorded during the delivery of the one therapy of the sequence of therapies;

generating septal and lateral displacement curves based on the regional displacement magnitudes calculated for the septal and lateral regions, respectively, the septal and lateral displacement curves each being a regional displacement magnitude verses time plot over the cardiac cycle;

providing a common magnitude reference to the septal and lateral displacement curves;

averaging each of the septal and lateral displacement curves over a predetermined number of cardiac cycles;

smoothing the averaged septal and lateral displacement curves; and

computing a relative phase Φ

_LS representing a phase relationship between the septal and lateral displacement curves based on a frequency analysis, the relative phase Φ

_LS indicative of one of the degrees of ventricular asynchrony associated with the one therapy of the sequence of therapies.

21. The method of claim 20, wherein computing the relative phase Φ

_LS comprises;

computing a septal displacement phase angle Φ

_S and a lateral displacement phase angle Φ

_L, based on Fourier transforms of the averaged septal and lateral displacement curves, respectively; and

computing the relative phase;

Φ

_LS=Φ

_L−

Φ

_S.

22. The method of claim 20, wherein delivering the sequence of therapies comprises delivering a sequence of cardiac pacing therapies each using one of a sequence of predetermined atrioventricular delays (AVDs).

23. The method of claim 22, wherein determining the approximately optimal therapy comprises determining an approximately optimal AVD based on the sequence of predetermined AVDs.

24. The method of claim 23, wherein selecting the approximately optimal AVD comprises selecting an AVD associated with an approximately maximum ventricular resynchronization from the sequence of predetermined AVDs.

25. The method of claim 24, wherein selecting the approximately optimal AVD comprises selecting an AVD associated with a minimum |Φ

_LS|.

26. The method of claim 23, wherein selecting the approximately optimal AVD comprises selecting an AVD associated with an approximately optimal unloading of stress on the lateral wall from the sequence of predetermined AVDs.

27. The method of claim 26, wherein selecting the approximately optimal AVD comprises selecting an AVD associated with a negative Φ

_LS.

28. The method of claim 20, wherein delivering the sequence of therapies comprises delivering a sequence of cardiac pacing therapies each using one of a sequence of predetermined pacing sites or pacing site combinations, the pacing sites or pacing site combinations each defined by parameters specifying one or more sites and timing of pacing pulse delivery with respect to each of the one or more sites.

29. The method of claim 28, wherein determining an approximately optimal therapy comprises determining an approximately optimal pacing site or pacing site combination based on the sequence of predetermined pacing sites or pacing site combinations.

30. The method of claim 29, wherein determining the approximately optimal pacing site or pacing site combination includes, if an approximately optimal pacing site combination is determined, determining a relative order of pacing pulse delivery to each site of the pacing site combination and an inter-site pacing delay between sites of the pacing site combination.

31. The method of claim 29, wherein selecting the approximately optimal pacing site or pacing site combination comprises selecting a pacing site or pacing site combination associated with an approximately maximum ventricular resynchronization from the sequence of predetermined pacing sites or pacing site combinations.

32. The method of claim 31, wherein selecting the approximately optimal pacing site or pacing site combination comprises selecting a pacing site or pacing site combination associated with a minimum |Φ

_LS|.

33. The method of claim 29, wherein selecting the approximately optimal pacing site or pacing site combination comprises selecting a pacing site or pacing site combination associated with an approximately optimal unloading of stress on a ventricular wall from the sequence of predetermined pacing sites or pacing site combinations.

34. The method of claim 33, wherein selecting the approximately optimal AVD comprises selecting an pacing site or pacing site combination associated with a negative Φ

_LS.