Coffee is the second most valuable traded commodity worldwide. Brazil is the world’s largest coffee producer, responsible for one third of the world production. A coffee plot exhibits high and low production in alternated years, a characteristic so called biennial yield. High yield is generally a result of suitable conditions of foliar biomass. Moreover, in high production years one plot tends to lose more leaves than it does in low production years. In both cases some correlation between coffee yield and leaf biomass can be deduced which can be monitored through time series of vegetation indices derived from satellite imagery. In Brazil, a comprehensive, spatially distributed study assessing this relationship has not yet been done. The objective of this study was to assess possible correlations between coffee yield and MODIS derived vegetation indices in the Brazilian largest coffee-exporting province. We assessed EVI and NDVI MODIS products over the period between 2002 and 2009 in the south of Minas Gerais State whose production accounts for about one third of the Brazilian coffee production. Landsat images were used to obtain a reference map of coffee areas and to identify MODIS 250 m pure pixels overlapping homogeneous coffee crops. Only MODIS pixels with 100% coffee were included in the analysis. A wavelet-based filter was used to smooth EVI and NDVI time profiles. Correlations were observed between variations on yield of coffee plots and variations on vegetation indices for pixels overlapping the same coffee plots. The vegetation index metrics best correlated to yield were the amplitude and the minimum values over the growing season. The best correlations were obtained between variation on yield and variation on vegetation indices the previous year (R = 0.74 for minEVI metric and R = 0.68 for minNDVI metric). Although correlations were not enough to estimate coffee yield exclusively from vegetation indices, trends properly reflect the biennial bearing effect on coffee yield.

O objetivo deste trabalho foi avaliar uma nova metodologia para mapeamento da cultura da soja no Estado de Mato Grosso, por meio de imagens Modis e de diferentes abordagens de classificação de imagens. Foram utilizadas imagens diárias e imagens de 16 dias. As imagens diárias foram diretamente classificadas pelo algoritmo Isoseg. As duas séries de imagens de 16 dias, referentes ao ciclo total e à metade do ciclo da cultura da soja, foram transformadas pela análise de componentes principais (ACP), antes de serem classificadas. Dados de referência, obtidos por interpretação visual de imagens do sensor TM/Landsat-5, foram utilizados para a avaliação da exatidão das classificações. Os melhores resultados foram obtidos pela classificação das imagens do ciclo total da soja, transformadas pela ACP: índice global de 0,83 e Kappa de 0,63. A melhor classificação de imagens diárias mostrou índice global de 0,80 e Kappa de 0,55. AACP aplicada às imagens do ciclo total da soja permitiu o mapeamento das áreas de soja com índices de exatidão melhores do que os obtidos pela classificação derivada das imagens de data única.

The objective of this study was to assess correlations between coffee yield and MODIS derived vegetation indices. We assessed EVI and NDVI MODIS products in the south of Minas Gerais state whose production accounts for about one third of the Brazilian coffee production. Correlations were observed between variations on yield of coffee plots and variations on vegetation indices for pixels overlapping the same coffee plots. The vegetation index metrics best correlated to yield were the amplitude and the minimum values over the growing season. The best correlations were obtained between variation on vegetation indices and variation in yield the following year. In general EVI presented better results than those of NDVI for all vegetation index metrics evaluated. Although correlations were not enough to estimate coffee yield exclusively from vegetation indices, trends properly reflect the biennial bearing effect on coffee yield.