MAPEAMENTO DO TEOR CLOROFILA TOTAL EM ÁREAS COM A CULTURA DO

MILHO A PARTIR DE IMAGENS MULTIESPECTRAIS TOMADAS POR AERONAVES

REMOTAMENTE PILOTADAS
1
George Deroco Martins 1, Denner Borges Rezende 2, Odair José Marques 3, Rafaela Souza Martins4,
Laura Cristina Moura Xavier5,Carlos Aberto Matias de Abreu Junior6,Rodrigo Bezerra de Araújo
Gallis7, Ricardo Luís Barbosa8, Alice Pedro Bom Paes9
1
Faculdade de Engenharia Civil (FECIV), km 1, LMG-746, Monte Carmelo-MG. E-mail: [email protected];
2
Instituto de Ciências Agrárias (ICIAG), km 1, LMG-746, Monte Carmelo-MG. E-mail: [email protected];
3
Instituto de Ciências Agrárias (ICIAG), km 1, LMG-746, Monte Carmelo-MG. E-mail: [email protected];
4
Instituto de Ciências Agrárias (ICIAG), km 1, LMG-746, Monte Carmelo-MG. E-mail: [email protected];
5
Instituto de Instituto de Geografia (IG), km 1, LMG-746, Monte Carmelo-MG. E-mail:[email protected]
6
Instituto de Geografia (IG), km 1, LMG-746, Monte Carmelo-MG . E-mail: [email protected];
7
Instituto de Geografia (IG), km 1, LMG-746, Monte Carmelo-MG . E-mail: [email protected];
8
Instituto de Geografia (IG), km 1, LMG-746, Monte Carmelo-MG . E-mail: [email protected];
9
Instituto de Geografia (IG), km 1, LMG-746, Monte Carmelo-MG . E-mail: [email protected].

previously limited the low space and time resolutions of the

RESUMO orbiting satellites. Thus, the objective of this work was to
evaluate the potential of a low cost airborne multispectral
Em atividades relacionadas à agricultura de precisão, camera (MAPIR 3) in estimating the total chlorophyll
sensoriamento remoto tornou-se uma ferramenta content of an experimental area cultivated with maize. In
imprescindível no monitoramento agrícola, principalmente this context, from a prediction model composed of the NDVI
por propiciar o acompanhamento e o desenvolvimento das index, it was possible to estimate the total chlorophyll
culturas agrícolas ao longo do tempo, identificação de content of the maize crop with an accuracy of over 87%.
áreas afetadas por fitopatógenos e na estimativa de
variáveis agronômicas. No cenário atual, a inserção de Key words — corn, chlorolphyll, RPA, multiespectral
imagens multiespectrais tomadas por aeronaves câmera, prediction model
remotamente pilotadas possibilitou desenvolver uma série
de aplicações em agricultura, as quais, antes estavam 1. INTRODUÇÃO
limitadas às baixas resoluções espaciais e temporais dos
satélites orbitais. Assim, objetivou-se com este trabalho Nos últimos anos, o uso do sensoriamento remoto
avaliar o potencial de uma câmera multiespectral tornou-se imprescindível no monitoramento agrícola,
aerotransportada de baixo custo (MAPIR 3) em estimar o principalmente por propiciar o acompanhamento e o
teor de clorofila total de uma área experimental cultivada desenvolvimento das culturas agrícolas ao longo do tempo,
com milho. Neste contexto, a partir de um modelo de previsão de safra, determinação de novas áreas de cultivo,
predição composto pelo índice NDVI foi possível estimar o determinação da quantidade e da localização da área
teor de clorofila total da cultura do milho com uma cultivada, entre outras possibilidades, como, uma forma
acurácia superior a 87%. alternativa para identificação de áreas afetadas por
fitopatógenos [1].
Palavras-chave —milho, clorofila, ARP,camera O monitoramento remoto em função da resposta
multiespectral, modelo de predição. espectral da vegetação em condições de estresse é
particularmente importante na agricultura de precisão, uma
ABSTRACT vez que os custos relacionados com a intensa amostragem
de material vegetal (folhas) e material do solo ou no manejo
In activities related to precision agriculture, remote sensing agrícola para toda área de cultivo podem ser minimizados
has become an essential tool in agricultural monitoring, pelo mapeamento prévio de variáveis agronômicas de
mainly for providing monitoring and development of the interesse.
crop over time, identification of agricultural areas affected A agricultura de precisão considera exatamente esse
by pathogens and in the estimation of agricultural variables. cenário: a presença, distribuição e intensidade de uma
In the present scenario, the insertion of multispectral condição específica que causa redução no rendimento da
images taken by remotely piloted aircraft made it possible to cultura agrícola, devem ser identificadas e receber
develop a series of applications in agriculture, which had tratamento diferenciado. Assim, o uso de dados

