An optimal sampling pattern for variable density sampling of a continuous signal uses a statistical knowledge of the signal to determine an autocorrelation matrix from which a basis set is identified. Sampling is performed at locations determined from an eigenvector matrix, and the sampled output provides coefficients for the basis set. The reconstructed signal output is a summation of the multiplication of the coefficients and the basis set.