Rapid whole-brain electric field mapping in transcranial magnetic stimulation using deep learning

A review of algorithms and software for real-time electric field modeling techniques for transcranial magnetic stimulation

Transcranial magnetic stimulation (TMS) is a device-based neuromodulation technique increasingly used to treat brain diseases. Electric field (E-field) modeling is an important technique in several TMS clinical applications, including the precision stimulation of brain targets with accurate stimulation density for the treatment of mental disorders and the localization of brain function areas for neurosurgical planning. Classical methods for E-field modeling usually take a long computation time. Fast algorithms are usually developed with significantly lower spatial resolutions that reduce the prediction accuracy and limit their usage in real-time or near real-time TMS applications. This review paper discusses several modern algorithms for real-time or near real-time TMS E-field modeling and their advantages and limitations. The reviewed methods include techniques such as basis representation techniques and deep neural-network-based methods. This paper also provides a review of software tools that can integrate E-field modeling with navigated TMS, including a recent software for real-time navigated E-field mapping based on deep neural-network models.

ViTab Transformer Framework for Predicting Induced Electric Field and Focality in Transcranial Magnetic Stimulation

Transcranial magnetic stimulation is an electromagnetic induction-based non-invasive therapeutic technique for neurological diseases. For finding new clinical applications and enhancing the efficacy of TMS in existing neurological disorders, the current study focuses on a deep learning-based prediction model as an alternative to time-consuming electromagnetic (EM) simulation software. The main bottleneck of the existing prediction models is to consider very few input parameters of a standard coil such as coil type and coil position for predicting an output of electric field value. To overcome this limitation, a transformer-based prediction model titled as ViTab transformer is developed in this work to predict electric field (E-max), focality or area of stmulation (S-half), and volume of stimulation (V-half) by considering several input parameters such as sources of MRI images, types of coils, coil position, rate of change of current, brain tissues conductivity, and coil distance from the scalp. The proposed framework consists of a vision and a tab transformer to handle both image and tabular-type data. The prediction performance of the offered model is evaluated in terms of coefficient determination, R2 score, for E-max, V-half, and S-half in the testing phase. The obtained result in terms of R2 score for E-max, V-half, and S-half are found 0.97, 0.87, and 0.90 respectively. The results indicate that the suggested ViTab transformer model can predict electric field as well as focality more accurately than the current state-of-the-art methods. The reduced computational time, as well as efficient prediction accuracy, resembles that ViTab transformer can assist the neuroscientist and neurosurgeon prior to providing superior TMS treatment in near future.

Fast and accurate computational E-field dosimetry for group-level transcranial magnetic stimulation targeting

Transcranial magnetic stimulation (TMS) is used to study brain function and treat mental health disorders. During TMS, a coil placed on the scalp induces an E-field in the brain that modulates its activity. TMS is known to stimulate regions that are exposed to a large E-field. Clinical TMS protocols prescribe a coil placement based on scalp landmarks. There are inter-individual variations in brain anatomy that result in variations in the TMS-induced E-field at the targeted region and its outcome. These variations across individuals could in principle be minimized by developing a large database of head subjects and determining scalp landmarks that maximize E-field at the targeted brain region while minimizing its variation using computational methods. However, this approach requires repeated execution of a computational method to determine the E-field induced in the brain for a large number of subjects and coil placements. We developed a probabilistic matrix decomposition-based approach for rapidly evaluating the E-field induced during TMS for a large number of coil placements due to a pre-defined coil model. Our approach can determine the E-field induced in over 1 Million coil placements in 9.5 h, in contrast, to over 5 years using a brute-force approach. After the initial set-up stage, the E-field can be predicted over the whole brain within 2-3 ms and to 2% accuracy. We tested our approach in over 200 subjects and achieved an error of <2% in most and <3.5% in all subjects. We will present several examples of bench-marking analysis for our tool in terms of accuracy and speed. Furthermore, we will show the methods' applicability for group-level optimization of coil placement for illustration purposes only. The software implementation link is provided in the appendix.

In-vivoverified anatomically aware deep learning for real-time electric field simulation

Objective.Transcranial magnetic stimulation (TMS) has emerged as a prominent non-invasive technique for modulating brain function and treating mental disorders. By generating a high-precision magnetically evoked electric field (E-field) using a TMS coil, it enables targeted stimulation of specific brain regions. However, current computational methods employed for E-field simulations necessitate extensive preprocessing and simulation time, limiting their fast applications in the determining the optimal coil placement.Approach.We present an attentional deep learning network to simulate E-fields. This network takes individual magnetic resonance images and coil configurations as inputs, firstly transforming the images into explicit brain tissues and subsequently generating the local E-field distribution near the target brain region. Main results. Relative to the previous deep-learning simulation method, the presented method reduced the mean relative error in simulated E-field strength of gray matter by 21.1%, and increased the correlation between regional E-field strengths and corresponding electrophysiological responses by 35.0% when applied into another dataset.In-vivoTMS experiments further revealed that the optimal coil placements derived from presented method exhibit comparable stimulation performance on motor evoked potentials to those obtained using computational methods. The simplified preprocessing and increased simulation efficiency result in a significant reduction in the overall time cost of traditional TMS coil placement optimization, from several hours to mere minutes.Significance.The precision and efficiency of presented simulation method hold promise for its application in determining individualized coil placements in clinical practice, paving the way for personalized TMS treatments.

