Embodiments of the present disclosure generally relate to neurostimulation (NS) systems, and more particularly to model-based programming guidance for implantation of spinal cord stimulation (SCS) systems.
NS systems are devices that generate electrical pulses and deliver the pulses to nervous tissue to treat a variety of disorders. For example, SCS has been used to treat chronic and intractable pain. Another example is deep brain stimulation, which has been used to treat movement disorders such as Parkinson's disease and affective disorders such as depression. While a precise understanding of the interaction between the applied electrical energy and the nervous tissue is not fully appreciated, it is known that application of electrical pulses depolarize neurons and generate propagating action potentials into certain regions or areas of nerve tissue. The propagating action potentials effectively mask certain types of pain transmitted from regions, increase the production of neurotransmitters, or the like. For example, applying electrical energy to the spinal cord associated with regions of the body afflicted with chronic pain can induce “paresthesia” (a subjective sensation of numbness or tingling) in the afflicted bodily regions. Inducing this artificial sensation replaces the feeling of pain in the body areas effectively masking the transmission of non-acute pain sensations to the brain.
Computational modeling of SCS, through coupled three dimensional (3-D) electrical field and nerve fiber kinetic models, can provide a tool for assessing the effectiveness of the SCS and/or a placement of the NS system within the patient. However, current modeling approaches commonly require commercial software packages that involve computationally-intensive steps, such as: obtain magnetic resonance imaging (MRI) of the patient; perform tissue segmentation on the medical images to create a 3-D spinal cord (SC) geometrical model; position the implanted leads within the SC model; specify stimulation contacts on the lead and set boundaries to contact-voltage/current condition; mesh the models; and in two stages solve for the electrical fields and activation regions in the dorsal column (DC) and dorsal root (DR) of the SC; and determine stimulation thresholds and activated dermatomal fiber zones. The process requires multiple software packages and specialized personnel to perform the tasks, which make conventional modeling approach difficult in the clinical setting.
However, some of the SCS systems available are not MRI-compatible, requiring MRI images to be taken prior to implant of the SCS system, and other modalities (e.g., X-rays, computed tomography (CT) scan) are needed to determine SCS lead position after implant. Moreover, detailed SC anatomy is difficult to ascertain with clinical MRI sequences, with dermatomal fiber tracts from such MRI images being difficult to visualize. Further, solving the computational model with the conventional approach is time-consuming, making this difficult to use in the clinical setting during an office visit or SCS implant. A need exists to overcome the shortcomings of traditional modeling methods.
In accordance with one embodiment, a method for modeling patient-specific spinal cord stimulation (SCS). The method may include acquiring impedance and evoked compound action potential (ECAP) signals from a lead positioned proximate to a spinal cord (SC). The lead may include at least one electrode. The method may include determining a patient-specific anatomical model based on the impedance and ECAP signals, and transforming a dorsal column (DC) map template based on a DC boundary of the patient-specific anatomical model. Optionally, the method may include acquiring additional impedance and ECAP signals. Each pair of impedance and ECAP signals may be acquired while the patient is in different patient postures or positions, such as sitting, standing, supine, or the like. Additionally or alternatively, the method may include detecting neural tissue damage based on the impedance and ECAP signals.
Further, the method may include extruding the patient-specific anatomical model along an SC axis to create a three dimensional model within a structural grid with material index assigned to each element. The geometry and location of the distal SCS leads are mapped into the 3D grids with their material properties assigned. The method may include solving a fully coupled extracellular and intracellular domain (e.g., a Bidomain Model) for electromagnetic fields in the extracellular domain and electrical propagation along neurons in the intracellular domain including dorsal column (DC), dorsal root (DR), or dorsal root ganglion (DRG).
Additionally, the method may include receiving patient responses at one or more pre-selected stimulation configurations and/or patient positions. The patient response may correspond to coverage zones on a body map as indicated on user interface at specific stimulation pulses for each of the pre-selected stimulation configurations.
Furthermore, the method may include using predetermined or user defined stimulation configurations for emitting a stimulation pulse from the at least one electrode, measuring a stimulation return signal in response to the stimulation pulse, and iteratively repeating the emitting and measuring operation of the method for each electrode to form a solution matrix corresponding to the coupled extracellular and intracellular domains.
In an embodiment, a system for modeling patient-specific spinal cord stimulation (SCS). The system may include a lead positioned proximate to a spinal cord (SC). The lead may include at least one electrode configured to acquire impedance and evoked compound action potential (ECAP) signals. The system may also include a system in communication with the lead. The system may include a memory device, a processor and a display. The system may be configured to determine a patient-specific anatomical model based on the impedance and ECAP signals, transform a dorsal column (DC) map template based on a DC boundary of the patient-specific anatomical model, and map the transformed DC map template to the patient-specific anatomical model. The system may also include algorithms to solve extracellular and intracellular domain electrical fields and propagation along neurons. The system may also include the user interfaces to collect patient responses and compare with model solutions.
In an embodiment, a method for differentiating spinal cord (SC) damage. The method may include emitting a stimulation waveform from at least one electrode of a lead. The lead being positioned proximate to a SC. The method may include acquiring impedance and evoked action compound action potential (ECAP) signals. The impedance and ECAP signals are based on the stimulation waveform. The method may further include selecting a first impedance and ECAP measurement and a second impedance and ECAP measurement from the impedance and ECAP signals. The second impedance and ECAP measurement is temporally separated from the first impedance and ECAP measurement. The method may include detecting SC tissue damage based on a difference between the second ECAP measurement and the first ECAP measurement, and between the second impedance measurement and first impedance measurement, and adjusting at least one therapy parameter to change the stimulation waveform or adjusting the position of the lead based on detection of the SC tissue damage.
While multiple embodiments are described, still other embodiments of the described subject matter will become apparent to those skilled in the art from the following detailed description and drawings, which show and describe illustrative embodiments of disclosed inventive subject matter. As will be realized, the inventive subject matter is capable of modifications in various aspects, all without departing from the spirit and scope of the described subject matter. Accordingly, the drawings and detailed description are to be regarded as illustrative in nature and not restrictive.
Embodiments described herein include a patient-specific spinal cord stimulation (SCS) programming guidance system 100 that may be used during SCS implantation, during office visits, and the like. The system 100 may be used by a clinician and/or patient to determine and/or select optimal SCS settings (e.g., amplitude, duration, frequency, type of stimulation pulses, selection of electrode configurations, or the like) that target a region of interest.
Optionally, the touchscreen 124 may be integrated with the display 122. The keyboard 126 (e.g., a typewriter keyboard 136) allows the user to enter data to the displayed fields, as well as interface with the RF subsystem 130. Furthermore, custom keys 128, for example, may turn on/off the system 100. The printer 112 prints copies of reports 140 for a physician to review or to be placed in a patient file, and the speaker 110 provides an audible warning (e.g., sounds and tones 142) to the user. The parallel I/O circuit 118 interfaces with a parallel port 144. The serial I/O circuit 120 interfaces with a serial port 146. The floppy drive 116 accepts diskettes 148. Optionally, the serial I/O port may be coupled to a USB port or other interface capable of communicating with a USB device such as a memory stick. The CD-ROM drive 114 accepts CD-ROMs 150.
The CPU 102 typically includes a microprocessor, a microcontroller, or equivalent control circuitry, designed specifically to control interfacing with the system 100 and with the NS system 200. The CPU 102 may include RAM or ROM memory, logic and timing circuitry, state machine circuitry, and I/O circuitry to interface with the NS system 200. The display 122 (e.g., may be connected to the video display 132). The display 122 displays various information related to the processes described herein. The touchscreen 124 may display graphic information relating to the NS system 200 (e.g., stimulation levels, stimulation waveforms, ECAP measurements) and include a graphical user interface.
The system 100 includes components 52-58 that may employ structures or aspects of various embodiments (e.g., systems and/or methods) discussed herein. In various embodiments, certain components (or operations) may be added, certain components may be combined, certain components may be performed simultaneously, certain components may be performed concurrently, certain components may be split into multiple components, certain components may be performed in a different order, or certain components may be re-performed in an iterative fashion. The components 52-58 may be a part of the CPU 102. Additionally or alternatively, the components 52-58 may be algorithms or instructions performed by the CPU 102 stored in memory (e.g., ROM 104, RAM 106, hard drive 108). Optionally, the components 52-58 may be separate modules in communication with the CPU 102. Optionally, one or more of the components 52-58 may be located external to the system 102. The system 102 may receive data from the components 52-58 via the RF subsystem 130.
Generally, the component 52 may be configured to generate lead position geometrical models, which illustrate a position of one or more electrodes on a lead with respect to a spinal cord (SC) tissue of interest based on measured impedance and evoked compound action potential (ECAP) signals. For example, the impedance signal is received by one or more electrodes, which is used to determine a distance between a select electrode and the SC tissue of interest. The geometrical models may additionally define anatomical parameters of interest (e.g., a thickness of cerebral spinal fluid) and landmark locations (e.g., position of the electrodes with respect to the SC tissue of interest) of the SC at select vertebral levels of interest based on the impedance and ECAP signals. For example, based on a structure of the ECAP signal (e.g., slope, peak to peak amplitude, peak latency, peak duration) received by one or more electrodes may be used to determine thickness of the cerebral spinal fluid. The component 52 may also be configured to use lookup tables (e.g., databases) of impedance and ECAP information from computational simulations, which are compared with the measured impedance and ECAP signals to determine a distance from the electrodes to dura/cerebral spinal fluid (CSF) and CSF thickness. Based on the measured impedance and ECAP signals, the component 52 may adjust a SC model template to match the patient's anatomy to generate the lead position geometric models.
The component 54 may be configured to transform (e.g., morph) a dermatomal map template of the SC dorsal column (DC) to fit a derived patient specific DC model based on an anatomy of a patient. Optionally, the transforming of the dermatomal map template by the component 54 may include meshing the dermatomal zone map template and deforming the mesh to adjust to the outer boundary shape of the patient's DC anatomy, such as, using Delaunay triangulation. The outer boundary shape of the patient's DC anatomy may be determined based on the anatomical parameters of interest and landmark locations determined from the lead position geometric model. Additionally or alternatively, the component 54 may include shifting pixels of the dermatomal zone map template according to a distance the nodes moved from the original mesh to the deformed mesh during the transformation operation. Optionally, variations may be added to the transforming or pixel movement to account for uncertainties in the dermatomal zone size and boundary locations.
Based on the lead position geometric model generated from the component 52 and the derived patient specific DC model from the component 54, the component 56 may be configured to generate a two dimensional (2D) and/or three dimensional (3D) SC model. The SC model may be used to compute electrical fields and neural transmembrane potential in the DC, dorsal root (DR), or dorsal root ganglion (DRG). Further the SC model may be used to determine the activation regions and activation of mapped dermatomal zones for corresponding electrodes. The SC model may include a fully coupled extracellular and intracellular domain for electromagnetic fields in the extracellular domain and electrical propagation along neurons in the intracellular domain. A technical effect of the system 100 eliminates the use of commercial simulation software using conventional finite element analysis (FEA) and finite difference (FD) methods. An additional technical effect of the system 100 may be to perform a discretization from a 2D SC model to a 3D SC model with the SCS lead that is much faster than other meshing methods used in FEA.
The component 58 may be configured to apply pre-selected or selected electrode stimulation configurations received from the clinician and/or patient. Additionally or alternatively, the component 58 may be configured to test and refine the SC model based on feedback by the patient. Optionally, the component 58 may include a user-interface, which may include a graphical user interface, for entering stimulation parameters with a display of the corresponding activation regions and activated dermatomal zones.
The RF subsystem 130 includes a central processing unit (CPU) 152 in electrical communication with RF circuitry 154, which may communicate with both memory 156 and an analog out circuit 158. The analog out circuit 158 includes communication circuits to communicate with analog outputs 164. The system 100 may wirelessly communicate with the NS system 200 using a telemetry system. Additionally or alternatively, the system 100 may wirelessly communicate with the NS system 200 utilize wireless protocols, such as Bluetooth, Bluetooth low energy, WiFi, MICS, and the like. Alternatively, a hardwired connection may be used to connect the system 100 to the NS system 200.
Optionally, the system 100 may transmit the stimulation database request to an implantable pulse generator (IPG) 250 (
The NS system 200 includes the IPG 250 that is adapted to generate electrical pulses for application to tissue of a patient. The IPG 250 typically comprises a metallic housing or can 259 that encloses a controller 251, pulse generating circuitry 252, a charging coil 253, a battery 254, a far-field and/or near field communication circuitry 255, battery charging circuitry 256, switching circuitry 257, sensing circuitry 258, memory 261, and the like. The controller 251 typically includes a microcontroller or other suitable processor for controlling the various other components of the device. Software code may be stored in memory 261 of the IPG 250 or integrated with the controller 251 for execution by the microcontroller or processor to control the various components of the device.
The IPG 250 may comprise a separate or an attached extension component 270. If the extension component 270 is a separate component, the extension component 270 may connect with a “header” portion of the IPG 250 as is known in the art. If the extension component 270 is integrated with the IPG 250, internal electrical connections may be made through respective conductive components. Within the IPG 250, electrical pulses are generated by the pulse generating circuitry 252 and are provided to the switching circuitry 257. The switching circuitry 257 connects to outputs of the IPG 250. Electrical connectors (e.g., “Bal-Seal” connectors) within the connector portion 271 of the extension component 270 or within the IPG header may be employed to conduct various stimulation pulses. The terminals of one or more leads 210 are inserted within the connector portion 271 or within the IPG header for electrical connection with respective connectors. Thereby, the pulses originating from the IPG 250 are provided to the one or more leads 210. The pulses are then conducted through the conductors of the lead 210 and applied to tissue of a patient via one or more electrodes (e.g., array of electrodes 211). Any suitable known or later developed design may be employed for connector portion 272.
The lead 210 is connected to a flat, thin, paddle structure 214 and connect thereto in a general longitudinal alignment with the length of the paddle structure 214. The paddle structure 214 may be formed from a medical grade, substantially inert material, for example, polyurethane, silicone, or the like. A front surface or face 212 of the paddle structure 214 is depicted in
Each of the electrodes 211 are mutually separated by non-conducting or insulative material of the paddle, which electrically isolate each electrode 211 from adjacent electrodes 211. The non-conducting material may include one or more insulative materials and/or biocompatible materials to allow the paddle structure 214 and lead 210 to be implantable within the patient. Non-limiting examples of such materials include polyimide, polyetheretherketone (PEEK), polyethylene terephthalate (PET) film (also known as polyester or Mylar), polytetrafluoroethylene (PTFE) (e.g., Teflon), or parylene coating, polyether bloc amides, polyurethane.
The electrodes 211 may be formed of non-corrosive, highly conductive material. For example, stainless steel, MP35N, platinum, platinum alloys, or the like. The electrodes 211 may be set to function as cathodes, anodes or set to a high-impedance state for a given pulse according to the pulses generated from the IPG 250. The electrodes 211 may be configured to emit the pulses in an outward radial direction proximate to or within a stimulation target. The electrodes 211 may also be configured to acquire electrical potential measurements (e.g., voltage, current) or electrical signals for the sensory circuit 258, such as evoked compound activation potentials (ECAP) emitted from the stimulation target. ECAP signals may be generated by neuronal transmembrane currents of neurons activated following or in response to a stimulation pulse from one or more of the electrodes 211.
Optionally, the IPG 250 may have more than one lead 210 connected via the connector portion 271 of the extension component 270 or within the IPG header. Additionally or alternatively, the electrodes 211 of each lead 210 may be configured separately to emit current pulses or measure electrical signals emitted from and/or proximate the stimulation target.
It should be noted that in other embodiments the electrodes 211 may be in various other formations or structures. For example, the electrodes 211 may be in the shape of a ring such that each electrode 211 continuously covers the circumference of the exterior surface of the lead 210 to form a percutaneous lead structure. Each of the ring electrodes 211 are separated by non-conducting rings, which electrically isolate each electrode 211 from an adjacent electrode 211. In another example, the electrodes 211 may be in the shape of a split or non-continuous ring such that the pulse may be directed in an outward radial direction adjacent to the electrodes 211. Further examples of a fabrication process of the electrodes 211 is disclosed in U.S. patent application Ser. No. 12/895,096, entitled, “METHOD OF FABRICATING STIMULATION LEAD FOR APPLYING ELECTRICAL STIMULATION TO TISSUE OF A PATIENT,” which is expressly incorporated herein by reference.
The lead 210 may comprise a lead body 272 of insulative material about a plurality of conductors within the material that extend from a proximal end of lead 210 (proximate to the IPG 250) to its distal end (proximate to the paddle structure 214). The conductors electrically couple a plurality of the electrodes 211 to a plurality of terminals (not shown) of the lead 210. The terminals are adapted to receive electrical pulses and the electrodes 211 are adapted to apply the pulses to the stimulation target of the patient. Also, sensing of physiological signals may occur through the electrodes 211, the conductors, and the terminals. It should be noted that although the paddle structure 214 of the lead 210 is depicted with a five by four array of electrodes 211, in other embodiments, the lead 210 may be connected to any suitable number of electrodes 211 (e.g., an array with more electrodes 211 than shown in
Although not required for all embodiments, the lead body 272 of the lead 210 may be fabricated to flex and elongate upon implantation or advancing within the tissue (e.g., nervous tissue) of the patient towards the stimulation target and movements of the patient during or after implantation. By fabricating the lead body 272, according to some embodiments, the lead body 272 or a portion thereof is capable of elastic elongation under relatively low stretching forces. Also, after removal of the stretching force, the lead body 272 may be capable of resuming its original length and profile. For example, the lead body may stretch 10%, 20%, 25%, 35%, or even 50% at forces of about 0.5, 1.0, and/or 2.0 pounds of stretching force. Fabrication techniques and material characteristics for “body compliant” leads are disclosed in greater detail in U.S. Provisional Patent Application No. 60/788,518, entitled “Lead Body Manufacturing,” which is expressly incorporated herein by reference.
For implementation of the components within the IPG 250, a processor and associated charge control circuitry for an IPG is described in U.S. Pat. No. 7,571,007, entitled “SYSTEMS AND METHODS FOR USE IN PULSE GENERATION,” which is expressly incorporated herein by reference. Circuitry for recharging a rechargeable battery (e.g., battery charging circuitry 256) of an IPG (e.g., the IPG 250) using inductive coupling and external charging circuits are described in U.S. Pat. No. 7,212,110, entitled “IMPLANTABLE DEVICE AND SYSTEM FOR WIRELESS COMMUNICATION,” which is expressly incorporated herein by reference.
An example and discussion of “constant current” pulse generating circuitry (e.g., pulse generating circuitry 252) is provided in U.S. Patent Publication No. 2006/0170486 entitled “PULSE GENERATOR HAVING AN EFFICIENT FRACTIONAL VOLTAGE CONVERTER AND METHOD OF USE,” which is expressly incorporated herein by reference. One or multiple sets of such circuitry may be provided within the IPG 250. Different pulses on different electrodes 211 may be generated using a single set of the pulse generating circuitry 252 using consecutively generated pulses according to a “multi-stimset program” as is known in the art. Complex pulse parameters may be employed such as those described in U.S. Pat. No. 7,228,179, entitled “Method and apparatus for providing complex tissue stimulation patterns,” and International Patent Publication Number WO 2001/093953 A1, entitled “NEUROMODULATION THERAPY SYSTEM,” which are expressly incorporated herein by reference. Alternatively, multiple sets of such circuitry may be employed to provide pulse patterns (e.g., tonic stimulation waveform, burst stimulation waveform) that include generated and delivered stimulation pulses through various electrodes of one or more leads 211 as is also known in the art. Various sets of parameters may define the pulse characteristics and pulse timing for the pulses applied to the various electrodes 211 as is known in the art. Although constant current pulse generating circuitry is contemplated for some embodiments, any other suitable type of pulse generating circuitry may be employed such as constant voltage pulse generating circuitry.
The sensing circuitry 258 may measure an electric potential (e.g., voltage, current) over time of the stimulation target or proximate tissue through at least one of the electrodes 211 that is proximate to the stimulation target. The electric potential (EP) measurements may correspond to the ECAP signal generated by the stimulation target in response to pulses emitted from the electrodes 211. For example, the sensing circuitry 258 may measure an ECAP signal from an Aβ sensory fiber or neural tissue of the SC processed from the EP sensed from one or more of the electrodes 211 on the lead 210. The sensing circuitry 258 may include amplifiers, filters, analog to digital converters, memory storage devices (e.g., RAM, ROM), digital signal processors, and/or the like. Optionally, the sensing circuitry 258 may store the EP in the memory 261.
The system 100 may be implemented to charge/recharge the battery 254 of the IPG 250 (although a separate recharging device could alternatively be employed), to access the memory 261, and to program the IPG 250 on the pulse specifications while implanted within the patient. Although, in alternative embodiments separate programmer devices may be employed for charging and/or programming the NS system 200. The system 100 may be a processor-based system that possesses wireless communication capabilities. Software may be stored within a non-transitory memory of the system 100, which may be executed by the processor to control the various operations of the system 100. Optionally, a “wand” 265 may be electrically connected to the system 100 through suitable electrical connectors (not shown). The electrical connectors may be electrically connected to a telemetry component 266 (e.g., inductor coil, RF transceiver) at the distal end of wand 265 through respective wires (not shown) allowing bi-directional communication with the IPG 250.
The user may initiate communication with the IPG 250 by placing the wand 265 proximate to the NS system 200. Preferably, the placement of the wand 265 allows the telemetry system of the wand 265 to be aligned with the far-field and/or near field communication circuitry 255 of the IPG 250. The system 100 preferably provides one or more user interfaces 268 (e.g., graphical user interface, display, touchscreen 124, keyboard 126, mouse, buttons, or the like) allowing the user to operate (e.g., adjust the pulse settings) the IPG 250. The system 100 may be controlled by the user (e.g., doctor, clinician) through the user interface 268 allowing the user to interact with the IPG 250. The user interface 268 may permit the user to move electrical stimulation along and/or across one or more of the lead(s) 210 using different electrode 211 combinations, for example, as described in U.S. Patent Application Publication No. 2009/0326608, entitled “METHOD OF ELECTRICALLY STIMULATING TISSUE OF A PATIENT BY SHIFTING A LOCUS OF STIMULATION AND SYSTEM EMPLOYING THE SAME,” which is expressly incorporated herein by reference. Optionally, the user interface 268 may permit the user to designate which electrodes 211 are to stimulate (e.g., emit current pulses, in an anode state, in a cathode state) the stimulation target, to measure the ECAP or impedance (e.g., connecting to the sensing circuitry 158) resulting from the current pulses, remain inactive (e.g., floating), or the like. Additionally or alternatively, the system 100 may access or download the electrical measurements from the memory 261 acquired by the sensing circuitry 258.
Also, the system 100 may permit operation of the IPG 250 according to one or more spinal cord stimulation (SCS) programs or therapies to treat the patient. Each SCS program may include one or more sets of stimulation parameters of the pulse including pulse amplitude, stimulation level, pulse width, pulse frequency or inter-pulse period, pulse repetition parameter (e.g., number of times for a given pulse to be repeated for respective stimset during execution of program), biphasic pulses, monophasic pulses, etc forming a drive signal or stimulation waveform. The IPG 250 may modify its internal parameters in response to the control signals from the system 100 to vary the stimulation characteristics of the stimulation pulses transmitted through the lead 210 to the tissue of the patient. NS systems, stimsets, and multi-stimset programs are discussed in PCT Publication No. WO 01/93953, entitled “NEUROMODULATION THERAPY SYSTEM,” and U.S. Pat. No. 7,228,179, entitled “METHOD AND APPARATUS FOR PROVIDING COMPLEX TISSUE STIMULATION PATTERNS,” which are expressly incorporated herein by reference.
In various embodiments, certain steps (or operations) may be omitted or added, certain steps may be combined, certain steps may be performed simultaneously, certain steps may be performed concurrently, certain steps may be split into multiple steps, certain steps may be performed in a different order, or certain steps or series of steps may be re-performed in an iterative fashion. Furthermore, it is noted that the following is just one possible method of determining a patient-specific anatomical model based on impedance and evoked compound action potential. It should be noted, other methods may be used, in accordance with embodiments herein.
One or more methods may (i) acquire impedance and evoked compound action potential (ECAP) signals from a lead positioned proximate to a spinal cord (SC) and (ii) determine a patient-specific anatomical model based on the impedance and ECAP signals.
Beginning at 302, a lead 402 (shown in
The lead 402 is connected via a lead body 410 to an IPG 412 (e.g., the IPG 250). Optionally, the lead 402 may be positioned at a selected vertebral level, which may be used to select an SC model template, as further described herein.
At 304, a drive signal is emitted from at least one electrode 511 of the lead 402. For example, the drive signal may be generated from the IPG 412, such as from generating circuitry (e.g., the generating circuitry 252), and conducted to at least one of the electrodes 511 via switching circuitry (e.g., the switching circuitry 257) and the lead 402. The drive signal may represent a current pulse (e.g., a monophasic pulse) or a series of current pulses (e.g., a biphasic pulse, tri-phasic pulses), a sinusoidal waveform, a burst waveform, and/or the like which are emitted from at least one of the electrodes 511 with a predetermined amplitude and pulse width. Additionally or alternatively the drive signal may be a voltage pulse with a predetermined amplitude and pulse width. The drive signal is used by the IPG (e.g., IPG 412, IPG 250) and/or the system 100 to determine anatomical parameters of interest (e.g., CSF thickness) and landmark locations (e.g., position of the electrodes 511 with respect to SC tissue of interest) of the SC based on impedance (at 306) and ECAP signals (at 308) resulting from the drive signal. Additionally or alternatively, in connection with
Optionally, the drive signal may be a series of pulses. At least one subset of the series of pulses may be shaped (e.g., pulse width, amplitude, frequency) to facilitate measurement of the impedance signal. Additionally or alternatively, at least another subset of the series of pulses may be shaped to facilitate measurement of the ECAP signal. It should be noted that in at least one embodiment the drive signal may be shaped to facility measurement of both the impedance signal and the ECAP signal.
At 306, an impedance signal is acquired from the lead based on the drive signal. The magnitude of the impedance signal is based on electrical properties, such as the conductivity, of the SC tissue proximate to the electrode(s) emitting the drive signal and/or acquiring the impedance signal. For example, if the proximate SC tissue is primarily epidural fat the impedance signal will have a high magnitude due to the low conductivity of epidural fat. In another example, if the proximate SC tissue is primarily CSF the impedance signal will have a low magnitude due to the high conductivity of the CSF. The impedance signal may be measured and stored in the memory 261. Based on the impedance signals, the system 100 and/or controller 251 may determine the electrical conductivity of the CSF and/or epidural fat layer of the SC tissue.
The impedance signal may correspond to a discrete impedance value or a calculated impedance over time during the drive signal. For example, the drive signal may be a current pulse emitted from two or more electrodes 511 (e.g., one electrode is in a cathode state, one electrode is in an anode state). Alternatively, the drive signal may be emitted from at least one of the electrodes 511 and the housing (e.g., the can 259). A voltage may be measured from the two or more electrodes 511 delivering the drive signal by the sensing circuitry 258 and/or the controller 251 comparing the two voltage potentials of the two or more electrodes 511 during the drive signal. Additionally or alternatively, the voltage may be measured by one or more alternative electrodes 511 (e.g., not emitting the drive pule) during the drive signal and received by the sensing circuitry 258. Using the measured voltage and the drive signal, the controller 251 can determine the impedance signal (variable Zc), using Equation 1, by dividing the measured voltage (variable Vmeasured) by the known stimulation current (variable Istim) for the drive signal.
Z
c
=V
measured
/I
stim (Equation 1)
In another example, the drive signal may be a voltage pulse emitted from two or more electrodes 511. Alternatively, the drive signal may be emitted from at least one of the electrodes 511 and the housing (e.g., the can 259). A measured current may be measured from the at least one electrode 511 emitting the drive signal or by an alternative electrode 511 (e.g., not emitting the drive signal) during the drive signal, and received by the sensing circuitry 258. Using the measured current and the known stimulation voltage value for the drive signal, the controller 251 can determine the impedance signal (variable Zc) using, Equation 2, by dividing the stimulation voltage (variable Vstim) of the drive signal by the measured current (variable Imeasured).
Z
c
=V
stim
/I
measured (Equation 2)
It should be noted in various embodiments, one or more of the electrodes 511 may be used for emitting a stimulation waveform (e.g., the drive signal) and measuring of the impedance signal (or ECAP signal at 308 generated by the corresponding neurons of the stimulation target). Additionally or alternatively, a subset of the electrodes 511 may be used for emitting the drive signal and another subset of the electrodes 511 may be used to measure the impedance signal (or ECAP signal at 308).
Returning to
At 310, a distance (e.g., 702-710) between the at least one electrode and a SC tissue of interest is determined based on the impedance signal. The determination at 310 is further described in connection with
The electrodes 511a-e are each a distance 702-710, respectively, away from the dura layer 514 surrounding the CSF 508. At different points in time (e.g., t0, t1, t2, t3, t4) a corresponding electrode 511a-e emits the drive signal (e.g., a monophasic current pulse), and a corresponding electrode 511f-j senses an electrical response (e.g., sensed voltage, sensed current) in response to the drive signal for the impedance signal. It, should be noted that in other embodiments, the electrode sensing the electrical response and emitting the drive signal may be interchanged (e.g., electrode 511f emits the drive signal, electrode 511a senses the electrical response), located in different columns on the electrode array 511, or measured from the same electrodes (e.g., electrode 511a-e) emitting the drive signal. Optionally, more than one electrode may sense the electrical response. Additionally or alternatively, the same electrode may sense the electrical response for more than one electrode emitting the drive signal. Optionally, more than one electrode may emit the drive signal (e.g., biphasic pulse, tri-phasic pulse).
The impedance signal is plotted on the graphical illustration 800. A vertical axis 802 represents magnitude (e.g., ohms), and the horizontal axis 804 is time. The impedance signal 806 corresponds to a drive signal emitted from the electrode 511a, the impedance signal 807 corresponds to a drive signal emitted from the electrode 511b, the impedance signal 808 corresponds to a drive signal emitted from the electrode 511c, the impedance signal 809 corresponds to a drive signal emitted from the electrode 511d, the impedance signal 810 corresponds to a drive signal emitted from the electrode 511e.
The magnitude of the impedance signal 806-810 corresponds to a distance (e.g., 702-710) from the emitting electrode (e.g., 511a-e) to the dura layer 514 or CSF 508. For example, larger distances (e.g., 702, 710) yield a larger impedance (e.g., 806, 810) because the lead is further from the highly conductive CSF 508. Alternatively, relatively smaller distances (e.g., 706) yield a comparatively smaller impedance (e.g., 808) since the lead is closer to the highly conductive CSF 508.
The impedance signals 806-810 may be compared by the system 100 and/or the controller 251 to an impedance distance database to determine the distances 702-710. The impedance distance database may be stored in memory (e.g., memory 261, ROM 104, RAM 106, a hard drive 108) and based on traditional three dimensional electromagnetic computation simulations or priori information of experimentally measured impedance measurements and electrode positioned across a patient population. The impedance distance database includes a collection of candidate impedances with associated distances. The impedance distance database may be used as a look up table to match a corresponding impedance signal magnitude with a distance. Additionally, the values of the distances and/or corresponding impedance signals magnitudes of the impedance distance database may be dependent on the structure of the lead 402 and the electrodes 511. For example, the impedance distance database may have different sets of values of the distances. A first set of distance values correspond to impedance signals for the paddle structure 510, and a second set of distance values for another structure where the electrodes 511 are in the shape of a ring such that each electrode 511 continuously covers the circumference of the exterior surface of the lead 402. It should be noted that the impedance distance database may have more than two sets of values of the distances with corresponding impedance signal magnitudes for various other electrode 511 and/or lead 402 structures as described herein.
In at least one embodiment, for each impedance signal 806-810, the system 100 and/or the controller 251 may find a matching impedance signal magnitude in the impedance distance database (or interpolation with closest data points) with a corresponding distance from an electrode contact to the dura layer 514 or CSF 508. For example, the system 100 may compare the impedance signal magnitude for the impedance signal 808 to the impedance distance database. The impedance distance database may list the impedance magnitude with a corresponding distance of 0.2 mm. Based on the listed distance, the system 100 may determine that the distance 706 is 0.2 mm from the electrode 511c to the SC tissue of interest (e.g., the dura layer 514).
Based on the combination of distances 704-710, the system 100 may determine an overall position of the paddle based on y and x axes 701 and 703 centered proximate to the paddle 510 on the SC tissue of interest. Optionally, the x and y axes 703, 701 may be centered on one of the electrodes, for example, the center column electrode, the electrode 511c, on the paddle structure. The y axis 701 corresponds to a dorsal/ventral position (e.g., dorsal position is based away from the x axis 703, ventral position is based towards the x axis 703), for example, based from the distances 702-710. The x axis 703 corresponds to a medial/lateral position (e.g., medial position is on the y axis 701, lateral position is based away from the y axis 701 in either a left or right direction along the x axis 703).
For example, the system 100 and/or the controller 251 may determine the medial/lateral position of the paddle structure 510 based on an interrelation of the magnitude of the impedance signals or morphology (e.g., changes in slope, peaks, rising/falling edges) of a line graph (e.g., an impedance curve 1220) formed from the impedance signals as further described in connection to
Based on the falling edge 1022 and the interrelation of the magnitudes of the impedance signals 1006-1007 and 1008-1010 being approximately the same, respectively, the system 100 and/or the controller 251 may determine that the falling edge 1022 corresponds to a portion of the paddle structure 510, such as electrodes 511a-b, is located remotely from the SC tissue of interest. Since the falling edge 1022 decreases in magnitude in the direction of the arrow 1024, the system 100 and/or the controller 251 can determine that the paddle structure 510 is shifted laterally.
Additionally, based on the combination of distances of the electrodes 511 from the SC tissue of interest, the system 100 and/or the controller 251 can determine an angle of the paddle structure 510 relative to the SC tissue of interest. For example, the system 100 and/or controller 251 may determine that the paddle structure 510 is positioned at an angle based on a slope of an impedance curve 1220, as described further in connection to
Returning to
The system 100 and/or the controller 251 may compare one or more ECAP signals features representing the morphology of the ECAP signal 802 such as a slope, a maximum and/or minimum amplitude, time delay from onset of stimulation to peak ECAP amplitude, or the like, with an ECAP database to determine the thickness characteristic of the SC tissue of interest at or near the corresponding electrode (e.g., the electrode 711c-d). Based on the combination of thickness measurements for each electrode, the system 100 and/or the controller 251 can determine an overall thickness of the SC tissue of interest. The ECAP database may be stored in memory (e.g., memory 261, ROM 104, RAM 106, a hard drive 108) and based on traditional three dimensional electromagnetic computation simulations or priori information of experimentally measured ECAP signals and electrode positioned across a patient population. The ECAP database includes a collection of candidate ECAP signal features with associated thickness characteristics of the SC tissue of interest. The ECAP database may be used as a look up table to match one or more ECAP signal features with thickness characteristic of SC tissue of interest.
In at least one embodiment, a peak to peak 1320 measurement from the ECAP signal 1302 may be measured by the system 100 and/or the controller 251, and compared with the ECAP database to find a matching peak to peak ECAP signal in the ECAP database (or interpolation with closest data points) with a corresponding CSF thickness. For example, the system 100 may compare the peak to peak 1320 to the ECAP database. The ECAP database may list a peak to peak matching the peak to peak 1320 with a corresponding CSF thickness of 2.8 mm. Based on the listed CSF thickness, the system 100 may determine that the CSF thickness 730 is 2.8 mm.
At 314, the method 300 adjusts an SC model template 1400 based on the distance and thickness characteristic.
After the selection of the SC model template 1400, the system 100 and/or the controller 251 may adjust the SC model template 1400 based on the distance and thickness characteristic operations determined at 310 and 312.
For example, the SC model template 1400 has a CSF thicknesses 1406 based on a position of a lead 1412. The thickness characteristic measurements determined at 312 correspond to CSF thicknesses 1506. The system 100 may determine a ratio along the x and y axes 1402, 1404 based on the CSF thicknesses 1506 compared with the CSF thicknesses 1406. The ratio is used to adjust the SC model template 1400 to form the adjusted SC model template 1500. The system 100 may further adjust the position (e.g., medial/lateral, distance from the dura layer, angle relative to the SC tissue) of the lead 1412 in the SC model template 1400. The position of the lead 1412 is adjusted based on the distance characteristics determined at 310 with respect to the adjusted SC model template 1500 to determine the adjusted position of the lead 1412.
Optionally, the SC model template 1400 may also be adjusted based on one or more medical images acquired by a medical imaging system (e.g., magnetic resonance imaging, computed tomography, X-ray, or the like). The one or more medical images of the patient may be imported by the method 300 for determination of patient anatomy and electrode position. For example, segmentation of pre-operative MRI could provide detailed information about SC tissue geometric structure. In another example, once the lead 402 is implanted within the patient, post-operative CT or X-ray may be used to acquire positions of the electrodes 511. In another example, for an MRI compatible SCS system, both SC anatomy and location of the electrodes 511 may be obtained from post-operative MRI.
The medical images may also be used for image-based calibration of the impedance distance and/or ECAP databases. For example, the system 100 and/or controller 251 may use the impedance signal for investigating SC tissue properties, such as the electrical conductivity of the CSF layer or of the epidural fat, to be used in building and/or updating the impedance distance database. In another example, the impedance distance database and/or ECAP database may be based on nominal values of tissue conductivity based on priori information. Based on the medical images, the impedance distance and ECAP signals corroborated with positioned information of the electrodes 511 and/or SC tissue properties are used for minor adjustment of conductivity away from these nominal values.
Additionally or alternatively, additional impedance and ECAP signals may be acquired. Each set of impedance and ECAP signals may correspond to different patient postures or positions, such as sitting, standing, supine, or the like. For example, a first set of impedance signal(s) and ECAP signal(s) may be acquired by the sensing circuitry 258 while the patient is in a first position, such as in a sitting position. When the patient adjusts to a second position, such as a standing position, the second set of impedance and ECAP signals may be acquired and so on for all desired postures. It should be noted that in other embodiments more than two different patient positions may be used. Optionally, the impedance and ECAP signals corresponding to multiple positions of the patient may be stored in a position measurement database on memory (e.g., memory 261, ROM 104, RAM 106, hard drive 108). Corresponding models (e.g., patient-specific anatomical models) at different postures may be built with associated changes in electrode location and SC anatomy. Through stimulation configurations with each patient-specific anatomical model, SCS programming guidance with each posture may be established and stored in the NS system 200 and/or the SCS programming guidance system 100.
A technical effect of acquiring impedance and ECAP signals at multiple positions allow the patient-specific anatomical model to be configured for different postures for SCS programming and stored in memory (e.g., ROM 104, RAM 106, hard drive 108). Allowing, patients to implement the desired SCS programming when using a graphical user interface of the system 100 when the posture or position of the patient changes.
A technical effect of the method 1600 is a transformed DC map template (e.g., a transformed DC map template 2000 of
One or more methods may (i) generate a meshed DC map template and (ii) translate the mesh nodes of the meshed DC map template based on a patient-specific anatomical model.
Beginning at 1602, the method 1600 automatically traces outer DC boundaries 1502 and 1704 of a DC map template (e.g., dermatomal zone map template is shown in
At 1604, the method 1600 transforms the boundaries of the DC map template and the patient-specific anatomical model. The system 100 may translate and/or scale the DC boundaries 1502 and 1704 of the adjusted SC model template 1500 and the DC map template 1700 so each DC boundary 1502 and 1704 are overlapping. For example, the minimum and maximum value of the DC boundaries 1502 and 1704 along an x and y axes.
At 1606, the method 1600 generates a meshed DC map template 1800.
At 1608, the method 1600 translates mesh nodes of the meshed DC map template based on the patient-specific anatomical model. The locations of the mesh nodes 1802 are shifted, as shown in
Once the descriptors are defined, the system 100 may adjust the background triangulation to incorporate the desired modification of the DC boundary 1704 to that of the DC boundary 1502 of the adjusted SC model template 1500 to form a transformed DC map template 2000 as shown in
Optionally, variations may be added during the translation operation at 1608 to account for uncertainties in dermatomal zone sizes and boundary locations. For example, artificial noise may be added during the movement such that a total distance moved by each mesh node 1802 may be adjusted by a small value using a random number generator centered on zero within the system 100. Additionally or alternatively, variations may be added by adjusting the spacing of the mesh nodes 1802 or spacing of the sample points.
At 1610, the method 1600 calculates a shifted distance of the mesh nodes 1802. For example, the system 100 may calculate distance based on a position of one or more mesh nodes 1802 from before the translation operation at 1608 to a position of the one or more mesh nodes 1802 after the translation operation at 1608. Optionally, the method 1600 may include applying an image filter, such as an averaging image filter, based on a user-specified pixel size. The user specified pixel size may be received by the system 100 through the user interface. Additionally or alternatively, the user specified pixel size may be predetermined and stored in memory (e.g., ROM 104, RAM 106, hard drive 108). The averaging image filter may be used to blend pixel locations of the transformed DC map template 2000.
Beginning at 2101, the method 2100 maps the transformed DC map template 2000 to the patient-specific anatomical model. For example, the transformed DC map template 2000 and the adjusted SC model template 1500 may be one or more 2D geometries based on an X-Y plane (e.g., defined by the axes 703 and 701) based on a row of electrodes (e.g., 812) of the electrode array 511. The system 100 may define a coordinate origin (e.g., the coordinate origin 1403) aligned for both the transformed DC map template 2000 and the adjusted SC model template 1500 located at the distal end (e.g., the distal end 216) of the lead 402 at the Z axis (e.g., the axis 516), centered on the adjusted SC model template 1500 (e.g., for the x and y direction). Optionally, a file may be generated by the system 100 containing all of the dimensions of the electrodes 511 (e.g., size of the electrodes, spacing between electrodes, electrode layout design, and the like).
At 2102, the method 2100 extrudes the patient-specific anatomical model along the Z axis defined by the SC to create a three dimensional (3D) grid. The system 100 may extrude along the Z axis (e.g., the axis 516) defined by the SC along the lead 402 (or paddle structure 510) to create the 3D patient-specific anatomical structure. Optionally, the 3D patient-specific anatomical structure may be bounded by a 3D rectangular box to enclose the electrode array(s) 511 and SC tissue. The anatomical contours of the adjusted SC model template 1500 in the X-Y plane may be mapped into a 2D grid in the X-Y plane with pre-determined dx and dy increments. A material index (e.g. based on electrical conductivity) for the various SC tissues may be assigned by the system 100 for each grid point. The 2D grids with material index are duplicated along the Z axis to create a 3D grid within the rectangular box of the model. Optionally, the electrodes 511 and paddle structure 510 or lead 402 is mapped into the 3D grids with a corresponding material index.
At 2104, the method 2100 emitting a stimulation pulse from the at least one electrode (e.g., the electrode 511a). The stimulation pulse may be a simulated unit pulse (e.g., constant current or voltage) from at least one electrode within the 3D patient-specific anatomical model.
At 2106, the method 2100 measures an electric field generated in response to the stimulation pulse based on the 3D patient-specific anatomical model. For example, the electric field generated may be a simulated response based from a bidomain model (e.g., extracellular and intracellular domain) as corresponding to a position of the at least one electrode within the 3D patient-specific anatomical model.
At 2108, the method iteratively repeats the emitting and measuring operations at 2104 and 2106 for each electrode of the 3D patient-specific anatomical model to form a solution matrix corresponding to a bidomain model. The solution matrix may be formed by the system such that the matrix deposition/factoring would allow post processing for any combinations of stimulation contact configurations based on the superposition principle. For example, based on the superposition principle of the measured stimulation return signals for each individual electrodes 511 (e.g., 511a, 511b, 511c, 511d, 511e) within a selected combination is the sum of the stimulation return signals. A technical effect of the solution matrix is that the need to solve for each contact configuration tested is not needed, and therefore, greatly speeds up the calculation for the activated dermatomal zones.
Optionally, the method 2100 may include adjusting the patient-specific anatomical model based on the stimulation return signals (e.g., actual dermatomal zone activations) in response to emitted stimulation signals from at least one of the electrodes 2208. The stimulation return signal may include sensed ECAP signals 2206 made across a subset of electrodes 2208 (e.g., the electrodes 511 within the box) of the lead 402, and may be used to determine the activated dermatomal zones within the DC.
Additionally or alternatively, a cross impedance may be determined by the system 100 and/or controller 251 based on the sensed ECAP signals 2206. For example, the cross impedance (variable Zc) may be equal to the ECAP signal (variable Vsense) divided by stimulation current of the SCS signal (variable Istim) as noted in Equation 3. This cross impedance may be used concurrently with the impedance signal measured at 306. Optionally the cross impedance may be used instead of the impedance signal measured at 306, for example, if the electrodes 511 have a high polarized impedance.
Z
c
=V
sense
/I
stim (Equation 3)
The technical effect of the solution matrix used concurrently with the 3D patient-specific anatomical model forms the bidomain model. The bidomain model may be a finite element analysis (FEA) solution coupled with the biophysical cellular model of neural elements to determine the effect of SCS on the transmembrane potentials of nerve fibers in the SC. This allows the clinician to detect and locate sites of neural activation resulting from the SCS emitted from one or more of the electrodes 511. The biophysical cellular model may be a predetermined template included within the 3D patient-specific anatomical model used to define electrical and geometrical properties of a neuron for an applied stimulation from at least one of the electrodes 511. For example, the biophysical model may determine the response of a neuron based on an electrical characteristic of the applied stimulation pulse (e.g., amplitude, frequency, pulse width). The biophysical cellular model may include transmembrane ion channels, transmembrane capacitance, intracellular conductance, applied stimulation fields or the like.
In at least one embodiment, the graphical user interface of the system 100 may include graphical icons, scroll bars, buttons, and the like which may receive or detect user or touch inputs 134 for the system 100 when selections are made by the user. For example, the graphical user interface may allow the user to test the neural response to various pre-selected electrode configurations for input validation.
Optionally, the GUI 2400 may include use of a predetermined or user-defined trolling algorithm that steps through a selected set of electrode configurations. In at least one embodiment a GUI 2500, as shown in
Additionally or alternatively, a GUI 2600 of the system 100, as shown in
Additionally or alternatively, in connection with a method 2650 shown in
Neural tissue damage may occur during implantation of the lead 402 into the patient, or stimulation of the neurons of the stimulation target in response to the drive signal, and/or the like. For example, pressure applied to the SC by the lead 402 during implantation damages the corresponding neural tissue. Alternatively, stimulation may cause neuronal tissue damage through over-activation of neurons or generation of potentially harmful electrochemical reactions. Damaged neurons may not generate an ECAP in response to the drive signal, and thereby not contributing to the ECAP signal reducing the magnitude of the ECAP signal (e.g., the peak to peak magnitude 1320 of
Substrate changes correspond to changes in a composition of tissue surrounding the CSF (e.g., the CSF 508 of
In various embodiments, certain steps (or operations) may be omitted or added, certain steps may be combined, certain steps may be performed simultaneously, certain steps may be performed concurrently, certain steps may be split into multiple steps, certain steps may be performed in a different order, or certain steps or series of steps may be re-performed in an iterative fashion. It should be noted, other methods may be used, in accordance with embodiments herein.
One or more methods may (i) emit a stimulation waveform from at least one electrode of a lead, (ii) acquire impedance and evoked compound action potential (ECAP) signals, (iii) select a first impedance and ECAP measurement and a second impedance and ECAP measurements from the impedance and ECAP signals, (iv) detect spinal cord (SC) tissue damage based on a difference between the second ECAP measurement and the first ECAP measurement, and (v) adjust at least one stimulation parameter to change the stimulation waveform or adjust the location of the lead 402 based on detection of the SC tissue damage.
Beginning at 2652, the lead 402 (shown in
At 2654, the lead acquires impedance and ECAP signals. The impedance and ECAP signals are generated in response to the stimulation waveform. As further described in connection to 306 of
The stimulation waveform 2702 may be emitted from one or more electrodes 511 as shown in
The impedance signal 2748 is formed by impedance measurements 2750-2758 plotted with respect to a vertical axis 2718 representing a magnitude (e.g., ohms) acquired during the pulses 2704-2712 of the stimulation waveform 2702. Additionally or alternatively, in other embodiments the impedance signal 2748 may be measured continuously by the sensing circuitry 258 and/or controller 251.
The impedance signal 2748 may be acquired similar to and/or the same as the impedance signal acquired at 306 of
The ECAP measurements 2760 may correspond to a plurality of peak to peak measurements 2762-2770 plotted with respect to a vertical axis 2720 representing magnitude. The ECAP measurements 2760 may be derived from ECAP signals measured at one or more of the electrodes 511 during the stimulation waveform 2702 and received by the sensing circuitry 258.
For example, the electrode 511h may measure the EP during each of the pulses 2704-2712. The EP may be generated from neurons of the stimulation target in response to the pulses 2704-2712 emitted by one or more electrodes 511. The measured EP measurements form an ECAP signal for each pulse 2704-2712, such as the ECAP signal 1302 and/or 1303 shown in
The controller 251 may determine the peak to peak measurements 2762-2770 based on a difference in magnitude (e.g., EP) between the maximum peak (e.g., 1317, 1319) and the minimum peak (e.g., 1316, 1318) for each ECAP signal corresponding to the pulses 2704-2714.
Returning to
For example, the selected impedance and ECAP measurements may correspond to two different pulses within the stimulation waveform 2702. The controller 251 may select the impedance measurement 2752 and the peak to peak ECAP measurement 2764 corresponding to the pulse 2706 as the first impedance and ECAP measurement, and the impedance measurement 2754 and the peak to peak measurement 2766 as the second impedance and ECAP measurement. It should be noted in other embodiments the ECAP measurements may correspond to amplitudes of the maximum peak (e.g., 1317, 1319), amplitudes of the minimum peak (e.g., 1316, 1318), and/or the like of the ECAP signal.
At 2660, the controller 251 may detect damage to neural tissue of SC, a substrate change, or movement of the lead 402 (
The impedance column 2802 corresponds to differences determined by the controller 251 between the first impedance measurement (e.g., the impedance measurement 2752) and the second impedance measurement (e.g., the impedance measurement 2754) selected at 2658 of
For example, the controller 251 may determine whether the impedance increases, decreases, and/or does not change from the first impedance measurement to the second impedance measurement based on a predetermined impedance threshold. The predetermined impedance threshold may be stored in the memory 261, and correspond to a minimum amount of change between the two selected impedance measurements that will be determined and/or verified by the controller 251 as an increase and/or decrease in impedance.
For example, the predetermined impedance threshold may correspond to an amount of change between the impedance measurements that is greater than electrical noise (e.g., additive noise, electromagnetic interference, signal-to-noise ratio) generated by the operation of the IPG 250. The predetermined impedance threshold may be a percentage corresponding to an amount of change between the second and the first impedance measurement (e.g., three percent change, five percent change). In another example, the predetermined impedance threshold may correspond to a rate of change (e.g., slope, differential) determined from the first and second impedance measurements. The controller 251 may compare a difference or delta between the second impedance measurement (e.g., the impedance measurement 2754) and the first impedance measurement (e.g., the impedance measurement 2752) with the predetermined impedance threshold to determine whether the impedance has changed (e.g. increased, decreased) from the first impedance measurement to the second impedance measurement.
For example, the controller 251 may select the impedance measurement 2752 and the impedance measurement 2754. The controller 251 may calculate a difference 2770 of, for example a two percent drop, between the impedance measurements 2752-2754. The controller 251 may compare the difference 2770 with the predetermined impedance threshold of, for example three percent. It should be noted that in other embodiments the predetermined impedance threshold may be greater than or less than three percent. If the magnitude of the difference 2770 is greater than the threshold, the controller 251 may determine that the impedance changed (e.g., increased, decreased) between the first and second impedance measurements. Alternatively, if the magnitude of the difference 2770 is lower than the threshold, the controller 251 may determine that the impedance did not change. For example, the controller 251 may determine that since the difference 2770 (e.g., a two prevent) is below the predetermined impedance threshold (e.g., three percent) the impedance measurements did not change.
Based on a determination that the impedance has not changed between impedance measurements, the controller 251 may determine that the substrate surrounding the lead 402 has not changed. For example, as indicated along row 2808 in the decision table 2800, a substrate change is shown based on a change in impedance. Since the controller 251 determined that the impedance measurements did not change, the controller 251 may determine that a substrate change surrounding the lead 402 has not occurred.
The ECAP amplitude column 2804 corresponds to differences determined by the controller 251 between the first ECAP measurement (e.g., the peak to peak measurement 2764) and the second ECAP measurement (e.g., the peak to peak measurement 2766) selected at 2658 of
For example, the controller 251 may determine whether the ECAP measurement increases, decreases, and/or does not change from the first ECAP measurement to the second ECAP measurement based on a predetermined ECAP amplitude threshold. The predetermined ECAP amplitude threshold may be stored in the memory 261, and correspond to an amount of change between two ECAP measurements that can be determined and/or verified by the controller 251 as an increase and/or decrease in ECAP amplitude.
For example, the predetermined ECAP amplitude threshold may correspond to an amount of change between the ECAP measurements that is greater than electrical noise (e.g., additive noise, electromagnetic interference, signal-to-noise ratio) generated by the operation of the IPG 250. The predetermined ECAP amplitude threshold may be a percentage corresponding to an amount of change between the second and the first ECAP measurements (e.g., three percent change, five percent change). In another example, the predetermined ECAP amplitude threshold may correspond to a rate of change (e.g., slope, differential) determined from the first and the second ECAP measurements. The controller 251 may compare a difference or delta between the second ECAP (e.g., the peak to peak measurement 2766) and the first ECAP measurement (e.g., the peak to peak measurement 2764) with the predetermined ECAP amplitude threshold to determine whether the ECAP amplitude has changed from the first ECAP measurement to the second ECAP measurement.
For example, the controller 251 may select the peak to peak measurement 2764 and the peak to peak measurement 2766. The controller 251 may calculate a difference 2774 (e.g., at 2658 of
Based on a determination that the ECAP amplitudes have not changed between ECAP measurements, the controller 251 may determine that the lead 402 has not moved relative to the SC tissue between the ECAP measurements and the SC tissue is not damaged. For example, as indicated along rows 2810-2820 in the decision table 2800, a movement in the lead 402 and damage to the SC tissue is shown based on a change in ECAP amplitude. Since the controller 251 determined that the ECAP measurements did not change, the controller 251 may determine that the lead 402 did not move between the ECAP measurements and/or the SC tissue is not damaged.
In various embodiments, the controller 251 may compare the ECAP measurements to a critical threshold. The critical threshold may be stored in the memory 261. The critical threshold may correspond to a change between the two ECAP measurements indicating neural tissue damage of the SC. The critical threshold is greater than the ECAP amplitude threshold. For example, the critical threshold may be a percentage corresponding to an amount of change between the first and the second ECAP measurement (e.g., thirty percent change, forty percent change, fifty percent change). In another example, the predetermined EACP amplitude threshold may correspond to a rate of change (e.g., slope, differential) determined from the first and second ECAP measurement.
Beginning at 2902, the controller 251 determines whether the ECAP amplitude has decreased between the first ECAP measurement and the second ECAP measurement more than the critical threshold. For example, the controller 251 may select at 2656 (
Alternatively, if the magnitude of the difference 2774 was greater than the critical threshold, the controller 251 may determine, at 2912, neural tissue damage of the SC at and/or proximate to the stimulation target due to overstimulation.
At 2904, the controller 251 determines whether the ECAP amplitude has increased between the first ECAP measurement and the second ECAP measurement. The controller 251 may compare the difference 2776 determined at 2658 with the predetermined ECAP amplitude threshold to determine whether an increase of the ECAP measurement corresponds to lead 402 movement. For example, the difference 2776 corresponds to an increase in the second ECAP measurement (e.g., the peak to peak measurement 2768) of twenty percent. The controller 251 may compare the difference 2776 with the predetermined ECAP amplitude threshold of, for example, three percent. Since the difference 2776 is greater than the predetermined ECAP amplitude threshold, the controller 251 may determine that the ECAP amplitude increased.
Alternatively, if the difference 2776 is lower than the predetermined ECAP amplitude threshold, the controller 251 may determine that the ECAP amplitude did not increase and that the lead 402 did not move towards the SC tissue and/or stimulation target.
At 2913, the controller 251 determines whether the ECAP amplitude decreased above the critical threshold within a time period. The time period may be based on a number of subsequent pulses of the stimulation waveform 2702, relative to the peak to peak measurement used in 2904. For example, the time period may correspond to when the successive pulse 2712 occurs relative to the pulse 2710 which corresponds to the peak to peak measurement 2768. Additionally or alternatively, the time period may be based on a number of a subsequent stimulation waveforms, relative to the peak to peak measurement used in 2904. In various embodiments, the time period may be based on time length, such as seconds, minutes, and/or hours.
Optionally, a length of the time period may be based on a procedure being performed on the patient. For example, during implantation of the lead 210 into the patient, the controller 251 may have the time period be a short period (e.g., several seconds, less than minutes). When the implantation is complete, the controller 251 may have the time period be longer relative to the implantation procedure (e.g., several hours, days).
After the time period, the controller 251 may determines whether the ECAP amplitude has decreased between the ECAP measurement used in 2904, the peak to peak measurement 2768, and a second ECAP measurement after the time period. For example, the time period may be based on the subsequent pulse, the pulse 2712. The controller 251 may select the corresponding peak to peak measurement 2770 as the second ECAP measurement. The controller 251 may compare the peak to peak measurements 2768-2770 to determine a difference 2778, for example a thirty-five percent decrease, between the ECAP measurements (e.g., 2768-2770). For example, the peak to peak measurement 2770 is thirty-five percent lower than the peak to peak measurement 2768. The controller 251 may compare the difference 2778 with the predetermined critical threshold, for example of thirty percent. Since, the difference 2774 is lower than the critical threshold, the controller 251 may determine, at 2915, that neural tissue damage of the SC due to pressure from the lead 402 is detected.
Alternatively, if the controller 251 determines at 2913 that the ECAP amplitude did not decrease above the critical threshold after the time period, then at 2914, the controller 251 may determine that the ECAP amplitude increase at 2904 is due to lead 402 movement toward the SC tissue.
At 2906, the controller 251 determines whether the ECAP amplitude has decreased. The controller 251 may compare a difference between the first ECAP measurement and the second ECAP measurement determined at 2658 with the predetermined ECAP amplitude threshold to determine whether a decrease from the first ECAP measurement to the second ECAP measurement corresponds to lead 402 moving away from the SC and/or stimulation. If the difference is greater than the predetermined ECAP amplitude threshold, the controller 251 at 2916 may determine that the ECAP amplitude decrease is due to lead 402 movement away from the SC tissue and/or stimulation target. Alternatively, if the difference is lower than the predetermined ECAP amplitude threshold, the controller 251 may determine that the lead 402 did not move away from the SC tissue and/or stimulation target.
At 2908, the controller 251 determines whether the impedance has changed (e.g., increased, decreased) between the first impedance measurement and the second impedance measurement above the predetermined impedance threshold. For example, the controller 251 may select at 2656 (
Returning to
For example, if the controller 251 detects the lead 402 has moved away from the SC tissue and/or stimulation target, the controller 251 may increase the amplitude of the pulses 2704-2712 forming the stimulation waveform 2702, increase the duration of pulses 2704-2712 of the stimulation waveform 2702, increase the number of pulses 2704-2712 forming the stimulation waveform 2702, and/or the like.
In another example, if the controller 251 detects the lead 402 has moved towards the SC tissue, the controller 251 may decrease the amplitude of the pulses 2704-2712 forming the stimulation waveform 2702, decrease the duration of pulses 2704-2712 of the stimulation waveform 2702, decrease the number of pulses 2704-2712 forming the stimulation waveform 2702, and/or the like.
Optionally, the controller 251 may have the IPG 250 stop delivering the stimulation waveform 2702. For example, if the controller 251 detects SC tissue damage, the controller 251 may have the IPG 250 stop delivering pulses 2704-2712 forming the stimulation waveform 2702 to the electrodes 511.
Additionally or alternatively, the controller 251 may transmit to an external device, such as the system 100, via the far-field and/or near field communication circuitry 255 when the controller 251 detects damage to the SC tissue, a substrate change, or movement of the lead, at 2662. For example, the controller 251 may transmit an alert signal to the system 100 when the substrate proximate to the lead 402 has changed and/or neural damage of the SC tissue is detected.
The controller 251, the CPU 102, and the CPU 152 may include any processor-based or microprocessor-based system including systems using microcontrollers, reduced instruction set computers (RISC), application specific integrated circuits (ASICs), field-programmable gate arrays (FPGAs), logic circuits, and any other circuit or processor capable of executing the functions described herein. Additionally or alternatively, the controller 251, the CPU 102, and the CPU 152 may represent circuit modules that may be implemented as hardware with associated instructions (for example, software stored in a tangible and non-transitory computer readable storage medium, such as a computer hard drive, ROM, RAM, or the like) that perform the operations described herein. The above examples are exemplary only, and are thus not intended to limit in any way the definition and/or meaning of the term “controller.” The controller 251, the CPU 102, and the CPU 152 may execute a set of instructions that are stored in one or more storage elements, in order to process data. The storage elements may also store data or other information as desired or needed. The storage element may be in the form of an information source or a physical memory element within the controller 251, the CPU 102, and the CPU 152. The set of instructions may include various commands that instruct the controller 251, the CPU 102, and the CPU 152 to perform specific operations such as the methods and processes of the various embodiments of the subject matter described herein. The set of instructions may be in the form of a software program. The software may be in various forms such as system software or application software. Further, the software may be in the form of a collection of separate programs or modules, a program module within a larger program or a portion of a program module. The software also may include modular programming in the form of object-oriented programming. The processing of input data by the processing machine may be in response to user commands, or in response to results of previous processing, or in response to a request made by another processing machine.
It is to be understood that the subject matter described herein is not limited in its application to the details of construction and the arrangement of components set forth in the description herein or illustrated in the drawings hereof. The subject matter described herein is capable of other embodiments and of being practiced or of being carried out in various ways. Also, it is to be understood that the phraseology and terminology used herein is for the purpose of description and should not be regarded as limiting. The use of “including,” “comprising,” or “having” and variations thereof herein is meant to encompass the items listed thereafter and equivalents thereof as well as additional items.
It is to be understood that the above description is intended to be illustrative, and not restrictive. For example, the above-described embodiments (and/or aspects thereof) may be used in combination with each other. In addition, many modifications may be made to adapt a particular situation or material to the teachings of the invention without departing from its scope. While the dimensions, types of materials and coatings described herein are intended to define the parameters of the invention, they are by no means limiting and are exemplary embodiments. Many other embodiments will be apparent to those of skill in the art upon reviewing the above description. The scope of the invention should, therefore, be determined with reference to the appended claims, along with the full scope of equivalents to which such claims are entitled. In the appended claims, the terms “including” and “in which” are used as the plain-English equivalents of the respective terms “comprising” and “wherein.” Moreover, in the following claims, the terms “first,” “second,” and “third,” etc. are used merely as labels, and are not intended to impose numerical requirements on their objects. Further, the limitations of the following claims are not written in means—plus-function format and are not intended to be interpreted based on 35 U.S.C. §112(f), unless and until such claim limitations expressly use the phrase “means for” followed by a statement of function void of further structure.
The present application relates to and claims priority from U.S. provisional application Ser. No. 62/088,451 filed Dec. 5, 2014, entitled “SPINAL CORD STIMULATION GUIDANCE SYSTEM AND METHOD OF USE,” which is hereby expressly incorporated by reference in its entirety.
Number | Date | Country | |
---|---|---|---|
62088451 | Dec 2014 | US |