We introduce an inversion algorithm for tomographic images of layered media. The algorithm is based on a multiscattering series expansion of the Green function that, unlike the Born series, converges unconditionally. Our inversion algorithm obtains images of the medium that improves iteratively as we use more and more terms in the multiscattering series. We present the derivation of the multiscattering series, formulate the inversion algorithm and demonstrate its performance through numerical experiments