Next Article in Journal
Novel Methodology to Recover Road Surface Height Maps from Illuminated Scene through Convolutional Neural Networks
Next Article in Special Issue
The Role of Artificial Intelligence in Decoding Speech from EEG Signals: A Scoping Review
Previous Article in Journal
Elevator Car Vibration Signal Denoising Method Based on CEEMD and Bilateral Filtering
Previous Article in Special Issue
A Comprehensive Review of Endogenous EEG-Based BCIs for Dynamic Device Control
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Estimation of ANT-DBS Electrodes on Target Positioning Based on a New PerceptTM PC LFP Signal Analysis

1
INESC TEC—Instituto de Engenharia de Sistemas e Computadores, Tecnologia e Ciência, Faculty of Engineering, University of Porto, 4200-465 Porto, Portugal
2
Neurophysiology Unit, Neurology Department, Centro Hospitalar Universitário de São João, 4200-319 Porto, Portugal
3
Neurosurgery Department, Centro Hospitalar Universitário de São João, 4200-319 Porto, Portugal
*
Author to whom correspondence should be addressed.
Sensors 2022, 22(17), 6601; https://doi.org/10.3390/s22176601
Submission received: 1 July 2022 / Revised: 27 August 2022 / Accepted: 28 August 2022 / Published: 1 September 2022
(This article belongs to the Special Issue Brain Activity Monitoring and Measurement)

Abstract

:
Deep brain stimulation of the Anterior Nucleus of the Thalamus (ANT-DBS) is an effective therapy in epilepsy. Poorer surgical outcomes are related to deviations of the lead from the ANT-target. The target identification relies on the visualization of anatomical structures by medical imaging, which presents some disadvantages. This study aims to research whether ANT-LFPs recorded with the PerceptTM PC neurostimulator can be an asset in the identification of the DBS-target. For this purpose, 17 features were extracted from LFPs recorded from a single patient, who stayed at an Epilepsy Monitoring Unit for a 5-day period. Features were then integrated into two machine learning (ML)-based methodologies, according to different LFP bipolar montages: Pass1 (nonadjacent channels) and Pass2 (adjacent channels). We obtained an accuracy of 76.6% for the Pass1-classifier and 83.33% for the Pass2-classifier in distinguishing locations completely inserted in the target and completely outside. Then, both classifiers were used to predict the target percentage of all combinations, and we found that contacts 3 (left hemisphere) and 2 and 3 (right hemisphere) presented higher signatures of the ANT-target, which agreed with the medical images. This result opens a new window of opportunity for the use of LFPs in the guidance of DBS target identification.

1. Introduction

Epilepsy is the second most common neurological disease, affecting around 0.5% of the population worldwide [1]. Treatment with adequately chosen antiseizure drugs allows for adequate control in up to 60% of patients, but the remainder have persisting and frequently disabling seizures (drug-resistant epilepsy, DRE) [2]. For a minority of those patients, resective surgery is a safe and effective option, but a substantial number of patients do not achieve long-lasting remission [3]. Moreover, many patients with DRE are not suitable candidates for cortical resections due to overlapping eloquent cortex or multifocal seizure onsets. Neurostimulation is often considered in these patients, including vagal nerve stimulation or deep brain stimulation (DBS). For adults with focal DRE, the DBS of the anterior nucleus of the thalamus (ANT) is an established approved treatment [4].
The ANT has long been considered a suitable DBS anatomical target for patients with refractory epilepsy, due to its therapeutic action on the interruption of seizure spread [5,6,7,8]. The ANT is in the superior region of the thalamus and is separated from the rest of the thalamus by the anterior medullary lamina. It consists of three subnuclei designated as anteroventral (AV), anterodorsal (AD) and anteromedial (AM) nuclei [9]. After the first ANT-DBS surgery performed by Cooper and Upton in 1980 [10,11,12], several studies have reported the therapeutic efficacy as well as the safety of this procedure [13,14,15,16,17,18]. The larger-scale multicenter trial, the SANTE trial, has reported a seizure reduction of 56–69% [19,20].
It is generally accepted that stimulation of the ANT is more effective compared to stimulation outside this region, since lead locations deviating from the target region (the AV nucleus) have been related with poorer surgical outcomes [21,22]. Furthermore, it is known that changes in stimulation parameters (current voltage, pulse width, frequency and stimulation cycling) produce a minimal effect compared with the choice of the active contact selection [21,22]. This highlights the importance of an accurate identification of the precise location of stimulation.
The identification of the DBS target relies on the visualization of anatomical structures in magnetic resonance imaging (MRI), which can be performed indirectly or directly [21]. In the indirect methodology, the target is defined in the brain atlas using common landmarks, such as the anterior and posterior commissures; the direct method can be performed using 3T MRI techniques [23], which allow for the direct visualization of white-matter structures involving the ANT, such as the external medullary lamina (EML), the internal medullary lamina (IML) and the mammillothalamic tract (MMT) [23,24]. However, both methods present some disadvantages: the indirect method is limited due to anatomical variations of the target structure in the stereotactic space between individuals [23,24]; the direct method, on the other hand, requires advanced imaging techniques, which may be limited in some DBS centers [23].
Several studies have investigated the potential of complementary methods in the guidance of the target identification, such as the analysis of single-unite (SU) recordings during DBS [21,24,25,26,27]. Hodaie et al. were the first to describe the electrophysiological properties of SU signals from the ANT of anesthetized patients during DBS along transventricular trajectory. They found bursting activity characterized as low-threshold calcium spikes (LTS), which were mostly observed during sleep [25]. Later, Schaper et al., found the same pattern using an extraventricular trajectory to the ANT [27]. The same result was also reported by Mottonen et al. [28]. Considering both trajectories, the same authors showed that SU signals were able to distinguish between the ANT from the ventral anterior nucleus [28]. In all these studies, an increase in the firing rate was found at the entrance of the ANT, and a decrease was found when exiting this structure. Despite these advances, it is still unknown whether there is a relation between the neuronal firing properties of the ANT and the clinical DBS outcome in epilepsy [27,28]. On the other hand, LTS bursts were not specific to ANT, since they were also found in other structures, such as the circularis and dorsomedial nucleus of the thalamus [27]. This highlights the need to find new methods that can complement the existing techniques for DBS targeting.
Local field potentials (LFPs) can be recorded by DBS leads and reflect the synchronous pre- and postsynaptic activity of neural populations [29,30,31,32,33,34]. Unlike SU recordings, LFPs can detect focal network rhythms and can be recorded both intra and extraoperatively [35]. Temporary recordings of LFPs can be performed by externalized DBS leads (e.g., [36]), but this recording type presents disadvantages, such as the tome restriction of recording, the influence of microlesions provoked by the oedema around the lead [37] and the increased risk of infection [38]. Recently, the Medtronic company (Medtronic Inc. (Dublin, Ireland)) launched a new neurostimulator device, the PerceptTM PC, which can record LFPs at the same time it stimulates, whether the patient is in the hospital or not (UC202013078EE©Medtronic2020). Compared to the previous systems (e.g., the Active TM PC+S), the PerceptTM PC system has longer battery life and is also able to stream data continuously in real time, as well as to correlate data with patient logged events.
To understand the potential role of DBS-LFPs for DBS targeting, a literature survey in PubMed and Web of Science was conducted, following the keywords (“PerceptTM PC” [tiab] or “ActivaTM PC+S” [tiab] and “Local Field Potentials”). A total of 12 human-based studies were found [39,40,41,42,43,44,45,46,47,48,49,50]. These studies are summarized in Table 1. Most of them research and report biomarkers of motor fluctuations in Parkinson’s disease (PD). These biomarkers can be found by extracting features from LFPs. The most common type of features can be categorized into spectral (features extracted in the frequency domain reflecting potential fluctuations associated with neural activity); morphological (features describing the signal morphology); statistical (features describing signal variance distribution) and multivariate (features that capture correlations between channels) [51]. For PD, LFP power density was the most frequent feature extracted since beta activity has been seen to correlate with PD symptoms. Regarding epilepsy, no study was found. Moreover, none of these studies focus of the study of LFPs in the optimization of DBS targeting.
In this work, we aimed to study whether ANT-LFPs recorded extraoperatively can be an asset in the identification of DBS target structures, to complement structural imaging approaches, making its joint usage a swifter and more reliable procedure. For this purpose, LFP signals were recorded from a single epilepsy patient, who stayed at an epilepsy monitoring unit (EMU) for 5 days, for simultaneous video-electroencephalography (vEEG) and Percept PC-LFP recordings. From time-domain LFP signals recorded during periods with no stimulation, 17 features were extracted and then integrated into two machine learning (ML)-based methodologies, developed to identify which led contacts (in both hemispheres) show a “signal-signature” of the ANT target. This paper aims to present the processing and classification methodology and discuss the results and their capability to support neurosurgeons to use LFPs an surrogate marker of the AV nucleus position. To the best of our knowledge, we present the first-ever study to use multichannel LFP signals, collected months after electrodes’ implantation (already out of the influence of the microlesion effect), to identify the best-targeted contacts of the implanted electrodes.

