Numerical analysis of the factorization method for EIT with a piecewise constant uncertain background

We extend the factorization method for electrical impedance tomography to the case of background featuring uncertainty. We first describe the algorithm for the known but inhomogeneous backgrounds and indicate expected accuracy from the inversion method through some numerical tests. Then we develop three methodologies to apply the factorization method to the more difficult case of a piecewise constant but uncertain background. The first one is based on a recovery of the background through an optimization scheme and is well adapted to relatively low-dimensional random variables describing the background. The second one is based on a weighted combination of the indicator functions provided by the factorization method for different realizations of the random variables describing the uncertain background. We show through numerical experiments that this procedure is well suited to the case where many realizations of the measurement operators are available. The third strategy is a variant of the previous one when measurements for the inclusion-free background are available. In this case, a single pair of measurements is sufficient to achieve comparable accuracy to the deterministic case.