Methods for identifying neuronal spikes

• US 7,957,793 B2
• Filed: 12/22/2005
• Issued: 06/07/2011
• Est. Priority Date: 12/22/2004
• Status: Active Grant

• Current Assignee: Wisconsin Alumni Research Foundation (University of Wisconsin System)
• Original Assignee: Wisconsin Alumni Research Foundation (University of Wisconsin System)
• Inventors: Huang, He, Gale, John T., Montgomery, Erwin B. Jr.
• Primary Examiner(s): Marmor, II; Charles A
• Assistant Examiner(s): Natnithithadha; Navin
• Application Number:
  US11/315,808
• Publication Number:
  US 20060167369A1
• Time in Patent Office: 1,993 Days
• Field of Search: 600544-547
• US Class Current: 600/544
• CPC Class Codes: A61B 5/24 Detecting, measuring or recording bioelectric or biomagnetic signals of the body or parts thereof
• CPC Class Codes: A61B 5/7217 of noise originating from a therapeutic or surgical apparatus, e.g. from a pacemaker

Abstract

A method for identifying neuronal spikes (extracellular action potentials) is described wherein measured microelectrode readings from tissue are reviewed to identify spikes (successive readings having prolonged rises and/or falls). The frequency of such spikes as a function of their amplitude assumes a bimodal distribution wherein higher amplitude spikes represent neuronal spikes (signal) and lower amplitude spikes represent noise, and thus the higher amplitude spikes can be assumed to be neuronal spikes. Neuronal spikes from the same neuron can then be assumed to have substantially the same waveform shape and period, with the only significant difference between them being the scaling of their amplitudes (i.e., the amplitudes of spikes from the same neuron tend to be proportionate at any given time along their period). Thus, by testing identified neuronal spikes for matching timing and for proportional amplitudes, the neuronal spikes may further be identified as coming from the same or different neurons.

Claims

1. A method of identifying neuronal spikes comprising the steps of:

a. obtaining amplitude readings sampled over time from an electrode situated within tissue;

b. identifying series of successive readings having at least one of

(1) continuously increasing amplitude, and

(2) continuously decreasing amplitude, each such series of readings representing a spike candidate;

c. identifying spikes within the spike candidates, the spikes being chosen from the spike candidates having greatest amplitude; and

d. performing at least one of the steps of;

(1) placing a therapeutic electrode within the tissue, and

(2) delivering electrical stimulation from a therapeutic electrode to the tissue, in dependence on the identified spikes.

2. The method of claim 1 wherein the step of identifying spikes within the spike candidates includes the step of choosing as spikes at least some of the spike candidates having amplitudes greater than the most frequently occurring amplitudes of the spike candidates.

3. The method of claim 2 wherein the spike candidates chosen as spikes have amplitudes at least twice as great as the most frequently occurring amplitudes of the spike candidates.

4. The method of claim 1 wherein each series of successive readings includes at least three readings.

5. The method of claim 1 wherein the step of identifying spikes within the spike candidates includes the step of:

a. identifying the lowest amplitude about which the greatest number of spike candidates occur;

b. setting a spike threshold amplitude greater than this lowest amplitude;

c. identifying as spikes at least some of the spike candidates having amplitudes greater than the spike threshold amplitude.

6. The method of claim 5 further comprising the step of assigning a start time to each identified spike, each spike having a start time corresponding to a zero amplitude for the spike occurring prior to the time when the spike attains the spike threshold amplitude.

7. The method of claim 6 where the start time corresponds to the second zero amplitude occurring prior to the time when the spike attains the spike threshold amplitude.

8. The method of claim 1 wherein the step of identifying spikes within the spike candidates further includes the steps of:

a. identifying the two largest groups of spike candidates wherein;

(1) each group has at least substantially similar amplitude; and

(2) one group has larger amplitude than the other;

b. choosing as spikes at least some of the spike candidates in the group having larger amplitude.

9. The method of claim 1 wherein the step of identifying spikes within the spike candidates further includes the steps of:

a. identifying two groups of spike candidates wherein;

(1) each group has high incidences of spike candidates with at least substantially similar amplitudes, and

(2) one group has a median amplitude at least twice as great as the median amplitude of the other group; and

b. choosing as spikes at least some of the spike candidates in the group having the median amplitude at least twice as great as the median amplitude of the other group.

10. The method of claim 1 further comprising the step of grouping the spikes, each group having:

(1) at least substantially identical time between the amplitude maxima and amplitude minima of the spikes therein, and

(2) proportionate amplitude maxima, and also proportionate amplitude minima, between the spikes therein.

11. The method of claim 10 wherein each group also has at least one of:

(1) amplitude maxima,

(2) amplitude maxima,

