We describe a new approach to the index reconstruction of three-dimensional optical systems with rotational symmetry, which is based on sampling ray paths that lie in the sagittal plane. Since the observed rays are distorted by the optical system itself, they cannot be used directly for index reconstruction. We present an iterative procedure to compute the true ray paths and then to find the index distribution. The utility of the method is verified on the model problem.