This paper concerns the computation of nonlinear crest distributions for irregular Stokes waves, and a numerical algorithm based on the Fast Fourier Transform (FFT) technique has been developed for carrying out the nonlinear computations. In order to further improve the computational efficiency, a new Transformed Rayleigh procedure is first proposed as another alternative for computing the nonlinear wave crest height distributions, and the corresponding computer code has also been developed. In the proposed Transformed Rayleigh procedure, the transformation model is chosen to be a monotonic exponential function, calibrated such that the first three moments of the transformed model match the moments of the true process. The numerical algorithm based on the FFT technique and the proposed Transformed Rayleigh procedure have been applied to calculating the wave crest distributions of a sea state with a Bretschneider spectrum and a sea state with the surface elevation data measured at the Poseidon platform. It is demonstrated in these two cases that the numerical algorithm based on the FFT technique and the proposed Transformed Rayleigh procedure can offer better predictions than those from using the empirical wave crest distribution models. Meanwhile, it is found that our proposed Transformed Rayleigh procedure can compute nonlinear crest distributions more than 25 times faster than the numerical algorithm based on the FFT technique.