Gustavo Sato dos Santos

Multifrequency Electrical Impedance Tomography Using Spectral Constraints

Multifrequency electrical impedance tomography (MFEIT) exploits the dependence of tissue impedance on frequency to recover an image of conductivity. MFEIT could provide emergency diagnosis of pathologies such as acute stroke, brain injury and breast cancer. We present a method for performing MFEIT using spectral constraints. Boundary voltage data is employed directly to reconstruct the volume fraction distribution of component tissues using a nonlinear method. Given that the reconstructed parameter is frequency independent, this approach allows for the simultaneous use of all multifrequency data, thus reducing the degrees of freedom of the reconstruction problem. Furthermore, this method allows for the use of frequency difference data in a nonlinear reconstruction algorithm. Results from empirical phantom measurements suggest that our fraction reconstruction method points to a new direction for the development of multifrequency EIT algorithms in the case that the spectral constraints are known, and may provide a unifying framework for static EIT imaging.

A method for reconstructing tomographic images of evoked neural activity with electrical impedance tomography using intracranial planar arrays

A method is presented for reconstructing images of fast neural evoked activity in rat cerebral cortex recorded with electrical impedance tomography (EIT) and a 6 × 5 mm(2) epicortical planar 30 electrode array. A finite element model of the rat brain and inverse solution with Tikhonov regularization were optimized in order to improve spatial resolution and accuracy. The optimized FEM mesh had 7 M tetrahedral elements, with finer resolution (0.05 mm) near the electrodes. A novel noise-based image processing technique based on t-test significance improved depth localization accuracy from 0.5 to 0.1 mm. With the improvements, a simulated perturbation 0.5 mm in diameter could be localized in a region 4 × 5 mm(2) under the centre of the array to a depth of 1.4 mm, thus covering all six layers of the cerebral cortex with an accuracy of <0.1 mm. Simulated deep brain hippocampal or thalamic activity could be localized with an accuracy of 0.5 mm with a 256 electrode array covering the brain. Parallel studies have achieved a temporal resolution of 2 ms for imaging fast neural activity by EIT during evoked activity; this encourages the view that fast neural EIT can now resolve the propagation of depolarization-related fast impedance changes in cerebral cortex and deeper in the brain with a resolution equal or greater to the dimension of a cortical column.

Imaging fast electrical activity in the brain with electrical impedance tomography

Imaging of neuronal depolarization in the brain is a major goal in neuroscience, but no technique currently exists that could image neural activity over milliseconds throughout the whole brain. Electrical impedance tomography (EIT) is an emerging medical imaging technique which can produce tomographic images of impedance changes with non-invasive surface electrodes. We report EIT imaging of impedance changes in rat somatosensory cerebral cortex with a resolution of 2 ms and < 200 μm during evoked potentials using epicortical arrays with 30 electrodes. Images were validated with local field potential recordings and current source-sink density analysis. Our results demonstrate that EIT can image neural activity in a volume 7 × 5 × 2 mm in somatosensory cerebral cortex with reduced invasiveness, greater resolution and imaging volume than other methods. Modeling indicates similar resolutions are feasible throughout the entire brain so this technique, uniquely, has the potential to image functional connectivity of cortical and subcortical structures.

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.

Interventional Photoacoustic Imaging of the Human Placenta with Ultrasonic Tracking for Minimally Invasive Fetal Surgeries

Image guidance plays a central role in minimally invasive fetal surgery such as photocoagulation of inter-twin placental anastomosing vessels to treat twin-to-twin transfusion syndrome (TTTS). Fetoscopic guidance provides insufficient sensitivity for imaging the vasculature that lies beneath the fetal placental surface due to strong light scattering in biological tissues. Incomplete photocoagulation of anastamoses is associated with postoperative complications and higher perinatal mortality. In this study, we investigated the use of multi-spectral photoacoustic (PA) imaging for better visualization of the placental vasculature. Excitation light was delivered with an optical fiber with dimensions that are compatible with the working channel of a fetoscope. Imaging was performed on an ex vivo normal term human placenta collected at Caesarean section birth. The photoacoustically-generated ultrasound signals were received by an external clinical linear array ultrasound imaging probe. A vein under illumination on the fetal placenta surface was visualized with PA imaging, and good correspondence was obtained between the measured PA spectrum and the optical absorption spectrum of deoxygenated blood. The delivery fiber had an attached fiber optic ultrasound sensor positioned directly adjacent to it, so that its spatial position could be tracked by receiving transmissions from the ultrasound imaging probe. This study provides strong indications that PA imaging in combination with ultrasonic tracking could be useful for detecting the human placental vasculature during minimally invasive fetal surgery.

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.