Real-Time Computation of Brain E-Field for Enhanced Transcranial Magnetic Stimulation Neuronavigation and Optimization

Transcranial Magnetic Stimulation (TMS) coil placement and pulse waveform current are often chosen to achieve a specified E-field dose on targeted brain regions. TMS neuronavigation could be improved by including real-time accurate distributions of the E-field dose on the cortex. We introduce a method and develop software for computing brain E-field distributions in real-time enabling easy integration into neuronavigation and with the same accuracy as 1st -order finite element method (FEM) solvers. Initially, a spanning basis set (< 400) of E-fields generated by white noise magnetic currents on a surface separating the head and permissible coil placements are orthogonalized to generate the modes. Subsequently, Reciprocity and Huygens' principles are utilized to compute fields induced by the modes on a surface separating the head and coil by FEM, which are used in conjunction with online (real-time) computed primary fields on the separating surface to evaluate the mode expansion. We conducted a comparative analysis of E-fields computed by FEM and in real-time for eight subjects, utilizing two head model types (SimNIBS's 'headreco' and 'mri2mesh' pipeline), three coil types (circular, double-cone, and Figure-8), and 1000 coil placements (48,000 simulations). The real-time computation for any coil placement is within 4 milliseconds (ms), for 400 modes, and requires less than 4 GB of memory on a GPU. Our solver is capable of computing E-fields within 4 ms, making it a practical approach for integrating E-field information into the neuronavigation systems without imposing a significant overhead on frame generation (20 and 50 frames per second within 50 and 20 ms, respectively).

[The influence of tissue conductivity on the calculation of electric field in the transcranial magnetic stimulation head model]

In transcranial magnetic stimulation (TMS), the conductivity of brain tissue is obtained by using diffusion tensor imaging (DTI) data processing. However, the specific impact of different processing methods on the induced electric field in the tissue has not been thoroughly studied. In this paper, we first used magnetic resonance image (MRI) data to create a three-dimensional head model, and then estimated the conductivity of gray matter (GM) and white matter (WM) using four conductivity models, namely scalar (SC), direct mapping (DM), volume normalization (VN) and average conductivity (MC), respectively. Isotropic empirical conductivity values were used for the conductivity of other tissues such as the scalp, skull, and cerebrospinal fluid (CSF), and then the TMS simulations were performed when the coil was parallel and perpendicular to the gyrus of the target. When the coil was perpendicular to the gyrus where the target was located, it was easy to get the maximum electric field in the head model. The maximum electric field in the DM model was 45.66% higher than that in the SC model. The results showed that the conductivity component along the electric field direction of which conductivity model was smaller in TMS, the induced electric field in the corresponding domain corresponding to the conductivity model was larger. This study has guiding significance for TMS precise stimulation.

Fast And Accurate Population Level Transcranial Magnetic Stimulation via Low-Rank Probabilistic Matrix Decomposition (PMD)

Transcranial magnetic stimulation (TMS) is used to study brain function and treat mental health disorders. During TMS, a coil placed on the scalp induces an E-field in the brain that modulates its activity. TMS is known to stimulate regions that are exposed to a large E-field. Clinical TMS protocols prescribe a coil placement based on scalp landmarks. There are inter-individual variations in brain anatomy that result in variations in the TMS-induced E-field at the targeted region and its outcome. These variations across individuals could in principle be minimized by developing a large database of head subjects and determining scalp landmarks that maximize E-field at the targeted brain region while minimizing its variation using computational methods. However, this approach requires repeated execution of a computational method to determine the E-field induced in the brain for a large number of subjects and coil placements. We developed a probabilistic matrix decomposition-based approach for rapidly evaluating the E-field induced during TMS for a large number of coil placements. Our approach can determine the E-field induced in over 1 Million coil placements in 9.5 hours, in contrast, to over 5 years using a bruteforce approach. After the initial set-up stage, the E-field can be predicted over the whole brain within 2-3 milliseconds and to 2% accuracy. We tested our approach in over 200 subjects and achieved an error of < 2% in most and < 3.5% in all subjects. We will present several examples of bench-marking analysis for our tool in terms of accuracy and speed across and its applicability for population level optimization of coil placement.

Highlights

A method for practical E-field informed population-level TMS coil placement strategies is developed.This algorithm enables the determination of E-field informed optimal coil placement in seconds enabling it’s use for close-loop and on-the-fly reconfiguration of TMS.After the initial set-up stage of less than 10 hours, the E-field can be predicted for any coil placement across the whole brain within in 2-3 milliseconds.