Influência da direção de semeadura do milho nas variáveis reflectância e índice de vegetação verde normalizado

Influence of the row direction of corn on the reflectance and green normalized difference vegetation index

METEOROLOGIA E CLIMATOLOGIA AGRÍCOLA

Influência da direção de semeadura do milho nas variáveis reflectância e índice de vegetação verde normalizado

Influence of the row direction of corn on the reflectance and green normalized difference vegetation index

Eduardo G. de Souza^I; Gracielle A. Orlando^II; Miguel A. Uribe-Opazo^III; Edson A. da Silva^IV; Peter Scharf^V

^ICCET/UNIOESTE. CP 711, CEP 85814-110, Cascavel, PR. Fone: (45) 220-3199. Email: [e-mail]

^IIUNIOESTE. E-mail: [e-mail]

^IIICCET/UNIOESTE. E-mail: [e-mail]

^IVCCET/UNIOESTE. E-mail: [e-mail]

^VUniversity of Missouri. 214 E Waters Hall, Columbia, MO 65211, USA. E-mail: scharfp@missouri

RESUMO

Com este trabalho, objetivou-se estudar a influência da direção de semeadura do milho (Zea Mays L.) na sua reflectância na banda 550 nm (r₅₅₀) e no Índice de Vegetação Verde da Diferença Normalizado - IVVDN, visando fornecer informações necessárias à viabilização de um sistema de aplicação localizada de nitrogênio (N) em tempo real. O sistema de semeadura utilizado foi o direto, administrando-se os nutrientes de acordo com análise do solo e com uma taxa de N de 160 kg ha^-1. As linhas de semeadura foram localizadas na direção Norte-Sul e Leste-Oeste. Conclui-se que o milho semeado na direção Leste-Oeste apresentou valores para reflectância r₅₅₀ e índice de vegetação IVVDN menores e menos dependente da hora do dia que o milho semeado na direção Norte-Sul. O IVVDN mostrou-se mais apropriado para uso em sistemas de aplicação de fertilizantes a taxas variáveis em tempo real utilizando-se técnicas de sensoriamento remoto, por se mostrar menos sensível á direção de semeadura.

Palavras-chave:Zea mays L., IVVDN, sensoriamento remoto

ABSTRACT

The objective of this work was to study the influence of the row direction of corn (Zea mays) on its spectral properties measured in the field, including the green normalized difference vegetation index (Green NDVI). Spectral properties of corn are known to be related to the nitrogen (N) status of the corn and may form the basis for real-time site-specific rates of N application. Row direction effects on spectral properties may be substantial before canopy closure. If so, it will be necessary to account for these effects in order to accurately predict the N need of corn from spectral measurements. Corn was planted directly, with no tillage, in North-South (NS) and East-West (EW) oriented rows. Nitrogen was applied as ammonium nitrate at a rate of 160 kg N ha^-1. Reflectance was measured with spectral radiometers placed 30 cm above individual plants. Reflectance and green NDVI were bly dependent on time of day for corn in NS rows, and much less so for corn in EW rows. Reflectance and green NDVI were lower for corn in EW rows than corn in NS rows, but lack of replication makes it difficult to be certain that this was due to row direction. Green NDVI was less sensitive than reflectance to row direction and therefore may be more appropriate for use in a real-time variable-rate nitrogen application system.

Key words:Zea mays L., NDVI, remote sensing

INTRODUÇÃO

A técnica de Agricultura de Precisão (AP) consiste na aplicação localizada dos insumos agrícolas em função da necessidade específica local e se fundamenta na existência da variabilidade espacial dos fatores produtivos e, portanto, da própria quantidade produzida pela cultura (Balastreire et al., 1997). O desenvolvimento de tecnologias de aplicação de fertilizantes, a taxas variáveis, mostrou a necessidade de melhores métodos para se estimar as variabilidades dos fatores produtivos dentro dos talhões.

A medição da reflectância espectral é considerada promissora abordagem não invasiva para o sensoriamento instantâneo da deficiência de N em milho (Tumbo et al., 2000). O problema é que esta característica das plantas não é bem conhecida; ela possibilita a prática de aplicação de fertilizantes a taxas variáveis em tempo real e é influenciada por fatores como ângulo zenital, estádio vegetativo da planta e nebulosidade do dia, entre outros.