1338
radiométricos, como ferramentas de planejamento, controle experimentais da cultura de milho, sob tratamentos
e análise de dados agrícolas, prevê então, o monitoramento inseticidas para controle de Spodoptera frugiperda.
das variações espaciais e espectrais dos fatores que limitam
a produção, orientando no processo de tomada de decisão, 2. MÉTODOS
na aplicação de insumos em áreas específicas e no manejo
diferenciado das culturas no campo de produção [2]. Para estimar os TCTs na cultura do milho, esta pesquisa
A dimensão espacial de variáveis agrícolas está foi divida nas etapas: (1) montagem da área experimental:
relacionada diretamente com a capacidade de amostragem talhões de milho sob diferentes formas de manejo; (2)
dos sensores e, indiretamente, com o nível de aquisição de mensuração in situ do TCT; (3) aquisição e processamento
dados. Dados medidos in situ ou em laboratório, com campo de imagens multiespectrais; (4) análise da acurácia dos
de visada restrito, conseguem capturar condições específicas modelos de estimativa; (5) geração dos mapas.
da folha, enquanto que sensores imageadores, instalados em
plataformas remotas, determinam uma área de amostragem 2.1. Área de Estudo
maior na superfície terrestre. Por incorporar em um único
elemento de resolução, a resposta espectral de outros alvos, O experimento foi instalado, no período da safrinha de
como o solo, além daquela específica da planta ou plantas, 2018, no Campo Demonstrativo e Experimental da
imagens aéreas podem indicar a dispersão espacial das Universidade Federal de Uberlândia, Campus Monte
diferentes proporções de planta e solo, dentro da cultura Carmelo. A área de estudo, localiza-sena região do Alto
agrícola. Ou seja, a partir da caracterização espectral e Paranaíba, estado de Minas Gerais, apresenta altitude
espacial da cultura sadia e infectada por fitopatógenos, por aproximada de 908 m. O solo da área é caracterizado como
exemplo, dados multiescala de sensoriamento remoto LATOSSOLO Vermelho eutroférrico de textura argilosa
permitem identificar as áreas infectadas e inferir graus de [5].
sanidade da cultura [1].
As imagens multiespectrais do sensoriamento remoto
obtidos em multiescala, com vistas ao monitoramento e 2.2. Instalação e condução do experimento
mapeamento de culturas agrícolas encontram-se em
constante avanço. Isso ocorre devido à grande quantidade de O delineamento foi realizado em blocos com tratamentos
dados obtidos por intermédio de plataformas terrestres, atualizados, contendo dois híbridos de milho (um
convencional SHS4080® da Santa Helena Sementes e um
aeronaves e orbitais, disponibilizando produtos
transgênico Fórmula Viptera® da Syngenta Seeds) e seis
radiométricos oriundos de fontes diversas [3].
tratamentos inseticidas para controle de S. frugiperda
Atualmente, uma série de novos sensores multiespectrais
instalados em plataformas espaciais incluem bandas (lagarta-do-cartucho), com quatro repetições, totalizando em
espectrais específicas e apropriadas para monitorar 48 parcelas. Cada parcela tinha dimensões de 6,0 x 2,0 m
(12 m2), contendo quatro linhas de plantas, com
atividades agrícolas. Sensores como RapidEye e o
espaçamento de 0,5 m nas entrelinhas, sendo a área útil
WorldView 2 e 3 dispõem de uma banda espectral no
composta pelas duas linhas centrais com dimensões de 5,0 x
vermelho limítrofe, a qual é capaz de detectar variações
1,0 (5 m2).
sutis de clorofila nas folhas e fornecer informações sobre
estrutura e a condição sanitária da vegetação [4]. Com intuito de simular condições diversas do TCT ao
Quanto aos tipos de plataformas aéreas, as Aeronaves longo do talhão, o manejo e tratos culturais preconizados
para a cultura foram aplicados, variando apenas a aplicação
Remotamente Pilotadas (ARPs) apresentam-se como uma
de inseticidas.
ferramenta promissora em aplicações no setor agrícola e em
missões de reconhecimento, principalmente pela redução de
custo e do tamanho dos equipamentos colocados a bordo. 2.3 Mensuração da clorofila
Além do mais, a alternativa de solução baseada em
imagens tomadas de plataformas aéreas constitui uma O TCT foi mensurado em 60 pontos georreferenciados
com uma única campanha de campo no estádio V8 (oito
solução que pode ser viabilizada, caso seja baseada na
folhas da planta estão completamente expandidas).
integração de diferentes sensores instalados em plataformas
As medidas foram realizadas na última folha expandida
ARPs, as quais permitem maior flexibilidade na realização
da planta com um clorofilômetro SPAD 502 plus, da Konica
das tomadas de imagens principalmente no caso do
monitoramento de culturas agrícolas, sendo uma excelente Minolta, com precisão de ± 1,0 unidade SPAD (SPAD para
ferramenta para a Agricultura de Precisão. o valor entre 0,0 e 60,0 em temperatura/umidade normais).
Assim, a proposta deste trabalho foi avaliar o potencial
de imagens multiespectrais tomadas por ARPs em estimar 2.4 Aquisição e processamento de imagens
variáveis agronômicas. Para tanto, a variável agronômica multiespectrais
estimada foi o teor de clorofila total (TCT) em talhões
Concomitante as medidas do TCT, as imagens
multiespectrais foram adquiridas a partir de uma câmera

