Estimates of particle-size made by operators in the field and laboratory represent a vast and relatively untapped data archive. The wide spatial distribution of particle-size estimates makes them ideal for constructing geological models and soil maps. This study uses a large data set from the Netherlands (n = 4837) containing both operator estimates of particle-size and complete particle-size distributions measured by laser granulometry. Operator estimates are inaccurate and imprecise relative to measured laser data; only 16.68% of samples were successfully classified using the Dutch classification scheme for fine sediment. Operator estimates of sediment particle-size encompass the same range of percentage values as particle-size distributions measured by laser analysis. However, the distributions measured by laser analysis show that most of the sand percentage values lie between 0 and 1, so the majority of the variability in the data is lost because operator estimates are made to the nearest 1% at best, and more frequently to the nearest 5%. Operator estimates made by three technicians trained by the Geological Survey of the Netherlands are found not to be influenced by bias, rather they exhibit very similar levels of accuracy and precision. This study compares five different methods of modelling complete particle-size distributions from sparse data: (i) a four-part Pearson's probability distribution function, (i) a log-linear interpolation, (iii) a logit-linear interpolation, (iv) a logistic probability distribution function and (v) a logit constrained cubic-spline (logit-CCS) interpolation. The logit-CCS interpolation performed best across all the samples used, although the performance of all models was very similar for normal Gaussian, skewed and peaked distributions. Predictions for bimodal distributions using the Pearson's, logit-linear and logistic models are markedly less accurate than both log-linear and logit-CCS interpolation models. Although the logit-CCS interpolation model produces the best predictions of continuous particle-size distributions, the low accuracy and precision of operator estimates does not warrant the use of such a complex algorithm. Given this, it is suggested that a standard log-linear interpolation is the most effective means of modelling complete particle-size distributions from sparse data. Interpolation-based models are recommended over probability distribution functions because they allow for a greater degree of flexibility and will always honour the available input data.