We investigate the application of multifrequency electrical impedance tomography (MFEIT) to imaging the brain in stroke patients. The use of MFEIT could enable early diagnosis and thrombolysis of ischaemic stroke, and therefore improve the outcome of treatment. Recent advances in the imaging methodology suggest that the use of spectral constraints could allow for the reconstruction of a one-shot image. We performed a simulation study to investigate the feasibility of imaging stroke in a head model with realistic conductivities. We introduced increasing levels of modelling errors to test the robustness of the method to the most common sources of artefact. We considered the case of errors in the electrode placement, spectral constraints, and contact impedance. The results indicate that errors in the position and shape of the electrodes can affect image quality, although our imaging method was successful in identifying tissues with sufficiently distinct spectra.

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.

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.

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.

Photoacoustic tomography can, in principle, provide quantitatively accurate, high-resolution, images of chromophore distributions in 3D in vivo. However, achieving this goal requires not only dealing with the optical fluence-related spatial and spectral distortion but also having access to high quality, calibrated, measurements and using image reconstruction algorithms free from inaccurate assumptions. Furthermore, accurate knowledge of experimental parameters, such as the positions of the ultrasound detectors and the illumination pattern, is necessary for the reconstruction step. A meticulous and rigorous experimental phantom study was conducted to show that highly-resolved 3D estimation of chromophore distributions can be achieved: a crucial step towards in vivo implementation. The phantom consisted of four 580 µm diameter tubes with different ratios of copper sulphate and nickel sulphate as hemoglobin analogues, submersed in a background medium of intralipid and india ink. The optical absorption, scattering, photostability, and Grüneisen parameter were characterised for all components independently. A V-shaped imaging scanner enabled 3D imaging with the high resolution, high sensitivity, and wide bandwidth characteristic of Fabry-Pérot ultrasound sensors, but without the limited-view disadvantage of single-plane scanners. The optical beam profile and position were determined experimentally. Nine wavelengths between 750 and 1110 nm were used. The images of the chromophore concentrations were obtained using a model-based, two-step, procedure, that did not require image segmentation. First, the acoustic reconstruction was solved with an iterative time-reversal algorithm to obtain images of the initial acoustic pressure at each of the nine wavelengths for an 18×17×13 mm 3 volume with 50µm voxels. Then, 3D high resolution estimates of the chromophore concentrations were obtained by using a diffusion model of light transport in an iterative nonlinear optimisation scheme. Among the lessons to be drawn from this study, one is fundamental: in order to obtain accurate estimates of chromophores (or their ratios) it is not only necessary to model the light fluence accurately, but it is just as crucial to obtain accurate estimates of the initial acoustic pressure distributions, and to account for variations in the thermoelastic efficiency (Grüneisen parameter).

We propose a combined reconstruction-classification method for simultaneously recovering absorption and scattering in turbid media from images of absorbed optical energy. This method exploits knowledge that optical parameters are determined by a limited number of classes to iteratively improve their estimate. Numerical experiments show that the proposed approach allows for accurate recovery of absorption and scattering in 2 and 3 dimensions, and delivers superior image quality with respect to traditional reconstruction-only approaches.

We propose a novel multispectral reconstruction-classification method for simultaneously recovering absorption and scattering coefficients from images of absorbed optical energy. In contrast with pre-existing chromophore reconstruction methods, this approach does not require prior knowledge of the characteristic spectra of the absorbers, which is not always available. Numerical experiments performed on anatomically realistic 3D phantoms show that this approach allows for improved recovery of both the optical absorption and scattering with respect to reconstruction-only methods, and accurate classification of chromophores of clinical interest.

No abstract