A DeepONet for Inverting the Neumann-to-Dirichlet Operator in Electrical Impedance Tomography: An Approximation Theoretic Perspective and Numerical Results

In this work, we consider the non-invasive medical imaging modality of Electrical Impedance Tomography (EIT), where the goal is to recover the conductivity in a medium from boundary current-to-voltage measurements, i.e., the Neumann-to-Dirichlet (N–t–D) operator. We formulate this inverse problem as an operator-learning task, where the aim is to approximate the implicitly defined map from N–t–D operators to admissible conductivities. To this end, we employ a Deep Operator Network (DeepONet) architecture, thereby extending operator learning beyond the classical function-to-function setting to the more challenging operator-to-function regime. We establish a universal approximation theorem that guarantees that such operator-to-function maps can be approximated arbitrarily well by DeepONets. Furthermore, we provide a computational implementation of our approach and compare it against the iteratively regularized Gauss–Newton (IRGN) method. Our results show that the proposed framework yields accurate and robust reconstructions, outperforms the baseline, and demonstrates strong generalization. To our knowledge, this is the first work that combines rigorous approximation-theoretic guarantees with DeepONet-based inversion for EIT, thereby opening a principled and interpretable pathway for use of DeepONets in such inverse problems.