This paper proposes an iterative technique to reconstruct the source term in transport equations, which account for scattering effects, from boundary measurements. In the two-dimensional setting, the full outgoing distribution in the phase space ( position and direction) needs to be measured. In three space dimensions, we show that measurements for angles that are orthogonal to a given direction are sufficient. In both cases, the derivation is based on a perturbation of the inversion of the two-dimensional attenuated Radon transform and requires that (the anisotropic part of) scattering be sufficiently small. We present an explicit iterative procedure, which converges to the source term we want to reconstruct. Applications of the inversion procedure include optical molecular imaging, an increasingly popular medical imaging modality.