Spatial interpolation methods are widely used in various fields and have been studied by several scholars with one or a few specific sampling datasets that do not reflect the complexity of the spatial characteristics and lead to conclusions that cannot be widely applied. In this paper, three factors that affect the accuracy of interpolation have been considered, i.e., sampling density, sampling mode, and sampling location. We studied the inverse distance weighted (IDW), regular spline (RS), and ordinary kriging (OK) interpolation methods using 162 DEM datasets considering six sampling densities, nine terrain complexities, and three sampling modes. The experimental results show that, in selective sampling and combined sampling, the maximum absolute errors of interpolation methods rapidly increase and the estimated values are overestimated. In regular-grid sampling, the RS method has the highest interpolation accuracy, and IDW has the lowest interpolation accuracy. However, in both selective and combined sampling, the accuracy of the IDW method is significantly improved and the RS method performs worse. The OK method does not significantly change between the three sampling modes. The following conclusion can be obtained from the above analysis: the combined sampling mode is recommended for sampling, and more sampling points should be added in the ridges, OPEN ACCESS ISPRS Int. J. Geo-Inf. 2015, 4 2743 valleys, and other complex terrain. The IDW method should not be used in the regular-grid sampling mode, but it has good performance in the selective sampling mode and combined sampling mode. However, the RS method shows the opposite phenomenon. The sampling dataset should be analyzed before using the OK method, which can select suitable models based on the analysis results of the sampling dataset.