Based on the power spectral density (PSD) function of stochastic irregularities of the standard grade road and by means of inverse fast Fouerier transform (IFFT) based on discretized PSD sampling, an equivalent sample of stochastic road surface model in time domain was built. A one-dimensional model of stochastic road was developed into a 2D model of stochastic road surface. Through computer simulation practice based on the MATlab, a 2D sample of stochastic road surface in time domain was regenerated. Furthermore, given the sample data, the PSD was estimated and then compared with the theoretical 2D PSD Equation deduced from the one-dimensional PSD expreesion so as to prove the effectiveness and accuracy of the time-domain model regeneration of 2D stochastic road surface by means of IFFT method. The 2D stochastic road surface model directly provided basic road excitation input data for virtual prototyping (VP) and virtual proving ground (VPG) technology.