The solution of the radiation transfer equation for the Earth's atmosphere needs to account for the re ectivity of the ground. When using the spherical harmonics method, the solution for this term involves an integral with a particular measure that presents numerical challenges. We are interested in computing a high order Gauss quadrature rule for this measure. We show that the two classical al...