In this paper an expansion method, based on Legendre or any orthogonal polynomials, is developed to find numerical solutions of two-dimensional linear Fredholm integral equations. We estimate the error of the method, and present some numerical examples to demonstrate the accuracy of the method.