Shannon's sampling theory is based on the reconstruction of bandlimited signals which requires infinite number of uniform time samples. Indeed, one can only have finite number of samples for numerical implementation. In this paper, as a dual of the bandlimited reconstruction, a solution for time-limited signal reconstruction from nonuniform samples is proposed. The system model we present is based on the idea that time-limited signals which are also nearly bandlimited can be well approximated by a low-dimensional subspace. This can be done by using prolate spheroidal wave functions as the basis. The order of the projection on this basis is obtained by means of the time-frequency dimension of the signal, especially in the case of non-stationary signals. The reconstruction requires the estimation of the nonuniform sampling times by means of an annihilating filter. We obtain the reconstruction para ters by solving a linear system of equations and show that our finite-dimensional 1 is not ill-conditioned. The practical aspects of our method including the dimensionality reduction are demonstrated by processing synthetic as well as real signals. (C) 2009 Elsevier B.V. All rights reserved.