System for interpolating surface potential values for use in calculating current density

• US 5,247,436 A
• Filed: 09/20/1991
• Issued: 09/21/1993
• Est. Priority Date: 08/14/1987
• Status: Expired due to Term

• Current Assignee: Sterling Trustees, L.L.C.
• Original Assignee: Micro-Tek Incorporated
• Inventors: Stone, Charles L. Jr.
• Primary Examiner(s): Envall, Jr., Roy N.
• Assistant Examiner(s): Bodendorf, Andrew
• Application Number:
  US07/763,862
• Time in Patent Office: 732 Days
• Field of Search: 364/413.05, 128/731
• US Class Current: 600/372
• CPC Class Codes: A61B 5/291 for electroencephalography [EEG]
• CPC Class Codes: A61B 5/375 using biofeedback
• CPC Class Codes: G06F 2218/22 Source localisation; Inverse modelling

Abstract

A system for mapping and/or displaying spatiotemporal features of underlying event-related neural activity involves measuring the event-related evoked potential values at a limited number of points over the region of interest, deriving a first grid of potential values at points over the area of interest representing at least about a five fold increase in the number of points by a Kriging-type of spatial statistical interpolation, and using the first grid to derive for the region of interest a second grid of points of current density values by a Laplacian conversion. A plurality of second grids of current density values spaced apart in time are derived and displayed on a monitor in cartoon-type fashion to provide a cartoon-type display of the features of the event-related evoked response over the region of interest varying with time. Various other displays are possible for highlighting specific features.

Claims

1. A process for deriving the current density at points over an area of interest as a measure of the underlying neural activity comprising:

the step of measuring physiological potentials at a limited number of points over the area of interest by means of electrodes to yield electrical signals indicative thereof;

the step of using the electrical signals to estimate the potential at each of an increased number of points over the area of interest by Kriging interpolation of the measured potentials; and

the step of deriving from said estimated potentials a current density value at each of said points.

2. The process of claim 1 in which the Kriging interpolation is ordinary Kriging.

3. The process of claim 1 in which the Kriging interpolation is ordinary Kriging with a power function variogram model.

4. The process of claim 3 in which the variogram is used with a power function raised to between 1.4 and 1.8.

5. The process of claim 1 in which the Kriging interpolation is universal Kriging utilizing information available as to the underlying neural activity in the area of interest.

6. The process of claim 1 wherein the step of Kriging interpolation and the step of deriving current density values are calculated in a single step.

7. The process of claim 6 wherein the single step is performed by a microprocessor.

8. The process of claim 1 wherein the estimating step and the deriving step is performed by a microprocessor.

9. The process of each of claims 1, 2, 3, 4, and 5 in which the current density is calculated from a Laplacian transformation of said estimated point potentials.

10. The process of each of claims 1, 2, 3, 4 and 5 in which the step of deriving the current density involves a five point difference operator to approximate a second spatial derivative at each point.

11. The process of each of claims 1, 2, 3, 4, and 5 in which said Kriging interpolation generates at least about a five fold increase in the number of points.

12. Apparatus for displaying the features of underlying neural activity over an area of interest comprising:

means for measuring physiological potentials at a plurality of points and for generating electrical signals indicative thereof;

a first electronic means for processing said signals to derive a grid of estimated potential values, corresponding to the area of interest, of increased density by a Kriging interpolation program;

a second means for processing said grid of estimated potential values to derive a corresponding grid of current density values; and

means for displaying in perspective a three dimensional image representing said grid of current density values and displaying thereby the current density features of underlying neural activity over the area of interest.

13. The apparatus of claim 12 in which the electronic means for processing said grid of estimated potential values uses a Laplacian transformation to derive a corresponding grid of current density values.

14. The apparatus of claim 12 in which the electronic means for processing said grid of estimated potential values uses a five point difference operator to derive a corresponding grid of current density values.

15. The apparatus of claim 12 in which the interpolation is ordinary Kriging.

16. The apparatus of claim 12 in which the interpolation is universal Kriging.

17. The apparatus of claim 12 in which the Kriging interpolation is ordinary Kriging with a power function variogram model.

18. The apparatus of claim 17 in which the variogram is used with a power function raised to between 1.4 and 1.8.

51. The apparatus of claim 12, 27 or 44 wherein the said first electronic means for deriving a grid of estimated potential values and the said second electronic means for deriving said grid of current density values comprises a microprocessor.

19. Apparatus for displaying features of underlying neural activity over an area of interest comprising:

an electrode array positioned over the area of interest for measuring potentials at a corresponding number of points over such area and means for generating electrical signals indicative thereof;

signal processing means for providing interpolated values by a Kriging interpolation program of potentials at points intermediate between said corresponding numbers of points for providing an increased number of potential values;

