Polynomial approximation method for tomographic reconstruction of three-dimensional refractive index fields with limited data

We propose a polynomial approximation method (PAM) for reconstruction of three-dimensional refractive index fields by interferometric tomography using limited data. Based on the assumption that the fields to be reconstructed are usually smooth and can be decomposed into a finite order of (orthogonal) polynomials, a set of linear equations can be constructed using both the measured projection data and the Radon transform of the basis functions. By solving these equations, the least-squares solutions of expansion coefficients can be obtained and then substituted back to yield the desired fields. Numerical results have demonstrated that the proposed method is fast, robust to noise and can achieve satisfactory results for refractive index fields with limited projection views and large opaque objects.