scispace - formally typeset
Search or ask a question
Author

Patricia Brunner

Bio: Patricia Brunner is an academic researcher from Graz University of Technology. The author has contributed to research in topics: Magnetic induction tomography & Inverse problem. The author has an hindex of 8, co-authored 16 publications receiving 309 citations.

Papers
More filters
Journal ArticleDOI
TL;DR: An inverse solver based on the Gauss-Newton-one-step method for differential imaging is developed, and four different regularization schemes are implemented and tested, demonstrating the feasibility of difference imaging with MIT.
Abstract: Magnetic induction tomography (MIT) of biological tissue is used to reconstruct the changes in the complex conductivity distribution inside an object under investigation. The measurement principle is based on determining the perturbation ?B of a primary alternating magnetic field B0, which is coupled from an array of excitation coils to the object under investigation. The corresponding voltages ?V and V0 induced in a receiver coil carry the information about the passive electrical properties (i.e. conductivity, permittivity and permeability). The reconstruction of the conductivity distribution requires the solution of a 3D inverse eddy current problem. As in EIT the inverse problem is ill-posed and on this account some regularization scheme has to be applied. We developed an inverse solver based on the Gauss?Newton-one-step method for differential imaging, and we implemented and tested four different regularization schemes: the first and second approaches employ a classical smoothness criterion using the unit matrix and a differential matrix of first order as the regularization matrix. The third method is based on variance uniformization, and the fourth method is based on the truncated singular value decomposition. Reconstructions were carried out with synthetic measurement data generated with a spherical perturbation at different locations within a conducting cylinder. Data were generated on a different mesh and 1% random noise was added. The model contained 16 excitation coils and 32 receiver coils which could be combined pairwise to give 16 planar gradiometers. With 32 receiver coils all regularization methods yield fairly good 3D-images of the modelled changes of the conductivity distribution, and prove the feasibility of difference imaging with MIT. The reconstructed perturbations appear at the right location, and their size is in the expected range. With 16 planar gradiometers an additional spurious feature appears mirrored with respect to the median plane with negative sign. This demonstrates that a symmetrical arrangement with one ring of planar gradiometers cannot distinguish between a positive conductivity change at the true location and a negative conductivity change at the mirrored location.

103 citations

