A direct D-bar reconstruction algorithm for recovering a complex conductivity in 2D

A direct reconstruction algorithm for complex conductivities in W-2,W-infinity(Omega), where Omega is a bounded, simply connected Lipschitz domain in R-2, is presented. The framework is based on the uniqueness proof by Francini (2000 Inverse Problems 6 107-19), but equations relating the Dirichlet-to-Neumann to the scattering transform and the exponentially growing solutions are not present in that work, and are derived here. The algorithm constitutes the first D-bar method for the reconstruction of conductivities and permittivities in two dimensions. Reconstructions of numerically simulated chest phantoms with discontinuities at the organ boundaries are included.

Effect of cooling and re-warming on cerebral and whole body electrical impedance

Cerebral electrical impedance is useful for the detection of cerebral edema following hypoxia in newborn infants. Thus it may be useful for determining neurological outcome or monitoring treatment. Hypothermia is a promising new therapy currently undergoing trials, but will alter impedance measurements. This study aimed to define the relationship between temperature and both cerebral and whole body electrical impedance, and to derive correction factors for adjustment of impedance measurements during hypothermia. In eight anaesthetized 1-2 day old piglets rectal, tympanic and scalp temperatures were monitored continuously. Following baseline readings at a rectal temperature of 39degreesC, piglets were cooled to 32degreesC. Four piglets were re-warmed. Cerebral and whole body impedance were measured at each 0.5degreesC as rectal temperature decreased. There was a strong linear relationship between both cerebral and whole body impedance and each of the temperatures measured. There was no difference in the relationship between impedance and rectal, tympanic or scalp temperatures. The relationship for impedance and rectal temperature was the same during cooling and re-warming. Using the correction factors derived it will be possible to accurately monitor cerebral and whole body fluid distribution during hypothermic treatment.

Ground penetrating radar (GPR) and electrical resistivity tomography (ERT) survey on the neotectonic activity of the Padul-Nigüelas Fault Zone (PNFZ), Southern Spain

The Padul-Nigüelas Fault Zone (PNFZ) is situated at the south-western mountain front of the Sierra Nevada (Spain) in an extensive regime and belongs to the internal zone of the Betic Cordilleras. The aim of this study is a collection of new evidence for neotectonic activity of the fault zone with classical geological field work and modern geophysical methods, such as ground penetrating radar (GPR). Among an apparently existing bed rock fault scarp with triangular facets, other evidences, such as deeply incised valleys and faults in the colluvial wedges, are present in the PNFZ. The preliminary results of our recent field work have shown that the synsedimentary faults within the colluvial sediments seem to propagate basinwards and the bed rock fault is only exhumed due to erosion for the studied segment (west of Marchena). We will use further GPR data and geomorphologic indices to gather further evidences of neotectonic activity of the PNFZ.

A Comparative Analysis Between Genetic Algorithms and Complex Nonlinear Least Squares on Electrical Impedance Characterization

This paper compares the performance of the complex nonlinear least squares algorithm implemented in the LEVM/LEVMW software with the performance of a genetic algorithm in the characterization of an electrical impedance of known topology. The effect of the number of measured frequency points and of measurement uncertainty on the estimation of circuit parameters is presented. The analysis is performed on the equivalent circuit impedance of a humidity sensor.

Persistence of ventilatory defect after resolution of pulmonary interstitial emphysema in a preterm baby

Pulmonary interstitial emphysema is a common complication of mechanical ventilation in preterm babies. We report a case of severe unilateral pulmonary interstitial emphysema in a premature newborn, treated with high-frequency oscillatory ventilation, lateral decubitus positioning and selective intubation. After complete radiological resolution of the pulmonary emphysema in the left lung, the patient was studied by electrical impedance tomography and a marked reduction of ventilation was identified in the left lung despite radiological resolution of the cysts. This finding indicates that functional abnormalities may persist for longer periods after radiologic resolution of such lesions.

Pulmonary lesion induced by low and high positive end-expiratory pressure levels during protective ventilation in experimental acute lung injury

Objective: To investigate the effects of low and high levels of positive end-expiratory pressure (PEEP), without recruitment maneuvers, during lung protective ventilation in an experimental model of acute lung injury (ALI). Design: Prospective, randomized, and controlled experimental study. Setting: University research laboratory. Subjects: Wistar rats were randomly assigned to control (C) [saline (0.1 ml), intraperitoneally] and ALI [paraquat (15 mg/kg), intra peritoneally] groups. Measurements and Main Results: After 24 hours, each group was further randomized into four groups (six rats each) at different PEEP levels = 1.5, 3, 4.5, or 6 cm H(2)O and ventilated with a constant tidal volume (6 mL/kg) and open thorax. Lung mechanics [static elastance (Est, L) and viscoelastic pressure (Delta P2, L)] and arterial blood gases were measured before (Pre) and at the end of 1-hour mechanical ventilation (Post). Pulmonary histology (light and electron microscopy) and type III procollagen (PCIII) messenger RNA (mRNA) expression were measured after 1 hour of mechanical ventilation. In ALI group, low and high PEEP levels induced a greater percentage of increase in Est, L (44% and 50%) and Delta P2, L (56% and 36%) in Post values related to Pre. Low PEEP yielded alveolar collapse whereas high PEEP caused overdistension and atelectasis, with both levels worsening oxygenation and increasing PCIII mRNA expression. Conclusions: In the present nonrecruited ALI model, protective mechanical ventilation with lower and higher PEEP levels than required for better oxygenation increased Est, L and Delta P2, L, the amount of atelectasis, and PCIII mRNA expression. PEEP selection titrated for a minimum elastance and maximum oxygenation may prevent lung injury while deviation from these settings may be harmful. (Crit Care Med 2009; 37:1011-1017)