2. Materials and Methods

2.1. Subject

A 54-year-old right-handed patient, diagnosed with focal refractory epilepsy of hypoxic-ischemic etiology with bilateral perisylvian ulegryc lesions, was admitted to the epilepsy monitoring unit (EMU) of the University Hospital S. João, Porto, Portugal, for simultaneous vEEG and ANT-DBS recordings, over 5 days (20–24 July 2020). This monitoring took place approximately one month after a bilateral ANT-DBS implantation. The patient’s seizures were markedly active by sleep and were characterized by bilateral spasms or tonic posturing, sometimes follows by hyperkinetic movements, and left or bilateral clonic activity. Ictal patterns in scalp EEG were in the midline frontocentral regions or obscured by muscle artefacts; interictally, rare epileptiform discharges were seen over the same region. Background activities awake and asleep were normal. This study complies with the requirements of medical confidentially according to the Declaration of Helsinki and data protection. The patient’s data were pseudonymized to not reveal names, initials, birthdays or other private information.

2.2. DBS Lead Contact Localization

The DBS leads were automatically detected by neurosurgeons with the “Lead localization” module of a commercially available software (Elements, BrainLab, Munich, Germany), using preoperative MRIs (fast gray matter acquisition T1 inversion recovery (FGATIR) 3T MRI and fluid attenuated inversion recovery (FLAIR) T2 MRI) and postoperative computed tomography (CT) (Figure 1). Each DBS lead included four contacts (0, 1, 2 and 3) and was implanted via transventricular trajectory. The ANT target of each hemisphere was directly identified by the visualization of the MMT. By inspecting the image, the most inserted contact in the left target is the 3L; and in the right, the 2R.

2.3. LFP Recording

Time-domain LFPs were recorded using the BrainSense Survey mode of the PerceptTM PC system, which is used to give a broad spatial overview of LFP signals during stimulation from which clinicians usually select a frequency band of interest for each patient for chronic recordings (UC202013078EE©Medtronic2020). The Survey mode acquired data in 6 electrode channels (sampling frequency of 250 Hz) in bipolar reference, according to two types of combinations, i.e., Passes: Pass1 (channels 0–3, 1–3, 0–2), which are stimulation-compatible pairs; and Pass2 (channels 0–1, 1–2, 2–3), which are immediately adjacent pairs. A total of 11 BrainSense Surveys were performed during the patient’s EMU stay, which corresponds to a total of 231 s of time-domain LFPs, recorded in each Pass. Pass1 and Pass2 data were then divided in segments of 5 s, making a total of 46 segments for each Pass.

2.4. Dataset Selection

In this work, we aimed to assess whether LFP signals are capable of distinguishing between anatomical structures underlying lead contacts from which LFPs were recorded. For this purpose, we analyzed and compared LFPs recorded from contacts that were mostly inserted into the ANT target and contacts farthest from this target, using an ML-based methodology. Considering Figure 1, contacts 3L and 3R were closer to the target and contacts 0L and 0R were farthest. Since LFP signals were recorded using a bipolar reference, signal spatial locations had to be considered as the midpoints between pairs of contacts [52], to assign each recorded signal to the respective anatomical location. From Pass1 and Pass2 reference combinations, it was possible to identify 5 LFP spatial locations at each hemisphere, indicated as A, B, C, D and E (Figure 2B). The corresponding channels are summarized in Table 2, as well as the proportion of the ANT target underlying each contact and a qualitative measure created to express the level of involvement of each location in the target, based on medical imaging analysis.
This patient was implanted with the Medtronic electrode model 3389, which contains contacts with a length of 1.5 mm separated from each other by 0.5 mm (Figure 2A). Depending on the type of bipolar reference used, the distance between each contact pair to its spatial location (d) varies for Pass2, d = 1 mm for all possible combinations; and for Pass1, d = 2 mm for 0–2 and 1–3 pairs, or d = 3 mm for the 0–3 pair. Due to these differences, we considered the Pass1 and Pass2 signals to be of a different nature. Therefore, we developed two ML-based models for differentiating on-target (signals with physical locations completely inserted in the target: +++) and off-target signals (signals with physical locations completely outside the target: −). The general methodology used is in Figure 3.
For the Pass1 model, we assigned LFP signals with spatial locations DRight as on-target and with BLeft as off-target; for the Pass2 model, LFPs with spatial locations ERight were assigned as on-target and with ALeft as off-target. For each model, corresponding 5 s-segmented signals were used for feature extraction.

