Spherical harmonic transform plays an important role in research in shape description. Current computation methods involve expensive voxelization, and are prone to numerical errors associated with the size of the voxels. This paper describes a fast and accurate technique for computing spherical harmonic coecients directly from the description of the mesh.