The Wigner-Kirkwood series expansion of equilibrium quantum statistics in the phase space formalism is generalized to include relativistic systems. For partition functions and all thermodynamic functions the series yields quantum corrections in terms of powers of h2 including systematic relativistic corrections given by modified Hankel functions 9@nc2/kT). An application to the symmetric relativistic quantum Toda oscillator is sketched.