Segundo Wood et al. (1993) sistemas tecnificados de produção de milho requerem grande quantidade de N, mas o excesso à necessidade nutricional da planta pode lixiviar e atingir o lençol freático, contribuindo para a eutrofização de lagos e fluxos de água. O aumento da concentração de nitratos no solo e na água tem exigido melhores técnicas de administração de nitrogênio, particularmente em áreas de produção de milho.

A clorofila é responsável pela coloração verde das folhas e tem mostrado alta correlação com a concentração de N das folhas, especialmente quando ele é deficiente deste nas folhas (Wolfe et al., 1988; Lohry, 1989; Schepers et al., 1992 citados por Blackmer et al., 1994). Thomas & Oerther (1972) citados por Blackmer et al. (1994), mostraram que o nível de N nas folhas poderia ser rapidamente estimado pela quantidade de luz refletida no comprimento de onda de 550 nm. Neste comprimento de onda há menos absorção de energia que na região do azul (de 400 a 500 nm) e vermelho (de 600 a 700 nm).

Davis (1957) mostrou que a reflectância das plantas oscila com a elevação do sol. Seus valores variaram de 22% ao meio dia para aproximadamente 43% ao amanhecer e 48% no pôr-do-sol. Jackson et al. (1979) coordenaram um estudo sobre a dependência da reflectância espectral do trigo com a configuração da cultura (ou seja, a direção das linhas de semeadura) e a posição do sol. O trigo foi semeado em três disposições: em uma configuração ele foi semeado com um denso estande sem direção preferencial; nas outras duas configurações, o trigo foi semeado em linhas, espaçadas 30 cm, uma na direção Norte-Sul (NS) e outra na direção Leste-Oeste (LO). Eles calcularam as reflectâncias para as bandas verde (500-600 nm), vermelho (600-700 nm) e infravermelho próximo (700-800 nm), entre outras. As medidas foram efetuadas quando as plantas estavam com aproximadamente 35 cm de altura. A influência da posição do sol na reflectância foi pequena para os solos nus e para o estande denso; no entanto, para as configurações em linha a reflectância mostrou-se altamente dependente do ângulo zenital e direção de semeadura. Os autores afirmaram, ainda, que a dependência da reflectância na banda visível ao ângulo zenital pode ser explicada com base na alta absorção das folhas verdes.

Em culturas em crescimento, as plantas absorvem a maior parte da radiação visível e sombreiam diferentes proporções do solo e de outras partes da planta, dependendo do ângulo zenital, da direção de semeadura e da altura da planta. Quando a direção de semeadura é a NS e se tem um baixo índice de área foliar, a maior parte do solo é sombreada pela manhã e a tarde, e menos sombreado ao meio dia solar; assim, a reflectância do dossel em qualquer banda atinge seu mínimo com o sol do meio dia. Para linhas de semeadura LO a fração do solo iluminada muda menos que para a direção NS e a reflectância apresenta menores mudanças. A banda do vermelho mostrou-se mais sensível que a banda verde, em relação ao ângulo zenital, já a banda do infravermelho, por sua alta transmitância através das folhas verdes, precisa de três a quatro vezes mais camadas de folhas para alcançar a mesma reflectância da luz visível. Para um alto ângulo zenital a luz penetra na cultura de tal forma que muitas folhas são atingidas, logo a reflectância é alta. Perto do meio dia solar a quantidade de folhas atingida diretamente é menor, ocasionando uma reflectância também menor no infravermelho. Verificou-se, ainda, que a direção de semeadura tem menor efeito na reflectância do infravermelho que na banda visível. O decréscimo da reflectância no infravermelho e o aumento na banda visível com o aumento da elevação solar indicam que índices que utilizem a relação da reflectância no infravermelho com a reflectância no verde ou no vermelho, deverão ser acentuadamente influenciados pelo ângulo zenital. Os autores concluíram que a configuração da cultura é um dos fatores mais determinantes na reflectância espectral do trigo e que os índices de vegetação a serem utilizados em estudos temporais deverão levar em consideração a elevação solar. Conse-qüentemente, a elevação solar ou o ângulo zenital deveriam ser especificados para o desenvolvimento de modelos apropriados de correção a serem utilizados no processamento de imagens orbitais e sub-orbitais adquiridas em diferentes datas.

