  We describe the first steps in the extension of the Symanzik O($a$)
improvement program for Wilson-type quark actions to anisotropic lattices, with
a temporal lattice spacing smaller than the spatial one. This provides a fully
relativistic and computationally efficient framework for the study of heavy
quarks. We illustrate our method with accurate results for the quenched
charmonium spectrum.
