The orthogonal invariant descriptors are among the best region based shape descriptors which are used in many image processing and pattern recognition applications. The Chebyshev Harmonic Fourier Moments (CHFMs) are one of such invariant descriptor. They, however, suffer from high time complexity and numerical instability at high orders of moment. In this paper, we propose a fast method based on the recursive computation of radial kernel function of CHFMs which not only reduces time complexity but also improves their numerical stability.