Scalable Gaussian Processes with Latent Kronecker Structure