Sampling lies at the heart of signal processing. The celebrated Shan-non - Nyquist theorem states that in order to reconstruct a continuous or discrete time signal from uniform samples one must sample at a rate twice the highest frequency present in the signal. Numerous signals and images of interest, however, are not even approximately bandlimited. While much progress has happened in recent years, reconstruction from sub-Nyquist samples still hinges on the use of random / incoherent (aggregate) sampling patterns, instead of uniform or regular sampling, which is far more simple, practical, and natural in many applications. In this work, we study regular sampling and reconstruction of three- or higher-dimensional signals (tensors). We prove that exact tensor reconstruction from regular samples is feasible under mild conditions on the rank of the tensor. Furthermore we cast the functional magnetic resonance imaging (fMRI) acceleration task as a regular tensor sampling problem and provide an algorithmic framework that effectively handles the reconstruction task. Experiments based on synthetic data and real fMRI data showcase the effectiveness of our approach.