Terrestrial laser scanning can provide high-resolution, two-dimensional sampling of soil surface roughness. While previous studies demonstrated the usefulness of these roughness measurements in geophysical applications, questions about the number of required scans and their resolution were not investigated thoroughly. Here, we suggest a method to generate digital elevation models, while preserving the surfaces stochastic properties at high frequencies and additionally providing an estimate of their spatial resolution. We also study the impact of the number and positions of scans on roughness indices estimates. An experiment over a smooth and isotropic soil plot accompanies the analysis, where scanning results are compared to results from active triangulation. The roughness measurement conditions for ideal sampling are revisited and updated for diffraction-limited sampling valid for close-range laser scanning over smooth and isotropic soil roughness. Our results show that terrestrial laser scanning can be readily used for roughness assessment on scales larger than 5 cm, while for smaller scales, special processing is required to mitigate the effect of the laser beam footprint. Interestingly, classical roughness parametrization (correlation length, root mean square height (RMSh)) was not sensitive to these effects. Furthermore, comparing the classical roughness parametrization between one- and four-scan setups shows that the one-scan data can replace the four-scan setup with a relative loss of accuracy below 1% for ranges up to 3 m and incidence angles no larger than 50, while two opposite scans can replace it over the whole plot. The incidence angle limit for the spectral slope is even stronger and is 40. These findings are valid for scanning over smooth and isotropic soil roughness.