RESUMO Perdas na produção agrícola, causadas principalmente por temperaturas extremas, chuvas e secas, têm exigido a adaptação da agricultura brasileira às mudanças climáticas em curso. Objetivou-se avaliar, por meio de uma revisão integrativa, os efeitos das práticas agrícolas estabelecidas por duas políticas públicas brasileiras (Plano ABC e Planaveg) na resiliência e adaptação dos sistemas agrícolas do país, considerando-se os efeitos no clima local e eventos extremos. Todas as coleções da biblioteca da Embrapa foram utilizadas como fontes principais e o Google Acadêmico como fonte secundária. A busca incluiu publicações de 1998 a abril de 2022 e resultou em 334 dados, extraídos de 51 referências, abrangendo pelo menos um fator e uma prática. A revisão mostrou que as práticas propostas podem contribuir para a resiliência e adaptação dos setores agropecuário e forestal do Brasil às mudanças climáticas nas principais regiões produtoras. Estudos indicam que a incorporação de componente arbóreo nos sistemas de produção pode proporcionar um ambiente mais resiliente a secas, aumento da umidade do ar, redução do impacto da temperatura e dos riscos de incêndio, entre outros benefícios. A adoção das práticas do Plano ABC e Planaveg pelo setor agropecuário é uma opção viável, estratégica e bem alinhada com uma economia de baixo carbono, ao mesmo tempo em que aumenta a resiliência aos impactos presentes e futuros das mudanças climáticas.

ABSTRACT Losses in agricultural production, mainly caused by extreme temperature, rainfall and drought, have demanded the adaptation of the Brazilian agriculture to the ongoing climate change. This study aimed to evaluate, through an integrative review, the effects of agricultural practices established by two Brazilian public policies (ABC Plan and Planaveg) on the resilience and adaptation of the country’s agricultural systems, considering the effects on local climate and extreme events. All the Embrapa library collections were used as main sources, and the Google Scholar was used as a secondary source. The search included publications from 1998 to April 2022 and resulted in 334 data, extracted from 51 references, covering at least one factor and one practice. The review showed that the proposed practices may contribute to the resilience and adaptation of the Brazilian agricultural and forestry sectors to climate change in the main production regions. Studies have indicated that the incorporation of the arboreal component into production systems may provide an environment more resilient to droughts, increased air humidity, reduced temperature impact and reduced fire risks, among other benefits. The adoption of the ABC Plan and Planaveg practices by the agricultural sector is an option feasible, strategic and well-aligned with a low-carbon economy, while it increases the resilience to present and future climate change impacts.

ABSTRACT: The aim of this study was to evaluate the performance of pedotransfer functions (PTFs) available in the literature to estimate soil bulk density (ρb) in different regions of Brazil, using different metrics. The predictive capacity of 25 PTFs was evaluated using the mean absolute error (MAE), mean error (ME), root mean squared error (RMSE), coefficient of determination (R2) and the regression error characteristic (REC) curve. The models performed differently when comparing observed and estimated ρb values. In general, the PTFs showed a performance close to the mean value of the bulk density data, considered as the simplest possible estimation of an attribute and used as a parameter to compare the performance of existing models (null model). The models developed by Benites et al. (2007) (BEN-C) and by Manrique and Jones (1991) (M&J-B) presented the best results. The separation of data into two layers according to depth (0-10 cm and 10-30 cm) demonstrated better performance in the 10-30 cm layer. The REC curve allowed for a simple and visual evaluation of the PTFs.

Resumo A modelagem por imputação múltipla com cadeias encadeadas é uma área de importância crescente cujos modelos e métodos têm sido muitas vezes desenvolvidos para lidar com aplicações específicas. No presente estudo, o modelo de imputação múltipla foi utilizado para estimar os dados de precipitação diária. Foram utilizados registros de precipitação de várias estações meteorológicas, obtidos a partir do sistema Agritempo para duas zonas climaticamente homogéneas. Os valores de precipitação obtidos para duas datas fixas (02 de Maio de 2005 e 20 de Janeiro de 2005) usando o modelo de imputação múltipla foram comparados com as técnicas geoestatisticas de krigagem ordinária e cokrigagem ordinária com a altitude como variável auxiliar. O modelo de imputação múltipla foi mais de 16% melhor para a primeira Zona e mais de 23% melhor para a segunda zona, comparadas com as estimativas de precipitação obtidas pelas técnicas geoestatísticas. O modelo provou ser uma técnica versátil, com resultados coerentes com as condições de diferentes zonas e épocas.