Aortic blood pressure measured via EIT: investigation of different measurement settings.

Electrical impedance tomography (EIT) allows the measurement of intra-thoracic impedance changes related to cardiovascular activity. As a safe and low-cost imaging modality, EIT is an appealing candidate for non-invasive and continuous haemodynamic monitoring. EIT has recently been shown to allow the assessment of aortic blood pressure via the estimation of the aortic pulse arrival time (PAT). However, finding the aortic signal within EIT image sequences is a challenging task: the signal has a small amplitude and is difficult to locate due to the small size of the aorta and the inherent low spatial resolution of EIT. In order to most reliably detect the aortic signal, our objective was to understand the effect of EIT measurement settings (electrode belt placement, reconstruction algorithm). This paper investigates the influence of three transversal belt placements and two commonly-used difference reconstruction algorithms (Gauss-Newton and GREIT) on the measurement of aortic signals in view of aortic blood pressure estimation via EIT. A magnetic resonance imaging based three-dimensional finite element model of the haemodynamic bio-impedance properties of the human thorax was created. Two simulation experiments were performed with the aim to (1) evaluate the timing error in aortic PAT estimation and (2) quantify the strength of the aortic signal in each pixel of the EIT image sequences. Both experiments reveal better performance for images reconstructed with Gauss-Newton (with a noise figure of 0.5 or above) and a belt placement at the height of the heart or higher. According to the noise-free scenarios simulated, the uncertainty in the analysis of the aortic EIT signal is expected to induce blood pressure errors of at least ± 1.4 mmHg.

Whither lung EIT: Where are we, where do we want to go and what do we need to get there?

Breathing moves volumes of electrically insulating air into and out of the lungs, producing conductivity changes which can be seen by electrical impedance tomography (EIT). It has thus been apparent, since the early days of EIT research, that imaging of ventilation could become a key clinical application of EIT. In this paper, we review the current state and future prospects for lung EIT, by a synthesis of the presentations of the authors at the 'special lung sessions' of the annual biomedical EIT conferences in 2009-2011. We argue that lung EIT research has arrived at an important transition. It is now clear that valid and reproducible physiological information is available from EIT lung images. We must now ask the question: How can these data be used to help improve patient outcomes? To answer this question, we develop a classification of possible clinical scenarios in which EIT could play an important role, and we identify clinical and experimental research programmes and engineering developments required to turn EIT into a clinically useful tool for lung monitoring.

Parallel modeling of the electric field distribution in the brain

The term "Brain Imaging" identi�es a set of techniques to analyze the structure and/or functional behavior of the brain in normal and/or pathological situations. These techniques are largely used in the study of brain activity. In addition to clinical usage, analysis of brain activity is gaining popularity in others recent �fields, i.e. Brain Computer Interfaces (BCI) and the study of cognitive processes. In this context, usage of classical solutions (e.g. f MRI, PET-CT) could be unfeasible, due to their low temporal resolution, high cost and limited portability. For these reasons alternative low cost techniques are object of research, typically based on simple recording hardware and on intensive data elaboration process. Typical examples are ElectroEncephaloGraphy (EEG) and Electrical Impedance Tomography (EIT), where electric potential at the patient's scalp is recorded by high impedance electrodes. In EEG potentials are directly generated from neuronal activity, while in EIT by the injection of small currents at the scalp. To retrieve meaningful insights on brain activity from measurements, EIT and EEG relies on detailed knowledge of the underlying electrical properties of the body. This is obtained from numerical models of the electric �field distribution therein. The inhomogeneous and anisotropic electric properties of human tissues make accurate modeling and simulation very challenging, leading to a tradeo�ff between physical accuracy and technical feasibility, which currently severely limits the capabilities of these techniques. Moreover elaboration of data recorded requires usage of regularization techniques computationally intensive, which influences the application with heavy temporal constraints (such as BCI). This work focuses on the parallel implementation of a work-flow for EEG and EIT data processing. The resulting software is accelerated using multi-core GPUs, in order to provide solution in reasonable times and address requirements of real-time BCI systems, without over-simplifying the complexity and accuracy of the head models.

Heterogeneous multicore systems for signal processing