A Reconstruction-Classification Method for Multifrequency Electrical Impedance Tomography

Multifrequency Electrical Impedance Tomography is an imaging technique which distinguishes biological tissues by their unique conductivity spectrum. Recent results suggest that the use of spectral constraints can significantly improve image quality. We present a combined reconstruction-classification method for estimating the spectra of individual tissues, whilst simultaneously reconstructing the conductivity. The advantage of this method is that a priori knowledge of the spectra is not required to be exact in that the constraints are updated at each step of the reconstruction. In this paper, we investigate the robustness of the proposed method to errors in the initial guess of the tissue spectra, and look at the effect of introducing spatial smoothing. We formalize and validate a frequency-difference variant of reconstruction-classification, and compare the use of absolute and frequency-difference data in the case of a phantom experiment.

Investigation of potential artefactual changes in measurements of impedance changes during evoked activity: implications to electrical impedance tomography of brain function

Abstract Electrical impedance tomography (EIT) could provide images of fast neural activity in the adult human brain with a resolution of 1 ms and 1 mm by imaging impedance changes which occur as ion channels open during neuronal depolarization. The largest changes occur at dc and decrease rapidly over 100 Hz. Evoked potentials occur in this bandwidth and may cause artefactual apparent impedance changes if altered by the impedance measuring current. These were characterized during the compound action potential in the walking leg nerves of Cancer pagurus, placed on Ag/AgCl hook electrodes, to identify how to avoid artefactual changes during brain EIT. Artefact-free impedance changes (δZ) decreased with frequency from −0.045 ± 0.01% at 225 Hz to −0.02 ± 0.01% at 1025 Hz (mean ± 1 SD, n = 24 in 12 nerves) which matched changes predicted by a finite element model. Artefactual δZ reached c.300% and 50% of the genuine membrane impedance change at 225 Hz and 600 Hz respectively but decreased with frequency of the applied current and was negligible above 1 kHz. The proportional amplitude (δZ (%)) of the artefact did not vary significantly with the amplitude of injected current of 5–20 µA pp. but decreased significantly from −0.09 ± 0.024 to −0.03 ± 0.023% with phase of 0 to 45°. For fast neural EIT of evoked activity in the brain, artefacts may arise with applied current of >10 µA. Independence of δZ with respect to phase but not the amplitude of applied current controls for them; they can be minimized by randomizing the phase of the applied measuring current and excluded by recording at >1 kHz.

A Registration Approach to Endoscopic Laser Speckle Contrast Imaging for Intrauterine Visualisation of Placental Vessels

Intrauterine interventions such as twin-to-twin transfusion syndrome procedure require accurate mapping of the fetal placental vasculature to ensure complete photocoagulation of vascular anastomoses. However, surgeons are currently limited to fetoscopy and external ultrasound imaging, which are unable to accurately identify all vessels especially those that are narrow and at the periphery. Laser speckle contrast imaging LSCI is an optical method for imaging blood flow that is emerging as an intraoperative tool for neurosurgery. Here we explore the application of LSCI to minimally invasive fetal surgery, with an endoscopic LSCI system based on a 2.7-mm-diameter fetoscope. We establish using an optical phantom that it can image flow in 1-mm-diameter vessels as far as 4 mm below the surface. We demonstrate that a spatiotemporal algorithm produces the clearest images of vessels within 200 ms, and that speckle contrast images can be accurately registered using groupwise registration to correct for significant motion of target or probe. When tested on a perfused term ex vivo human placenta, our endoscopic LSCI system revealed small capillaries not evident in the fetoscopic images.

Comparison of Different Quadratic Regularization for Electrical Impedance Tomography

Electrical Impedance Tomography (EIT) is an ill-conditioned problem in which regularization is used to calculate a stable and accurate solution by incorporating some form of prior knowledge into the solution. A number of regularization terms have been proposed for EIT since they influenced the performance of EIT a lot. However, it is difficult to evaluate the performance of different regularization terms due to the influences of optimization and regularization parameters methods. This paper compares four widely used quadratic regularization methods, including Markov Random Field (MRF), Noser, Laplace and Tikhonov, with L curve and Newton- Krylov subspace method for both simulation and phantom experiments. The results demonstrate that the MRF not only generates clear images in all trials but also have moderate imaging speed, fitting the general applications, and Tikhonov approach is fastest but with strong artifacts, fitting the application with high time resolution.