In modern agroecology, one of the most pressing problems is the problem of spatial data mapping. The development of information technology opens up a wide range of approaches for solving this problem. One of these approaches is based on the use of geostatistical methods. This study was carried out with the aim of developing ideas about the applicability of the ordinary kriging method for predicting the spatial distribution of the agro-ecological indicator with identifying the boundaries of in-field heterogeneity according to remote sensing data. For the model computational experiment, aerial photographs of the agricultural field in the red and near infrared ranges were used, which made it possible to obtain sets of uniformly distributed values of the vegetative index NDVI that were randomly generated. The high spatial resolution of the images allowed us to analyze the observational data for the studied agricultural field.