This paper presents applications of boundary element methods to electrical impedance tomography. An algorithm for imaging the interior of a domain that consists of regions of constant conductivity is developed, that makes use of a simpler parametrization of the shapes of the regions to achieve efficiency. Numerical results from tests of this algorithm on synthetic data are presented, and show that the method is quite promising.