2.5. Feature Extraction

A total of 17 features from different feature types were extracted (Table 3). These features were selected according to their consensual significance in the epilepsy state of the art [51,53]. All features were extracted using MATLAB R2016b. Statistical features were extracted using the EEG Feature Extraction Toolbox, while the remainder were extracted using MATLAB custom-made functions. As multivariate features, the node strength (the sum of weights links connected to the node) was computed of each adjacency matrix estimated with metrics present in Table 3. Node strength values were computed using the MATLAB Brain Connectivity Toolbox.
All features were ranked for each dataset, using the MATLAB rankfeatures function, which used an independent evaluation criterion for binary classification (t-test pooled variation estimation), to assess which features were more discriminative for the two considered classes (Figure 4). Since there was no agreement on the most relevant features for each dataset (set of features extracted from each in- and off-target spatial location, according to each Pass), a dimensionality reduction procedure through a projection method was performed.

2.6. Feature Reduction

While a high number of features may help to discriminate classes, there is a greater potential for overfitting, especially in small datasets. Moreover, irrelevant features can also blur the boundaries between classes [54]. For dimensionality reduction, principal component analysis (PCA) was used, a feature projection method that uses orthogonal transformations for transforming a set of observations of possibly correlated variables into a set of linearly uncorrelated variables (principal components) [7]. The PCA was performed using the pca MATLAB function. The minimum number of principal components (PC) in each Pass-based methodology was chosen as the minimum number that kept 99% of the variance of the entire dataset [55]. The dataset transformed in the reduced space was then used for classification.

2.7. Classification

The classification was performed according to two methods: Pass1 method (Pass1 LFP signals) and Pass2 method (Pass2 LFP signals). For each dataset in the reduced space, we assigned each sample with a class label 0 (off-target) and 1 (in-target). Then, each dataset was split into 80% for training and 20% for test. The training set was also split into 80% for a second set of training (training set 2) and 20% for cross-validation (CV dataset). In training set 2, three classification models were trained using standard parameters (Table 4) and tested in the CV dataset. The best classification model (with higher accuracy (acc)) was then retrained in the training set 2 dataset and tested in the CV dataset with varying model parameters, to identify the best parameter (the parameter for which the classifier had higher accuracy).
The Pass1 dataset was reduced into 13 principal components and neural networks with a layer size of 2 were considered the best model and parameter to discriminate classes (acc = 64.29% and acc = 100%, respectively). The Pass2 dataset was reduced into 11 principal components and the kNN with a squared Euclidean distance (seuclidean) were considered the best model and classification model (acc = 94.44% for bot classifications). These models were used to train and test the training and the test sets of each dataset. To assess and compare the classification performance of each method in the discrimination of signals coming from in- and off-target physical locations, we computed the accuracy (number of corrected classifications divided by the total number of classifications), sensitivity (describes the effectiveness of the classifier to classify correctly in-target samples) and specificity (describes the effectiveness of the classifier to classify correctly off-target samples). The classification was performed using the MATLAB machine learning package [56].

2.8. Prediction

The best model of each methodology was used to predict the Target % of the remaining spatial locations relative to each pass, as well as those used to build each method (Table 5). For this purpose, LFP signals related to each location were 5 s-segmented and features were extracted for each segment (Table 3), and then reduced using the number of principal components obtained for each methodology (13 for the Pass1 model and 11 for the Pass2 model). Our goal was to make use of the discriminative power of each classifier to distinguish between in- and off-target LFP signal physical locations, and consequently predict which electrodes were more inserted in the target, which may be the best choices for chronic stimulation.

3. Results

The Pass1 classifier (NN with 2 hidden layers) distinguished LFPs from DRight (in-target) and BLeft (off-target) spatial locations with an accuracy of 76.62%, sensitivity of 81.82% and specificity of 71.43%. The Pass2 classifier (kNN with a seuclidean distance) distinguished LFPs from ERight (in-target) and ALeft (off-target) physical locations with an accuracy of 83.33%, sensitivity of 100% and specificity of 66.6%. Predictions performed for each classifier are summarized in Table 6.

4. Discussion

Both classifiers (Pass1 and Pass2) classified the two class groups (in-target and out-target signals) with an average accuracy of 80%. Specificity, which quantifies the effectiveness of classifiers to identify correctly negative class labels (i.e., non-target signals), was also similar for both classifiers (70%, in average). Sensitivity, which classifies the effectiveness of classifiers to classify correctly positive class labels, i.e., target signals, was higher for the Pass2 classifier than Pass 1 (100% against 81.82%). Therefore, we considered that the Pass2 classifier had high performance, which was expected since spatial locations in Pass2 signals are closest to LFP channels.
The Pass1 classifier identified spatial locations completely outside the target (−) as 40% of the target; physical locations marginally in the target (+) as 50%; physical locations partially but not completely inserted in the target (++) as 50–70%. The spatial location completely inserted in the target (+++) was used as a classifier class label. Since the classifier results increased with the degree of involvement of each spatial location to the target, the Pass2 classifier predictions seem to agree with the medical imaging results (Figure 1).
The Pass2 classifier identified spatial locations completely outside the target (−) as 50% of the target; spatial locations marginally in the target (+) as 60%; spatial locations partially but not completely inserted in the target (++) as 50%. The spatial location completely inserted in the target (+++) was also used as a classifier label. These predictions seem to be worse than those obtained with the Pass1 classifier, since DLeft (+) was classified to be more target percentage than CRight (++).
The target percentage of spatial locations used as class labels for each classifier were predicted by the opposite classifier (Table 5). Pass1 classified ALeft as 50% (which agrees with other predictions made by his classifier) and ERight as 60% (higher target percentage obtained than other (++) spatial locations); The Pass2 classifier classified BLeft as 60% (higher target percentage than other (+) spatial locations) and DRight as 60% (the same order of magnitude as other (++) predictions). Therefore, the Pass1 classifier seems to predict better Pass2 spatial locations than the opposite situation. This result was expected since spatial locations of nonadjacent channels include those for adjacent channels. Note that we considered Pass1 and Pass2 signals from different natures, since distance between channels and their midpoints vary according to the bipolar montage used. However, the only way to obtain a quantitative evaluation of the target prediction of spatial locations used to build each classifier is to use the opposite classifier. This is a limitation of the proposed methodology.
The gamma band power (30–100 Hz) was considered the most significant feature to discriminate between target and non-target signals in both montages and for the considered dataset (Figure 4). For the Pass2 classifier, the alpha band power (8–13 Hz) was equally relevant, whereas for the Pass1 classifier, all spectral features were considered discriminative. Considering morphological features, in the Pass2 classifier only the average mean had discriminative power; mean peaks and minimum values had middle discriminative power and the maximum value was the less discriminative. For the Pass1 classifier, morphological features were medially discriminative. Considering statistical features, variance and kurtosis seem to have good discriminative power for the Pass2 classifier, whereas for the Pass1 classifier only kurtosis appeared to have middle-range discriminative power. Considering multivariate features, the phase locking value had good discriminative power and cross-correlation was medially discriminative power in both classifiers. The phase lag index also discriminated well the two classes for the Pass2 classifier, and coherence had middle discriminative power. For the Pass1 classifier, phase lag index and coherence were not good features to discriminate between two classes.