1339
multiespectral MAPIR 3, contendo as bandas do intervalo classes em intervalos iguais, apresentando a magnitude do
espectral médio do verde (550 nm), vermelho (660 nm) e TCT em ordem crescente de cinco unidades de valores de
infravermelho próximo (IVP) (850 nm). SPAD.
A câmera foi instalada na posição nadir em um ARP,
drone Phantom 4 Pro e, o voo foi realizado às 12 h com uma
altitude de 165 m e e um configuração da resolução espacial 3. RESULTADOS E DISCUSSÕES
das imagens em 12,5 cm. No momento da aquisição de
imagens, foram instaladas na superfície próxima do talhão, A Tabela 1 apresenta a correção obtida entre as bandas
placas de papel fotográfico nas cores: azul, verde, vermelho, da câmera, índice de vegetação e o TCT.
branco e preto para a posterior calibração radiométricas das
imagens adquiridas. Tabela 1 – Correlação entre bandas - índice de vegetação
O mosaico de imagens foi gerado no software Pixel4D. e o teor de clorofila total
Na etapa de fototriangulação, não foram utilizados pontos Banda Correlação p-valor
de apoio materializados no terreno, para formação do Verde (R550) 0,45 0,001
mosaico, apenas utilizou-se as coordenadas geográficas do Vermelho (R660) 0,40 0,006
centro perspectivo das imagens e os ângulos de orientação IVP (R850) 0,55 0
exterior proveniente de um sistema inercial embaracado na Índice
própria câmera MAPIR 3. A calibração radiométrica do NDVI 0,70 0
mosaico foi realizada de forma automática com o software
Mapir Camera Control. A partir da Tabela 1 observa-se que as bandas da câmara
e índice de vegetação apresentaram uma correlação
2.5 Definição de modelos predição para estimativa do significativa com o TCT, onde o maior p-valor foi de 0,006
teor de clorofila para a banda do vermelho. A maior correlação registrada foi
para o índice NDVI, 0,70, e a menor com a banda do
Para geração dos modelos empíricos, primeiramente, vermelho, 0,40.
foram avaliadas a correlação entre o TCT e as bandas A alta correlação entre o índice NDVI e o TCT pode ser
multiespectrais e o índice de vegetação NDVI [6]. atribuída à sensibilidade do índice com a quantidade de
O critério para escolha das bandas e dos índices de nitrogênio na folha. Isto se justifica ao fato de que o NDVI e
vegetação para composição dos modelos foi definido a partir as leituras do clorofilômetro se tornam mais sensíveis às
dos canais que apresentaram uma significativa e maior taxas de nitrogênio em cada estádio fenológico da planta.
correlação com o TCT. Deste modo, definiu-se como Constatações semelhantes a estes resultados observam-se
significativa à correlação superior a 50% e um nível de em [7], onde foi avaliada a correlação do NDVI com o TCT
significância de 5%. em todas as fases do desenvolvimento fenológico do milho.
Os modelos predição foram definidos a partir de Neste estudo, o autor constatou que correlação do índice
equações de regressão lineares simples entre valores com o nitrogênio estimado por medidas de SPAD aumentou
radiométricos medidos na imagem e a concentração do em estádios mais avançados de desenvolvimento da cultura,
TCT. Para tanto se considerou a radiância e o TCT de 50 possivelmente devido ao maior acúmulo de nitrogênio pela
amostras (definidas aleatoriamente) entre os 60 pontos de planta.
coletas contidos dentro da área (Figura 1). A Tabela 2 apresenta os modelos de predição do teor de
A validação da precisão dos modelos de predição foi clorofila total gerados, os RMSE/RMSE (%) e os
realizada estimando-se o TCT em 10 (definidos correspondentes coeficientes de determinação (R²).
aleatoriamente) dos 60 pontos amostrados in situ. A partir
dos dados de discrepância entre o TCT observado e Tabela 2 – Modelos de estimativa do teor de clorofila e
estimado, foi calculado o erro médio quadrático (RMSE) da análise da acurácia. (onde R representa banda da câmara
amostra. O RMSE e erro (RMSE(%)) foi calculado a partir MAPIR 3)
do modelo da equação 1 e 2. Banda Modelo RMSE/ RMSE (%) R²
Verde (R550) 3,9+0,0007* R550 6,18/13,1% 21
∑
√ (1) Vermelho (R660) 17,8+0,0014*R660 6,38/13,5% 15
IVP (R850) 21,1+0,0017*R850 6,01/13,0% 27
Índice
NDVI 34,5+225*NDVI 5,86/12,5% 46
Onde, xi e xmed representam o valor estimado e medido in
situ do TCT, respectivamente, e n é o número de amostras.
A partir da Tabela 2, observa-se que embora haja uma
Ao final, a partir do modelo que obteve o menor RMSE
diferença significativa entre o maior (46%) e o menor
foi gerada uma carta contendo o mapeamento do TCT sobre
coeficiente de determinação (21%), os índices de acerto dos
os talhões milho. O mapa coroplético foi divido em cinco

