Semina: Ciências Agrárias 1676-546X: Issn: Semina - Agrarias@
Semina: Ciências Agrárias 1676-546X: Issn: Semina - Agrarias@
Semina: Ciências Agrárias 1676-546X: Issn: Semina - Agrarias@
ISSN: 1676-546X
[email protected]
Universidade Estadual de Londrina
Brasil
Groff, Edson Cristiano; Nanni, Marcos Rafael; Pinheiro Povh, Fabrício; Cezar, Everson
Características agronômicas associadas com índices de vegetação medidos por sensores
ativos de dossel na cultura da soja
Semina: Ciências Agrárias, vol. 34, núm. 2, marzo-abril, 2013, pp. 517-525
Universidade Estadual de Londrina
Londrina, Brasil
Resumo
A soja sempre se destacou pela sua importância econômica no cenário agrícola nacional tendo sido
observado aumento do seu cultivo em todo o território brasileiro. Com o uso crescente da agricultura de
precisão (AP) é de suma importância entender a variabilidade das variáveis agronômicas. Atualmente,
sensores espectrais em nível terrestre vêm sendo utilizados para estabelecer relações entre as respostas
espectrais e parâmetros de crescimento da cultura. Isso é possível por meio da leitura dos índices de
vegetação, como o Índice de Vegetação da Diferença Normalizada utilizando a banda verde (GNDVI).
Este estudo foi conduzido entre outubro de 2008 e abril de 2009, em Ponta Grossa – PR, com o objetivo
de avaliar a resposta espectral da cultura da soja por meio do GNDVI e sua correlação com algumas
variáveis agronômicas, como: produtividade de grãos, produtividade de massa seca, densidade de plantas
e altura de plantas. Leituras do GNDVI foram realizadas em seis estádios fenológicos (EF) utilizando o
sensor óptico ativo Crop Circle. Análises de regressão foram utilizadas para determinar quais variáveis
agronômicas estabeleciam relação com o GNDVI. O GNDVI na cultura da soja estabeleceu correlações
signi cativas apenas com a produtividade de grãos e a produtividade de MS no EF V9. Assim, para se
estabelecer a correlação do GNDVI com essas variáveis recomenda-se que as leituras sejam realizadas
nesse EF.
Palavras-chave: Índice de vegetação, resposta espectral e sensoriamento remoto
Abstract
Soybean has always stood out for its economic importance in the national agricultural scenario with
an increase in its cultivation throughout the Brazilian territory. With the increasing use of precision
agriculture (PA) is of paramount importance to understand the variability of agronomic variables.
Currently, spectral sensors in the ground-level have been used to establish relationships between spectral
response and crop growth parameters. This is possible through the reading of vegetation indices such
as Normalized Difference Vegetation Index using the green band (GNDVI). This study was conducted
between October to 2008 and April to 2009, in Ponta Grossa – PR, with the objective of evaluating the
spectral response of soybean through GNDVI and its correlation with some agronomic variables, such
as: productivity of grain, productivity of dry mass, plant density and plant height. GNDVI readings
were taken on six phenological stages using Crop Circle active optical sensor. Regression analyzes were
used to determine which agronomic variables established relationship with the GNDVI. The GNDVI
1
Engº Agrº, MSc. em Agronomia pela Universidade Estadual de Maringá, UEM, Maringá, PR. E-mail: [email protected]
2
Prof. Dr. Deptº Agronomia, UEM/CCA, Maringá, PR. E-mail: [email protected]
3
Engº Agrº Dr. Coordenador na área de Agricultura de Precisão da Fundação ABC, Castro, PR. E-mail: fabricio@fundacaoabc.
org.br
4
Discente de Doutorado em Agronomia, UEM, PR. E-mail: [email protected]
*
Autor para correspondência
Recebido para publicação 12/09/11 Aprovado em 21/12/12
517
Semina: Ciências Agrárias, Londrina, v. 34, n. 2, p. 517-526, mar./abr. 2013
Groff, E. C. et al.
in soybean established signi cant correlations only with grain and MS productivity in EF V9. Thus, to
establish the correlation of GNDVI with these variables it is recommended that the readings are carried
out in this EF.
Key words: Vegetation index, spectral responses, remote sensing
Implantação e manejo da cultura da soja Figura 1. Sensor óptico ativo acoplado à barra de
pulverização, coletor de dados e receptor GPS.
A cultura da soja foi implantada, em uma área
de 7,50 ha, utilizando-se a cultivar BRS 133. A
semeadura foi realizada no dia 18 de outubro de
2008, em sistema de semeadura direta, sobre a
resteva de aveia preta, previamente dessecada,
utilizando-se o espaçamento entre linhas de 0,42 m
e população inicial de 306.000 plantas ha-1.
No estabelecimento da cultura, foram retiradas
amostras de solo para a realização da análise
química conforme resultados: pH (CaCl2): 5,0;
fósforo (resina): 27 mg.dm-3; H + Al: 6,2 cmol.
dm-3; Al: 0,3 cmol.dm-3; K: 0,2 cmol.dm-3; Ca: 4,9
cmol.dm-3; Mg 2,4 cmol.dm-3; S.B.: 7,5 cmol.dm-3; Fonte: Elaboração dos autores.
C.T.C.: 13,7 cmol.dm-3; M.O.: 54 mg.dm-3. E com
base nos resultados, foi realizada na semeadura,
As leituras foram realizadas com passadas a
a adubação com 300 kg ha-1 do adubo formulado
cada 9 m, respeitando-se as linhas de pulverização.
NPK 00-20-20 e quarenta dias após a semeadura
A velocidade de deslocamento do trator foi ajustada
efetuou-se a adubação de cobertura com 100 kg ha-1
a m de gerar um ponto a cada 1 m totalizando,
de cloreto de potássio seguindo-se a recomendação
aproximadamente, 4.000 pontos na área total.
proposta pela EMBRAPA 2001.
Estes pontos foram interpolados pelo método da
média simples, com resolução espacial de 5 x 5
Avaliação do índice de vegetação da diferença m, utilizando-se o programa Spring 4.3.3. Com
normalizada (GNDVI) base na grade regular de pontos, gerados pela
interpolação, amostrou-se, em torno de cada ponto
Foram realizadas seis leituras do GNDVI em
georreferenciado, nove pontos para a leitura do
diferentes estádios fenológicos (EF): (12/12/08 EF
GNDVI. Conforme especi cação do fabricante,
V5, 24/12/08 EF V9, 06/01/09 EF R1, 15/01/09 EF
as leituras foram realizadas entre 0,80 e 1,20 m
R3, 29/01/09 EF R6 e 26/03/09 EF R8), utilizando-
de altura do dossel (ápice da cultura), pois, neste
se um sensor óptico ativo, modelo Crop Circle®
intervalo, o sensor garante a qualidade da leitura
ACS-210 (Holland Scienti c), instalado na barra de
(HOLLAND SCIENTIFIC, 2004).
um pulverizador hidráulico e, acoplados no trator,
o coletor de dados modelo GeoSCOUT GLS-420
(com cartão de memória para armazenamento de Determinação das variáveis agronômicas
dados), o GPS modelo SR20 da marca Leica, de
Produtividade de grãos
precisão métrica (Figura 1) e a antena externa. O
sensor fornece os valores de GNDVI, por meio da As avaliações das variáveis agronômicas
medição da re ectância do dossel nos comprimentos seguiram os pontos georreferenciados distribuídos
de onda de 590 nm (visível) e 880 nm (infravermelho aleatoriamente, conforme a Figura 2.
próximo).
519
Semina: Ciências Agrárias, Londrina, v. 34, n. 2, p. 517-526, mar./abr. 2013
Groff, E. C. et al.
520
Semina: Ciências Agrárias, Londrina, v. 34, n. 2, p. 517-526, mar./abr. 2013
Características agronômicas associadas com índices de vegetação medidos por sensores ativos de dossel...
Variáveis
Equações r2 p
Dependente Independente
GNDVI V5 y = 2,8929x +1,418** 0,16 0,00
GNDVI V9 y = 7,5610x – 1,840** 0,59 0,00
GNDVI R1 y = 4,5236x – 0,234* 0,09 0,03
Produtividade
GNDVI R3 y = 5,8861x – 1,406 0,04 ns
GNDVI R6 y = 6,4586x – 1,8316* 0,10 0,02
GNDVI R8 y = 1,2237x + 2,476 0,01 ns
ns: não-signi cativo. * e ** Signi cativo a 5 e 1% de probabilidade, respectivamente.
Fonte: Elaboração dos autores.
No EF V5, o baixo coe ciente de determinação o GNDVI. Esta observação está de acordo com
pode ter resultado da grande interferência do solo Myneni et al. (2002) que analisaram o NDVI para
nesta fase. Conforme Dalmolin et al. (2005), vários diferentes dosséis, e veri caram que a saturação não
elementos que compõem a superfície do solo (matéria permite a diferenciação de zonas que caracterizem a
orgânica, textura, composição mineralógica, variação do dossel e que possam ser correlacionadas
rugosidade e umidade) podem afetar a re ectância com a produtividade.
do dossel e, por consequência, o GNDVI. Esse
No EF V9, o coe ciente de determinação
conjunto de elementos é denominado de “superfície
observado foi de 0,59, pois nesse estádio a cultura
de fundo” (HUETE, 1985). Esta superfície de
encontrava-se no pleno desenvolvimento e sem
fundo, presente na entre linha da cultura quando a
a presença de ores que poderiam in uenciar
mesma ainda não fechou completamente, interfere
a resposta espectral do dossel. Assim, para se
na resposta espectral do dossel que, por sua vez,
estabelecer a correlação das medidas de GNDVI
afeta a determinação do GNDVI.
com a produtividade de grãos, recomenda-se que as
A partir do EF R1, os baixos valores encontrados leituras sejam realizadas neste estádio fenológico.
nos coe cientes de determinação podem resultar Esse é um fator que limita o uso deste índice
do pleno desenvolvimento da cultura, que faz com para a estimativa da produtividade, pois eventos
que a re ectância torne-se homogênea saturando que ocorram posteriormente a este EF não serão
521
Semina: Ciências Agrárias, Londrina, v. 34, n. 2, p. 517-526, mar./abr. 2013
Groff, E. C. et al.
quanti cados, como, p.e., o ataque de pragas, as uma vez que os valores de IAF variaram de 2,0 a
doenças e outros danos ocasionados por fenômenos 4,5, ou seja, aquém da saturação do NDVI para o
climáticos. milho. Macedo e Rudorff (2003) observaram que o
NDVI explicou 54% da variação na produtividade
Batista; Rudorff; Oviedo (1988) avaliaram
de grãos do milho safrinha, e os valores de NDVI
duas cultivares de soja e encontraram coe ciente
variaram de 0,60 a 0,90 e a produtividade de 0,5 a
de determinação superior (0,76). Ma et al. (2001)
5,0 t ha-1.
também observaram coe ciente de determinação
superior, de 0,80 no EF R5. Porém Holzapfel
(2007), avaliando a correlação entre o NDVI e a
GNDVI versus produtividade de massa seca
produtividade de grãos de canola, obteve coe cientes
de determinação, em diferentes EF, variando entre De acordo com os coe cientes de determinação
0,00 e 0,63, para o ano de 2005 e de 0,25 a 0,53 para (r2), observou-se melhor correlação entre a
o ano de 2006, mostrando que as relações variam de produtividade de MS e o GNDVI somente para
um ano para o outro no mesmo EF. o EF V9 (Tabela 2), assim como observado para
a produtividade de grãos (Tabela 1). No entanto,
Rudorff et al. (2003) avaliaram o NDVI na
os valores de correlação foram menores. Como,
estimativa da produtividade de milho e observaram
de certa forma, a produtividade de grãos depende
que não houve relação direta entre o NDVI e a
da produtividade de MS, era de se esperar que
produtividade de grãos para produtividades acima
ocorresse o mesmo comportamento. Santos Júnior
de 5 t ha-1. Estes autores atribuíram a limitação
et al. (2001) observaram, para a cultura da soja,
do uso do NDVI, na estimativa de produtividades
relação da produtividade de massa com o NDVI.
superiores a 5 t ha-1, ao elevado índice de área foliar
(IAF). Entretanto, isto deve ser visto com cuidado
Variáveis
Equações r2 p
Dependente Independente
GNDVI V5 y = 2,8638x + 3,6094 (*) 0,04 0,04
GNDVI V9 y = 12,485x – 1,0747 (**) 0,28 0,00
GNDVI R1 y = 7,8714x + 3,2841 0,04 ns
MS
GNDVI R3 y = 14,197x – 0,8292 0,02 ns
GNDVI R6 y = 20,688x – 5,3701 0,10 ns
GNDVI R8 y = 6,3323x + 7,0779 0,03 ns
ns: não-signi cativo. * e ** Signi cativo a 5 e 1% de probabilidade, respectivamente.
Fonte: Elaboração dos autores.
No entanto, Epstein et al. (2005), usando imagens de trigo, e obtiveram r2 de 0,74; essa correlação
de satélite para medição do NDVI e sua relação com foi conseguida apenas até a população de 1.000
a produtividade de massa de gramíneas, obtiveram plantas m-2, a partir dessa população, o NDVI não
coe ciente de determinação mais elevado (0,60) foi considerado um estimador da densidade. Ma et
que os obtidos nesse estudo. Brandão et al. (2007) al. (2001), estudando três densidades de plantas de
utilizaram imagens de satélite, provenientes do soja, em duas texturas de solo (arenosa e argilosa)
LANDSAT, para avaliarem a relação entre o NDVI e em dois anos consecutivos, concluíram que a
e a massa total, de quatro tipos de vegetação no densidade de plantas não apresenta relação com o
Estado do Ceará: vegetação semidensa, densa, NDVI. Já Ahmadi e Mollazade (2009), avaliando o
nativa e rala, e encontraram elevados coe cientes NDVI, obtido por imagens de satélite, encontraram
de determinação, respectivamente, 0,90, 0,99, 0,80 boa correlação para a densidade de plantas de soja
e 0,74. Lamb; Trotter; Schneider (2009) avaliaram (r2 = 0,87).
a massa de sorgo (Sorghum bicolor), usando sensor
ativo Crop Circle®, e obtiveram coe ciente de
determinação de 0,90. Hancock e Dougherty (2007), GNDVI versus altura das plantas
avaliando o NDVI, em alfafa, e sua relação com a Para a altura das plantas, os coe cientes de
MS produzida encontraram r2 de 0,58. determinação foram baixos (0,00005 a 0,06600)
Nos trabalhos anteriormente citados, para assim como observado para a densidade de plantas.
diversas culturas, encontrou-se r2 elevado, assim, Não houve relação entre o GNDVI e a altura das
pode-se a rmar que existe potencial para o uso plantas nos EF avaliados.
do NDVI na estimativa da produtividade de MS, Araújo (2004) avaliando o NDVI com
porém, no presente trabalho, para a cultura da soja, videogra a aérea na cultura da soja, obteve r2 de
isso não foi observado. 0,14 (p<0,01) o que demonstra a limitação do uso do
NDVI na estimativa da altura das plantas. Xavier et
al. (2006), ao avaliarem a altura de plantas de trigo,
GNDVI versus densidade de plantas
encontraram r2 de 0,68, na fase de espigamento da
Para todos os EF avaliados, os coe cientes de cultura. Já Motomiya; Molin; Chiavegato (2009),
determinação foram baixos (0,0036 a 0,0891) o avaliando a relação entre o NDVI e a altura de
que está de acordo com Araújo (2004), que obteve plantas de algodão, encontraram r2 de 0,74.
r2 de 0,0838 avaliando a relação entre a densidade
Isso pode estar relacionado ao espaçamento da
de plantas e o NDVI, obtidos por videogra a aérea
cultura, uma vez que para o algodão o espaçamento
para a cultura da soja.
entre linhas é maior e, assim, a cultura não fecha
De acordo com Mercante et al. (2009), baixos por completo, permitindo que a re ectância, maior
valores para a cultura da soja podem ocorrer em ou menor, seja caracterizada pelo NDVI, sendo
função do hábito de crescimento indeterminado e possível estabelecer relações entre este e a altura, o
da alta capacidade de compensar falhas no estande, que não ocorreu para a cultura da soja.
pela emissão de ramos laterais. Sendo assim, este
fator pode limitar a percepção das mudanças do
GNDVI, quando existe variação na densidade de Conclusões
plantas. O GNDVI na cultura da soja estabeleceu
Phillips et al. (2004) utilizaram os valores correlações signi cativas apenas com a
de NDVI, para estimar a densidade de plantas produtividade de grãos e a produtividade de MS no
523
Semina: Ciências Agrárias, Londrina, v. 34, n. 2, p. 517-526, mar./abr. 2013
Groff, E. C. et al.
524
Semina: Ciências Agrárias, Londrina, v. 34, n. 2, p. 517-526, mar./abr. 2013
Características agronômicas associadas com índices de vegetação medidos por sensores ativos de dossel...
LAMB, D. W.; TROTTER, M. G.; SCHNEIDER, NANNI, M. R.; DEMATTÊ, J. A. M.; FIORIO, P. R. Soil
D. A. Measuring and mapping crop vigour using an discrimination analysis by spectral response in the ground
active optical sensor in an ultra low level aircraft. In: level. Pesquisa Agropecuária Brasileira, Brasília, v. 39,
SYMPOSIUM ON PRECISION AGRICULTURE IN n. 10, p. 995-1006, 2004.
AUSTRALASIA, 13., 2009, Amidale. Proccedings...
PHILLIPS, S. B.; KEAHEY, D. A.; WARREN, J. G.;
Armidale: University of New England, 2009. p. 31-40.
MULLINS, G. L. Estimating winter wheat tiller density
MA, B. L.; DWYER, L. M.; COSTA, C.; COBER, E. using spectral re ectance sensors for early-spring,
R.; MORRISON, M. J. Early prediction of soybean variable-rate nitrogen applications. Agronomy Journal,
yield from canopy re ectance measurements. Agronomy Madison, v. 96, n. 3, p. 591-600, 2004.
Journal, Madison, v. 93, n. 6, p. 1227-1234, 2001.
RUDORFF, B. F. T.; MOREIRA, M. A.; ADAMI, M.;
MACEDO, M. A.; RUDORFF, B. F. T. Geotecnologias FREITAS, J. G. Resposta espectral do milho e sua
no seguro agrícola do milho safrinha. Agricultura em São relação com altos níveis de produtividade de grãos. In:
Paulo, São Paulo, v. 50, n. 1, p. 43-52, 2003. SIMPÓSIO BRASILEIRO DE SENSORIAMENTO
REMOTO, 11., Belo Horizonte. Anais... Belo Horizonte:
MACHADO, H. M. Determinação da biomassa de
INPE, 2003. p. 2597-2603.
cana-de-açúcar considerando a variação espacial de
dados espectrais do satélite LANDSAT 7 – ETM+. SANTOS JÚNIOR, R. F.; SANTOS, J. M.; RUDORFF,
2003. Dissertação (Mestrado em Engenharia Agrícola) B. F. T. MARCHIORATO, I. A. Espectrorradiometria de
– Faculdade de Engenharia Agrícola. Universidade campo no visível e infravermelho próximo para detecção
Estadual de Campinas, Campinas. de áreas infestadas com Meloidogyne javanica em soja.
In: CONGRESSO BRASILEIRO DE NEMATOLOGIA,
MERCANTE, E.; LAMPARELLI, R. A. C.; URIBE-
23., 2001, Marília. Anais... Marília: [s.n.], 2001. p. 62-62.
OPAZO, M. A.; ROCHA, J. V. Características espectrais
da soja ao longo do ciclo vegetativo com imagens landsat SHANAHAN, J. F.; SCHEPERS, J. S.; FRANCIS, D.
5/TM em área agrícola no Oeste do Paraná. Engenharia D.; VARVEL, G. E.; WILHELM, W.; TRINGE, J. M.;
Agrícola, Jaboticabal, v. 29, n. 2, p. 328-338, 2009. SCHLEMMER, M. R.; MAJOR, D. Use of Remote-
sensig imagery to estimate corn grain yield. Agronomy
MOREIRA, M. A. Fundamentos do sensoriamento
Faculty Publications, 2001. Paper 9. Disponível em:
remoto e metodologias de aplicação. 3. ed. São José dos
<http://digitalcommons.unl.edu/agronomyfacpub/9>.
Campos: Universidade Federal de Viçosa, 2007.
Acesso em: 15 jun. 2012.
MOTOMIYA, A. V. A.; MOLIN, J. P.; CHIAVEGATO,
SUGAWARA, L. M.; RUDORFF, B. F.; ADAMI, M.
E. J. Utilização de sensor óptico ativo para detectar
Viabilidade de uso de imagens Landsat em mapeamento
de ciência foliar de nitrogênio em algodoeiro. Revista
de área cultivada com soja no Estado do Paraná. Pesquisa
Brasileira Engenharia Agrícola Ambiental, Campina
Agropecuária Brasileira, Brasília, v. 43, n. 12, p. 1777-
Grande, v. 13, n. 2, p. 137-145, 2009.
1783, dez. 2008.
MOTTA, J. L. G.; FONTANA, D. C. Veri cação da
XAVIER, A. C.; RUDORFF, B. F. T.; MOREIRA, M.
acurácia da estimativa de área cultivada com soja
A.; ALVARENGA, B. S.; FREITAS, J. D.; SALOMON,
através de classi cação digital em imagens Landsat. In:
M. V. Hyperspectral eld re ectance measurements to
SIMPÓSIO BRASILEIRO DE SENSORIAMENTO
estimate wheat grain yield and plant height. Scientia
REMOTO, 5., 2001, Foz do Iguaçu. Anais... Foz do
Agricola, Piracicaba, v. 63, n. 2, p. 130-138, 2006.
Iguaçu: INPE, 2001. p. 123-129.
MYNENI, R. B.; KNYAZIKHIN, Y.; PRIVETTE, J.
L.; GLASSY, J.; TIAN, Y.; WANG, Y.; SONG, X.;
ZHANG, Y.; SMITH, G. R.; LOTSCH, A.; FRIEDL,
M.; MORISETTE, J. T.; VOTAVA, P.; NEMANI, R. R.;
RUNNING, S.W. Global products of vegetation leaf area
and fraction absorbed PAR from year one of MODIS
data. Remote Sensing of Environment , New York, v. 28,
n. 2, p. 214-231, 2002.
525
Semina: Ciências Agrárias, Londrina, v. 34, n. 2, p. 517-526, mar./abr. 2013