5. Conclusions

In this work, we extracted 17 features from off-stimulation time-domain LFPs, recorded from one patient; and two ML-based methods, designed to identify which lead contacts presented a “signal signature” of the ANT target, the anteroventral (AV) nucleus. These methods were designed using different LFP bipolar montages: Pass1 (nonadjacent channels) and Pass2 (adjacent channels). After the best classification model for each dataset was identified, we obtained an accuracy of 76.62% for the Pass1 classifier and 83.33% for the Pass2 classifier. Then, we used both classifiers to predict the target percentage of all possible spatial locations, with the aim of identifying which contacts (in the left and right hemisphere) were closer or inserted to the target, and consequently, should be chosen for chronic stimulation.
We found that the spatial location ELeft had highest target percentage in the left hemisphere and spatial locations DRight and ERight in the right hemisphere. Converting spatial locations into channels (Table 2), we concluded that LFP signals from 2–3 L and 2–3 R and 1–3 R channels may present higher ANT target signatures. These results agree with those obtained by clinical images (Figure 1).
We concluded that the Pass1 and Pass2 classifiers were able to predict ANT target signatures for this patient. These results may be useful in support of the early selection of the best stimulation electrode, avoiding the trial-and-error process that usually takes place in the months after the implantation. To the best of our knowledge, we presented the first-ever study using multichannel LFPs collected months after the electrodes’ implantation (already out of the influence of the microlesion effect) in the guidance of the DBS target identification. However, further studies should be carried out in the future, with a larger number of patients, to validate this result.

Author Contributions

Conceptualization, E.M.L., R.R. and J.P.S.C.; methodology, E.M.L. and J.P.S.C.; software, E.M.L.; validation, R.R., M.R. and C.C.; formal analysis, E.M.L., R.R., M.R., C.C., D.D. and J.P.S.C.; investigation, E.M.L., R.R., M.R., C.C., D.D. and J.P.S.C.; resources, R.R., M.R. and C.C.; writing—original draft preparation, E.M.L.; writing—review and editing, E.M.L., R.R., M.R., C.C., D.D. and J.P.S.C.; supervision, J.P.S.C. All authors have read and agreed to the published version of the manuscript.

Funding

This work is financed by National Funds through the Portuguese funding agency, FCT–Fundação para a Ciência e Tecnologia, within the project LA/P/0063/2020. The first author is supported by a PhD grand from FCT—Fundação para a Ciêncua e Tecnologia (SFRH/BD/140273/2018).

Institutional Review Board Statement

The study was conducted in accordance with the Declaration of Helsinki and approved by the Institutional Review Board and Ethics Committee of Centro Hospitlar de São João/Faculdade de Medicina da Universidade do Porto.

Informed Consent Statement

Informed consent was obtained from all subjects involved in the study.

Data Availability Statement

Not applicable.

Acknowledgments

The first author thanks the support of FCT–Fundação para a Ciêncua e Tecnologia; Neurophysiology technicians at the Epilepsy Monitoring Unit: Ângela Santos, Ana Luísa Gonçalves, Sofia Marques, André Marques, Carlos Fernanrdes.

Conflicts of Interest

The authors declare no conflict of interest. The funders had no role in the design of the study; in the collection, analysis, or interpretation of data; in the writing of the manuscript; or in the decision to publish the results.

