Simulation-based feasibility study of monitoring of intracerebral hemorrhages and detection of secondary hemorrhages using electrical impedance tomography.

Purpose: We present a simulation-based feasibility study of electrical impedance tomography (EIT) for continuous bedside monitoring of intracerebral hemorrhages (ICH) and detection of secondary hemorrhages. Approach: We simulated EIT measurements for six different hemorrhage sizes at two different hemorrhage locations using an anatomically detailed computational head model. Using this dataset, we test the ICH monitoring and detection performance of our tailor-made, patient-specific stroke-monitoring algorithm that utilizes a novel combination of nonlinear region-of-interest difference imaging, parallel level sets regularization and a prior-conditioned least squares algorithm. We compare the results of our algorithm to the results of two reference algorithms, a total variation regularized absolute imaging algorithm and a linear difference imaging algorithm. Results: The tailor-made stroke-monitoring algorithm is capable of indicating smaller changes in the simulated hemorrhages than either of the reference algorithms, indicating better monitoring and detection performance. Conclusions: Our simulation results from the anatomically detailed head model indicate that EIT equipped with a patient-specific stroke-monitoring algorithm is a promising technology for the unmet clinical need of having a technology for continuous bedside monitoring of brain status of acute stroke patients.

Effects of atlas-based anatomy on modelled light transport in the neonatal head.

Objective.Diffuse optical tomography (DOT) provides a relatively convenient method for imaging haemodynamic changes related to neuronal activity on the cerebral cortex. Due to practical challenges in obtaining anatomical images of neonates, an anatomical framework is often created from an age-appropriate atlas model, which is individualized to the subject based on measurements of the head geometry. This work studies the approximation error arising from using an atlas instead of the neonate's own anatomical model.Approach.We consider numerical simulations of frequency-domain (FD) DOT using two approaches, Monte Carlo simulations and diffusion approximation via finite element method, and observe the variation in (1) the logarithm of amplitude and phase shift measurements, and (2) the corresponding inner head sensitivities (Jacobians), due to varying segmented anatomy. Varying segmentations are sampled by registering 165 atlas models from a neonatal database to the head geometry of one individual selected as the reference model. Prior to the registration, we refine the segmentation of the cerebrospinal fluid (CSF) by separating the CSF into two physiologically plausible layers.Main results.In absolute measurements, a considerable change in the grey matter or extracerebral tissue absorption coefficient was found detectable over the anatomical variation. In difference measurements, a small local 10%-increase in brain absorption was clearly detectable in the simulated measurements over the approximation error in the Jacobians, despite the wide range of brain maturation among the registered models.Significance.Individual-level atlas models could potentially be selected within several weeks in gestational age in DOT difference imaging, if an exactly age-appropriate atlas is not available. The approximation error method could potentially be implemented to improve the accuracy of atlas-based imaging. The presented CSF segmentation algorithm could be useful also in other model-based imaging modalities. The computation of FD Jacobians is now available in the widely-used Monte Carlo eXtreme software.

Ultrasound-modulated optical tomography: recovery of amplitude of vibration in the insonified region from boundary measurement of light correlation.

We address a certain inverse problem in ultrasound-modulated optical tomography: the recovery of the amplitude of vibration of scatterers [p(r)] in the ultrasound focal volume in a diffusive object from boundary measurement of the modulation depth (M) of the amplitude autocorrelation of light [φ(r,τ)] traversing through it. Since M is dependent on the stiffness of the material, this is the precursor to elasticity imaging. The propagation of φ(r,τ) is described by a diffusion equation from which we have derived a nonlinear perturbation equation connecting p(r) and refractive index modulation [Δn(r)] in the region of interest to M measured on the boundary. The nonlinear perturbation equation and its approximate linear counterpart are solved for the recovery of p(r). The numerical results reveal regions of different stiffness, proving that the present method recovers p(r) with reasonable quantitative accuracy and spatial resolution.