U.S. patent application number 13/264142 was filed with the patent office on 2012-09-27 for neurocranial electrostimulation models, systems, devices and methods.
This patent application is currently assigned to Research Foundation of the City University of New York. Invention is credited to Marom Bikson, Abhishek Datta, Jacek Dmochowski, Lucas C. Parra, Yuzhuo Su.
Application Number | 20120245653 13/264142 |
Document ID | / |
Family ID | 42983111 |
Filed Date | 2012-09-27 |
United States Patent
Application |
20120245653 |
Kind Code |
A1 |
Bikson; Marom ; et
al. |
September 27, 2012 |
NEUROCRANIAL ELECTROSTIMULATION MODELS, SYSTEMS, DEVICES AND
METHODS
Abstract
Methods, systems, and apparatus, including computer programs
encoded on a computer storage medium, for developing transcranial
electrical stimulation protocols are disclosed. In one aspect, a
method includes the actions of accepting an image model of target
tissue, obtaining a forward model having a first electrode
configuration and first electrical stimulation parameters based on
electrical stimulation of the target tissue, accepting electrode
configuration changes or electrical stimulation parameter changes
resulting in a second electrode configuration or second electrical
stimulation parameters, determining an optimized tissue model using
a least square methodology and based on the second electrode
configuration or second electrical stimulation parameter changes,
comparing the optimized tissue model with a desired outcome, and
providing a confirmation of the optimized model with the desired
outcome.
Inventors: |
Bikson; Marom; (Brooklyn,
NY) ; Datta; Abhishek; (New York, NY) ; Parra;
Lucas C.; (Brooklyn, NY) ; Dmochowski; Jacek;
(New York, NY) ; Su; Yuzhuo; (Houston,
TX) |
Assignee: |
Research Foundation of the City
University of New York
New York
NY
|
Family ID: |
42983111 |
Appl. No.: |
13/264142 |
Filed: |
April 13, 2010 |
PCT Filed: |
April 13, 2010 |
PCT NO: |
PCT/US2010/030943 |
371 Date: |
June 1, 2012 |
Related U.S. Patent Documents
|
|
|
|
|
|
Application
Number |
Filing Date |
Patent Number |
|
|
61168859 |
Apr 13, 2009 |
|
|
|
Current U.S.
Class: |
607/45 |
Current CPC
Class: |
A61N 1/0476 20130101;
A61N 1/0484 20130101; A61N 1/0529 20130101; G16H 30/20 20180101;
G16H 50/50 20180101; A61N 1/36025 20130101; A61N 1/0456
20130101 |
Class at
Publication: |
607/45 |
International
Class: |
A61N 1/36 20060101
A61N001/36 |
Claims
1. A method performed by data processing apparatus, the method
comprising: obtaining an image of target tissue; assigning to the
image tissue electrical conductance values; arranging a plurality
of stimulation electrodes around the target tissue; computing, from
the locations of electrodes and tissue electrical conductances, a
forward model of the response of the tissue to applied currents;
defining desired tissue response; optimizing one or more electrical
stimulation parameters using the forward model to obtain the
desired tissue response.
2. The method of claim 1 wherein the desired tissue response
includes field intensities or currents in a portion of the target
tissue and minimal stimulation of a second, different part of the
tissue.
3. The method of claim 1 wherein the desired tissue response
includes a change in the volume of tissue activated or a
physiological response of the target tissue.
4. The method of claim 1 wherein the desired tissue response
includes a change in the volume of tissue activated or a
physiological response of tissue or muscle that is separate and
different from the target tissue.
5. The method of claim 1 wherein the desired tissue response
includes strict constraints of maximum allowable currents of field
intensities at various tissue locations.
6. The method of claim 1 wherein the forward model is computed
using an finite-element model of the tissue properties.
7. The method of claim 1 wherein the electrical conductance are
non-isotropic and or non-uniform.
8. The method of claim 1 wherein the parameters altered include
changing the voltage, current, activation time, location, sequence
or number of electrodes.
9. The method of claim 1 where desired response is optimized with a
minimum number of electrodes.
10. The method of claim 1 wherein the step of optimizing a new
electrical stimulation pattern adjusts the results of the forward
model using least squares methodology and any of its derivative
forms such as constrained least squares, penalizes least squares,
ridge regression, elastic nets, etc.
11. The method of claim 7 wherein the step of optimizing a new
electrical stimulation pattern determines a volume of tissue
activated that is different than the volume of tissue activated to
determine the forward model.
12. The method of claim 1 wherein the image is derived from a
pre-existing image library or a target tissue specific image.
13. The method of claim 1 wherein the image is derived from a
fluoroscopic image, an MRI image, a CT image, or a combination of
imaging techniques.
14. The method of claim 1 wherein the target tissue is transcranial
tissue.
15. The method of claim 1 wherein the electrodes comprise at least
two or more electrodes.
16. The method of claim 1 wherein the electrodes comprise at least
10 or more electrodes.
17. The method of claim 1 wherein the electrodes comprise at least
100 or more electrodes.
18. The method of claim 1 wherein the electrodes comprise at least
200 or more electrodes.
19. The method of claim 1 wherein the electrodes comprise at least
256 or more electrodes.
20. The method of claim 1 wherein the plurality of electrodes are
place around the target tissue based on anatomical landmarks.
21. The method of claim 18 wherein the plurality of electrodes are
placed around the target tissue using the International 10-20
System.
22. The method of claim 1 wherein the plurality of electrodes are
placed around the target tissue in a pattern on the skin, below the
skin or within the target tissue.
23. The method of claim 1 wherein the electrical stimulation
applied is a direct current of 0 to 10 mA.
24. The method of claim 1 wherein the electrical stimulation
applied is an alternating current of 0-10 mA and 0 H-1 kHz.
25. The method of claim 1 wherein the electrical stimulation
applied is the same for each electrode in the plurality of
electrodes.
26. The method of claim 1 wherein the electrical stimulation
applied is different for each electrode in the plurality of
electrodes.
27. A computer storage medium encoded with a computer program, the
program comprising instructions that when executed by data
processing apparatus cause the data processing apparatus to perform
operations comprising: accepting an image model of target tissue;
obtaining a forward model having a first electrode configuration
and first electrical stimulation parameters based on electrical
stimulation of the target tissue; accepting electrode configuration
changes or electrical stimulation parameter changes resulting in a
second electrode configuration or second electrical stimulation
parameters; determining an optimized tissue model using a least
square methodology and based on the second electrode configuration
or second electrical stimulation parameter changes; comparing the
optimized tissue model with a desired outcome; and providing a
confirmation of the optimized model with the desired outcome.
28. The computer storage medium of claim 24 wherein the image model
of target tissue comprises a standard library image, a target
tissue specific image, a fluoroscopy image, an MRI image, a CT
image or a combination of images.
29. The computer storage medium of claim 24 wherein the forward
model is pre-calculated.
30. The computer storage medium of claim 24 wherein the forward
model is calculated utilizing a processor or processors in
communication with the storage medium.
31. The computer storage medium of claim 24 wherein the forward
model is a finite element model.
32. The computer storage medium of claim 24 wherein the first and
second electrode configurations are positions of electrodes around
or in the target tissue.
33. The computer storage medium of claim 24 wherein the first and
second electrical stimulation parameters comprise voltage,
amperage, current, duration, timing, or sequence of the electrical
stimulation.
34. The computer storage medium of claim 24 wherein the optimized
tissue model is determined using a linear least squares
methodology.
35. The computer storage medium of claim 24 wherein the desired
outcome comprises a change in the volume of tissue activated, a
physiological response of the target tissue, a physiological
response of tissue or muscle separate and different from the target
tissue.
36. The computer storage medium of claim 24 wherein the target
tissue is transcranial tissue.
Description
CROSS-REFERENCE TO RELATED APPLICATIONS
[0001] This application claims the benefit under 35 U.S.C.
.sctn.119(e) of U.S. Patent Application No. 61/168,859, entitled
"Neurocranial Electrostimulation Models, systems, Devices and
Methods," filed Apr. 13, 2009, which is incorporated herein by
reference in its entirety.
BACKGROUND
[0002] This specification generally relates to a method and
apparatus of electrostimulation of body tissue, including methods
of transcranial electrostimulation. For example, plastic changes in
brain function can be safely induced in humans by low-intensity
electrical stimulation through scalp electrodes. Such electrical
stimulation is known as transcranial electrostimulation (TES).
These changes can be potentially used for therapeutic or
performance enhancing applications. Currently available devices are
rudimentary in the sense that they do not target brain modulation
and do not apply insights from biophysical studies to functionally
target brain function. Moreover, current methods of TES require
costly and time consuming computations in order to alter treatment
methods.
SUMMARY
[0003] This specification describes systems and methods relating to
electrostimulation of body tissue including methods relating to
electrostimulation of transcranial tissue.
[0004] In general, one innovative aspect of the subject matter
described in this specification can be embodied in methods that
include the actions of obtaining a model image of target tissue,
developing a forward model of the electric field induced in the
target tissue based on electrode configurations and tissue
properties, and determining from this forward model the optimal
stimulation parameters using a least squares approach to obtain a
desired stimulation outcome. Other embodiments of this aspect
include corresponding systems, apparatus, and computer programs,
configured to perform the actions of the methods, encoded on
computer storage devices.
[0005] Embodiments and aspects can include a computer program for
optimization of a forward model having been obtained using finite
element model (FEM) computations by changing electrode
configurations or electrostimulation parameters utilizing: matrix
calculations, linear matrix addition, transposition, or inversion;
optimization algorithms; least-squares fits or least square
minimization; source localizations; local minimums; constraints
based on the number of electrodes; constraints based on the total
current; constraints based on the current in a specific electrode;
constraints that the currents in specific sets of electrodes must
add to zero; current control; and/or voltage control.
[0006] An aspect of the subject matter described in this
specification can be embodied in methods that include the actions
of: obtaining an image of target tissue; assigning to the image
tissue electrical conductance values; arranging a plurality of
stimulation electrodes around the target tissue; computing, from
the locations of electrodes and tissue electrical conductances, a
forward model of the response of the tissue to applied currents;
defining desired tissue response; optimizing one or more electrical
stimulation parameters using the forward model to obtain the
desired tissue response.
[0007] These and other embodiments can each optionally include one
or more of the following features. The desired tissue response
includes field intensities or currents in a portion of the target
tissue and minimal stimulation of a second, different part of the
tissue. The desired tissue response includes a change in the volume
of tissue activated or a physiological response of the target
tissue. The desired tissue response includes a change in the volume
of tissue activated or a physiological response of tissue or muscle
that is separate and different from the target tissue. The desired
tissue response includes strict constraints of maximum allowable
currents of field intensities at various tissue locations. The
forward model is computed using an finite-element model of the
tissue properties. The electrical conductance are non-isotropic and
or non-uniform. The parameters altered include changing the
voltage, current, activation time, location, sequence or number of
electrodes. The desired response is optimized with a minimum number
of electrodes. The step of optimizing a new electrical stimulation
pattern adjusts the results of the forward model using least
squares methodology and any of its derivative forms such as
constrained least squares, penalizes least squares, ridge
regression, elastic nets, etc. The step of optimizing a new
electrical stimulation pattern determines a volume of tissue
activated that is different than the volume of tissue activated to
determine the forward model. The image is derived from a
pre-existing image library or a target tissue specific image. The
image is derived from a fluoroscopic image, an MRI image, a CT
image, or a combination of imaging techniques. The target tissue is
transcranial tissue. The electrodes comprise at least two or more
electrodes. The electrodes comprise at least 10 or more electrodes.
The electrodes comprise at least 100 or more electrodes. The
electrodes comprise at least 200 or more electrodes. The electrodes
comprise at least 256 or more electrodes. The plurality of
electrodes are place around the target tissue based on anatomical
landmarks. The plurality of electrodes are placed around the target
tissue using the International 10-20 System. The plurality of
electrodes are placed around the target tissue in a pattern on the
skin, below the skin or within the target tissue. The electrical
stimulation applied is a direct current of 0 to 10 mA. The
electrical stimulation applied is an alternating current of 0-10 mA
and 0 H-1 kHz. The electrical stimulation applied is the same for
each electrode in the plurality of electrodes. The electrical
stimulation applied is different for each electrode in the
plurality of electrodes.
[0008] In general, another aspect of the subject matter described
in this specification can be embodied in methods that include the
actions of accepting an image model of target tissue; obtaining a
forward model having a first electrode configuration and first
electrical stimulation parameters based on electrical stimulation
of the target tissue; accepting electrode configuration changes or
electrical stimulation parameter changes resulting in a second
electrode configuration or second electrical stimulation
parameters; determining an optimized tissue model using a least
square methodology and based on the second electrode configuration
or second electrical stimulation parameter changes; comparing the
optimized tissue model with a desired outcome; and providing a
confirmation of the optimized model with the desired outcome.
[0009] Particular embodiments of the subject matter described in
this specification can be implemented so as to realize one or more
of the following advantages. Electrode configuration or placement
of individual electrodes in the electrode configuration can
facilitate: control the volume of (brain) tissue activated by
neurocranial electrostimulation; determination of the volume of
influence during stimulation; determination of the
neurophysiological outcome of stimulation; determination of the
clinical or behavioral outcome of stimulation; targeting of a
specific brain region such as the cortex, a cortical region, the
hippocampus, deep brain structures, axons, axons of passage,
pre-frontal cortex, motor region, or sensory regions; treatment of
a neurological or psychiatric disease; prevention of tissue damage
or prevent cognitive side-effects; determination of the hazards of
electrostimulation; accommodation of individual factors, for
example to optimize treatment based on patient specific anatomical
features; control of stimulation dose; triggering a specific
desired response; reduction of stimulation artifact; and creation
of a sham stimulation or establishing a control stimulation
condition.
[0010] Particular embodiments of the subject matter described in
this specification can be implemented so as to realize one or more
of the following additional advantages. Specific anode/cathode
relationships can be established between electrodes, such as
neurocranial electrodes, to stimulate the tissue near these
electrodes, such as brain tissue near the neurocranial electrodes.
Electrodes and other hardware including cranial caps and circuitry
can be fabricated based on an optimized model. Clinicians can
determine a course of treatment or even if a treatment should
proceed based on an optimized model without repeated trial and
error stimulation protocols on an actual patient. Electrical
stimulation systems for patient treatment can be programmed in
accordance with electrode configurations and stimulation parameters
identified using the optimized stimulation model. Clinicians can
see multiple optimized models without having to run multiple FEM
analyses; the cost of stimulation modeling, as reflected in
necessary resources (computers), both manual and computational
time, supplementary on-site expertise and technical help needed,
financial cost, and ability to deploy the optimization in a wide
range of environment are reduced. Specifically, 1) the "wait time"
for stimulation optimization is reduced and optimization can be
done locally by a clinician, for example on a lap-top; 2) the
resulted optimized stimulation parameters are superior than those
achieved with other (iterative methods) and moreover can include
additional practical constraints as needed, for example minimizing
skin irritation and damage or the number of electrodes (robustness,
simplicity, and accessibility of technology/therapy and complex
stimulators can be avoided) The system and methods outlined in this
advantage thus increases the safety and efficacy of Neurocranial
stimulation.
[0011] The details of one or more embodiments of the subject matter
described in this specification are set forth in the accompanying
drawings and the description below. Other features, aspects, and
advantages of the subject matter will become apparent from the
description, the drawings, and the claims.
BRIEF DESCRIPTION OF THE DRAWINGS
[0012] FIG. 1 is an illustration an exemplary image of target
tissue.
[0013] FIG. 2 is an illustration of an exemplary image of target
tissue.
[0014] FIG. 3 illustrates an exemplary electrode configuration
surrounding target tissue.
[0015] FIG. 4 illustrates an exemplary electrode configuration
incorporated into a wearable device.
[0016] FIG. 5 illustrates a process for determining a forward FEM
model of electrical stimulation of target tissue.
[0017] FIG. 6 illustrates a process for determining an optimized
model of electrical stimulation of target tissue according to the
present invention.
[0018] FIG. 7 illustrates a process for determining an optimized
model of electrical stimulation of target tissue according to the
present invention.
[0019] FIG. 8 illustrates a process for determining an optimized
model of electrical stimulation of target tissue according to the
present invention.
[0020] Like reference numbers and designations in the various
drawings indicate like elements.
DETAILED DESCRIPTION
[0021] When applying electrostimulation to body tissue, for example
neurocranial electrostimulation, decisions must be made regarding
the configuration of electrodes and the parameters of the
electrical inputs to the electrodes in order to obtain a desired
result from the electrostimulation. That is clinicians must
determine (1) the positioning of electrodes (how many and where
does each electrode go around the targeted tissue), and (2) how
much current (or voltage) to apply to each electrode, the duration
of the current (or voltage), and the sequence of the activation of
electrodes with the specified current (or voltage) for the
identified duration. These decisions regarding electrode
configuration and stimulation parameters fundamentally affect the
outcome of the electrical stimulation on the specific target
tissue. One advantage of electrical stimulation is the ability to
customize therapy for diseases, anatomical targets, and
individuals--but at the same time the many permutations of
stimulation configurations (electrode locations, current intensity)
make determining an optimal, or even good, configuration
challenging. This applies for positioning of electrodes in an
implant, but for Neurocranial stimulation using arrays of surface
electrode, of potentially high number, there are even more
potential permutations. Moreover, electrode configuration and
stimulation parameters have typically been determined using
multiple trial and error procedures on the patient along with
lengthy and costly computer analyses and modeling that are specific
to each electrode configuration and stimulation parameter. Computer
modeling using forward Finite Element Modeling ("FEM") analyses
have been run for "generic" or standard library tissue images and
data, such as generic head models. Forward FEM models have also
been run based on individual head data such as measurements and MRI
anatomical scans. Such FEM models reduce patient trial and error
procedures but involve expansive computer simulations that are
specific to each electrode configuration and stimulation parameter
associated with each tissue image/model. This is costly and time
consuming. Moreover, there are many permutations of potential
stimulation configuration resulting in an intractable number of
stimulation that would have be run iteratively--in this way an
ideal configuration is unlikely to be found. Implementations of the
present invention eliminate the need to repeated patient trial and
error procedures or multiple matrix or FEM calculations, thereby
reducing the time and cost of optimizing electrode configuration
and stimulation parameters to achieve a desired outcome from an
electrostimulation procedure. Moreover, the implementation in the
present invention yields the best electrode configuration and
currents with or without additional constraints as practically
relevant to the clinical application.
[0022] Aspects of the present invention relate to the speed of
optimization and the feasibility of optimization of a matrix
solution of a particular electrode configuration having a
particular stimulation parameter (current or voltage intensity). In
one aspect, a series of models are run that calculate the electric
field response in the target tissue when one electrode of a
plurality of electrodes in an electrode configuration is stimulated
at a given current or voltage intensity. In one embodiment, the
reference electrode in each case may be an additional fixed
electrode or set fixed point in the tissue. In another embodiment,
the solutions are from a set of pairs. This model is repeated for
each electrode position in the plurality of electrodes in the
electrode configuration surrounding the target tissue, (e.g., 100
electrode positions=100 simulations). The same stimulation
parameter (e.g., voltage or current intensity) is used for each of
run electrode or distinct stimulation parameters are run for each
electrode. Regardless, generally only a single stimulation
parameter is run for each electrode. Modeling each electrode in the
electrode configuration in this manner results in a "forward-model"
that is specific to each electrode configuration using a particular
voltage or current intensity. In some aspects, the electrode
positions can reflect locations on the subject's heads (e.g. 100
electrodes positions around the subjects head). The forward model
can then be optimized.
[0023] Once this series (e.g. 100 simulations) of models are all
solved and the solutions saved as the forward model, predictions
can be made as to what electrical fields will develop in the target
tissue when any combination of electrodes is activates at any
intensity. This is optimizing the forward model and produces the
"optimized model." Moreover, optimizing the forward model according
to aspects of the present invention can be done rapidly without the
time consuming and costly process of developing a separate matrix
or forward model for each desired intensity. Aspects of the
invention use pre-solved solutions for each electrode (at one
intensity per electrode) to predict what will happen if any
combination of the electrodes are activated with any combination of
electrode intensities. As such, a clinician wanting to know the
outcome of a specific electrode configuration (locations and
intensities), can use aspects herein, to immediately predict the
electric fields developed in the target tissue in response to the
stimulation specified in the optimized model.
[0024] Optimization in accordance with aspects of the present
invention can also allow determination of the optimal electrode
configuration (electrode number, position of each electrode,
current at each electrode) based on an outcome specified by the
clinician, such as production of a desired electrical field at a
desired tissue location. For example the clinician may want to
target one part of the brain without affecting neighboring regions
of the brain. Rather than the clinician trying the infinite number
of combinations possible, the optimization algorithm uses the
forward model solved solution series, and automatically calculates
the optical electrode configuration.
[0025] Developing the Model Image:
[0026] Obtaining a two-dimensional or 3-dimensional image of the
target tissue, such as a cranial image, and ascertaining tissue
properties including electrical characteristics of the tissue such
as conductivities, resistivities, permittivities, capacitances,
impedances, or applied energies or combinations thereof, allows for
development of an "image model" of the target tissue. Image models
can be developed from generic library tissue images, fluoroscopy
images, MRI images, CT scans or images, other imaging modalities
known in the art, or a combination of images. Tissue properties can
be developed from library information, testing of the target
tissue, or approximated from image data, as described in U.S.
Patent Application Pub. No. US 2007/0043268 A1, incorporated herein
by reference in its entirety.
[0027] A relationship of tissue resistivity to MRI gray scale that
can be correlated to tissue types can be expressed by the
formula:
R(V)=K(1-v).sup.E+D,
where [0028] R=Resistivity; [0029] V=Numeric value of MRI data*;
*The V value can be either simple MRI data values or combined
values from multiple MRIs or multiple types of MRIs. Exemplary
values include K=1600, E=4 and D=65. [0030] K=Multiplier value;
[0031] E=Exponent; and [0032] D=Density value.
[0033] Anisotropies/directionalities can be inferred from the
anatomy or determined based on the MRI data, or a combination
thereof. A direct determination is accomplished by diffusion tensor
MRI (DT-MRI, or DTI). The indirect is accomplished by inferring the
direction of fibers, specifically nerve fibers, by the general
anatomy. DT-MRI data are sometimes called Anisotrophic MRIs.
[0034] FIG. 1 illustrates MRIs of three different types: T1, T2 and
PD. Below each MRI is a gray scale. The gray scale for the T1 MRI
appears to resolve three peaks which may correspond to three
distinct tissue types having three different resistivities. The
gray scale for T2 shows one, or possibly two, peaks, and the gray
scale for PD shown one peak at a different resistivity than T2 or
T1. By utilizing information from different MRI types, it is
possible to enhance gray scale segmentation.
[0035] FIG. 2 illustrates a MRI and a plot of resistivities of
tissues showing multiple resolved peaks achieved by gray scale
differentiation of tissues of different resistivities. The gray
scale for the MRI shown in FIG. 1 can resolve multiple peaks
corresponding to various tissue types including compact bone,
cancellous bone, white matter, soft tissue, gray matter, skin,
blood and cerebral spinal fluid. Other resolvable tissues may
include cancerous tissue, inflammatory tissue and ischemic tissue,
as well as eye fluid. By having enhanced resolution of tissues, it
is possible to assign more correctly the vector resistivities or
other electrical values to brain or other body tissues, and thereby
developing a more accurate image model.
[0036] Electrode Configuration:
[0037] The electrode configuration is placed around the target
tissue. The electrode configuration includes two or more electrodes
(e.g., three or more electrodes, four or more electrodes, 8 or more
electrodes, 10 or more electrodes, 100 or more electrodes, 200 or
more electrodes, 256 or more electrodes). The electrodes deliver an
electrical stimulation having a desired current or voltage
intensity. For examples currents between about 0.1 mA and about 100
mA can be delivered to the target tissue via the electrodes in the
electrode configuration. For example, voltages of between about
0.1V and about 100V can be used. Not all electrodes in the
configuration need to be active. That is, in certain embodiments,
some electrodes may have zero current, 0 volts while others may
deliver the desired current or voltage intensity. FIG. 3
illustrates an exemplary electrode configuration around target
cranial tissue, where electrodes 310 are placed around the subjects
head 315. Electrodes 310 can be placed randomly or according to
established systems, such as the International 10-20 system.
Electrodes can be placed around the target tissue based on
anatomical landmarks. Electrodes 310 can be placed on the surface
of the skin, under the skin, or within the target tissue.
Electrodes 310 can be incorporated into a wearable device, such as
the skull cap 430 illustrated in FIG. 4.
[0038] Developing the Forward Model
[0039] Development of the forward model involves determining the
electric field response in the target tissue for each electrode in
the electrode configuration at a specified stimulation parameter
such as current or voltage intensity. The electric field response
of the target tissue can be described as the volume of tissue
activated in response to the electrical stimulation.
[0040] There are different ways to determine the "volume of tissue
activated". At the most basic level this includes induced brain
voltage (rarely used), electric fields (same as current density),
and electric field derivative (same as activating function). On a
more complex level, the calculations of induced voltage/field in
the brain are than passed through another set of equations that
more accurately predict brain activation. These second stage
filters can to include 1) at the simplest a threshold (e.g.
electric field greater than x is full activation, less than x is no
activation); 2) the electric field may be broken down into
direction where the direction themselves are based on the cortical
geometry (for example the electric field perpendicular to the
cortical surface) or sub-cortical anatomy (for example direction of
axons); 3) the effects of electric fields on neurons may be
directly computational determined (e.g. a (detailed) model so a
cortical neuron). On a yet more complex level, system level and
control theory analysis may be used to predict effects of cognition
and behavior. Any combination of this "volume of tissue activated"
solutions can form a "prior-set" of solutions wherein a set of
solutions for each electrode in the electrode configuration forms
the forward model.
[0041] The volume of tissue activated can be determined based on a
spheres model, an individualized spheres model, and anatomical
model based on gross measurements, and anatomical model based on
individual images such as MRI. The above "volume of tissue
activated" can be determined for a sub-set of potential stimulation
configurations. Where a stimulation configuration used a set of one
or more electrodes, in a specific arrangement, and a certain level
of current delivered to each electrode. The number of independent
configurations are tested and the "volume of tissue activated" is
determined for each electrode configuration. This forms a
pre-determined set of solutions. There is a time cost associated
with determined the volume of tissue activated.
[0042] Finite Element Modeling (FEM) as a method for determining
the volume of tissue activated can be used with the image model and
the specified electrode configuration to produce the forward model.
FEM techniques are described in U.S. Patent Application Pub. Nos.
US 2007/0043268 A1 and US 2006/0017749 A1, both of which are
incorporated herein by reference in their entirety. FEM analysis
may be done using any number of commercial computer programs, such
as FEMLAB by Comsol Pty. Ltd. of Burlington, Mass.
[0043] FIG. 5 is a flow chart illustrating generally one example of
a method of using a model to calculate a volume of activation, as
discussed above. Portions of the method may be embodied in any
machine-accessible medium carrying instructions for executing acts
included in the method. Such a method applies to deep brain
stimulation (DBS) or any other electrical tissue stimulation. At
300, imaging data of an anatomic volume of a patient is obtained.
In one example, this includes obtaining imaging data of a patient's
brain using an imaging modality, such as computed tomography (CT)
or magnetic resonance (MR) imaging modalities, for example, or
another suitable imaging modality. The anatomic volume need not be
all or part of the patient's brain, but could be all or part of any
other anatomic structure.
[0044] At 302, in one example, diffusion tensor imaging (DTI) data
is obtained (this may occur at 300, such as where a DTI MR imaging
modality is used at 300). In one example, the DTI data is obtained
from the same patient being analyzed. Alternatively, "atlas" or
library DTI data is obtained from at least one other patient. If
atlas DTI data from another patient is used, it is typically
spatially scaled to correspond to the anatomic size and shape of
the patient being analyzed. In one example, the atlas DTI data is
based on a composite from more than one other patient. The
composite atlas DTI data typically spatially scales DTI data from
the different patients before combining into the composite DTI
atlas. The atlas DTI data avoids the need to obtain DTI data from
the particular patient being analyzed. This is useful, for example,
when a non-DTI imaging modality is used at 300.
[0045] At 304, a tissue conductivity model, or image model as
discussed above, is created for all or part of the anatomic volume.
The tissue conductivity model is typically a non-uniform spatial
distribution. Such a model more accurately represents inhomogeneous
and anisotropic characteristics of the tissue anatomy. For example,
the conductivity of brain tissue varies from one brain region to
another. Moreover, conductivity of the nervous system is
preferential to a particular direction that is also dependent on
the particular location in the brain. In one example, a non-uniform
tissue conductivity model is created by transforming the DTI data
into conductivity data, such as by using linear transform
techniques known in the art.
[0046] It should be noted that it is not required to obtain
non-uniform tissue conductivity data using DTI. There exist several
alternatives to using DTI based approximations for the anisotropic
and inhomogeneous tissue properties for the patient specific finite
element volume conductor model. One example technique would be a
simple designation of a white matter and a grey matter conductivity
tensor, as discussed above. These two universal conductivity
tensors could then be applied to the nodes of the FEM mesh using
co-registration with the anatomical MRI. In this manner, the
individual voxels of the MRI data are designated as either white
matter or grey matter using post-processing image analysis. Then,
each such voxel is assigned a conductivity dependent on whether it
was classified as white matter or grey matter, which white matter
voxels having a different conductivity value than grey matter
voxels. A second example technique would define individual
conductivity tensors for designated brain regions (e.g., nuclei,
sub-nuclei, fiber tracts, etc.). This method would allow for a more
detailed representation of the tissue electrical properties than
the first example technique. The conductivity tensor of each
designated brain region is defined, in one example, using explicit
experimental tissue impedance results and anatomical information
provided by a human brain atlas. In this technique, the anatomical
MRI is sub-divided into different designated brain regions on a
voxel-by-voxel basis using post-processing image analysis. The
appropriate conductivity tensors for each designated brain region
is then co-registered with the nodes of the FEM mesh.
[0047] At 306, a finite element model (FEM) is created using the
conductivity data obtained at 304. In one example, the FEM model
uses a default boundary condition that is appropriate for a typical
electrode contact morphology. However, in another example, the FEM
model includes an electrode-specific boundary condition that is
tailored to the morphology of a particular electrode contact or
contacts to be used in the DBS or other procedure. The FEM model
provides for non-uniform conductivity in the tissue, such as by
using a DTI-derived other conductivity value at each node in the
FEM mesh. The FEM model may include aspects that are not obtained
from the DTI-derived data. In one such example, the FEM mesh models
a thin encapsulation sheath about the electrode lead body, as
discussed above, which is not derived from the DTI data.
[0048] At 308, in one example, the FEM is solved for the electric
potential distribution or the second difference (.DELTA.sup.2V) of
the electric potential distribution, such as by using FEM solver
software. In one example, the FEM is solved for a normalized
stimulation amplitude of 1V or 0.1 mA. In another example, for a
different electric stimulation amplitude, the resulting electric
potential distribution (or second difference of the electric
potential distribution) is multiplied by a scale ratio of the
different electric stimulation amplitude to the normalized electric
stimulation amplitude.
[0049] Scoring of Solutions Based on VOA
[0050] In one embodiment as outlined below, solutions are optimized
to reduce the error between a desired electric field distribution
and the electric field distribution achieved by applying the
optimized current to electrodes. In other embodiments additional
steps are taken to SCORE solutions based on additional analysis
including volume of activation (VOA). Scoring may be integrated
into the optimization process as described below, or used as an
additional processing, data display, and user interface step. The
prior art solved and scored each electrode configuration and
current individually and iteratively without optimization at
tremendous computational and time cost.
[0051] At 310, a volume of tissue activation (VOA) or other volume
of influence is calculated, in one example, using the second
difference of the electric potential distribution. The VOA
represents the region in which any neurons therein are expected to
typically be activated, that is, they are expected to generate
propagating action potentials at the stimulus frequency in response
to the electrical stimulation delivered at the stimulation
electrode contact. Conversely, neurons outside the VOA are expected
to typically remain unactivated in response to the electrical
stimulation. Neurons in the volume of influence are expected to be
modulated (neuromodulatoin) in a manner related to the strength of
regional electric field or some function there-of. In one example,
a particular threshold value of the second difference of the
electric potential distribution defines the boundary surface of the
VOA.
[0052] As discussed above, the particular threshold value defining
the boundary of the VOA is determined as follows. First, model
neuronal elements are positioned relative to the electrode using
known neuroanatomical information about specific fiber pathways and
nuclei of interest near the electrode. These generalized positions
of the model neuronal elements are then refined, such as by using
explicit "patient-specific" information provided in the DTI or
anatomical MR imaging data. For example, the DTI imaging data
describes the inhomogeneous and anisotropic tissue properties near
the electrode. In this example, such DTI imaging data is used to
explicitly define one or more axonal trajectories, if needed, or to
help define nuclear boundaries specified in the anatomical MRI.
[0053] A model of these neurons is then created. In one example,
the neurons are modeled using an axon model, which is a simplified
form of a neuron model. An example of an axon model is described in
Cameron C. McIntyre et al., "Modeling the Excitability of Mammalian
Nerve Fibers Influence of Afterpotentials on the Recovery Cycle,"
J. Neurophysiology, Vol. 87, February 2002, pp. 995-1006, which is
incorporated by reference herein in its entirety, including its
disclosure of axon models. In another example, a more generalized
neuronal model is used, an example of which is described in Cameron
C. McIntyre et al., "Cellular Effects of Deep Brain Stimulation:
Model-Based Analysis of Activation and Inhibition," J.
Neurophysiology, Vol. 91, April 2004, pp. 1457-1469, which is
incorporated by reference herein in its entirety, including its
disclosure of neuronal models. The neuron model describes how the
neurons will respond to an applied electric field, that is, whether
the neuron will fire and whether the neurons will generate a
propagating action potential.
[0054] In one example, using this neuron model to simulate how the
neurons (located as determined from the DTI-derived conductivity
data, in one example) behave, the threshold value of the second
difference of electric field that will result in such propagating
action potentials is calculated. The stimulating influence of the
electric field is applied to the model neurons to define a
threshold value. This threshold value is then used to define the
boundary of the VOA in the non-uniform conductivity tissue, as
discussed above.
[0055] It should be noted that calculation of explicit threshold
criteria for each patient is not required. For example, in a more
generalized situation, threshold criteria will have already been
determined using the detailed neuron models under a wide variety of
different stimulation conditions. Once these threshold criteria
have been determined, they need not be re-determined for each
subsequent patient.
[0056] It should also be noted that using a threshold criteria upon
the second difference of the potential distribution in the tissue
medium is a simplified technique for quickly determining a VOA or
other volume of influence. The intermediate step of using the
second difference of the potential distribution is not required. In
an alternate example, the FEM model is directly coupled to a
detailed neuron model, such as a multi-compartment neuron model
that is oriented and positioned in the FEM model to represent at
least one actual nerve pathway in the anatomic volume.
[0057] At 312, the calculated VOA region is displayed, such as on a
computer monitor. In one example, the VOA is displayed superimposed
on the displayed imaging data or a volumetric representation
derived from such imaging data. In another example, an anatomic
boundary or other representation of an anatomic structure is
superimposed on the VOA and imaging data or the like. The anatomic
boundary data is typically obtained from an atlas of brain anatomy
data, which can be scaled for the particular patient, as discussed
above. Alternatively, the anatomic representation is extracted from
the imaging data for the patient being analyzed. In one example,
the anatomic representation is a line depicting one or more
boundaries between particular nucleus structures or other regions
of the brain, such as the STN, IC, or ZI illustrated above in FIG.
1B.
[0058] In any case, by viewing a representation emphasizing one or
more brain regions displayed together with the VOA, the user can
then determine whether a particular anatomic region falls within or
outside of the modeled VOA. The user may want a particular anatomic
region to be affected by the DBS, in which case that region should
fall within the modeled VOA. Alternatively, the user may want a
particular region to be unaffected by the DBS, such as to avoid
certain unwanted DBS stimulation side effects, as discussed above.
This evaluation of whether the VOA is properly located can
alternatively be performed by, or assisted by, a computer
algorithm.
[0059] For example, the computer algorithm can evaluate various
VOA's against either or both of the following input criteria: (a)
one or more regions in which activation is desired; or (b) one or
more regions in which activation should be avoided. In one example,
at 314, the computer algorithm creates a score of how such
candidate VOAs map against desired and undesired regions. In one
example, the score is computed by counting how many VOA voxels map
to the one or more regions in which activation is desired, then
counting how many VOA voxels map to the one or more regions in
which activation is undesired, and subtracting the second quantity
from the first to yield the score. In another example, these two
quantities may be weighted differently such as, for example, when
avoiding activation of certain regions is more important than
obtaining activation of other regions (or vice-versa). In yet
another example, these two quantities may be used as separate
scores.
[0060] At 316, the score can be displayed to the user to help the
user select a particular VOA (represented by a particular electrode
location and parameter settings). Alternatively, the algorithm can
also automatically select the target electrode location and
parameter settings that provide the best score for the given input
criteria.
[0061] In one example, the VOA is displayed on a computer display
monitor of an image-guided surgical (IGS) workstation, such as the
StealthStation.RTM. from the Surgical Navigation Technologies, Inc.
(SNT) subsidiary of Medtronic, Inc., for example. The VOA can be
displayed on the IGS workstation monitor with at least one of the
imaging data representing the anatomic volume, the target electrode
location, a burr hole or other anatomic entry point, a trajectory
between the anatomic entry point and the target electrode location,
or an actual electrode location.
[0062] In one IGS workstation example, the displayed VOA
corresponds to a target electrode location.
[0063] After the electrode is positioned at the target location,
there remains the challenging task of adjusting the DBS stimulation
parameters (e.g., the particular electrode contact(s) of a
plurality of electrode contacts disposed on the same DBS leadwire,
pulse amplitude, pulsewidth, electrode "polarity" (i.e., monopolar
or bipolar electrode return path), electrode pulse polarity (i.e.,
positive or negative), frequency, etc.). In one example, the IGS
workstation or a DBS pulse generator programmer includes the
above-described VOA methods to assist the user in selecting an
appropriate combination of DBS stimulation parameters, such as by
using the scoring techniques discussed above. In essence, for each
electrode position, a forward model can be calculated. As discussed
above, this is time consuming and costly.
[0064] FIG. 6 is a flow chart illustrating generally one example of
a method of using an FEM model to calculate a volume of activation,
as discussed above, and using the volume of tissue activation to
select a particular electrode configuration. Portions of the method
may be embodied in any machine-accessible medium carrying
instructions for executing acts included in the method. Such a
method applies to selecting an electrode configuration for deep
brain stimulation (DBS) or for any other electrical tissue
stimulation. At 500, a set of N candidate electrodes are defined,
where N is an integer greater than 1. Defining the candidate
electrode configurations typically includes providing information
about the size, shape, or arrangement of electrode contacts on a
leadwire. Such information is typically in a form in which it can
be used as input to a finite element model (FEM). At 502, a FEM is
created for each candidate electrode position. The FEMs typically
use non-uniform conductivity model of a desired region of interest.
At 504, each FEM is solved for a second difference in the electric
potential distribution. At 506, a volume of activation (VOA) is
computed for each candidate electrode morphology from its
corresponding second difference in the electric potential
distribution. The boundary of the VOA is typically obtained from a
threshold value that is based on a neuron or axon model, as
discussed above. At 508, the VOAs are scored, as discussed above,
or otherwise evaluated to select one or more electrode locations
that exhibit a desired VOA, or a VOA that is deemed more desirable
than the VOA of one or more other electrode morphologies. At 510,
at least one electrode is manufactured using the selected at least
one electrode morphology.
[0065] Optimizing the Forward Model:
[0066] The above image collection and modeling techniques require
time intensive and costly computations to determine the volume of
tissue activated for a particular electrode configuration using
specific stimulation parameters. Moreover, changes to electrode
positions and or changes to stimulation parameters, such as current
or voltage intensities, require recomputation to verify the new
volume of tissue activated. For example, to solve and score a
multitude of solution as described above is very costly
computationally and may not result in an optimal configuration.
[0067] It has been found that changes to the forward model, and
particularly changes to current or voltage intensities can be
optimized or predicted using linear approximations. This avoids the
need for re-running costly computational forward models and greatly
enhances the speed to which electrical stimulation procedures can
be modeled, outcomes predicted, and treatments prescribed.
[0068] In an aspect, the governing equation for the optimization
process is given by the simple linear relationship:
E=MI (1)
[0069] where:
[0070] E is a vector that represents the field magnitudes at
multiple locations in the tissue.
[0071] I is a vector that represents the currents at multiple
electrodes.
[0072] M is a matrix that linearly relates the currents I to the
fields E and is called a "forward model".
[0073] The goal of optimization is to choose an optimal current
vector I* with a net-zero to current such that desired field
magnitudes E* are achieved, while maintaining maximum current
limits:
I * = argmin I E * - MI 2 subject to constraint I < I ma x and
sum ( I ) = 0 ( 2 ) ##EQU00001##
[0074] Other constraints such as "minimum number of non-zero
currents", "maximum allowable field intensities" or "maximum
allowable current densities" can be incorporated into this approach
using additional optimization criteria (such as an L1-norm penalty
term) or linear boundary constraints.
[0075] The forward model M can be computed from an 3D distribution
of electrical conductances and the locations of the electrodes and
locations for which the field is to be computed. The values in
matrix M are the solution to the quasi-static Laplace equation
(simplification of Maxwell equations) under Dirichlet boundary
conditions which govern the relationship between static currents I,
ie. the boundary conditions, and the resulting electric fields E,
i.e. the solutions, in a purely resistive material. Each column in
matrix M represents the solutions for a unit current in a given
electrode pair (all other electrode carrying no current). Each row
in matrix M represents a different location in the tissue. If N
electrodes are used and a single electrode is used as a common
reference for all pairs, then matrix M has N-1 columns. The number
of rows of matrix M scales with the number of 3D locations in the
tissue for which one would like to specify field magnitudes. To be
more specific, there is one such matrix for each field orientation:
M.sub.x for lateral-direction field E.sub.x, M.sub.y for
depth-direction field E.sub.y and M.sub.z for vertical-direction
field E.sub.z:
[0076] M=[M.sub.x.sup.T, M.sub.y.sup.T, M.sub.y.sup.T].sup.T, and
E=[E.sub.x.sup.TE.sub.y.sup.TE.sub.y.sup.T].sup.T. Thus the number
of rows is precisely 3 times the number of tissue locations of
interest.
[0077] The 3D current distribution can be obtained, for instance,
from a segmentation of a 3D image of the tissue such as MRI, CT,
DTI, etc. Such a segmentation can be reformulated in a
finite-element model (FEM), which is then used to provide an
efficient numerical solution to the Laplace equation. Commercial
software such as Abacus or Comsole is available to perform these
numerical FEM computation.
[0078] In an example, the optimization problem given in equation
(2) can be formulated on the basis of the linear relationship (1)
and that matrix M does not have to be recomputed using
computationally intensive numerical approaches for every possible
choice of I. Thus the optimization can be solved strictly using
efficient least-squares optimization procedures and do not require
an expensive trial-and-error search for I.
[0079] FIG. 7 is a flow chart illustrating generally one example of
optimizing a model to determine or predict changes to a stimulation
procedure, as discussed above. Portions of the method may be
embodied in any machine accessible medium carrying instructions for
executing acts included in the method. Such a method applies to
electrical stimulation of bodily tissue, including transcranial
tissue. At 610, an image of target tissue is obtained, as discussed
above. The image can be obtained from generic library tissue
images, fluoroscopy images, MRI images, CT scans or images, or a
combination of images. The target tissue can be any bodily tissue
including transcranial tissue.
[0080] At 620 conductance values are determined for the target
tissue. Conductance values can be a uniform or a non-uniform
spatial distribution. Conductance values can be determined from
pre-existing data of generic tissue images, target tissue specific
data derived from testing, or determined from grey scale MRI images
of the target data as described above.
[0081] At 630, electrodes are configured about the target tissue.
For example, electrodes are distributed around a patient's head
according to the International 10-20 System. The number of
electrodes can range from 2 or more, 3 or more, 10 or more, 100 or
more 200 or more, and 256 or more. A current and voltage intensity
is determined. The applied current can be in the range of about 0.1
mA to about 100 mA. The voltage can be in the range of about 1V to
about 100V. In aspects, the same current or voltage intensity is
used for each electrode in the electrode configuration.
[0082] At 640, electrical stimulation is applied to each electrode,
using the same current and voltage intensity. For each electrode
stimulation, the volume of tissue activated is calculated.
Cataloguing each volume of tissue activated for each electrode
stimulation results in the development of the forward model at
650.
[0083] At 660 a desired tissue response is defined. For example, a
particular volume of tissue activated may be desired. Other
examples include maximizing one are of tissue activation while
minimizing other areas of tissue activation. A particular
physiological response may be desired.
[0084] At 670, the forward model is optimized, as described above,
to predict the electrode configuration and current or voltage
intensities required to produce the desired tissue response.
[0085] FIG. 8 is a flow chart illustrating generally one example of
optimizing a model to determine or predict changes to a stimulation
procedure, as discussed above. Portions of the method may be
embodied in any machine accessible medium carrying instructions for
executing acts included in the method. Such a method applies to
electrical stimulation of bodily tissue, including transcranial
tissue.
[0086] At 710, an image model, as described above is accepted. The
image model includes a 2-D or 3-D image of the target tissue along
with certain electrical characteristics of the target tissue. For
example, the image model can be of transcranial tissue and can
include tissue conductance values. It will be noted that the image
model need not be developed locally but can be down loaded from a
file. The image model can be generic image model or a patient
specific image model.
[0087] At 720 Finite Element Model (FEM) of the image model is
obtained. The FEM is based on a particular electrode configuration,
such as electrodes arranged according to the International 10-20
System. The number of electrodes can range from 2 or more, 3 or
more, 10 or more, 100 or more 200 or more, and 256 or more. The FEM
is also based on a particular stimulation parameter (e.g., current
and voltage intensities). A current and voltage intensity is
determined. The applied current can be in the range of about 0.1 mA
to about 100 mA. The voltage can be in the range of about 1 V to
about 100V. In aspects, the same current or voltage intensity is
used for each electrode in the electrode configuration. It will be
noted that the FEM does not need to be developed locally but can be
downloaded from a file.
[0088] At 730 a second electrode configuration (e.g., movement of
one or more electrodes in the first electrode configuration) and/or
a second stimulation parameter (e.g., a change in the current or
voltage intensity from the first stimulation parameter) is
accepted.
[0089] At 740 the FEM is optimized using the optimization formula
described above to determine the tissue response from the changes
in electrode configuration or stimulation parameters. In an
example, a desired tissue response is determined and the electrode
configuration and stimulation parameter predicted utilizing the
optimization formula described above. At 750 the optimized model
can be compared with a desired outcome or physiological response.
In an example, the optimized model can be compared against certain
constraints, such as maximum voltages or currents at one electrode
or all electrodes, or maximum electrical fields at any tissue
location.
[0090] At 760, the electrode configuration and stimulation
parameters can be adjusted based on the optimized model and the
desired outcome or known constraints. The desired outcome can be
adjusted based on the electrode configuration and stimulation
parameters used to optimize the model.
[0091] At 770, confirmation of the optimized model and the desired
outcome, along with the electrode configurations and stimulation
parameters can be provided to a user via a user interface, a
print-out or other means known in the art.
Additional Embodiments
[0092] Aspects of the present invention use a matrix set of
solutions to determine the "volume of tissue activated" or "region
of influence" for a separated electrode configuration (one in which
the volume of tissue activated has not been determined)--this is a
new electrode configuration of interest. This is done by combining
the set of solutions in a manner that minimizes error. The
solutions may, for example, be each weighted linearly such that
there is a weight associated with each solution, and the each
solution multiplied by this weight added. For example, a series of
solutions of electric fields in the brain. The weight may be
determined such that the electrical current injection by the
summation of each solution multiplied by its respective weight
equals the electrical current injection by the new electrode
configuration of interest. The weight may be determined by
minimizing the error between the combined electrical current
injected by the set of solution multiplied by the weights and the
electrical current injection by the configuration of interest. The
weight may be determined by other means including random selection
and testing or based in a linear summation or prior knowledge about
the electrode configurations. The time cost of determining these
weights, and then by multiple each prior solution bits respective
weights, and then addition of this new multiplied set of solutions,
is less than the time cost of running a forward model (e.g., a
finite element model) of the new configuration. The computational
cost of determining these weights, then multiplying each prior
solution by it's respective weights, and then addition of this new
multiplied set of solutions, is less than the computational cost of
running a forward model of the new configuration. The complexity
cost of determining these weights, then multiplying each prior
solution bits respective weights, and then addition of this new
multiplied set of solutions, is less than the complexity cost of
running a forward model of the new configuration.
[0093] In another aspect, the weights applied to each
prior-solution of volume of tissue activation are determined by an
optimization algorithm. The optimization algorithm determines the
optimal weight that may be applied to the prior-set of solutions to
achieve a desired outcome. The optimization algorithm achieves this
at a time and computational cost that is less than trying every
possible electrode configuration. The optimization algorithm
compares the set of prior solutions (e.g., prior volume of tissue
activation) with the solution corresponding to the desired clinical
outcome (e.g., new volume of tissue activation) and then determines
the weight necessary to make these as "close as possible." The
criteria for "as close as possible" are determined by pre-set
conditions of the operator and can include minimizing the
difference between the clinically desired "volume of tissue
activation" and the volume of tissue activation by a configuration
determined by the algorithm.
[0094] In one implementation the following steps are taken: The
electric fields generated (forward model FEM) from monopole
stimulation with a single electrode with 1 mA (or 1V) is calculated
using a head model. This prior-solution is done for 128 electrode
positions. The return electrode is at a fixed position for all
electrode stimulations. This results in a set of 128
prior-solutions of electric fields in the head. This can be done
"off-line" with a higher time and computational cost. The user
interface allows a clinician to select to activate any combination
of these 128 electrodes with any level of current (or voltage) at
the same time. This is done by the clinician indicating how much
current X to apply at each of the 128 electrodes. The position of
each of these electrodes corresponds to the position of electrodes
in the prior-solution. The current at each electrode does not have
to be the same and can be zero. (On-line) the electric field
generated by the clinician selected configuration is determined as
follows:
[0095] For each electrode selected by the clinician a level of
current X is specified. The prior solution (for 1 mA or 1 V)
associated with each electrode is multiplied by X. This new set of
solutions is then added together. The time and computational cost
associated with this operation is minimized. The results of this
operation are displayed to the clinicians. The results of this
operation can determine the conditions of stimulation.
[0096] In another implementation the following steps are taken: The
electric fields are generated (forward head model FEM) from
monopole stimulation with a single electrode with 1 mA (or 1V) and
calculated using a head model. This prior-solution is done for 128
electrode positions. The return electrode is at a fixed position
for all electrode stimulations. This results in a set of 128
prior-solutions of electric fields in the head. This can be done
"off-line" with a higher time and computational cost. The user
interface allows a clinician to select a desired electric field
distribution in the head. How much current to apply to each
electrode to achieve this desired electric field distribution in
the head is then calculated. This is done by comparing the set of
prior-solutions with the desired electric field distribution and
then determining how much current to apply to each electrode to
minimize the difference between the desired electric field
distribution and the electric field distribution induced by applied
said current to the electrodes. The results of this operation are
displayed to the clinicians. The results of this operation can
determine the conditions of stimulation.
[0097] In another example, software calculates the stimulation
efficacy and target specificity for a given stimulation protocol.
This software includes a 3-D representation of the head impedance
including the scalp, skull, CSF, and brain compartments. The 3-D
representation requires manually-assisted segmentation of
anatomical MRI scans. The user (clinician) can place two or more
electrodes on the segmented scalp surface and simulate brain
modulation and iteratively "explore" high density transcranial
electrostimulation configurations. In prior-art applications, for
each configuration change, the predicted brain modulation must be
"re-solved", which takes time. In addition, placing the electrodes
and determination of the current flows has been an ad-hock process.
The optimal configuration may ultimately depend on an integration
of patient, disease, and clinical experience related factors, and
it is thus useful for a clinician to be able to intuitively and
rapidly evaluate a range of configurations.
[0098] The high spatial specificity provided by in aspects of the
present invention offers the opportunity to design a stimulation
protocol that is tailored to the individual subject and the desired
target location. Specifically, the task for the targeting
methodologies disclosed herein is to determine the ideal electrode
locations and to compute the required current flows through each of
these electrodes; moreover the targeting processes should be
accessible to a clinician using a PC, desktop computer, or portable
computing device. Three technical challenges associated with the
targeting processes disclosed herein include: (1) Automatic patient
specific MRI segmentation; (2) Acceleration of the computational
time for predicting brain modulation in response to a given
stimulation configuration (3) Computation of the optimal electrode
locations and currents for a user selected brain target.
Forward Model
[0099] High resolution (gyri/sulci precise) MRI derived finite
element (FE) human head models were generated by segmenting grey
matter, white matter, CSF, skull, muscle, fatty tissue, eyes, blood
vessels, scalp etc. Each model may comprise >10 million elements
with >15 million degrees of freedom. The induced cortical
electric field/current density values are calculated and used to
predict regions of brain modulation.
Subject Specific
[0100] To obtain a subject specific prediction of current flows
existing software packages were used to automatically segment a
given subject's anatomical MRI. As significant effort by the
research community has already generated a number of publicly
available segmentation packages that focus on MRI segmentation
(Shattuck & Leahy, 2002, Dogdas 2005, see brainsuite.usc.edu).
These are adapted to integrate with the commercial FEM solver such
as SIMPLEWARE (SIMPLEWARE LTD. to import MRI scans, segmentation,
creating FE mesh, and COMSOL INC. for FEM computation). The result
of this integration work was a software package that takes an
anatomical MRI and computes the required forward models for each
electrode location.
[0101] Variability in anatomy is well characterized by head size,
gender, and age. Therefore, a set of standard forward models can be
computed and provided as part of the real-time clinician interface
to cover a range of subjects to be selected according to these
criteria. If a subject-specific MRI is available the corresponding
forward-model can be computed remotely using high-end
computers.
[0102] In another example, computationally costly segmentation and
FEM solving is run on a patient specific basis, while physician
targeting software can operate on real-time by virtue of using
pre-solved FEM solutions.
[0103] Another example relates to real-time optimization of
stimulation protocols. It has been found that current flows within
the brain are a linear superposition of the currents generated by
each electrode. Based on this observation the optimization problem
can be formulated as a constrained least-squares problem. Least
squares aims to minimize the square difference of a desired current
distribution (specified by the clinician) with the achievable
current distribution as predicted by the FEM. The method is
computationally efficient as the predicted current distribution can
be formulated as a linear combination of the "forward model" that
has been previously computed for each electrode. The forward model
computation using the FEM is computationally expensive and can be
performed remotely (along with segmentation).
[0104] Standard electrode locations to facilitate proper electrode
placement (10/20 system conventionally used in EEG). While this is
not typically required, one can increase specificity by confirming
electrode locations with conventional 3D electrode-localization
hardware (e.g. Fastrak, Patriot, see cortechsolutions.com). The
constraints in the optimization procedure implement the
requirements that currents at each electrode do not exceed some
threshold value (corresponding to skin irritation and perception
limits) as well as limits on maximum current flows within the
tissue (safety limits). Analogous optimization problems are adapted
with sensor arrays in acoustics (Parra, 2006). The computationally
expensive step of computing the FEM is performed in advance for a
set of standard electrode locations (, while the optimization of
the electrode currents to target a specific cortical area is
performed in real-time. A GUI allows the physician to define the
target location and inspect the expected optimal current
distribution.
Targeting Software
[0105] In another example, targeting software calculates the
stimulation efficacy and target specificity for a given stimulation
protocol. This software includes a 3-D representation of the head
impedance including the all significant anatomical compartments. In
the system the 3-D representation requires manually-assisted
segmentation of anatomical MRI scans. The user (experimenter) can
place two or more electrodes on the segmented scalp surface and
predict resulting current distributions on the brain. In the
current system, for each configuration change, the predicted
current distribution in the brain must be recomputed, which takes
time. In addition, placing the electrodes and deciding on the
adequacy of resulting current flows is currently an ad-hoc process.
The optimal configuration may ultimately depend on an integration
of various factors, including patient brain anatomy, desired target
area, safety and efficacy considerations. This software allows the
experimenter to intuitively and rapidly evaluate a range of
electrode configurations.
[0106] The high spatial specificity provided by HD-TES offers the
opportunity to design a stimulation protocol that is tailored to
the individual subject and the desired target location.
Specifically, the task for the HD-TES targeting software is to
determine the ideal electrode locations and to compute the required
current flows through each of these electrodes; moreover the
targeting software should be accessible to an experimenter using a
conventional PC. There are three technical challenges related to
this task: (1) Automatic patient specific MRI segmentation; (2)
Acceleration of the computational time for predicting brain
modulation in response to a given stimulation configuration (3)
Computation of the optimal electrode locations and currents for a
user selected brain target. These challenges can be addressed with
two software modules, one for computing the head model of the
current flows based on anatomical MRI, and the second for
optimizing the applied currents and electrode locations.
Subject Specific Forward Model.
[0107] To obtain a subject specific prediction of current flows,
existing software packages can be used to automatically segment a
given subject's anatomical MRI. As significant effort by the
research community has already generated a number of publicly
available segmentation packages that focus on MRI segmentation
(Shattuck & Leahy, 2002, Dogdas 2005, see brainsuite.usc.edu).
These packages can be adapted to integrate with a commercial FEM
solver (e.g., SIMPLEWARE LTD. to import MRI scans, segmentation,
creating FE mesh, and COMSOL INC. for FEM computation). The result
of this integration work is a software package that takes an
anatomical MRI and computes the required forward models for each
electrode location. Fully automated segmentation is a challenging
problem in the art. As such, the accuracy of automated tools are
compared to a hand-segmented MRI and evaluated to see if the
differences in automated vs hand-segmentation leads to
significantly different predicted cortical current
distributions.
[0108] Variability in anatomy is well characterized by head size,
gender, and age. Therefore, a set of standard forward models will
be generated and provided as part of the real-time targeting
module.
[0109] In an example, the computationally expensive step of
computing the FEM is performed in advance for a set of standard
electrode locations, while the optimization of the electrode
currents to target a specific cortical area is performed in
real-time. The GUI allows the experimenter to define the target
location and inspect the expected optimal current distribution in
real-time.
[0110] Embodiments of the subject matter and the operations
described in this specification can be implemented in digital
electronic circuitry, or in computer software, firmware, or
hardware, including the structures disclosed in this specification
and their structural equivalents, or in combinations of one or more
of them. Embodiments of the subject matter described in this
specification can be implemented as one or more computer programs,
i.e., one or more modules of computer program instructions, encoded
on computer storage medium for execution by, or to control the
operation of data processing apparatus. Alternatively or in
addition, the program instructions can be encoded on an
artificially-generated propagated signal, e.g., a machine-generated
electrical, optical, or electromagnetic signal, that is generated
to encode information for transmission to suitable receiver
apparatus for execution by a data processing apparatus. A computer
storage medium can be, or be included in, a computer-readable
storage device, a computer-readable storage substrate, a random or
serial access memory array or device, or a combination of one or
more of them. Moreover, while a computer storage medium is not a
propagated signal, a computer storage medium can be a source or
destination of computer program instructions encoded in an
artificially-generated propagated signal. The computer storage
medium can also be, or be included in, one or more separate
physical components or media (e.g., multiple CDs, disks, or other
storage devices).
[0111] The operations described in this specification can be
implemented as operations performed by a data processing apparatus
on data stored on one or more computer-readable storage devices or
received from other sources.
[0112] The term "data processing apparatus" encompasses all kinds
of apparatus, devices, and machines for processing data, including
by way of example a programmable processor, a computer, a system on
a chip, or multiple ones, or combinations, of the foregoing The
apparatus can include special purpose logic circuitry, e.g., an
FPGA (field programmable gate array) or an ASIC
(application-specific integrated circuit). The apparatus can also
include, in addition to hardware, code that creates an execution
environment for the computer program in question, e.g., code that
constitutes processor firmware, a protocol stack, a database
management system, an operating system, a cross-platform runtime
environment, a virtual machine, or a combination of one or more of
them. The apparatus and execution environment can realize various
different computing model infrastructures, such as web services,
distributed computing and grid computing infrastructures.
[0113] A computer program (also known as a program, software,
software application, script, or code) can be written in any form
of programming language, including compiled or interpreted
languages, declarative or procedural languages, and it can be
deployed in any form, including as a stand-alone program or as a
module, component, subroutine, object, or other unit suitable for
use in a computing environment. A computer program may, but need
not, correspond to a file in a file system. A program can be stored
in a portion of a file that holds other programs or data (e.g., one
or more scripts stored in a markup language document), in a single
file dedicated to the program in question, or in multiple
coordinated files (e.g., files that store one or more modules,
sub-programs, or portions of code). A computer program can be
deployed to be executed on one computer or on multiple computers
that are located at one site or distributed across multiple sites
and interconnected by a communication network.
[0114] The processes and logic flows described in this
specification can be performed by one or more programmable
processors executing one or more computer programs to perform
actions by operating on input data and generating output. The
processes and logic flows can also be performed by, and apparatus
can also be implemented as, special purpose logic circuitry, e.g.,
an FPGA (field programmable gate array) or an ASIC
(application-specific integrated circuit).
[0115] Processors suitable for the execution of a computer program
include, by way of example, both general and special purpose
microprocessors, and any one or more processors of any kind of
digital computer. Generally, a processor will receive instructions
and data from a read-only memory or a random access memory or both.
The essential elements of a computer are a processor for performing
actions in accordance with instructions and one or more memory
devices for storing instructions and data. Generally, a computer
will also include, or be operatively coupled to receive data from
or transfer data to, or both, one or more mass storage devices for
storing data, e.g., magnetic, magneto-optical disks, or optical
disks. However, a computer need not have such devices. Moreover, a
computer can be embedded in another device, e.g., a mobile
telephone, a personal digital assistant (PDA), a mobile audio or
video player, a game console, a Global Positioning System (GPS)
receiver, or a portable storage device (e.g., a universal serial
bus (USB) flash drive), to name just a few. Devices suitable for
storing computer program instructions and data include all forms of
non-volatile memory, media and memory devices, including by way of
example semiconductor memory devices, e.g., EPROM, EEPROM, and
flash memory devices; magnetic disks, e.g., internal hard disks or
removable disks; magneto-optical disks; and CD-ROM and DVD-ROM
disks. The processor and the memory can be supplemented by, or
incorporated in, special purpose logic circuitry.
[0116] To provide for interaction with a user, embodiments of the
subject matter described in this specification can be implemented
on a computer having a display device, e.g., a CRT (cathode ray
tube) or LCD (liquid crystal display) monitor, for displaying
information to the user and a keyboard and a pointing device, e.g.,
a mouse or a trackball, by which the user can provide input to the
computer. Other kinds of devices can be used to provide for
interaction with a user as well; for example, feedback provided to
the user can be any form of sensory feedback, e.g., visual
feedback, auditory feedback, or tactile feedback; and input from
the user can be received in any form, including acoustic, speech,
or tactile input. In addition, a computer can interact with a user
by sending documents to and receiving documents from a device that
is used by the user; for example, by sending web pages to a web
browser on a user's client device in response to requests received
from the web browser.
[0117] Embodiments of the subject matter described in this
specification can be implemented in a computing system that
includes a back-end component, e.g., as a data server, or that
includes a middleware component, e.g., an application server, or
that includes a front-end component, e.g., a client computer having
a graphical user interface or a Web browser through which a user
can interact with an implementation of the subject matter described
in this specification, or any combination of one or more such
back-end, middleware, or front-end components. The components of
the system can be interconnected by any form or medium of digital
data communication, e.g., a communication network. Examples of
communication networks include a local area network ("LAN") and a
wide area network ("WAN"), an inter-network (e.g., the Internet),
and peer-to-peer networks (e.g., ad hoc peer-to-peer networks).
[0118] The computing system can include clients and servers. A
client and server are generally remote from each other and
typically interact through a communication network. The
relationship of client and server arises by virtue of computer
programs running on the respective computers and having a
client-server relationship to each other. In some embodiments, a
server transmits data (e.g., an HTML page) to a client device
(e.g., for purposes of displaying data to and receiving user input
from a user interacting with the client device). Data generated at
the client device (e.g., a result of the user interaction) can be
received from the client device at the server.
[0119] While this specification contains many specific
implementation details, these should not be construed as
limitations on the scope of any inventions or of what may be
claimed, but rather as descriptions of features specific to
particular embodiments of particular inventions. Certain features
that are described in this specification in the context of separate
embodiments can also be implemented in combination in a single
embodiment. Conversely, various features that are described in the
context of a single embodiment can also be implemented in multiple
embodiments separately or in any suitable subcombination. Moreover,
although features may be described above as acting in certain
combinations and even initially claimed as such, one or more
features from a claimed combination can in some cases be excised
from the combination, and the claimed combination may be directed
to a subcombination or variation of a subcombination.
[0120] Similarly, while operations are depicted in the drawings in
a particular order, this should not be understood as requiring that
such operations be performed in the particular order shown or in
sequential order, or that all illustrated operations be performed,
to achieve desirable results. In certain circumstances,
multitasking and parallel processing may be advantageous. Moreover,
the separation of various system components in the embodiments
described above should not be understood as requiring such
separation in all embodiments, and it should be understood that the
described program components and systems can generally be
integrated together in a single software product or packaged into
multiple software products.
[0121] Thus, particular embodiments of the subject matter have been
described. Other embodiments are within the scope of the following
claims. In some cases, the actions recited in the claims can be
performed in a different order and still achieve desirable results.
In addition, the processes depicted in the accompanying figures do
not necessarily require the particular order shown, or sequential
order, to achieve desirable results. In certain implementations,
multitasking and parallel processing may be advantageous.
[0122] All references and publications disclosed or discussed
herein are incorporated by reference in their entirety.
* * * * *