This thesis explores the capabilities of heterogeneous multi-core systems, based on multiple Graphics Processing Units (GPUs) in a standard desktop framework. Multi-GPU accelerated desk side computers are an appealing alternative to other high performance computing (HPC) systems: being composed of commodity hardware components fabricated in large quantities, their price-performance ratio is unparalleled in the world of high performance computing. Essentially bringing “supercomputing to the masses”, this opens up new possibilities for application fields where investing in HPC resources had been considered unfeasible before. One of these is the field of bioelectrical imaging, a class of medical imaging technologies that occupy a low-cost niche next to million-dollar systems like functional Magnetic Resonance Imaging (fMRI). In the scope of this work, several computational challenges encountered in bioelectrical imaging are tackled with this new kind of computing resource, striving to help these methods approach their true potential. Specifically, the following main contributions were made: Firstly, a novel dual-GPU implementation of parallel triangular matrix inversion (TMI) is presented, addressing an crucial kernel in computation of multi-mesh head models of encephalographic (EEG) source localization. This includes not only a highly efficient implementation of the routine itself achieving excellent speedups versus an optimized CPU implementation, but also a novel GPU-friendly compressed storage scheme for triangular matrices. Secondly, a scalable multi-GPU solver for non-hermitian linear systems was implemented. It is integrated into a simulation environment for electrical impedance tomography (EIT) that requires frequent solution of complex systems with millions of unknowns, a task that this solution can perform within seconds. In terms of computational throughput, it outperforms not only an highly optimized multi-CPU reference, but related GPU-based work as well. Finally, a GPU-accelerated graphical EEG real-time source localization software was implemented. Thanks to acceleration, it can meet real-time requirements in unpreceeded anatomical detail running more complex localization algorithms. Additionally, a novel implementation to extract anatomical priors from static Magnetic Resonance (MR) scansions has been included.

Analog Signal Acquisition and Conditioning for Near-Field Capacitive Communication and Active Combined EEG-EIT Monitoring

The work of the present thesis is focused on the implementation of microelectronic voltage sensing devices, with the purpose of transmitting and extracting analog information between devices of different nature at short distances or upon contact. Initally, chip-to-chip communication has been studied, and circuitry for 3D capacitive coupling has been implemented. Such circuits allow the communication between dies fabricated in different technologies. Due to their novelty, they are not standardized and currently not supported by standard CAD tools. In order to overcome such burden, a novel approach for the characterization of such communicating links has been proposed. This results in shorter design times and increased accuracy. Communication between an integrated circuit (IC) and a probe card has been extensively studied as well. Today wafer probing is a costly test procedure with many drawbacks, which could be overcome by a different communication approach such as capacitive coupling. For this reason wireless wafer probing has been investigated as an alternative approach to standard on-contact wafer probing. Interfaces between integrated circuits and biological systems have also been investigated. Active electrodes for simultaneous electroencephalography (EEG) and electrical impedance tomography (EIT) have been implemented for the first time in a 0.35 um process. Number of wires has been minimized by sharing the analog outputs and supply on a single wire, thus implementing electrodes that require only 4 wires for their operation. Minimization of wires reduces the cable weight and thus limits the patient's discomfort. The physical channel for communication between an IC and a biological medium is represented by the electrode itself. As this is a very crucial point for biopotential acquisitions, large efforts have been carried in order to investigate the different electrode technologies and geometries and an electromagnetic model is presented in order to characterize the properties of the electrode to skin interface.

Influence of fluid and volume state on PaO(2) oscillations in mechanically ventilated pigs

ABSTRACT Varying pulmonary shunt fractions during the respiratory cycle cause oxygen oscillations during mechanical ventilation. In artificially damaged lungs, cyclical recruitment of atelectasis is responsible for varying shunt according to published evidence. We introduce a complimentary hypothesis that cyclically varying shunt in healthy lungs is caused by cyclical redistribution of pulmonary perfusion. Administration of crystalloid or colloid infusions would decrease oxygen oscillations if our hypothesis was right. Therefore, n = 14 mechanically ventilated healthy pigs were investigated in 2 groups: crystalloid (fluid) versus no-fluid administration. Additional volume interventions (colloid infusion, blood withdrawal) were carried out in each pig. Intra-aortal PaO(2) oscillations were recorded using fluorescence quenching technique. Phase shift of oxygen oscillations during altered inspiratory to expiratory (I:E) ventilation ratio and electrical impedance tomography (EIT) served as control methods to exclude that recruitment of atelectasis is responsible for oxygen oscillations. In hypovolemia relevant oxygen oscillations could be recorded. Fluid and volume state changed PaO(2) oscillations according to our hypothesis. Fluid administration led to a mean decline of 105.3 mmHg of the PaO(2) oscillations amplitude (P < 0.001). The difference of the amplitudes between colloid administration and blood withdrawal was 62.4 mmHg in pigs not having received fluids (P = 0.0059). Fluid and volume state also changed the oscillation phase during altered I:E ratio. EIT excluded changes of regional ventilation (i.e., recruitment of atelectasis) to be responsible for these oscillations. In healthy pigs, cyclical redistribution of pulmonary perfusion can explain the size of respiratory-dependent PaO(2) oscillations.