Índices de Vegetação Espectrais (IVE) são úteis para minimizar variações devido a fatores estranhos e maximizar a sensividade à variável de interesse, neste caso a concentração de clorofila na folha do milho (Daughtry et al., 2000). Esses índices são determinados a partir da baixa reflectância nos comprimentos de onda visíveis e da alta reflectância nos comprimentos de onda infravermelho-próximo que são característicos da assinatura espectral de uma vegetação verde. O número relatado de IVEs é elevado e eles se agrupam em três categorias: (1) índices intrínsecos, (2) índices ajustados ao solo, e (3) índices ajustados à atmosfera. Daughtry et al. (2000) concluíram que alguns IVE, como a relação infravermelho próximo/vermelho, minimizaram as contribuições da reflectância de fundo, enquanto outros, como a relação entre o infraver-melho próximo/verde responderam mais à concentração de clorofila da folha. Tumbo et al. (2000) encontraram uma forte correlação entre a relação de infravermelhos-próximo/verde e concentração de clorofila no estádio vegetativo V6 do milho, sob irradiância solar constante.

Objetivou-se, com este trabalho, estudar a influência da direção da semeadura do milho na reflectância na banda 550 nm e no índice verde de vegetação da diferença normalizada - IVVDN. O conhecimento da magnitude de influência permitirá encontrar-se correções matemáticas que possibilitarão a conversão da reflectância na correta taxa necessária de N.

MATERIAL E MÉTODOS

Os dados foram levantados em laboratórios e campos experimentais da Universidade do Missouri, Columbia, MO, EUA. O milho foi instalado através da técnica de semeadura direta, administrando-se os nutrientes de acordo com a análise do solo e a uma taxa de N de 160 kg ha^-1. As coordenadas geográficas locais são: Latitude 38° 54' N e Longitude 92° 12' W. As linhas de semeadura foram localizadas na direção Norte-Sul (NS) e Leste-Oeste (LO). O espaçamento entre plantas foi de 25 cm e entre linhas de 80 cm, proporcionando uma população aproximada de 50.000 plantas ha^-1.

As medidas foram tomadas no período de 03/09/2002 a 05/09/2002, estando o milho no estádio vegetativo V6. Quatro radiômetros foram montados a 25 cm acima do ápice das plantas e posicionados verticalmente para baixo, oferecendo um campo de visão de 28°, o que representa 12,5 cm de diâmetro de visão. Devido a trajetória do sol no mês de setembro e o adequado espaçamento entre radiômetros, evitaram-se sombras sobre as plantas que estavam sendo monitoradas. Os radiômetros 1 e 5 foram posicionados sobre as linhas orientadas na direção LO e os radiômetros 2 e 3 sobre as linhas orientadas na direção NS. Os radiômetros foram calibrados um contra os outros visando-se que diferenças nas reflectâncias calculadas não fossem originadas de diferenças entre os radiômetros mas, sim, devido as plantas sobre os seus campos de visão, tais como ao estado nutricional da planta e da geometria da radiação. Um radiômetro de campo de visão de 180° (radiômetro 4) foi posicionado voltado verticalmente para cima, com vistas a se estimar a irradiação global. Os dados espectrais da cultura do milho foram obtidos em intervalos de dois segundos e, posteriormente, agrupados em intervalos de dois minutos com vistas a se reduzir o volume de dados. Nos três dias de medição (03, 04 e 05/09/2002) levantaram-se dados das 8:00 às 18:00 h. A resposta espectral foi medida em duas bandas, cada uma com 10 nm de largura, centradas em 550 nm (verde) e 800 nm (infravermelho próximo). A reflectância do milho foi calculada por:

(1)

em que:

ri - reflectância na banda i (550 e 800 nm)

Rad(ji) - radiação medida no radiômetro j (1, 2, 3 e 5) e na banda i (550, 800 nm)