(3) the differences between amplitude maxima and minima, and

(4) root mean square amplitude,

matching to a predefined tolerance.

12. The method of claim 10 further comprising the step of generating, for each group of spikes, average values of at least some of the data of the spikes in the group.

13. The method of claim 1 wherein spike candidates are defined by series of successive readings having at least one of

a. continuously increasing amplitude, and

b. continuously decreasing amplitude,

wherein each such series extends over a period between approximately 0.075-0.175 ms.

14. A method of identifying neuronal spikes comprising the steps of:

a. obtaining amplitude readings sampled over time from an electrode situated within tissue;

b. identifying series of readings wherein a positive-amplitude peak is followed by a negative-amplitude valley, each such series of readings representing a spike candidate;

c. determining the amplitude distribution of the spike candidates;

d. identifying the amplitude within the distribution at which the highest proportion of spike candidates are clustered, this amplitude representing a noise peak;

e. selecting as spikes at least some of the spike candidates having amplitudes greater than the amplitude of the noise peak; and

d. performing at least one of the steps of;

(1) placing a therapeutic electrode within the tissue, and

(2) delivering electrical stimulation from a therapeutic electrode to the tissue, in dependence on the spikes.

15. The method of claim 14 wherein the spike candidates selected as spikes have amplitudes at least twice as great as the amplitude of the noise peak.

16. The method of claim 14:

a. further comprising the step of identifying the amplitude within the distribution at which the second highest proportion of spike candidates are clustered, this amplitude representing a signal peak; and

b. selecting as spikes spike candidates situated about the signal peak.

17. The method of claim 14 further comprising the step of assigning a start time to each spike, each spike having a start time corresponding to a zero amplitude for the spike occurring prior to the time of the spike'"'"'s positive-amplitude peak.

18. The method of claim 14 further comprising the step of assigning an end time to each spike, each spike having an end time corresponding to a zero amplitude for the spike occurring after to the time of the spike'"'"'s negative-amplitude peak.

19. The method of claim 14 further comprising the step of grouping the identified spikes, each group having:

a. at least substantially identical time between the amplitude maxima and amplitude minima of the identified spikes therein, and

b. at least one of;

(1) proportionate amplitude maxima between the spikes therein, and

(2) proportionate amplitude minima between the spikes therein.

20. The method of claim 19 further comprising the step, for each group, of calculating an average of at least some of the data contained within the spikes in the group.

21. The method of claim 19 further comprising the steps of:

a. identifying groups having at least a minimum number of spikes therein as containing signal spikes, and

b. identifying groups having less than the minimum number of spikes therein as containing noise spikes.

22. A method of identifying neuronal spikes comprising the steps of:

a. obtaining amplitude readings sampled over time from an electrode situated within tissue, the amplitude readings defining spikes;

b. grouping at least some of the spikes, each group having;

(1) at least substantially identical time between the amplitude maxima and amplitude minima of the spikes therein,

(2) proportionate amplitude maxima between the spikes therein, and

(3) proportionate amplitude minima between the spikes therein; and

d. performing at least one of the steps of;

(1) placing a therapeutic electrode within the tissue, and

(2) delivering electrical stimulation from a therapeutic electrode to the tissue,

in dependence on the grouped spikes.

23. The method of claim 22 further comprising the step of calculating, for each group, an average of at least some of the data of the spikes within the group.

24. The method of claim 22 further comprising the step of identifying each group as being:

a. a signal spike corresponding to a neuronal firing, or

b. a noise spike which does not correspond to a neuronal firing,

at least partially on the basis of the number of spikes within each group.

25. The method of claim 22 wherein the amplitude readings include signal spikes and noise spikes.

26. The method of claim 25 further comprising the step of identifying the signal spikes by:

a. identifying series of amplitude readings wherein a positive-amplitude peak is followed by a negative-amplitude valley, each such series of readings representing a spike candidate;

b. determining the amplitude distribution of the spike candidates;

c. identifying the amplitude within the distribution at which the highest proportion of spike candidates are clustered, this amplitude representing a noise peak;

d. selecting as spikes at least some of the spike candidates having amplitudes greater than the amplitude of the noise peak.

27. The method of claim 25 further comprising the step of identifying the signal spikes by:

a. identifying series of amplitude readings wherein a positive-amplitude peak is followed by a negative-amplitude valley, each such series of readings representing a spike candidate;

b. determining the amplitude distribution of the spike candidates;

c. identifying the two different amplitudes within the distribution at which the highest proportions of spike candidates are clustered, the lower of these amplitudes representing a noise peak and the higher of these amplitudes representing a signal peak;

d. selecting as spikes at least some of the spike candidates having amplitudes about the signal peak.