Journal ArticleDOI
TL;DR: A 14-channel multifrequency magnetic induction tomography system (MF-MIT) for biomedical applications and a FEM model and an eddy current solver are used to evaluate the experimental results and to reconstruct the images.
Abstract: We developed a 14-channel multifrequency magnetic induction tomography system (MF-MIT) for biomedical applications. The excitation field is produced by a single coil and 14 planar gradiometers are used for signal detection. The object under measurement was rotated (16 steps per turn) to obtain a full data set for image reconstruction. We make measurements at frequencies from 50 kHz to 1 MHz using a single frequency excitation signal or a multifrequency signal containing several frequencies in this range. We used two acquisition boards giving a total of eight synchronous channels at a sample rate of 5 MS s(-1) per channel. The real and imaginary parts of DeltaB/B(0) were calculated using coherent demodulation at all injected frequencies. Calibration, averaging and drift cancellation techniques were used before image reconstruction. A plastic tank filled with saline (D = 19 cm) and with conductive and/or paramagnetic perturbations was measured for calibration and test purposes. We used a FEM model and an eddy current solver to evaluate the experimental results and to reconstruct the images. Measured equivalent input noise voltage for each channel was 2 nV Hz(-1/2). Using coherent demodulation, with an integration time of 20 ms, the measured STD for the magnitude was 7 nV(rms) (close to the theoretical value only taking into account the amplifier's thermal noise). For long acquisition times the drift in the signal produced a bigger effect than the input noise (typical STD was 10 nV with a maximum of 35 nV at one channel) but this effect was reduced using a drift cancellation technique based on averaging. We were able to image a 2 S m(-1) agar sphere (D = 4 cm) inside the tank filled with saline of 1 S m(-1).

60 citations

Journal ArticleDOI
TL;DR: The non-linear inverse problem of MIT can be approximated by a linear one (depending only on the frequency), similar to in EIT, but with a limitation that these spectra do not provide absolute conductivities but preserve the shape of the spectrum.
Abstract: Magnetic induction tomography (MIT) of biological tissue is used for the reconstruction of the complex conductivity distribution kappa inside the object under investigation. It is based on the perturbation of an alternating magnetic field caused by the object and can be used in all applications of electrical impedance tomography (EIT) such as functional lung monitoring and assessment of tissue fluids. In contrast to EIT, MIT does not require electrodes and magnetic fields can also penetrate non-conducting barriers such as the skull. As in EIT, the reconstruction of absolute conductivity values is very difficult because of the method's sensitivity to numerical errors and noise. To overcome this problem, image reconstruction in EIT is often done differentially. Analogously, this concept has been adopted for MIT. Two different methods for differential imaging are applicable. The first one is state-differential, for example when the conductivity change between inspiration and expiration in the lung regions is being detected. The second one is frequency-differential, which is of high interest in motionless organs like the brain, where a state-differential method cannot be applied. An equation for frequency-differential MIT was derived taking into consideration the frequency dependence of the sensitivity matrix. This formula is valid if we can assume that only small conductivity changes occur. In this way, the non-linear inverse problem of MIT can be approximated by a linear one (depending only on the frequency), similar to in EIT. Keeping this limitation in mind, the conductivity changes between one or more reference frequencies and several measurement frequencies were reconstructed, yielding normalized conductivity spectra. Due to the differential character of the method, these spectra do not provide absolute conductivities but preserve the shape of the spectrum. The validity of the method was tested with artificial data generated with a spherical perturbation within a conducting cylinder as well as for real measurement data. The measurement data were obtained from a potato immersed in saline. The resulting spectra were compared with reference measurements and the preservation of the shape of the spectra was analyzed.

40 citations

Journal ArticleDOI
TL;DR: This article addresses the following questions: to which extent do different abdominal compartments contribute to the impedance? 2) How does the hydration state of tissues affect the data? 3) Can hydration and fat content be assessed independently?
Abstract: In a previous publication, it was demonstrated that the abdominal subcutaneous fat layer thickness (SFL) is strongly correlated with the abdominal electrical impedance when measured with a transversal tetrapolar electrode arrangement. This article addresses the following questions: 1) To which extent do different abdominal compartments contribute to the impedance? 2) How does the hydration state of tissues affect the data? 3) Can hydration and fat content be assessed independently? For simulating the measured data a hierarchical electrical model was built. The abdomen was subdivided into three compartments (subcutaneous fat, muscle, mesentery). The true anatomical structure of the compartment boundaries was modeled using finite-element modeling (FEM). Each compartment is described by an electrical tissue model parameterized in physiological terms. Assuming the same percent change of the fat fraction in the mesentery and the SFL the model predicts a change of 1,24 /spl Omega//mm change of the SFL compared to 1,1 /spl Omega//mm measured. 42% of the change stem from the SFL, 56% from the mesentery and 2% from changes of fat within the muscle compartment. A 1% increase of the extracellular water in the muscle is not discernible from a 1% decrease of the SFL. The measured data reflect not only the SFL but also the visceral fat. The tetrapolar electrode arrangement allows the measurement of the abdominal fat content only if the hydration remains constant.

29 citations

Journal ArticleDOI
TL;DR: The goals of this study were the analysis of the sensitivity for the PEIX, an estimate of the optimal electrode configuration and the determination of the required frequencies, which shows that eight to ten electrodes and six frequencies are theoretically sufficient to obtain a coefficient of variation.
Abstract: Currently only ionizing or invasive methods are used in clinical applications for the monitoring of extracellular lung water. Alternatively a method called focused conductivity spectroscopy (FCS) is suggested, which aims at reconstructing a pulmonary edema index (PEIX) by measuring the electrical conductivity of the region of interest (ROI) at several frequencies. In contrast to electrical impedance tomography (EIT) a minimum number of strategically placed electrodes is used. The goals of this study were the analysis of the sensitivity for the PEIX, an estimate of the optimal electrode configuration and the determination of the required frequencies. In order to calculate the solution of the FCS forward problem a realistic 3D model of a human torso was developed containing both lungs, the heart, the liver and the thorax musculature. The bioelectrical properties for each compartment were described with appropriate tissue models which relate the conductivity spectra to physiological parameters. The PEIX was defined as the interstitial volume fraction of the alveolar septa. Furthermore the model includes 48 electrodes subdivided into three layers. The optimal electrode configuration was selected by minimizing the number of electrodes, among certain subsets of these electrodes. The analysis shows that eight to ten electrodes and six frequencies are theoretically sufficient to obtain a coefficient of variation.