Rad4(i) - radiação medida no radiômetro 4 (voltado vertical-mente para cima) na banda i (550 e 800 nm)

Com base nos resultados reportados por Daughtry et al. (2000), Tumbo et al. (2000) e Souza et al. (2002), selecionou-se o índice de vegetação:

(2)

em que:

IVVDN - Índice verde de vegetação da diferença normalizada (Green Normalized Difference Vegetation Index, Green NDVI)

r₈₀₀ - radiação na banda 800 nm

r₅₅₀ - radiação na banda 550 nm

A análise estatística dos dados se constituiu de uma análise exploratória seguida de modelagem. Na elaboração dos modelos empregaram-se três métodos de seleção automática de variáveis: stepwise forward, stepwise backward e melhor subconjunto de variáveis. O coeficiente de determinação múltipla ajustado (R² ajustado) foi utilizado para mostrar os melhores modelos. Utilizou-se a média entre os dados dos radiômetros 2 e 3 para a direção NS e, para a direção LO, somente os dados do radiômetro 1, visto que o radiômetro 5 apresentou reflectâncias marcadamente menores quando comparadas aos outros, o que não era esperado, visto que todas as plantas estavam no mesmo estádio vegetativo e receberam a mesma quantidade de nutrientes. Como já comentado, este fato deve ter ocorrido por mau posicionamento do radiômetro sobre a planta.

Para elaborar os modelos da reflectância espectral do milho (r) e do índice de vegetação verde normalizado (IVVDN) selecionaram-se as variáveis: ângulo zenital (y, graus), hora do dia (t, hora) e direção de semeadura (dir), variável binária (0 para a direção LO e 1 para a direção NS). Essas variáveis foram selecionadas por se correlacionarem com a posição de Sol e o grau de iluminação da cultura, fatores influentes na reflectância (Davis, 1957; Jackson et al.; 1979; e Tumbo et al., 2000). Visando reduzir a esperada multicolinearidade pelo uso de modelos polinomiais, as variáveis quantitativas independentes foram transformadas para ângulo zenital centrado (y_c) e hora do dia centrada (t_c), através de

(3)

(4)

em que as constantes de centralização foram escolhidas para propiciar uma simetria dos dados, isto é, 45, como a metade do máximo ângulo zenital possível, e 13 como a hora aproximada do meio solar. Após uma análise exploratória dos dados optou-se por ajustar os dados ao seguinte modelo

(5)

em que:

y - variável resposta (r₅₅₀ ou IVVDN)

a, b, c, d, e, f, g, h, i - parâmetros do modelo a serem estimados

e - erro aleatório

RESULTADOS E DISCUSSÃO

A Figura 1 apresenta a dispersão dos dados da reflectância na banda 550 nm obtida com os quatro radiômetros para os três dias, em função do tempo. Verifica-se que, com exceção do radiômetro 5, os outros três radiômetros apresentaram resultados similares no período da tarde. Comparando-se os quatro radiômetros constata-se que, diferentemente dos radiômetros 1 e 5 (LO), as reflectâncias associadas aos radiômetros 2 e 3 (NS) não apresentaram visualmente diferenças entre si. Os radiômetros posicionados sobre milho semeado na direção LO mostraram reflectâncias médias maiores e menos dependentes da hora do dia que os radiômetros da direção NS, confirmando o encontrado por Jackson et al. (1979) de que a direção de semeadura do milho influencia na reflectância apresentada pelas plantas. A influência da hora do dia na reflectância foi maior no milho semeado na direção NS, pelo fato de que a irradiação atingiu diretamente as folhas pela manhã e principalmente elas, mas com o decorrer do dia a parcela dessa irradiação que atinge o solo foi aumentando, alcançando o seu máximo perto do meio dia. A influência da hora do dia na reflectância foi menor no milho semeado na direção LO, visto que a irradiação atinge igualmente solo e planta, praticamente todo o dia. Ressalte-se que é esperada uma diminuição da influência da direção de semeadura a medida em que a cultura cresce e maior parte do solo é coberta pelo dossel.