References

  1. Fisher, R.S.; Cross, J.H.; French, J.A.; Higurashi, N.; Hirsch, E.; Jansen, F.E.; Lagae, L.; Moshé, S.L.; Peltola, J.; Roulet Perez, E.; et al. Operational classification of seizure types by the International League Against Epilepsy: Position Paper of the ILAE Commission for Classification and Terminology. Epilepsia 2017, 58, 522–530. [Google Scholar] [CrossRef] [PubMed]
  2. Kwan, P.; Brodie, M.J. Early identification of refractory epilepsy. N. Engl. J. Med. 2000, 342, 314–319. [Google Scholar] [CrossRef] [PubMed]
  3. West, S.; Nolan, S.J.; Newton, R. Surgery for epilepsy: A systematic review of current evidence. Epileptic Disord. 2016, 18, 113–121. [Google Scholar] [CrossRef] [PubMed]
  4. Benbadis, S.R.; Geller, E.; Ryvlin, P.; Schachter, S.; Wheless, J.; Doyle, W.; Vale, F.L. Putting it all together: Options for intractable epilepsy: An updated algorithm on the use of epilepsy surgery and neurostimulation. Epilepsy Behav. 2018, 88, 33–38. [Google Scholar] [CrossRef]
  5. Mirski, M.A.; Ferrendelli, J.A. Interruption of the mammillothalamic tract prevents seizures in guinea pigs. Science 1984, 226, 72–74. [Google Scholar] [CrossRef]
  6. Mirski, M.A.; Ferrendelli, J.A. Anterior thalamic mediation of generalized pentylenetetrazol seizures. Brain Res. 1986, 399, 212–223. [Google Scholar] [CrossRef]
  7. Mirski, M.A.; Rossell, L.A.; Terry, J.B.; Fisher, R.S. Anticonvulsant effect of anterior thalamic high frequency electrical stimulationin the rat. Epilepsy Res. 1997, 28, 89–100. [Google Scholar] [CrossRef]
  8. Mirski, M.A.; McKeon, A.C.; Ferrendelli, J. Anterior thalamus and substantia nigra: Two distinct structures mediating experimental generalized seizures. Brain Res. 1986, 397, 377–380. [Google Scholar] [CrossRef]
  9. Gartner, L.P. A Textbook of Neuroanatomy; John Wiley & Sons: Hoboken, NJ, USA, 2009. [Google Scholar]
  10. Cooper, I.; Upton, A.; Amin, I. Reversibility of chronic neurologic deficits. Some effects of electrical stimulation of the thalamus and internal capsule in man. Stereotact. Funct. Neurosurg. 1980, 43, 244–258. [Google Scholar] [CrossRef]
  11. Cooper, I.; Upton, A.; Amin, I.; Garnett, S.; Brown, G.; Springman, M. Evoked metabolic responses in the limbic-striate system produced by stimulation of anterior thalamic nucleus in man. Int. J. Neurol. 1984, 18, 179–187. [Google Scholar]
  12. Upton, A.; Cooper, I.; Springman, M.; Amin, I. Suppression of seizures and psychosis of limbic system origin by chronic stimulation of anterior nucleus of the thalamus. Int. J. Neurol. 1985, 19, 223–230. [Google Scholar] [PubMed]
  13. Hodaie, M.; Wennberg, R.A.; Dostrovsky, J.O.; Lozano, A.M. Chronic anterior thalamus stimulation for intractable epilepsy. Epilepsia 2002, 43, 603–608. [Google Scholar] [CrossRef] [PubMed]
  14. Kerrigan, J.F.; Litt, B.; Fisher, R.S.; Cranstoun, S.; French, J.A.; Blum, D.E.; Dichter, M.; Shetter, A.; Baltuch, G.; Jaggi, J.; et al. Electrical stimulation of the anterior nucleus of the thalamus for the treatment of intractable epilepsy. Epilepsia 2004, 45, 346–354. [Google Scholar] [CrossRef]
  15. Andrade, D.; Zumsteg, D.; Hamani, C.; Hodaie, M.; Sarkissian, S.; Lozano, A.; Wennberg, R. Long-term follow-up of patients with thalamic deep brain stimulation for epilepsy. Neurology 2006, 66, 1571–1573. [Google Scholar] [CrossRef]
  16. Lee, K.; Jang, K.; Shon, Y. Chronic deep brain stimulation of subthalamic and anterior thalamic nuclei for controlling refractory partial epilepsy. In Advances in Functional and Reparative Neurosurgery; Springer: Berlin/Heidelberg, Germany, 2006; pp. 87–91. [Google Scholar]
  17. Osorio, I.; Overman, J.; Giftakis, J.; Wilkinson, S.B. High frequency thalamic stimulation for inoperable mesial temporal epilepsy. Epilepsia 2007, 48, 1561–1571. [Google Scholar] [CrossRef] [PubMed]
  18. Lee, K.J.; Shon, Y.M.; Cho, C.B. Long-term outcome of anterior thalamic nucleus stimulation for intractable epilepsy. Stereotact. Funct. Neurosurg. 2012, 90, 379–385. [Google Scholar] [CrossRef]
  19. Fisher, R.; Salanova, V.; Witt, T.; Worth, R.; Henry, T.; Gross, R.; Oommen, K.; Osorio, I.; Nazzaro, J.; Labar, D.; et al. Electrical stimulation of the anterior nucleus of thalamus for treatment of refractory epilepsy. Epilepsia 2010, 51, 899–908. [Google Scholar] [CrossRef] [PubMed]
  20. Salanova, V.; Witt, T.; Worth, R.; Henry, T.R.; Gross, R.E.; Nazzaro, J.M.; Labar, D.; Sperling, M.R.; Sharan, A.; Sandok, E.; et al. Long-term efficacy and safety of thalamic stimulation for drug-resistant partial epilepsy. Neurology 2015, 84, 1017–1025. [Google Scholar] [CrossRef]
  21. Lehtimäki, K.; Möttönen, T.; Järventausta, K.; Katisko, J.; Tähtinen, T.; Haapasalo, J.; Niskakangas, T.; Kiekara, T.; Öhman, J.; Peltola, J. Outcome based definition of the anterior thalamic deep brain stimulation target in refractory epilepsy. Brain Stimul. 2016, 9, 268–275. [Google Scholar] [CrossRef]
  22. Järvenpää, S.; Lehtimäki, K.; Rainesalo, S.; Möttönen, T.; Peltola, J. Improving the effectiveness of ANT DBS therapy for epilepsy with optimal current targeting. Epilepsia Open 2020, 5, 406–417. [Google Scholar] [CrossRef]
  23. Buentjen, L.; Kopitzki, K.; Schmitt, F.C.; Voges, J.; Tempelmann, C.; Kaufmann, J.; Kanowski, M. Direct targeting of the thalamicanteroventral nucleus for deep brain stimulation by T1-weighted magnetic resonance imaging at 3 T. Stereotact. Funct. Neurosurg. 2014, 92, 25–30. [Google Scholar] [CrossRef] [PubMed]
  24. Möttönen, T.; Katisko, J.; Haapasalo, J.; Tähtinen, T.; Kiekara, T.; Kähärä, V.; Peltola, J.; Öhman, J.; Lehtimäki, K. Defining the anterior nucleus of the thalamus (ANT) as a deep brain stimulation target in refractory epilepsy: Delineation using 3 T MRI and intraoperative microelectrode recording. NeuroImage Clin. 2015, 7, 823–829. [Google Scholar] [CrossRef] [PubMed]
  25. Hodaie, M.; Cordella, R.; Lozano, A.M.; Wennberg, R.; Dostrovsky, J.O. Bursting activity of neurons in the human anterior thalamic nucleus. Brain Res. 2006, 1115, 1–8. [Google Scholar] [CrossRef] [PubMed]
  26. Ohara, S.; Taghva, A.; Kim, J.; Lenz, F. Spontaneous low threshold spike bursting in awake humans is different in different lateral thalamic nuclei. Exp. Brain Res. 2007, 180, 281–288. [Google Scholar] [CrossRef]
  27. Schaper, F.L.; Zhao, Y.; Janssen, M.L.; Wagner, G.L.; Colon, A.J.; Hilkman, D.M.; Gommer, E.; Vlooswijk, M.C.; Hoogland, G.; Ackermans, L.; et al. Single-cell recordings to target the anterior nucleus of the thalamus in deep brain stimulation for patients with refractory epilepsy. Int. J. Neural Syst. 2019, 29, 1850012. [Google Scholar] [CrossRef] [PubMed]
  28. Möttönen, T.; Katisko, J.; Haapasalo, J.; Tähtinen, T.; Saastamoinen, A.; Peltola, J.; Öhman, J.; Lehtimäki, K. The correlation between intraoperative microelectrode recording and 3-Tesla MRI in patients undergoing ANT-DBS for refractory epilepsy. Stereotact. Funct. Neurosurg. 2016, 94, 86–92. [Google Scholar] [CrossRef]
  29. Levy, R.; Ashby, P.; Hutchison, W.D.; Lang, A.E.; Lozano, A.M.; Dostrovsky, J.O. Dependence of subthalamic nucleus oscillations on movement and dopamine in Parkinson’s disease. Brain 2002, 125, 1196–1209. [Google Scholar] [CrossRef]
  30. Priori, A.; Foffani, G.; Pesenti, A.; Tamma, F.; Bianchi, A.M.; Pellegrini, M.; Locatelli, M.; Moxon, K.; Villani, R. Rhythm-specific pharmacological modulation of subthalamic activity in Parkinson’s disease. Exp. Neurol. 2004, 189, 369–379. [Google Scholar] [CrossRef]
  31. Brown, P.; Williams, D. Basal ganglia local field potential activity: Character and functional significance in the human. Clin. Neurophysiol. 2005, 116, 2510–2519. [Google Scholar] [CrossRef]
  32. Kühn, A.A.; Trottenberg, T.; Kivi, A.; Kupsch, A.; Schneider, G.H.; Brown, P. The relationship between local field potential and neuronal discharge in the subthalamic nucleus of patients with Parkinson’s disease. Exp. Neurol. 2005, 194, 212–220. [Google Scholar] [CrossRef]
  33. Zaidel, A.; Spivak, A.; Grieb, B.; Bergman, H.; Israel, Z. Subthalamic span of β oscillations predicts deep brain stimulation efficacyfor patients with Parkinson’s disease. Brain 2010, 133, 2007–2021. [Google Scholar] [CrossRef] [PubMed]
  34. Marceglia, S.; Fumagalli, M.; Priori, A. What neurophysiological recordings tell us about cognitive and behavioral functions of the human subthalamic nucleus. Expert Rev. Neurother. 2011, 11, 139–149. [Google Scholar] [CrossRef]
  35. Rosa, M.; Marceglia, S.; Barbieri, S.; Priori, A. Local Field Potential and Deep Brain Stimulation (DBS). In Encyclopedia of Computational Neuroscience; Springer: New York, NY, USA, 2014. [Google Scholar]
  36. Kühn, A.A.; Kupsch, A.; Schneider, G.H.; Brown, P. Reduction in subthalamic 8–35 Hz oscillatory activity correlates with clinicalimprovement in Parkinson’s disease. Eur. J. Neurosci. 2006, 23, 1956–1960. [Google Scholar] [CrossRef]
  37. Mann, J.M.; Foote, K.D.; Garvan, C.W.; Fernandez, H.H.; Jacobson, C.E.; Rodriguez, R.L.; Haq, I.; Siddiqui, M.S.; Malaty, I.A.; Morishita, T.; et al. Brain penetration effects of microelectrodes and DBS leads in STN or GPi. J. Neurol. Neurosurg. Psychiatry 2009, 80, 794–798. [Google Scholar] [CrossRef] [PubMed]
  38. Constantoyannis, C.; Berk, C.; Honey, C.R.; Mendez, I.; Brownstone, R.M. Reducing hardware-related complications of deep brain stimulation. Can. J. Neurol. Sci. 2005, 32, 194–200. [Google Scholar] [CrossRef] [PubMed]
  39. Quinn, E.J.; Blumenfeld, Z.; Velisar, A.; Koop, M.M.; Shreve, L.A.; Trager, M.H.; Hill, B.C.; Kilbane, C.; Henderson, J.M.; Brontë- Stewart, H. Beta oscillations in freely moving Parkinson’s subjects are attenuated during deep brain stimulation. Mov. Disord. 2015, 30, 1750–1758. [Google Scholar] [CrossRef]
  40. Blumenfeld, Z.; Koop, M.M.; Prieto, T.E.; Shreve, L.A.; Velisar, A.; Quinn, E.J.; Trager, M.H.; Brontë-Stewart, H. Sixty-hertz stimulation improves bradykinesia and amplifies subthalamic low-frequency oscillations. Mov. Disord. 2017, 32, 80–88. [Google Scholar] [CrossRef]
  41. Neumann, W.J.; Staub-Bartelt, F.; Horn, A.; Schanda, J.; Schneider, G.H.; Brown, P.; Kühn, A.A. Long term correlation of subthalamic beta band activity with motor impairment in patients with Parkinson’s disease. Clin. Neurophysiol. 2017, 128, 2286–2291. [Google Scholar] [CrossRef]
  42. Syrkin-Nikolau, J.; Koop, M.M.; Prieto, T.; Anidi, C.; Afzal, M.F.; Velisar, A.; Blumenfeld, Z.; Martin, T.; Trager, M.; Bronte-Stewart, H. Subthalamic neural entropy is a feature of freezing of gait in freely moving people with Parkinson’s disease. Neurobiol. Dis. 2017, 108, 288–297. [Google Scholar] [CrossRef]
  43. Anidi, C.; O’Day, J.J.; Anderson, R.W.; Afzal, M.F.; Syrkin-Nikolau, J.; Velisar, A.; Bronte-Stewart, H.M. Neuromodulation targets pathological not physiological beta bursts during gait in Parkinson’s disease. Neurobiol. Dis. 2018, 120, 107–117. [Google Scholar] [CrossRef]
  44. Hell, F.; Plate, A.; Mehrkens, J.H.; Bötzel, K. Subthalamic oscillatory activity and connectivity during gait in Parkinson’s disease. NeuroImage: Clin. 2018, 19, 396–405. [Google Scholar] [CrossRef] [PubMed]
  45. Maling, N.; Lempka, S.F.; Blumenfeld, Z.; Bronte-Stewart, H.; McIntyre, C.C. Biophysical basis of subthalamic local field potentials recorded from deep brain stimulation electrodes. J. Neurophysiol. 2018, 120, 1932–1944. [Google Scholar] [CrossRef] [PubMed]
  46. Swann, N.C.; de Hemptinne, C.; Thompson, M.C.; Miocinovic, S.; Miller, A.M.; Gilron, R.; Ostrem, J.L.; Chizeck, H.J.; Starr, P.A. Adaptive deep brain stimulation for Parkinson’s disease using motor cortex sensing. J. Neural Eng. 2018, 15, 046006. [Google Scholar] [CrossRef]
  47. Ramirez Pasos, U.E.; Steigerwald, F.; Reich, M.M.; Matthies, C.; Volkmann, J.; Reese, R. Levodopa modulates functional connectivity in the upper beta band between subthalamic nucleus and muscle activity in tonic and phasic motor activity patterns in Parkinson’s disease. Front. Hum. Neurosci. 2019, 13, 223. [Google Scholar] [CrossRef] [PubMed]
  48. Yang, D.; Wang, Q.; Xu, C.; Fang, F.; Fan, J.; Li, L.; Du, Q.; Zhang, R.; Wang, Y.; Lin, Y.; et al. Transcranial direct current stimulation reduces seizure frequency in patients with refractory focal epilepsy: A randomized, double-blind, sham-controlled, and three-arm parallel multicenter study. Brain Stimul. 2020, 13, 109–116. [Google Scholar]
  49. Goyal, A.; Goetz, S.; Stanslaski, S.; Oh, Y.; Rusheen, A.E.; Klassen, B.; Miller, K.; Blaha, C.D.; Bennet, K.E.; Lee, K. The development of an implantable deep brain stimulation device with simultaneous chronic electrophysiological recording and stimulation in humans. Biosens. Bioelectron. 2021, 176, 112888. [Google Scholar] [CrossRef]
  50. Jimenez-Shahed, J. Device profile of the percept PC deep brain stimulation system for the treatment of Parkinson’s disease and related disorders. Expert Rev. Med. Devices 2021, 18, 319–332. [Google Scholar] [CrossRef]
  51. Osorio, I.; Zaveri, H.P.; Frei, M.G.; Arthurs, S. Epilepsy: The Intersection of Neurosciences, Biology, Mathematics, Engineering, and Physics; CRC Press: Boca Raton, FL, USA, 2016. [Google Scholar]
  52. Zaveri, H.P.; Duckrow, R.B.; Spencer, S.S. On the use of bipolar montages for time-series analysis of intracranial electroencephalograms. Clin. Neurophysiol. 2006, 117, 2102–2108. [Google Scholar] [CrossRef]
  53. Bastos, A.M.; Schoffelen, J.M. A tutorial review of functional connectivity analysis methods and their interpretational pitfalls. Front. Syst. Neurosci. 2016, 9, 175. [Google Scholar]
  54. Georgiy, R.; John, B.; Martha, J.; Richard, N. Patient-specific early seizure detection from scalp electroencephalogram. J. Clin. Neurophysiol. 2010, 27, 163–178. [Google Scholar]
  55. Wold, S.; Esbensen, K.; Geladi, P. Principal component analysis. Chemom. Intell. Lab. Syst. 1987, 2, 37–52. [Google Scholar]
  56. Machine Learning with Matlab. Available online: https://www.mathworks.com/matlabcentral/fileexchange/42744-machine-learning-with-matlab (accessed on 2 June 2022).
