James Avery

A Versatile and Reproducible Multi-Frequency Electrical Impedance Tomography System

A highly versatile Electrical Impedance Tomography (EIT) system, nicknamed the ScouseTom, has been developed. The system allows control over current amplitude, frequency, number of electrodes, injection protocol and data processing. Current is injected using a Keithley 6221 current source, and voltages are recorded with a 24-bit EEG system with minimum bandwidth of 3.2 kHz. Custom PCBs interface with a PC to control the measurement process, electrode addressing and triggering of external stimuli. The performance of the system was characterised using resistor phantoms to represent human scalp recordings, with an SNR of 77.5 dB, stable across a four hour recording and 20 Hz to 20 kHz. In studies of both haeomorrhage using scalp electrodes, and evoked activity using epicortical electrode mats in rats, it was possible to reconstruct images matching established literature at known areas of onset. Data collected using scalp electrode in humans matched known tissue impedance spectra and was stable over frequency. The experimental procedure is software controlled and is readily adaptable to new paradigms. Where possible, commercial or open-source components were used, to minimise the complexity in reproduction. The hardware designs and software for the system have been released under an open source licence, encouraging contributions and allowing for rapid replication.

Characterisation and imaging of cortical impedance changes during interictal and ictal activity in the anaesthetised rat

Epilepsy affects approximately 50 million people worldwide, and 20–30% of these cases are refractory to antiepileptic drugs. Many patients with intractable epilepsy can benefit from surgical resection of the tissue generating the seizures; however, difficulty in precisely localising seizure foci has limited the number of patients undergoing surgery as well as potentially lowered its effectiveness. Here we demonstrate a novel imaging method for monitoring rapid changes in cerebral tissue impedance occurring during interictal and ictal activity, and show that it can reveal the propagation of pathological activity in the cortex. Cortical impedance was recorded simultaneously to ECoG using a 30-contact electrode mat placed on the exposed cortex of anaesthetised rats, in which interictal spikes (IISs) and seizures were induced by cortical injection of 4-aminopyridine (4-AP), picrotoxin or penicillin. We characterised the tissue impedance responses during IISs and seizures, and imaged these responses in the cortex using Electrical Impedance Tomography (EIT). We found a fast, transient drop in impedance occurring as early as 12 ms prior to the IISs, followed by a steep rise in impedance within ~ 120 ms of the IIS. EIT images of these impedance changes showed that they were co-localised and centred at a depth of 1 mm in the cortex, and that they closely followed the activity propagation observed in the surface ECoG signals. The fast, pre-IIS impedance drop most likely reflects synchronised depolarisation in a localised network of neurons, and the post-IIS impedance increase reflects the subsequent shrinkage of extracellular space caused by the intense activity. EIT could also be used to picture a steady rise in tissue impedance during seizure activity, which has been previously described. Thus, our results demonstrate that EIT can detect and localise different physiological changes during interictal and ictal activity and, in conjunction with ECoG, may in future improve the localisation of seizure foci in the clinical setting.

Comparison of total variation algorithms for electrical impedance tomography

The applications of total variation (TV) algorithms for electrical impedance tomography (EIT) have been investigated. The use of the TV regularisation technique helps to preserve discontinuities in reconstruction, such as the boundaries of perturbations and sharp changes in conductivity, which are unintentionally smoothed by traditional l2 norm regularisation. However, the non-differentiability of TV regularisation has led to the use of different algorithms. Recent advances in TV algorithms such as the primal dual interior point method (PDIPM), the linearised alternating direction method of multipliers (LADMM) and the spilt Bregman (SB) method have all been demonstrated successful EIT applications, but no direct comparison of the techniques has been made. Their noise performance, spatial resolution and convergence rate applied to time difference EIT were studied in simulations on 2D cylindrical meshes with different noise levels, 2D cylindrical tank and 3D anatomically head-shaped phantoms containing vegetable material with complex conductivity. LADMM had the fastest calculation speed but worst resolution due to the exclusion of the second-derivative; PDIPM reconstructed the sharpest change in conductivity but with lower contrast than SB; SB had a faster convergence rate than PDIPM and the lowest image errors.

Correcting electrode modelling errors in EIT on realistic 3D head models.

Electrical impedance tomography (EIT) is a promising medical imaging technique which could aid differentiation of haemorrhagic from ischaemic stroke in an ambulance. One challenge in EIT is the ill-posed nature of the image reconstruction, i.e., that small measurement or modelling errors can result in large image artefacts. It is therefore important that reconstruction algorithms are improved with regard to stability to modelling errors. We identify that wrongly modelled electrode positions constitute one of the biggest sources of image artefacts in head EIT. Therefore, the use of the Fréchet derivative on the electrode boundaries in a realistic three-dimensional head model is investigated, in order to reconstruct electrode movements simultaneously to conductivity changes. We show a fast implementation and analyse the performance of electrode position reconstructions in time-difference and absolute imaging for simulated and experimental voltages. Reconstructing the electrode positions and conductivities simultaneously increased the image quality significantly in the presence of electrode movement.