O comportamento da dispersão dos dados da reflectância na banda 550 nm para o período de 03 a 05/09/2002, para os quatro radiômetros e em função do ângulo zenital está apresentado na Figura 2. A reflectância aumentou com o ângulo zenital, alcançando o máximo no início e no fim do dia. Além disto, a reflectância do milho semeado na direção LO se manteve mais constante com o ângulo zenital que o milho semeado na direção NS.

Na Tabela 1 tem-se a análise estatística descritiva dos dados de reflectância para os dias 03, 04 e 05/09/2002. Como já visto, o milho sob o radiômetro 5 apresenta, em todos os dias, menores reflectâncias quando comparadas com as dos outros radiômetros, o que não era esperado, visto que todas as plantas estavam no mesmo estádio vegetativo e receberam a mesma quantidade de nutrientes. Verifica-se também que os resultados associados aos radiômetros 1 e 5 (direção LO) apresentaram menor variabilidade, visto que mostraram os menores coeficien-tes de variação (CV). Percebe-se também que a direção de semeadura influenciou na reflectância, sendo que as plantas semeadas na direção NS indicaram, em média, valores maiores de reflectância (de 7,3 a 10,1%), quando comparadas com as semeadas na direção LO (de 3,7 a 6,9%); isto evidencia a interação entre a geometria da radiação e a arquitetura das plantas.

Outra forma de se avaliar a radiação refletida por uma cultura é a utilização de índices de vegetação. A Figura 3 apresenta a dispersão dos dados do índice verde de vegetação da diferença normalizada (IVVDN) referente aos quatro radiômetros para os três dias, em função do tempo. Comparando-se as Figuras 1 e 3, observa-se que o IVVDN apresentou menor dependência com a hora do dia que a reflectância na banda 550 nm. Este fato contraria Jackson et al. (1979) que encontraram acentuada influência da hora do dia em índices que usam a relação de infravermelho-próximo/verde, provavelmente por eles terem trabalhado com cultura do trigo. Com exceção do índice IVVDN 5 (radiômetro 5), os índices de vegetação da direção LO e NS apresentaram-se mais agrupados que as respectivas reflectâncias. Por outro lado, os radiômetros 2 e 3 (direção NS), assim como para a reflectância a 550 nm, mostraram também maiores valores de IVVDN quando comparados com os radiômetros 1 e 5 (direção LO). Os baixos valores de índice de vegetação observados pelo radiômetro 5 reiteram a suspeita de que houve problemas de posicionamento deste radiômetro sobre a planta.

A Tabela 2 apresenta a análise estatística descritiva para os valores do índice verde de vegetação da diferença normalizada. Observa-se que os radiômetros posicionados sobre milho semeado na direção NS apresentaram maiores valores de índice de vegetação e, ainda que as variabilidades dos IVVDN foram menores que as da reflectância r₅₅₀, visto que os coeficientes de variação também foram muito menores. Percebe-se que a direção de semeadura influenciou bem menos no IVVDN que na reflectância.

Desenvolveram-se modelos através do método de mínimos quadrados para descrever a dependência das variáveis respostas r₅₅₀ e IVVDN em relação às variáveis ângulo zenital centrado (y_c), hora do dia centrada (t_c) e direção de semeadura (dir). A explicação da variabilidade das variáveis dependentes através dos modelos, foi maior que 86%, visto que o R² encontrado foi 0,866 para o índice de vegetação IVVDN e de 0,901 para a reflectância r₅₅₀ (Tabela 3). Ressalte-se que a inclusão da variável binária dir (0 para a direção LO, 1 para a direção NS) possibilitou a construção de um único modelo para ambas as direções, para cada variável resposta.

A Figura 4 apresenta os intervalos de predição para a variável reflectância a nível de 95% de confiança. Os intervalos de predição, como era de se esperar, abrangem a maioria dos valores medidos. Vê-se, também, que os intervalos das direções NS e LO se sobrepõem entre 11:00 e 18:00 h. O intervalo da direção NS é muito mais largo que o LO, devido principalmente a pontos discrepantes.

