Consider a stationary spatio‐temporal random process and let be a sample from the process. Our object here is to predict, given the sample, for all t at the location so. To obtain the predictors, we define a sequence of discrete Fourier transforms using the observed time series. We consider these discrete Fourier transforms as a sample from the complex valued random variable . Assuming that the discrete Fourier transforms satisfy a complex stochastic partial differential equation of the Laplacian type with a scaling function that is a polynomial in the temporal spectral frequency ω, we obtain, in a closed form, expressions for the second‐order spatio‐temporal spectrum and the covariance function. The spectral density function obtained corresponds to a non‐separable random process. The optimal predictor of the discrete Fourier transform is in terms of the covariance functions. The estimation of the parameters of the spatio‐temporal covariance function is considered and is based on the recently introduced frequency variogram method. The methods given here can be extended to situations where the observations are corrupted by independent white noise. The methods are illustrated with a real data set.