Abstract Modeling by multiple enchained imputation is an area of growing importance. However, its models and methods are frequently developed for specific applications. In this study the model for multiple imputation was used to estimate daily rainfall data. Daily precipitation records from several meteorological stations were used, obtained from system AGRITEMPO for two homogenous climatic zones. The precipitation values obtained for two dates (Jan. 20th 2005 and May 2nd 2005) using the multiple imputation model were compared with geo-statistics techniques ordinary Kriging and Co-kriging with the altitude as an auxiliary variable. The multiple imputation model was 16% better for the first zone and over 23% for the second one, compared to the rainfall estimation obtained by geo-statistical techniques. The model proved to be a versatile technique, presenting coherent results with the conditions of different zones and times.

Resumo: O objetivo deste trabalho foi desenvolver e avaliar um método para estimar a produtividade de milho com uso de um número mínimo de parâmetros e de informações limitadas sobre o manejo da cultura. O método proposto estima rendimentos potenciais e atingíveis com base no nível tecnológico dos sistemas de produção e em modelos agrometeorológicos relativamente simples. A produtividade de milho foi estimada para as safras de 2000/2001 a 2007/2008, tendo-se considerado vários locais e regiões do Brasil, e comparada aos dados de produção reais de levantamentos oficiais. A produtividade estimada apresentou alta correlação com a observada (0,76≤R2<0,92; p<0,01), com eficiência do modelo (E1’) entre 0,45 e 0,73; erro médio relativo (MRE) entre -0,9 e 2,4%; e erro médio absoluto (MAE) inferior a 70 kg ha-1, de acordo com o nível tecnológico considerado. Com base nestes resultados, este modelo pode ser recomendado para estimativas de produtividade em todo o País, e contribuir para tornar este processo mais preciso e exato.

Abstract: The objective of this work was to develop and evaluate a method for estimating corn yield using a minimum number of parameters and limited information about crop management. The proposed method estimates potential and attainable yields based on the technological level of the production systems and on relatively simple agrometeorological models. Corn yield was estimated for the crop seasons from 2000/2001 to 2007/2008, considering several locations and regions in Brazil, and was compared with the actual yield data from official surveys. There was a high correlation between the estimated and observed yield (0.76≤R2<0.92; p<0.01), with model efficiency (E1’) ranging from 0.45 to 0.73; mean relative error (MRE) between -0.9 and 2.4%; and mean absolute error (MAE) of less than 70 kg ha-1, depending on the technological level adopted. Based on these results, the proposed yield model can be recommended to forecast yields all over the country, contributing to make this process more precise and accurate.

