A method of rapidly building radio environment map based on three-dimensional parabolic equation is studied, and the accuracy and applicability of several spatial interpolation methods are discussed, which provides an effective tool for network optimization and management. Firstly, the electric field distribution in large area and complex environment is calculated by using parabolic equation. Then, four typical research areas are selected with reference to the best statistical unit in geography, and Kriging, Inverse Distance Weighting(IDW), biharmonic spline interpolation and several interpolation methods based on triangulation are adopted to analyze the research area. It shows that several kinds of triangulation interpolation and biharmonic spline interpolation have the smallest RMSE and the shortest time in each region, but the degree of spatial fitting is the lowest; Kriging has the highest spatial fitting degree in each region, and the RMSE and time vary with the semi-variance model; IDW has a good spatial fitting and RMSE in each region, but it takes the longest time. Among several spatial interpolation methods, the Kriging of the exponential model can construct the radio environment map based on the three-dimensional parabolic equation with the highest accuracy and applicability.