Gaussian beam algorithms for rough surface underground imaging