O objetivo deste trabalho foi avaliar um sistema de estimativa de produtividade da cultura do arroz no Brasil, com base em modelos agrometeorológicos simples e no nível tecnológico dos sistemas produtivos. Este sistema de estimativa incorpora as bases conceituais propostas por Doorenbos & Kassam para produtividades potenciais e alcançáveis com ajustes empíricos para produtividade máxima e sensibilidade da cultura à deficiência hídrica, considerando-se cinco classes de produtividade de arroz. A produtividade de arroz foi estimada de 2000/2001 a 2007/2008 e comparada aos dados do IBGE. Análises de regressão entre as estimativas do modelo e dados do IBGE resultaram em coeficientes de determinação significativos, com menor dispersão na região Sul do que nas regiões Norte e Nordeste do país. O índice de eficiência do modelo (E1') variou de 0,01, nas classes de menor produtividade, a 0,45 nas classes de maior produtividade, e o erro médio absoluto variou de 58 a 250 kg ha‑1, respectivamente.

The objective of this work was to evaluate an estimation system for rice yield in Brazil, based on simple agrometeorological models and on the technological level of production systems. This estimation system incorporates the conceptual basis proposed by Doorenbos & Kassam for potential and attainable yields with empirical adjusts for maximum yield and crop sensitivity to water deficit, considering five categories of rice yield. Rice yield was estimated from 2000/2001 to 2007/2008, and compared to IBGE yield data. Regression analyses between model estimates and data from IBGE surveys resulted in significant coefficients of determination, with less dispersion in the South than in the North and Northeast regions of the country. Index of model efficiency (E1') ranged from 0.01 in the lower yield classes to 0.45 in higher ones, and mean absolute error ranged from 58 to 250 kg ha‑1, respectively.

O objetivo deste trabalho foi avaliar o efeito do aumento na temperatura, previsto pelo Painel Intergovernamental de Mudanças Climáticas (IPCC), sobre o risco climático para a produção de algodão no Brasil. O aumento previsto na temperatura fez com que a região Nordeste apresentasse a maior redução na área de baixo risco para cultivo de algodão, em comparação com as demais regiões do país. Essa área pode diminuir de 83 milhões de ha no ano de 2010 para aproximadamente 71 milhões de ha em 2040, o que significaria 15% de redução em 30 anos. As regiões Sudeste e Centro-Oeste apresentam pequena redução nas áreas aptas para o cultivo de algodão até 2040, enquanto a região Norte não teve redução nessas áreas. A elevação das temperaturas não é benéfica para a produção do algodão no Brasil porque a área de baixo risco para o cultivo econômico pode diminuir.

The objective of this work was to evaluate the effect of the temperature increase forecasted by the Intergovernmental Panel on Climate Change (IPCC) on agricultural zoning of cotton production in Brazil. The Northeastern region showed the highest decrease in the low-risk area for cotton cultivation due to the projected temperature increase. This area in the Brazilian Northeast may decrease from 83 million ha in 2010 to approximately 71 million ha in 2040, which means 15% reduction in 30 years. Southeastern and Center-Western regions had small decrease in areas suitable for cotton production until 2040, while the Northern region showed no reduction in these areas. Temperature increase will not benefit cotton cultivation in Brazil because dimension of low-risk areas for economic cotton production may decrease.

O objetivo deste trabalho foi quantificar a contribuição da variável auxiliar altitude, na estimativa da distribuição espacial da precipitação anual média no Estado de São Paulo. A estatística quadrado médio do erro (QME) foi usada em dois conjuntos de observações de precipitação anual média (1957 a 1997): o completo, com 1.027 observações, e o reduzido, com 445. Bolsões de precipitação foram perfeitamente definidos nos mapas de variabilidade espacial que utilizaram o conjunto completo de dados, e indicaram a existência de possíveis microclimas. O interpolador geoestatístico de krigagem ordinária apresentou desempenho 82 vezes mais preciso que o interpolador do inverso do quadrado da distância, quando o QME foi usado como critério de comparação para o conjunto de dados completo. Para o conjunto reduzido, essa magnitude foi de duas vezes. Os erros de estimação obtidos por krigagem ordinária foram menores no conjunto completo, enquanto os obtidos por cokrigagem ordinária foram menores no reduzido. Isso indica que esses interpoladores devem ser usados para determinação da distribuição espacial da precipitação anual média. O uso da altitude como variável auxiliar beneficia o interpolador de cokrigagem ordinária e define microrregiões mais uniformes quanto à distribuição espacial da precipitação anual média.

The objective of this work was to quantify the contribution of the auxiliary variable elevation in the estimate of the spatial distribution of annual average rainfall in the state of São Paulo, Brazil. The statistic mean square error (MSQR) was used in two observation sets for the mean annual rainfall (1957 to 1997): a full set, with 1,027 observations, and a reduced one, with 445. Areas of rainfall concentration were clearly defined in spatial variability maps that used the full data set, indicating possible microclimates. Ordinary kriging geostatistical interpolation had a performance 82 times more accurate than that of interpolation by the inverse square of the distance when the MSQR was used as a benchmark for the full set. For the reduced set, this magnitude was of two times. The estimation errors obtained by ordinary kriging were lower in the full set, whereas those obtained by ordinary cokriging were lower in the reduced one. This indicates that these interpolators should be used for determining the spatial distribution of annual average rainfall. The use of altitude as an auxiliary variable benefits the ordinary cokriging interpolator and defines more uniform micro regions as to the spatial distribution of annual average rainfall.

Um conceito de eficiência de produção agrícola foi derivado da Termodinâmica, para gerar um indicador para avaliação do desempenho de sistemas de produção agrícola e quantificação da interferência dos fatores determinantes desse desempenho. Em Termodinâmica, a eficiência de um processo é dada pela razão entre a energia recuperada e a energia utilizada. Para estabelecer essa relação em sistemas agrícolas, admitiu-se que a energia utilizada seja dada pela produtividade atingível, estimada por modelos de simulação desenvolvidos com base em variáveis ambientais. O método das zonas agroecológicas da FAO foi utilizado para a estimativa da produtividade atingível, e os dados do Instituto Brasileiro de Geografia e Estatística (IBGE), como produtividade observada. Avaliou-se a eficiência da produção de cana-de-açúcar no Estado de São Paulo, em duas safras, e calculou-se sua correlação com alguns fatores do meio físico que regulam a produção. Foi identificada forte relação da eficiência da produção agrícola com a aptidão agrícola do solo. Isto permitiu inferir sobre o efeito de fatores relativos ao agronegócio sobre a eficiência da produção agrícola. Foi possível identificar, ainda, que a radiação solar, a deficiência hídrica, a quantidade de chuva e a temperatura máxima são as principais variáveis climáticas que afetam a eficiência da produção agrícola.

A conceptual framework for crop production efficiency was derived using thermodynamic efficiency concept, in order to generate a tool for performance evaluation of agricultural systems and to quantify the interference of determining factors on this performance. In Thermodynamics, efficiency is the ratio between the output and input of energy. To establish this relationship in agricultural systems, it was assumed that the input energy is represented by the attainable crop yield, as predicted through simulation models based on environmental variables. The method of FAO's agroecological zones was applied to the assessment of the attainable sugarcane yield, while Instituto Brasileiro de Geografia e Estatística (IBGE) data were used as observed yield. Sugarcane efficiency production in São Paulo state was evaluated in two growing seasons, and its correlation with some physical factors that regulate production was calculated. A strong relationship was identified between crop production efficiency and soil aptitude. This allowed inferring the effect of agribusiness factors on crop production efficiency. The relationships between production efficiency and climatic variables were also quantified and indicated that solar radiation, annual rainfall, water deficiency, and maximum air temperature are the main factors affecting the sugarcane production efficiency.

Este trabalho teve como objetivo indicar algumas possíveis manipulações ecofisiológicas, com base em experimentos de arquitetura de plantas. Plantas de erva-mate, seringueira e videira foram caracterizadas no campo e no laboratório, e suas maquetes 3D foram reconstruídas com os programas V-Plants e PlantGLViewer. Indicações sobre a análise de crescimento e impactos ambientais, com a aplicação de diversas simulações, foram discutidas.

The objective of this work was to indicate some possible ecophysiological manipulations, based on experiments of plant architectural traits. Yerbamate, rubber tree and grapevine plants were measured in field and laboratory, and their 3D mockups were constructed using V-Plants and PlantGLViewer softwares. Indications about growth and environmental impact analysis were discussed, applying diverse simulations in modeling process.

Este trabalho teve por objetivo avaliar um sistema de previsão de safra de soja para o Brasil, baseado em modelos empíricos regionalizados para estimativa da produtividade, a partir de um banco de dados de área cultivada em escala municipal, e de um sistema de monitoramento agrometeorológico de abrangência nacional. Os modelos incorporam a base conceitual proposta por Doorenbos & Kassam, com ajustes empíricos para cada região do Brasil, considerando-se as diferenças quanto ao potencial produtivo das principais variedades e as peculiaridades dos sistemas de produção utilizados nas diferentes regiões. A base de informações de área cultivada de cada Estado foi constituída por dados provenientes do IBGE e da Conab. A produção de soja, entre as safras 2000/2001 e 2005/2006, foi estimada e comparada com os levantamentos da Conab. A análise estatística pelo teste t indica não haver diferença entre as estimativas e os dados oficiais. Bons ajustes (R²>0,87) foram obtidos para as produções regionalizadas, com desvios mais expressivos nos estados do Rio Grande do Sul, Santa Catarina, Mato Grosso do Sul, Maranhão, Piauí e Bahia. Para os dados consolidados em escala nacional, o maior desvio observado foi de 5,81%, na safra 2000/2001, e o menor de 0,62% na safra 2005/2006.

The aim of this work was to evaluate a system to forecast soybean crop yield for Brazil, based on regional empirical models to assess crop yield, with data from a national database of soybean cropped area in municipal scale, and from an agrometeorological monitoring system covering all Brazilian States. The models incorporated the conceptual bases proposed by Doorenbos & Kassam, with empirical adjustments for each region of Brazil, considering the diverse potential yield of most used varieties and the regional differences in the soybean production systems. Cultivated area database from each State was structured with data from IBGE and Conab. Soybean yield was estimated for harvests from 2000/2001 to 2005/2006 and compared to Conab surveys. Statistical analysis using Test t point out that there is no significant difference between estimates and official data. Good fittings were gotten for region grouped yield data (R²>0.87), with higher deviations for harvest assessment of Rio Grande do Sul, Santa Catarina, Mato Grosso do Sul, Maranhão, Piauí and Bahia. In national scale, the highest observed deviation was 5.81% for the 2000/2001 harvest, and the smaller one was 0.62% for the 2005/2006 yield.

O objetivo deste trabalho foi avaliar o impacto do aumento da resolução espacial e radiométrica da imagem pancromática do Ikonos-II na identificação de plantios de café (Coffea arabica), em comparação com as imagens do Landsat/ETM+. A área de estudo está localizada no Município de Pedregulho, SP, onde foram selecionados 50 talhões com plantios de café, e foram levantados dados referentes à altura, idade, espaçamento e variedade de cada talhão. As imagens permitiram a identificação de talhões com características diferentes em campo, tendo-se destacado a imagem do Ikonos-II, que apresentou melhor desempenho. Para os talhões com características iguais em campo, as imagens analisadas não se mostraram eficientes, independentemente do satélite utilizado. As correções atmosféricas e radiométricas, na imagem do Ikonos-II, não proporcionaram ganho efetivo nas análises realizadas. A maioria dos talhões identificados na imagem do Ikonos-II pode ser localizada na imagem do Landsat/ETM+ (68%). A correlação significativa entre a banda 4 do Landsat/ETM+ e o canal pancromático do Ikonos-II indica uma forma de ligação entre as imagens dos dois satélites.

The objective of this work was to assess the impact of the better spatial and radiometric resolutions of the Ikonos-II panchromatic image, for the identification of coffee (Coffea arabica) planting areas, in comparison with a Landsat/ETM+ image. The area of study is situated in the city of Pedregulho, State of São Paulo, Brazil, where 50 coffee fields were selected. Information about plants features, like height, age, spacing and variety were collected. Images allowed the identification of coffee areas with different field features, and the calibrated Ikonos-II image showed the best results. Considering the areas with similar field features, images from both satellites were not efficient in the coffee identification. The atmospheric and radiometric corrections applied on the Ikonos-II image did not improve the analyses results. More than half of the identified areas in the Ikonos-II image could be found in the Landsat/ETM+ image (68%). The significant correlation between Landsat/ETM+ band 4 and Ikonos-II panchromatic channel shows a link between both satellite images.

O objetivo deste trabalho foi identificar, no Brasil, regiões homogêneas quanto à distribuição de probabilidades de chuva e assim contribuir para estudos de riscos climáticos na agricultura. As regiões foram delimitadas mediante aplicação da análise de agrupamento hierárquica, com variáveis classificatórias definidas pela proporção de pêntadas secas e por medidas de posição, escala e forma das distribuições de freqüências da quantidade de chuva. A análise de agrupamento permitiu identificar 25 zonas pluviometricamente homogêneas em todo o território brasileiro.

The objective of this work was to identify homogeneous areas in Brazil according to the similarity of their rainfall probability distribution, in order to contribute to the studies of climatic risks in agriculture. These areas were delimited using the hierarchical cluster analysis, with discriminating variables defined by the proportion of dry pentads and by measures of location, scale and shape of the pluviometric frequency distributions. Cluster analysis was able to identify 25 rainfall homogeneous areas in the whole Brazilian territory.

A partir das indicações do último relatório do IPCC (International Pannel of Climatic Change), foram feitas várias simulações e avaliados os impactos que um aumento na temperatura média do ar de 1ºC, 3ºC e 5,8ºC e um incremento de 15% na precipitação pluvial teriam na potencialidade da cafeicultura brasileira, definida pelo atual zoneamento agroclimático do café (Coffea arábica L.) nos Estados de Goiás, Minas Gerais, São Paulo e Paraná. Os resultados indicaram uma redução de área apta para a cultura superior a 95% em Goiás, Minas Gerais e São Paulo, e de 75% no Paraná, no caso de um aumento na temperatura de 5,8ºC. Esses resultados são válidos se mantidas as atuais características genéticas e fisiológicas das cultivares de café arábica utilizadas no Brasil, que têm como limite de tolerância temperaturas médias anuais entre 18ºC e 23ºC.

According to the last report of the Intergovernmental Panel on Climate Change (IPCC), the global temperature is supposed to increase 1°C to 5.8°C and the rainfall 15% in the Tropical area. This paper analyses the effect that these possible scenarios would have in the agroclimatic zoning of the arabic coffee (Coffea arabica L.) main plantation areas in Brazil. The results indicated a reduction of suitable areas greater than 95% in Goiás, Minas Gerais and São Paulo and about 75% for Paraná in the case of a temperature increase of 5.8ºC. These results presume that all the physiological characteristics of the crop will be the same for the varieties analyzed and that the ideal climatic condition for economic development is mean annual temperatures between 18ºC and 23ºC.

A definição das melhores datas para plantio de milho, adotadas no Zoneamento de Riscos Climáticos do Estado de São Paulo, é baseada nas variações espaço-temporais do índice de satisfação das necessidades de água da planta (ISNA), sendo a média ponderada o método de interpolação utilizado na espacialização deste parâmetro. Este método não considera a propagação espacial do erro da estimativa, o que permite interpretações imprecisas sobre as melhores datas de plantio, principalmente nas situações de início e de final do ciclo. O objetivo deste trabalho foi comparar métodos de espacialização de valores numéricos dos índices agrometeorológicos e avaliar sua variação espacial. Foram utilizados os métodos da média ponderada, krigagem ordinária e krigagem por indicação. A krigagem por indicação foi o método mais apropriado para espacializar o índice ISNA e definir a melhor data de plantio do milho.

The definition of the best sowing dates for corn, adopted in the Zoning of Climatic Risks of the State of São Paulo, Brazil, was based on space-temporary variations of the water stress index (ISNA), using the pondered average as the interpolation method to spatialize this parameter. This method does not consider the spatial propagation of errors, leading to imprecise interpretations of the best sowing dates, mainly at the beginning and at the end of the cycle. The objective of this work was to compare methods of spatializing numeric values of agrometeorological indices and to evaluate their spatial space variation. The methods of the pondered means ordinary kriging and indication kriging were used. The indication kriging was the most suitable method for spatializing ISNA and to define the best sowing date for corn, in the State of São Paulo.

Para eliminar divergências na interpretação dos resultados e agilizar os atuais métodos de detecção de fraudes em café torrado e moído, foi estabelecido um método baseado na análise por imagem e fundamentado no princípio de que diferentes materiais de origem orgânica, como o pó de café, podem apresentar reflectâncias distintas em diferentes comprimentos de onda do espectro eletromagnético. Partiu-se da hipótese de que o pó de café adulterado, quando submetido a uma fonte artificial de iluminação, apresenta uma reflectância, nos canais vermelho (R), verde (G) e azul (B), diferente em relação à do pó de café não-adulterado. Após as etapas de limpeza, secagem e homogeneização, foram geradas imagens multiespectrais das amostras de café, por meio de uma lupa acoplada a uma câmara CCD (Charge Coupled Device). A quantificação das impurezas na amostra foi obtida utilizando-se curvas de calibração entre a área relativa obtida pela classificação supervisionada de imagens e a porcentagem de impurezas presentes nas amostras. Esse novo método permite agilidade da resposta, ausência de subjetividade nos resultados e não-destruição das amostras analisadas, e assegura um patamar mínimo de detecção de 95% das impurezas do produto.

In order to reduce divergences in the interpretation of the results obtained by current methods of fraud detection in powdered coffees, a new method to quantify coffee adulterants is presented in this study. This method is based on image analyses and on the principle that different organic materials found in ground coffees present distinct spectral signatures. The hypothesis is that the adultered coffees, when submitted to an artificial source of illumination, present different reflectance values in RGB (Red/Green/Blue) channels in comparison to those from pure coffees. After the cleaning, drying, and homogenization steps, multispectral images of coffee samples were generated by using a glass magnifier connected to a CCD (Charge Coupled Device) camera. The images were stored and classified by means of a digital image software. The spectral signatures of each component in the coffee sample were identified from their histograms in percentage, and the amounts (by weight) of adulterants in the coffee sample were obtained by calibration curves. The advantages of this method are the capability to analyse a great number of samples in a relatively short time period, the lack of subjectivity in the results, and the non destruction of the coffee samples. This method reaches a minimum accuracy of 95% in the quantification of adulteration in the powdered coffee.