1340
 modelos de predição obteve pouca variação, onde o maior concentração do teor de clorofila total com acurácia superior
erro foi para o modelo composto pela banda do vermelho à 87%.
(13,5%) e o menor erro para o modelo composto pela banda 3-Com mapeamento do teor de clorofila total foi
do IVP (12,5%). possível constatar que o NDVI derivado de uma câmara de
O fato dos coeficientes determinísticos não terem sido baixo custo é extremamente sensível a pequenas variações
maiores do que 50% justificam-se por serem compostos por de clorofila em produtos radiométricos de alta resolução
uma única banda e pelo fator multiesca, ou seja, há uma espacial.
imprecisão radiométrica aliada à distância do sensor/alvo e
pela resolução espacial da imagem. Em [8], os modelos para 6. REFERÊNCIAS
predição para o TCT no milho apresentaram coeficientes de
determinação próximas a 98%, porém, neste os modelos [1] MARTINS, George Deroco; GALO, Maria de Lourdes Bueno
foram compostos pelo NDVI derivado de médias Trindade; VIEIRA, Bruno Sergio. Detecting and Mapping Root-
espectrorradiométricas in situ. Knot Nematode Infection in Coffee Crop Using Remote Sensing
A Figura 1 apresenta o mapeamento do TCT Measurements. Ieee Journal Of Selected Topics In Applied
Earth Observations And Remote Sensing, [s.l.], v. 10, n. 12,
representados em valores de SPAD gerados a partir do p.5395-5403, dez. 2017. Institute of Electrical and Electronics
modelo de predição composto pelo índice NDVI. Engineers (IEEE)..

