Method and apparatus for predicting the efficacy of cardioversion

• US 5,730,144 A
• Filed: 07/10/1996
• Issued: 03/24/1998
• Est. Priority Date: 07/10/1996
• Status: Expired due to Fees

• Current Assignee: the united states of america as represented by the secretary of the navy
• Original Assignee: the united states of america as represented by the secretary of the navy
• Inventors: Thomas, James D., Grimm, Richard A., Chandra, Shalabh, Katz, Richard A.
• Primary Examiner(s): Kamm, William E.
• Assistant Examiner(s): LAYNO, CARL HERNANDZ
• Application Number:
  US08/687,099
• Time in Patent Office: 622 Days
• Field of Search: 128/695, 128/713, 607/6
• US Class Current: 600/526
• CPC Class Codes: A61B 8/06 Measuring blood flow measuring volume flow in general G01F, e.g. G01F1/66, G01F1/72; measuring speed of fluids in general G01P5/00
• CPC Class Codes: A61N 1/3925 Monitoring; Protecting
• CPC Class Codes: A61N 1/395 for treating atrial fibrillation

Abstract

An apparatus and method for predicting the efficacy of cardioversion as a dality for reverting a patient with atrial fibrillation to normal sinus rhythm. Blood flow through the patient'"'"'s atrium is measured, converted and processed using nonlinear or chaotic processing to obtain a differential radius signal. The number of excursions of the differential radius beyond a threshold value indicates whether cardioversion will be successful.

Claims

1. A method for determining the likelihood that cardioversion will cause a patient with atrial fibrillation to revert to normal sinus rhythm cardioversion including the steps of:

monitoring the flow of blood through the patient'"'"'s atrium;

generating a velocity signal representing the velocity of blood flow over a diagnostic time interval including a plurality of cardiac cycles;

generating a differential radius signal based upon chaotic and differential processing of the velocity signal over the diagnostic time interval; and

determining the number of instances that the differential radius signal exceeds a threshold value over the diagnostic time interval thereby to indicate whether the patient will revert to normal sinus rhythm by cardioversion.

2. A method as recited in claim 1 wherein said generation of said differential radius signal includes the step of digitizing the velocity signal to obtain a series of time samples at a sampling frequency.

3. A method as recited in claim 2 wherein said generation of the differential radius signal includes generating a time series vector representation for the series of time samples.

4. A method as recited in claim 3 wherein said generation of the differential radius signal includes:

generating an embedding delay value based upon the time series vector representation;

comparing the magnitudes of the terms of the vector representation at a given time and at a time delayed by the embedding delay value to obtain a chaotic radius; and

generating in response to each value of the chaotic radius the differential chaotic radius for the given time.

5. A method as recited in claim 4 wherein said step of determining an embedding delay includes:

generating a vector time delay interval in response to the data in the time series of samples; and

generating the time series vector representation based upon the value of the data in the time samples at intervals corresponding to the vector time delay interval.

6. A method as recited in claim 5 wherein said step of generating the vector time delay interval includes the step of obtaining average mutual information from the time series of samples.

7. A method as recited in claim 5 wherein said generation of the embedding delay value includes the step of obtaining global false nearest neighbor information from the time series vector representation.

8. A method as recited in claim 5 wherein said step of generating the differential chaotic radius includes comparing, at times corresponding to each vector time delay interval, the magnitude of the time series vector representation at that time and the magnitude of the time series vector representation at time offset by the embedding delay value.

9. A method as recited in claim 8 wherein said step of determining the threshold value includes selecting a value from a range of threshold values that define statistically significant data.

10. A method as recited in claim 5 wherein said step of generating the differential chaotic radius includes comparing, at times corresponding to each vector time delay interval, the magnitude of the time series vector representation at that time and the magnitude of the time series vector representation at time offset by the embedding delay value and wherein:

the sampling frequency is in the range of six to ten times the Nyquist sampling frequency for the velocity signal;

the vector time delay is an integer multiple of the sampling interval; and

the embedding function is an integer multiple of the vector time delay.

11. Apparatus for determining the likelihood cardioversion will cause a patient with atrial fibrillation to revert to normal sinus rhythm, said apparatus comprising:

means for monitoring the flow of blood through the patient'"'"'s atrium;

first generating means for generating a velocity signal representing the velocity of blood flow over a diagnostic time interval including a plurality of cardiac cycles;

second generating means for generating a differential radius signal based upon chaotic and differential processing of the velocity signal over the diagnostic time interval; and

threshold evaluation means for determining the number of instances that the differential radius signal exceeds a threshold value over the diagnostic time interval thereby to indicate whether the patient will revert to normal sinus rhythm by cardioversion.

12. Apparatus as recited in claim 11 wherein said second generating means includes means for digitizing the velocity signal to obtain a series of time samples at a sampling frequency.

13. Apparatus as recited in claim 12 wherein said second generating means includes means for generating a time series vector representation for the series of time samples.

14. Apparatus as recited in claim 13 wherein said second generating means includes:

means for generating an embedding delay value based upon the time series vector representation;

means for comparing the magnitudes of the terms of the vector representation at a given time and at a time delayed by the embedding delay value to obtain a chaotic radius; and

means for generating in response to each value of the chaotic radius the differential chaotic radius for the given time.

15. Apparatus as recited in claim 14 wherein said second generator means includes:

means for generating a vector time delay interval in response to the data in the time series of samples; and

time series vector representation generating means responding to the values of the data in the time samples at intervals corresponding to the vector time delay interval.

16. Apparatus as recited in claim 15 wherein said means for generating the vector time delay interval includes means for obtaining average mutual information from the time series of samples.

17. Apparatus as recited in claim 15 wherein said means for generating the embedding delay value includes means for obtaining global false nearest neighbor information from the time series vector representation.

18. Apparatus as recited in claim 15 wherein said second generating means includes means for comparing, at times corresponding to each vector time delay interval, the magnitude of the time series vector representation at that time and the magnitude of the time series vector representation at time offset by the embedding delay value.

19. Apparatus as recited in claim 18 wherein means for said threshold evaluation means includes means for selecting a value from a range of threshold values that define statistically significant data.

20. Apparatus as recited in claim 15 wherein said second generating means includes means for comparing, at times corresponding to each vector time delay interval, the magnitude of the time series vector representation at that time and the magnitude of the time series vector representation at time offset by the embedding delay value and wherein:

said digitizing means operates at a sampling frequency in the range of six to ten times the Nyquist sampling frequency for the velocity signal;

the vector time delay is an integer multiple of the sampling interval; and

the embedding function is an integer multiple of the vector time delay.