Characterising the frequency response of impedance changes during evoked physiological activity in the rat brain

OBJECTIVE Electrical impedance tomography (EIT) can image impedance changes associated with evoked physiological activity in the cerebral cortex using an array of epicortical electrodes. An impedance change is observed as the externally applied current, normally confined to the extracellular space is admitted into the conducting intracellular space during neuronal depolarisation. The response is largest at DC and decreases at higher frequencies due to capacitative transfer of current across the membrane. Biophysical modelling has shown that this effect becomes significant above 100 Hz. Recordings at DC, however, are contaminated by physiological endogenous evoked potentials. By moving to 1.7 kHz, images of somatosensory evoked responses have been produced down to 2 mm with a resolution of 2 ms and 200 μm. Hardware limitations have so far restricted impedance measurements to frequencies  <2 kHz. The purpose of this work was to establish the optimal frequency for extending EIT to image throughout the brain and to characterise the response at frequencies  >2 kHz using improved hardware. APPROACH Impedance changes were recorded during forepaw somatosensory stimulation in both cerebral cortex and the VPL nucleus of the thalamus in anaesthetised rats using applied currents of 1 kHz to 10 kHz. MAIN RESULTS In the cortex, impedance changed by -0.04 ± 0.02 % at 1 kHz, reached a peak of -0.13 ± 0.05 % at 1475 Hz and decreased to -0.05 ± 0.02 % at 10 kHz. At these frequencies, changes in the thalamus were -0.26 ± 0.1%, -0.4 ± 0.15 % and -0.08 ± 0.03 % respectively. The signal-to-noise ratio was also highest at 1475 Hz with values of -29.5 ± 8 and -31.6 ±10 recorded from the cortex and thalamus respectively. Signficance: This indicates that the optimal frequency for imaging cortical and thalamic evoked activity using fast neural EIT is 1475 Hz.

Feasibility of imaging epileptic seizure onset with EIT and depth electrodes

&NA; Imaging ictal and interictal activity with Electrical Impedance Tomography (EIT) using intracranial electrode mats has been demonstrated in animal models of epilepsy. In human epilepsy subjects undergoing presurgical evaluation, depth electrodes are often preferred. The purpose of this work was to evaluate the feasibility of using EIT to localise epileptogenic areas with intracranial electrodes in humans. The accuracy of localisation of the ictal onset zone was evaluated in computer simulations using 9M element FEM models derived from three subjects. 5 mm radius perturbations imitating a single seizure onset event were placed in several locations forming two groups: under depth electrode coverage and in the contralateral hemisphere. Simulations were made for impedance changes of 1% expected for neuronal depolarisation over milliseconds and 10% for cell swelling over seconds. Reconstructions were compared with EEG source modelling for a radially orientated dipole with respect to the closest EEG recording contact. The best accuracy of EIT was obtained using all depth and 32 scalp electrodes, greater than the equivalent accuracy with EEG inverse source modelling. The localisation error was 5.2 ± 1.8, 4.3 ± 0 and 46.2 ± 25.8 mm for perturbations within the volume enclosed by depth electrodes and 29.6 ± 38.7, 26.1 ± 36.2, 54.0 ± 26.2 mm for those without (EIT 1%, 10% change, EEG source modelling, n = 15 in 3 subjects, p < 0.01). As EIT was insensitive to source dipole orientation, all 15 perturbations within the volume enclosed by depth electrodes were localised, whereas the standard clinical method of visual inspection of EEG voltages, only localised 8 out of 15 cases. This suggests that adding EIT to SEEG measurements could be beneficial in localising the onset of seizures. HighlightsA new method is proposed to use EIT to localise and ictal activity in patients with depth electrodes.The location accuracy was improved with the best EIT protocol than EEG inverse source or SEEG detection in simulations.EIT was not sensitive to dipole orientation, while EEG detection varied with the field angle demonstrated in modelling.A combination of EIT and SEEG can potentially improve the diagnostic yield in epilepsy.

Optimisation of current injection protocol based on a region of interest.

OBJECTIVE Electrical impedance tomography has the potential to image fast neural activity associated with physiological or epileptic activity throughout the brain. These applications pose a particular challenge as expected voltage changes on the electrodes are less than 1% and geometrical constraints of the body under investigation mean that electrodes can not be evenly distributed around its boundary. Unlike other applications, however, information regarding the location of expected activity is typically available. An informative method for choosing current paths that maximise sensitivity to specific regions is desirable. APPROACH Two electrode addressing protocol generation methods based on current density vectors concentrated in a region of interest have been proposed. One focuses solely on maximising its magnitude while the other considers its distribution. The quality of reconstructed images using these protocols was assessed in a simulation study conducted in a human and rat mesh and compared to the protocol that maximises distance between injecting electrodes. MAIN RESULTS When implementing the protocol that focused on maximising magnitude, the current density concentrated in a region of interest increased by up to a factor of 3. When the distribution of the current was maximised, the spread of current density vectors increased by up to fivefold. For the small conductivity changes expected in the applications explored, image quality was best when implementing the protocol that maximised current density. The average image error when using this protocol was 7% better than when employing other protocols. SIGNIFICANCE We conclude that for fast neural EIT applications, the protocol that maximises current density is the best protocol to implement.