signal processing means for deriving a substantially corresponding number of current density values for the area of interest; and

means for displaying the current density values as a perspective three dimensional display of the current density values over the area of interest.

20. Apparatus in accordance with claim 19 in which said Kriging program utilizes ordinary Kriging.

21. Apparatus in accordance with claim 19 in which said Kriging program is a universal Kriging interpolation.

22. The apparatus of claim 19 in which the Kriging interpolation is ordinary Kriging with a power function variogram.

23. The apparatus in accordance with claim 22 in which the variogram is used with a power function raised to between 1.4 and 1.8.

24. The apparatus of claim 19 in which the Kriging interpolation is universal Kriging.

25. The apparatus of claim 19 in which the signal processing means derives the current density by performing a Laplacian transformation.

26. The apparatus of claim 19 in which the signal processing means derives the current density by using a five point difference operator.

27. Apparatus for displaying features of underlying neural activity at a desired point in an area of interest on a human or animal comprising:

a plurality of electrodes for positioning at points over the area of interest and means for deriving therefrom electrical signals as measures of the physiological potentials at such points;

a first electronic means to derive by a Kriging interpolation program a grid of estimated potential values of an area corresponding to the area of interest with a greater coordinate density than the number of original points;

a second electronic means for processing said grid of estimated potential values to derive a corresponding grid of current density values;

means for sampling the values of a selected point in the grid of estimated potential values and the grid of current density values corresponding to the desired point; and

means for displaying said sampled values of the estimated potential and of the current density at the desired point in successive sampling intervals in a plot against time.

28. Apparatus in accordance with claim 27 in which the plurality comprises fifteen electrodes.

29. The apparatus of claim 27 in which the Kriging interpolation is ordinary Kriging.

30. The apparatus of claim 27 in which the Kriging interpolation is ordinary Kriging with a power function variogram model.

31. The apparatus of claim 30 in which the variogram is used with a power function raised to between 1.4 and 1.8.

32. The apparatus of claim 27 in which the Kriging interpolation is universal Kriging.

33. The apparatus of claim 27 in which the signal processing means derives the current density by performing a Laplacian transformation of said grid points.

34. The apparatus of claim 27 in which the signal processing means derives the current density by using a five point difference operator to approximate a second spatial derivative at each grid point.

35. Apparatus for displaying features of underlying neural activity over an area of interest of a human or animal body comprising:

means for measuring physiological potentials at a plurality of points and for generating electrical signals indicative thereof;

a first signal processing means for deriving by a Kriging interpolation for the area of interest, potential values at a grid of points of a density greater than the measured points;

a second signal processing means for deriving current density values at said grid of points; and

means for displaying the current density values at points in the area of interest.

36. Apparatus in accordance with claim 35 in which the means for displaying current density values at points in the area of interest displays said values at successive sampling intervals whereby there is a dynamic display of the current density values with changing time.

37. The apparatus of claim 35 in which the Kriging interpolation is ordinary Kriging.

38. The apparatus of claim 35 in which the Kriging interpolation is ordinary Kriging with a power function variogram model.

39. The apparatus of claim 38 in which the variogram is used with a power function raised to between 1.4 and 1.8.

40. The apparatus of claim 35 in which the Kriging interpolation is universal Kriging.

41. The apparatus of claim 35 in which the signal processing means derives the current density by performing a Laplacian transformation of said grid points.

42. The apparatus of claim 35 in which the signal processing means derives the current density by using a five point difference operator.

43. The apparatus of claim 35 wherein the first signal processing means and the said second signal processing means comprises a microprocessor.

44. Apparatus for displaying features of underlying neural activity over an area of interest comprising:

means for measuring physiological potentials at a plurality of points and for generating electrical signals indicative thereof;

a first electronic means for processing said signals to derive a grid of estimated potential values, corresponding to the area of interest, of increased density by a Kriging interpolation program;

a second electronic means for processing said grid of estimated potential values to derive a corresponding grid of current density values; and

means for displaying in perspective a three dimensional image representing said grid of estimated potential values and displaying thereby potential features of underlying neural activity over the area of interest.

45. The apparatus of claim 44 in which the Kriging interpolation is ordinary Kriging.

46. The apparatus of claim 44 in which the Kriging interpolation is ordinary Kriging with a power function variogram model.

47. The apparatus of claim 46 in which the variogram is used with a power function raised to between 1.4 and 1.8.

48. The apparatus of claim 44 in which the Kriging interpolation is universal Kriging.

49. The apparatus of claim 44 in which the signal processing means derives the current density by performing a Laplacian transformation of said grid points.

50. The apparatus of claim 44 in which the signal processing means derives the current density by using a five point difference operator.