By using a laser profiler, the roughness of sowed and plowed surfaces was obtained. Through evaluation of the precision of fractal dimensions based on the Weierstrass–Mandelbort (W–M) function, we found the rescaled (R/S) analysis method and the structure spectral method were not suitable for the calculation of the fractal dimension on a soft terrain surface. Therefore, the fractal dimension, non-scale range and correlation coefficient for each kind of terrain were analyzed using the following fractal computational methods – (i) variate-difference method, (ii) power spectral density method, and (iii) root mean square method. The results showed that: fractal dimension of plowed terrain was large with small fluctuations, while its internal structure was complex. The power spectral density method was not robust enough to compute the fractal dimension of a soft terrain surface. The fractal dimension computed using the root mean square method was found to be more accurate for the soft terrain surface. Moreover, the correlation coefficient of linear regression when using the root mean square method was good and the range of non-scale variation was small.