26 citations


Cited by
More filters
Journal ArticleDOI
TL;DR: The Unified Form Language (UFL) as mentioned in this paper is a domain-specific language for representing weak formulations of partial differential equations with a view to numerical approximation, which has been used to effortlessly express finite element methods for complex systems of PDEs in near-mathematical notation.
Abstract: We present the Unified Form Language (UFL), which is a domain-specific language for representing weak formulations of partial differential equations with a view to numerical approximation. Features of UFL include support for variational forms and functionals, automatic differentiation of forms and expressions, arbitrary function space hierarchies for multifield problems, general differential operators and flexible tensor algebra. With these features, UFL has been used to effortlessly express finite element methods for complex systems of partial differential equations in near-mathematical notation, resulting in compact, intuitive and readable programs. We present in this work the language and its construction. An implementation of UFL is freely available as an open-source software library. The library generates abstract syntax tree representations of variational problems, which are used by other software libraries to generate concrete low-level implementations. Some application examples are presented and libraries that support UFL are highlighted.

338 citations

Posted Content
TL;DR: The Unified Form Language is presented, which is a domain-specific language for representing weak formulations of partial differential equations with a view to numerical approximation and generates abstract syntax tree representations of variational problems, which are used by other software libraries to generate concrete low-level implementations.
Abstract: We present the Unified Form Language (UFL), which is a domain-specific language for representing weak formulations of partial differential equations with a view to numerical approximation. Features of UFL include support for variational forms and functionals, automatic differentiation of forms and expressions, arbitrary function space hierarchies for multi-field problems, general differential operators and flexible tensor algebra. With these features, UFL has been used to effortlessly express finite element methods for complex systems of partial differential equations in near-mathematical notation, resulting in compact, intuitive and readable programs. We present in this work the language and its construction. An implementation of UFL is freely available as an open-source software library. The library generates abstract syntax tree representations of variational problems, which are used by other software libraries to generate concrete low-level implementations. Some application examples are presented and libraries that support UFL are highlighted.

218 citations

Journal ArticleDOI
Yuan Xu1, Bin He1
TL;DR: The promising pilot experimental results suggest the feasibility of the proposed MAT-MI approach, which promises to eliminate the shielding effects of other imaging modalities in which the current is applied directly with electrodes.
Abstract: We report our theoretical and experimental investigations on a new imaging modality, magnetoacoustic tomography with magnetic induction (MAT-MI) In MAT-MI, the sample is located in a static magnetic field and a time-varying (micros) magnetic field The time-varying magnetic field induces an eddy current in the sample Consequently, the sample will emit ultrasonic waves by the Lorentz force The ultrasonic signals are collected around the object to reconstruct images related to the electrical impedance distribution in the sample MAT-MI combines the good contrast of electrical impedance tomography with the good spatial resolution of sonography MAT-MI has two unique features due to the solenoid nature of the induced electrical field Firstly, MAT-MI could provide an explicit or simple quantitative reconstruction algorithm for the electrical impedance distribution Secondly, it promises to eliminate the shielding effects of other imaging modalities in which the current is applied directly with electrodes In the theoretical part, we provide formulae for both the forward and inverse problems of MAT-MI and estimate the signal amplitude in biological tissues In the experimental part, the experimental setup and methods are introduced and the signals and the image of a metal object by means of MAT-MI are presented The promising pilot experimental results suggest the feasibility of the proposed MAT-MI approach

218 citations

Journal ArticleDOI
TL;DR: Electrical impedance tomography is a promising functional tomography technology on the verge of its clinical application and potential applications in the future as well as limitations of EIT technology are described.
Abstract: Objective:This review outlines the basic principle, in addition to validated and upcoming clinical use of electrical impedance tomography (EIT). EIT generates functional tomograms of the thorax for detection of changes in regional lung aeration. These images allow an intraindividual comparison of ch

192 citations