Na Figura 5 observam-se os intervalos de predição a nível de 95% de confiança para o índice verde de vegetação da diferença normalizada. Os intervalos de predição são muito mais próximos que os encontrados para a variável r₅₅₀. Enquanto para a variável r₅₅₀ a direção de semeadura influencia na interação entre a variável dependente r₅₅₀ e a hora do dia, para a variável IVVDN, a direção de semeadura contribui principalmente com um valor aditivo à variável dependente r₅₅₀.

CONCLUSÕES

1. O direcionamento da semeadura do milho em relação à trajetória do Sol, influenciou a reflectância r₅₅₀ e o índice de vegetação IVVDN.

2. O milho semeado na direção Leste-Oeste apresentou valores de reflectância r₅₅₀ e índice de vegetação IVVDN menores e menos dependentes da hora do dia que o milho semeado na direção Norte-Sul.

3. A direção de semeadura influenciou de forma diferente as variáveis analisadas. Na variável r₅₅₀ a direção de semeadura resultou em uma interação significativa entre a variável dependente r₅₅₀ e a hora do dia. Para a variável IVVDN a direção de semeadura contribui principalmente com um valor aditivo à variável dependente IVVDN.

4. O índice verde de vegetação da diferença normalizada mostrou-se mais apropriado para uso em sistemas de aplicação de fertilizantes a taxas variáveis em tempo real, utilizando-se técnicas de sensoreamento remoto por se mostrar menos sensível à direção de semeadura.

AGRADECIMENTOS

Os autores agradecem o apoio financeiro propiciado pela University of Missouri, Columbia-MO, EUA, e pelo Conselho Nacional de Desenvolvimento Científico e Tecnológico CNPq. Agradecem, também, o apoio técnico disponibilizado por Larry Mueller e Bettina Coggeshall, nos trabalhos de campo.

LITERATURA CITADA

Protocolo 104 - 16/6/2003 – Aprovado em 6/2/2004

• Balastreire, L.A.; Elias, A.I.; Amaral, J.R. Agricultura de precisão: Mapeamento da produtividade da cultura de milho. Revista de Engenharia Rural, Piracicaba, v.8, p.97-111, 1997.
• Blackmer, T.M.; Schepers, J.S.; Varvel, G.E. Light reflectance compared with other nitrogen stress measurement in corn leaves. Agronomy Journal, Madison, v.86. p.934-938, 1994.
• Daughtry, C.S.T.; Walthall, C.L.; Kim, M.S.; Brown de Colstoun, E.; Mcmurtrey III, J.E. Estimating corn leaf chlorophyll concentration from leaf and canopy reflectance. Remote Sensing of Environment, New York, p.229-239, 2000.
• Davis, P.A. Exploring the atmosphere's first mile. Lettau, H.H.; Davidson, B. (eds.) New York: Proceedings of the Great Plains turbulence field program, Symposium Publications Division, Pergamon Press, 1957. p.377-383.
• Jackson, R.D.; Pinter, P.J.; Idso, S.B.; Reginato, R.J. Wheat spectral reflectance: Interactions between crop configuration, sun elevation, and azimuth angle. Applied Optics, v.18, n.3, p.730-33, 1979.
• Souza, E.G.; Scharf, P.; Sudduth, K. How corn reflectance and vegetation indices are influenced by sun position and atmospheric conditions. In: 2002 Mid-Central ASAE Meeting. Proceedings.... St. Joseph: ASAE. 2002. ASAE Paper No. MC02-06
• Tumbo, D.S.; Wagner, D.G.; Heinenmann, P.H. Real-world hyperspectral characteristics of corn plants under different chlorophyll levels. In: 2000 ASAE Annual International Meeting. Proceedings... Mineapolis: ASAE. 2000. ASAE Paper No. 001093
• Wood, C.W.; Reeves, D.W.; Himelrick, D.G. Relationships between chlorophyll meter readings and leaf chlorophyll concentration, N status, and crop yield: A review. Proceedings.... New Zealand: Agronomy Society., 1993. p.1-9.

Datas de Publicação

• Publicação nesta coleção
  28 Set 2004
• Data do Fascículo
  Abr 2004

Histórico

• Aceito
  06 Fev 2004
• Recebido
  16 Jun 2003