Estimation of soil particle size fractions (percentage of sand, clay and silt) with sparse data can be improved by taking into account the spatial correlation with other variables. Reflectance of the near infrared band (0.76-0.90μm) is used as an auxiliary variable in the prediction of soil particle size fractions on a 100m×100m grid in a 2200m×1300m field. The results of kriging using the textural information alone are compared with those of cokriging with the auxiliary variable. Cokriging gives a higher correlation coefficient and a lower mean squared error between estimates and measurements than kriging. The most significant improvement of the estimation is in areas of the field with sparse data for the estimated variables. The relative improvement is up to 90% in terms of reduced kriging variance and up to 33% for the mean squared error between actual measurements and estimated values. Only 17-25% of the original observations of texture are needed to obtain relatively accurate estimation when cokriging with reflectance data.