Frequency-dependent characterisation of impedance changes during epileptiform activity in a rat model of epilepsy

OBJECTIVE Electrical impedance tomography (EIT) can be used to image impedance changes associated with epileptiform activity and so holds therapeutic potential for improving presurgical localisation of the ictal onset zone in patients with treatment-resistant epilepsy. There are two principal impedance changes which occur during seizures that may be imaged with EIT: (a) a fast, transient impedance decrease over milliseconds due to hypersynchronous neuronal depolarisation in individual ictal discharges; and (b) a larger, slow impedance increase caused by cell swelling over the course of the seizure. The magnitude of these signals is highly dependent on the carrier frequency of applied current used for obtaining impedance measurements. The purpose of this work was to characterise the frequency response of the fast and slow impedance changes during epileptiform activity. APPROACH Seizures were induced in anaesthetised rats by electrically stimulating the cerebral cortex. During each seizure, impedance measurements were obtained by delivering 50 µA, through two electrodes on an epicortical array, at one of 20 frequencies in the 1-10 kHz range. Recordings were demodulated to determine the magnitude of fast and slow impedance responses at each frequency. MAIN RESULTS The fast impedance change during averaged ictal discharges reached a maximal amplitude and signal-to-noise ratio (SNR) of  -0.36%  ±  0.05% and 50.2  ±  11.3, respectively, at 1355 Hz. At this frequency, the slow impedance change had an amplitude of 4.61%  ±  1.32% and an SNR of 545  ±  125, which did not significantly change across frequency (p  >  0.01). SIGNIFICANCE We conclude that the optimal frequency for imaging epileptiform activity is 1355 Hz, which maximises the SNR of fast neural changes whilst enabling simultaneous measurement of slow changes. These findings will inform future investigations aimed at imaging epilepsy in subcortical brain structures, where SNR is considerably reduced, and those using parallel, multi-frequency EIT.

Feasibility of imaging evoked activity throughout the rat brain using electrical impedance tomography

Abstract Electrical Impedance Tomography (EIT) is an emerging technique which has been used to image evoked activity during whisker displacement in the cortex of an anaesthetised rat with a spatiotemporal resolution of 200 &mgr;m and 2 ms. The aim of this work was to extend EIT to image not only from the cortex but also from deeper structures active in somatosensory processing, specifically the ventral posterolateral (VPL) nucleus of the thalamus. The direct response in the cortex and VPL following 2 Hz forepaw stimulation were quantified using a 57‐channel epicortical electrode array and a 16‐channel depth electrode. Impedance changes of −0.16 ± 0.08% at 12.9 ± 1.4 ms and −0.41 ± 0.14% at 8.8Symbol1.9 ms were recorded from the cortex and VPL respectively. For imaging purposes, two 57‐channel epicortical electrode arrays were used with one placed on each hemisphere of the rat brain. Despite using parameters optimised toward measuring thalamic activity and undertaking extensive averaging, reconstructed activity was constrained to the cortical somatosensory forepaw region and no significant activity at a depth greater than 1.6 mm below the surface of the cortex could be reconstructed. An evaluation of the depth sensitivity of EIT was investigated in simulations using estimates of the conductivity change and noise levels derived from experiments. These indicate that EIT imaging with epicortical electrodes is limited to activity occurring 2.5 mm below the surface of the cortex. This depth includes the hippocampus and so EIT has the potential to image activity, such as epilepsy, originating from this structure. To image deeper activity, however, alternative methods such as the additional implementation of depth electrodes will be required to gain the necessary depth resolution. Symbol. No Caption available. HighlightsImpedance changes in the cortex and VPL have been characterised during forepaw stimulation.EIT was conducted with 114 epicortical electrodes placed over both hemispheres in the rat brain.Signal amplitude precluded imaging of ascending neural activity.Simulations show depth accuracy of EIT extends to superficial layers of the hippocampus.

Multi-frequency electrical impedance tomography and neuroimaging data in stroke patients

Electrical Impedance Tomography (EIT) is a non-invasive imaging technique, which has the potential to expedite the differentiation of ischaemic or haemorrhagic stroke, decreasing the time to treatment. Whilst demonstrated in simulation, there are currently no suitable imaging or classification methods which can be successfully applied to human stroke data. Development of these complex methods is hindered by a lack of quality Multi-Frequency EIT (MFEIT) data. To address this, MFEIT data were collected from 23 stroke patients, and 10 healthy volunteers, as part of a clinical trial in collaboration with the Hyper Acute Stroke Unit (HASU) at University College London Hospital (UCLH). Data were collected at 17 frequencies between 5 Hz and 2 kHz, with 31 current injections, yielding 930 measurements at each frequency. This dataset is the most comprehensive of its kind and enables combined analysis of MFEIT, Electroencephalography (EEG) and Computed Tomography (CT) or Magnetic Resonance Imaging (MRI) data in stroke patients, which can form the basis of future research into stroke classification.