Figure 1. DBS lead localization performed by the BrainLab Elements. Images above show the perspective where the electrodes have a higher level of contact with the brain structures.
Figure 1. DBS lead localization performed by the BrainLab Elements. Images above show the perspective where the electrodes have a higher level of contact with the brain structures.
Sensors 22 06601 g001
Figure 2. (A) Medtronic 3389 lead DBS model: each contact has a length of 1.5 mm, and the distance between two adjacent contacts is 0.5 mm; (B) LFP spatial locations of bipolar signals. Each spatial location corresponds to the midpoint between the considered channels.
Figure 2. (A) Medtronic 3389 lead DBS model: each contact has a length of 1.5 mm, and the distance between two adjacent contacts is 0.5 mm; (B) LFP spatial locations of bipolar signals. Each spatial location corresponds to the midpoint between the considered channels.
Sensors 22 06601 g002
Figure 3. Schematic representation of the general methodology used.
Figure 3. Schematic representation of the general methodology used.
Sensors 22 06601 g003
Figure 4. Ranking features results. Abbreviations—Spectral features: bpd (band power delta), bpt (band power theta), bpa (band power alpha), bpb (band power beta), bpg (band power gamma); Morphological features: am (absolute mean), mp (mean peaks), max (maximum), min (minimum); Multivariate features: crosscorr (cross-correlation), coh (coherence), plv (phase locking value), pli (phase lag index); Statistical features: var (variance), kurt (kurtosis), md (median), skew (skewness).
Figure 4. Ranking features results. Abbreviations—Spectral features: bpd (band power delta), bpt (band power theta), bpa (band power alpha), bpb (band power beta), bpg (band power gamma); Morphological features: am (absolute mean), mp (mean peaks), max (maximum), min (minimum); Multivariate features: crosscorr (cross-correlation), coh (coherence), plv (phase locking value), pli (phase lag index); Statistical features: var (variance), kurt (kurtosis), md (median), skew (skewness).
Sensors 22 06601 g004
Table 1. Local field potential studies. Abbreviations: x—not applicable; STN—subthalamic nucleus; ACC—anterior cingulate cortex.
Table 1. Local field potential studies. Abbreviations: x—not applicable; STN—subthalamic nucleus; ACC—anterior cingulate cortex.
ReferenceDescriptionNeurostimulatorAnatomical TargetNumber of Patients
Jimenz-Shahed et al. (2021) [50]Description of the sensing capabilities of the neurostimulatorPercept PCxx
Goyal et al. (2021) [49]Description of the sensing capabilities of the neurostimulatorPercept PCxx
Yang et al. (2020) [48]Study of which stimulation parameters provide pain relief without triggering after dischargesActiva PC+SACC3
Passos et al. (2019) [47]Study the impact of dopaminergic state and movement of beta functional connectivity between basal ganglia and lower motor neuronsActiva PC+SSTN8
Anidi et al. (2019) [43]Study whether STN-DBS affects beta burst dynamics and gait impairments differentially in freezers and non-freezersActiva PC+SSTN12
Hell et al. (2018) [44]Study of the STN function during gaitActiva PC+SSTN10
Maling et al. (2018) [45]Development of a patient-specific computational framework to analyze LFP recordingsActiva PC+SSTN1
Swann et al. (2018) [46]Study of the disfunction of disturbed neural networks in FLPsActiva PC+SSTN5
Neumann et al. (2017) [41]Study of long-term association of STN beta activity with parkinsonian motor signsActiva PC+SSTN15
Syrkin-Nikolau et al. (2017) [42]Study of STN neural features of freezers and non-freezersActiva PC+SSTN14
Blumenfeld et al. (2016) [40]Study whether the STN alpha/beta oscillation attenuation is causal to the bradykinesia improvementActiva PC+SSTN9
Quinn et al. (2015) [39]Study whether beta power is similar in different resting postures during forward walkingActiva PC+SSTN15
Table 2. Correspondence between bipolar LFP signal contacts, their spatial locations and the proportion of ANT target underlying each contact. Abbreviations—RH: right hemisphere; LH: left hemisphere; −: completely outside to the target; +: marginally in the target; ++: partially but not completely inserted in the target; +++: completely inserted in the target.
Table 2. Correspondence between bipolar LFP signal contacts, their spatial locations and the proportion of ANT target underlying each contact. Abbreviations—RH: right hemisphere; LH: left hemisphere; −: completely outside to the target; +: marginally in the target; ++: partially but not completely inserted in the target; +++: completely inserted in the target.
ContactsPhysical LocationMontageTarget % (LH)Target % (RH)
0–1APass2+
0–2BPass1+
1–2CPass2Pass2++
0–3CPass1Pass1++
1–3DPass1++++
2–3EPass2+++++
Table 3. List of extracted features.
Table 3. List of extracted features.
Feature TypeFeatureDescription
StatisticalVariance (var)Difference (spread) between the normalized squared sum of instantaneous values with the mean value
Skewness (skew)Distortion or asymmetry of the probability density function of the amplitude of time-series
Kurtosis (kurt)Sharpness of the probability density function of the amplitude of time-series
Median (md)Value separating the higher half from the lower half of a data sample
SpectralBand power delta (bpd)Power spectral density in the delta frequency band (1–3 Hz)
Band power theta (pbt)Power spectral density in the theta frequency band (4–7 Hz)
Band power alpha (pba)Power spectral density in the alpha frequency band (8–12 Hz)
Band power beta (bpb)Power spectral density in the beta frequency band (13–30 Hz)
Band power gamma (bpg)Power spectral density in the gamma frequency band (3–100 Hz)
MorphologicalAbsolute mean (am)Absolute average value of a data sample
Mean peaks (mp)Absolute average of the maximum values of a data sample
Maximum (max)Maximum value of a data sample
Minimum (min)Minimum value of a data sample
MultivariateCross-correlation (crosscorr)Measure of similarity of two time-series
Coherence (coh)Measures the causal relationship between two signals
Phase locking value (plv)Measure of the phase synchrony between two time-series
Phase lag index (pli)Evaluates the phase difference distribution across observations
Table 4. Classifiers and respective varying and standard model parameters.
Table 4. Classifiers and respective varying and standard model parameters.
ClassifierModel ParametersStandard
SVMKernel function: linear, radial, basis function (rbf), polynomialrbf
kNNDistance metrics: squared Euclidean (seuclidean), Euclidean, correlation, Spearmanseuclidean
NNHidden layer size (1:10)5
Table 5. Prediction of the target-% of lead DBS physical locations.
Table 5. Prediction of the target-% of lead DBS physical locations.
ModelSpatial Locations to Predict
Pass1CPass1-Left, CPass1-Right, BRight, DLeft, ALeft and ERight
Pass2CPass2-Left, CPass2-Right, ARight, ELeft, BLeft and DRight
Table 6. Pass1 and Pass2 classifiers ANT target percentage predictions. Pass1 classifier predicted ANT target percentages of Pass1 spatial locations (BRight, CPass1-Left and CPass1-Right and DLeft) and Pass2 spatial locations used as a Pass2 classifier class labels (ALeft and ERight); Pass2 classifier predicted ANT target percentages of Pass2 spatial locations (ARight, CPass2-Left and CPass2-Right and ELeft) and Pass1 spatial locations used as a Pass1-classifier class labels (BLeft and DRight). Gray cells indicate the predictions made by the opposite classifier. High prediction values for both hemispheres are written in bold. Abbreviations—LH: left hemisphere; RH: right hemisphere: completely outside the target; +: marginally in the target; ++: partially but not completely inserted in the target; +++: completely inserted in the target.
Table 6. Pass1 and Pass2 classifiers ANT target percentage predictions. Pass1 classifier predicted ANT target percentages of Pass1 spatial locations (BRight, CPass1-Left and CPass1-Right and DLeft) and Pass2 spatial locations used as a Pass2 classifier class labels (ALeft and ERight); Pass2 classifier predicted ANT target percentages of Pass2 spatial locations (ARight, CPass2-Left and CPass2-Right and ELeft) and Pass1 spatial locations used as a Pass1-classifier class labels (BLeft and DRight). Gray cells indicate the predictions made by the opposite classifier. High prediction values for both hemispheres are written in bold. Abbreviations—LH: left hemisphere; RH: right hemisphere: completely outside the target; +: marginally in the target; ++: partially but not completely inserted in the target; +++: completely inserted in the target.
Spatial LocationTarget (LH)Predictions (LH)Target (RH)Predictions (RH)
A54.35%+50.00%
B63.04%+56.52%
CPass141.30%++52.17%
CPass2-52.17%++54.35%
D+40.87%+++60.87%
E++69.57%+++60.87%
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Lopes, E.M.; Rego, R.; Rito, M.; Chamadoira, C.; Dias, D.; Cunha, J.P.S. Estimation of ANT-DBS Electrodes on Target Positioning Based on a New PerceptTM PC LFP Signal Analysis. Sensors 2022, 22, 6601. https://doi.org/10.3390/s22176601

AMA Style

Lopes EM, Rego R, Rito M, Chamadoira C, Dias D, Cunha JPS. Estimation of ANT-DBS Electrodes on Target Positioning Based on a New PerceptTM PC LFP Signal Analysis. Sensors. 2022; 22(17):6601. https://doi.org/10.3390/s22176601

Chicago/Turabian Style

Lopes, Elodie Múrias, Ricardo Rego, Manuel Rito, Clara Chamadoira, Duarte Dias, and João Paulo Silva Cunha. 2022. "Estimation of ANT-DBS Electrodes on Target Positioning Based on a New PerceptTM PC LFP Signal Analysis" Sensors 22, no. 17: 6601. https://doi.org/10.3390/s22176601

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop