Improving Coffee Yield Interpolation in the Presence of Outliers Using Multivariate Geostatistics and Satellite Data

Precision agriculture for coffee production requires spatial knowledge of crop yield. However, difficulties in implementation lie in low-sampled areas. In addition, the asynchronicity of this crop adds complexity to the modeling. It results in a diversity of phenological stages within a field and also continuous production of coffee over time. Big Data retrieved from remote sensing can be tested to improve spatial modeling. This research proposes to apply the Sentinel-2 vegetation index (NDVI) and the Sentinel-1 dual-polarization C-band Synthetic Aperture Radar (SAR) dataset as auxiliary variables in the multivariate geostatistical modeling of coffee yield characterized by the presence of outliers and assess improvement. A total of 66 coffee yield points were sampled from a 4 ha area in a quasi-regular grid located in southeastern Brazil. Ordinary kriging (OK) and block cokriging (BCOK) were applied. Overall, coupling coffee yield with the NDVI and/or SAR in BCOK interpolation improved the accuracy of spatial interpolation of coffee yield even in the presence of outliers. Incorporating Big Data for improving the modeling for low-sampled fields requires taking into account the difference in supports between different datasets since this difference can increase uncontrolled uncertainty. In this manner, we will consider, for future research, new tests with other covariates. This research has the potential to support precision agriculture applications as site-specific plant nutrient management.