[2] BECK, P. S. A et al., The feasibility of detecting trees

affected by Pine Wood Nematode using remote sensing. Europe
Comission. Joint Research Centre. Institute of Environment and
Sustainability, Italy. 2015. 34p.

[3] PONZONI, F.J.; SHIMABUKURO, Y.E. Sensoriamento

Remoto no Estudo da Vegetação. São José dos Campos, SP.
Editor Parêntese, 2015.

[4] OUMAR, Zakariyyaa; MUTANGA, Onisimo. Integrating

environmental variables and WorldView-2 image data to improve
the prediction and mapping of Thaumastocoris peregrinus (bronze
bug) damage in plantation forests. Isprs Journal Of
Photogrammetry And Remote Sensing, [s.l.], v. 87, p.39-46, jan.
2014.

[5] EMBRAPA. Empresa Brasileira de Pesquisa

Agropecuária. Sistema brasileiro de classificação de solos. 2ª.
Ed. Rio de Janeiro, Embrapa, 2006. 306p.
Figura 1. Carta do teor de clorofila total.
[6] ROUSE, J.W; HAAS, R.H; SCHELL, J.A; DEERING, DW,
onitoring vegetation system in the great plains with ERTS. In:
Na carta da Figura 1 observa-se que o modelo de
EARTH RESOURCES TECHNOLOGY SATELLITE
predição composto pelo NDVI foi sensível as variações de SYMPOSIUM, 3, 1973, Washington. Proceedings. Washington:
clorofila resultantes dos diferentes manejos inseticidas NASA, 1973.
aplicados sobre a área experimental.
Em alguns locais como na parte superior do talhão e na [7] POVH, Fabricio Pinheiro. Gestão da adubação nitrogenada em
porção inferior do talhão menor, nota-se alguns vazios milho utilizando sensoriamento remoto. Tese de Doutorado, [s.l.],
margeados por valores menores de SPAD (35-40), ou seja, p.1-108, 2011. Universidade de Sao Paulo Sistema Integrado de
locais em que não foi estimado o TCT. A esta ocorrência Bibliotecas - SIBiUSP.
justifica-se pela desfolha causada pela S. frugiperda.
[8] GASPAROTTO, Aline de Carvalho et al. Utilização do índice
de vegetação RNIR nas diferentes doses de nitrogênio em cultura
5. CONCLUSÕES agrícola. Anais Xvii Simpósio Brasileiro de Sensoriamento
Remoto - Sbsr, João Pessoa-pb, Brasil, 25 A 29 de Abril de
1-As bandas do visível apresentaram correlação 2015, Inpe, Pb, p.2435-2442, 2015.
significativa com índice de vegetação, porém com valores
não expressivos. Em contrapartida, a banda do
infravermelho e o índice NDVI apresentaram correlações
superiores a 0,50.
2- Mesmo com um coeficiente de determinação inferior
a 50%, a partir do índice NDVI foi possível estimar a

Powered by TCPDF (www.tcpdf.org)

1341