The tomographic reconstruction of supersonic flows faces two challenges. Firstly, techniques used in the past, such as the direct Fourier method (DFM) (Gottlieb and Gustafsson in On the direct Fourier method for computer tomography, 1998; Morton inTomographic imaging of supersonic flows, 1995) or various backprojection (Kak and Slaney in Principles of computerized tomographic imaging, vol. 33 in Classics in Applied Mathematics, 2001) techniques, have only been able to reconstruct areas of the flow which are upstream of any opaque objects, such as a model. Secondly, shock waves create sharp discontinuities in flow properties, which can be difficult to reconstruct both in position and in magnitude with limited data. This paper will present a reconstruction method, matrix inversion using ray-tracing and least squares conjugate gradient (MI-RLS), which uses geometric ray-tracing and a sparse matrix iterative solver Paige and Saunders in ACM Trans. Math. Softw. 8(1):43– 71, 1982) to overcome both of these challenges. It will be shown, through testing with a phantom object described in tomographic literature, that the results compare favourably to those produced by theDFMtechnique. Finally, the method will be used to reconstruct three-dimensional density fields from interferometric shock layer images, with good resolution.