Objective. Large structural brain changes, such as chronic stroke lesions, alter the current pathways throughout the patients' head and therefore have to be taken into account when performing transcranial direct current stimulation simulations. Approach. We implement, test and distribute the first MATLAB pipeline that automatically generates realistic and individualized volume conduction head models of chronic stroke patients, by combining the already existing software SimNIBS, for the mesh generation, and lesion identification with neighborhood data analysis, for the lesion identification. To highlight the impact of our pipeline, we investigated the sensitivity of the electric field distribution to the lesion location and lesion conductivity in 16 stroke patients' datasets. Main results. Our pipeline automatically generates 1 mm-resolution tetrahedral meshes including the lesion compartment in less than three hours. Moreover, for large lesions, we found a high sensitivity of the electric field distribution to the lesion conductivity value and location. Significance. This work facilitates optimizing electrode configurations with the goal to obtain more focal brain stimulations of the target volumes in rehabilitation for chronic stroke patients.

The inconsistent response to transcranial electric stimulation in the stroke population is attributed to, among other factors, unknown effects of stroke lesion conductivity on stimulation strength at the targeted brain areas. Volume conduction models are promising tools to determine optimal stimulation settings. However, stroke lesion conductivity is often not considered in these models as a source of inter-subject variability. The goal of this study is to propose a method that combines MRI, EEG, and transcranial stimulation to estimate the conductivity of cortical stroke lesions experimentally. In this simulation study, lesion conductivity was estimated from scalp potentials during transcranial electric stimulation in 12 chronic stroke patients. To do so, first, we determined the stimulation configuration where scalp potentials are maximally affected by the lesion. Then, we calculated scalp potentials in a model with a fixed lesion conductivity and a model with a randomly assigned conductivity. To estimate the lesion conductivity, we minimized the error between the two models by varying the conductivity in the second model. Finally, to reflect realistic experimental conditions, we test the effect rotation of measurement electrode orientation and the effect of the number of electrodes used. We found that the algorithm converged to the correct lesion conductivity value when noise on the electrode positions was absent for all lesions. Conductivity estimation error was below 5% with realistic electrode coregistration errors of 0.1° for lesions larger than 50 ml. Higher lesion conductivities and lesion volumes were associated with smaller estimation errors. In conclusion, this method can experimentally estimate stroke lesion conductivity, improving the accuracy of volume conductor models of stroke patients and potentially leading to more effective transcranial electric stimulation configurations for this population.