Deep Random Features for Scalable Interpolation of Spatiotemporal Data