Geotecnologias em Agricultura

Prévia do material em texto

UNIVERSIDADE DE SÃO PAULO - ESCOLA SUPERIOR DE AGRICULTURA "LUIZ DE QUEIROZ"
DEPARTAMENTO DE ENGENHARIA DE BIOSSISTEMAS
DISCIPLINA LEB450 - GEOTECNOLOGIAS APLICADAS ÀS CIÊNCIAS AGRÁRIAS
AULA PRÁTICA 13 – SENSORES EM AGRICULTURA DE PRECISÃO
Prof. José Paulo Molin e Luiz Gustavo Góes
1. Condutividade elétrica aparente do solo (CEa) - CEa 0-30 cm e CEa 0-90 cm
● Adicionar os dados de CEa em formato shapefile.
● Exibir a tabela de atributos da ‘CEa’.
● Mudar a simbologia: Propriedade > Simbologia> Graduado > ‘0-30’ (CEa profundidade
a 30 cm) > 5 classes, Quartil > Aplicar.
● Mostrar as classes e valores gerados em cada classe no mapa.
Realizar a interpolação de pontos pelo método IDW e gerar o mapa de superfície da CEa.
● Inserir o shapefile ‘Limite_Talhão’ na pasta ‘Condutividade Elétrica’
● Para interpolação: na Caixa de Ferramentas de Processamento, digitar “Interpolação” e
selecionar ‘Interpolação IDW’
Camada vetorial: ‘CEa’> atributo de interpolação: ‘0-30’; selecionar o ‘+’ > Tipo: pontos
> Distância para coeficiente P: 0.5 (É a potência utilizada na equação de interpolação.
Potência igual a 2 seria o inverso do quadrado da distância, que é o mais usado
comumente) > Extensão > Calcular a partir da camada > ‘Limite_Talhão’ > tamanho do
pixel em X e Y: 10.00 m > Interpolado > Salvar no arquivo > pasta
‘Dados\Condutividade Elétrica’, nome: “CEa_interpolado” > Executar.
● Recortar a superfície gerada na interpolação, com base no limite do talhão. Para isto:
Selecionar ‘Raster’ > Extrair > Recortar pela camada de máscara > Camada de entrada:
‘CEa_interpolado’ > camada máscara: ‘Limite_Talhão’ > Recortado (mascara) > nome:
“CEa_mapa” > Pasta (Dados\Condutividade Elétrica’) > Executar.
● Mudar a simbologia: Selecionar o ‘CEa_mapa’ > Simbologia > Tipo de renderização:
banda simples falsa-cor > Interpolar: Linear > Modo: Quartil, 5 classes > Classificar > OK.
CEa é uma forma de estimar textura do solo, resumidamente quanto maior a água presente
no solo, maior a CEa, quem governa água no solo é a textura, portanto, quanto maior a
condutividade, infere-se que o teor de argila é maior.
2. Refletância Crop Circle (em café) - bandas (red, red edge e NIR) e os IV (NDVI e NDRE)
● Características do Crop-Circle e bandas necessárias para os cálculos de NDVI e NDRE.
● Selecionar ‘Gerenciador de Fonte de dados’ > Texto delimitado > Nome do arquivo
(‘dados\Reflectancia\ Crop Circle.csv’) > Formato de arquivo ‘Delimitadores
personalizados ‘selecionar ponto e vírgula’> Definição de geometria: Coordenadas de
pontos > Campo X: Longitude, Campo Y: Latitude > Geometria SRC: EPSG 32983 >
Adicionar.
Interpolar os pontos a partir do atributo: NDVI
● Inserir ‘limite_talhão’ > pasta ‘Dados/Refletancia’ > ‘limite_talhão’
● Para interpolação: na Caixa de Ferramentas de Processamento, digitar “Interpolação” e
selecionar ‘Interpolação IDW’ > Camada vetorial: ‘Crop Circle’> atributo de
interpolação: ‘NDVI’; selecionar o ‘+’ > Tipo: pontos > Distância para coeficiente P: 2.0
(É a potência utilizada na equação de interpolação. Potência igual a 2 seria o inverso do
quadrado da distância, que é o mais usado comumente.) > Extensão > Calcular a partir
da camada > ‘Limite_Talhão’ > tamanho do pixel em X e Y: 10.00 m > Interpolado >
Salvar no arquivo > pasta ‘Dados\Refletancia, nome: “Crop Circle_interpolado” >
Executar.
● Recortar as imagens para o ‘limite_talhão’: Raster > Extrair > Recorte raster pela
camada de máscara.
Obs. Por se tratar de dados brutos existem outliers, portanto, valores de NDVI fora da
faixa -1 a 1, a legenda apresenta valores a partir ou igual a 0,8111 pois se trata de uma
área de café, geralmente nos carreadores e estradas são onde apresentaram os
menores valores de NDVI
● ‘Crop Circle_interpolado: Propriedades > Simbologia > Banda Simples falsa-cor >
Interpolar: Método discreto > Modo: Quartil, 5 classes > OK.
3. Produtividade (milho) - produtividade (VRYIELDMAS) e umidade dos grãos (moisture)
● Inserir a camada ‘Produtividade_Milho_MT_2016_S680.shp’ da pasta
“Dados/Produtividade”.
● Definir projeção: EPSG: 4326 (WGS 84)
● Abrir tabela de atributos da camada ‘Produtividade_Milho_MT_2016_S680.shp’ >
alternar modo de edição > abrir calculadora de campo > Criar um novo campo >
Nome do campo de saída: ‘X’ > Tipo de campo de saída: ‘número decimal (real) >
comprimento: ‘10’ e Precisão ‘5’ > na expressão, digite $x > OK > faça o mesmo
para a coordenada Y > digitando na expressão $y> OK > alternar edição e salvar.
Visualização de dados de produtividade e umidade de grãos
● Botão direto do mouse na camada ‘Produtividade_Milho_MT_2016_S680.shp’ >
propriedades > simbologia > alterar de símbolo simples para graduado com valor
VRYIELDMAS (produtividade) > modo: quartil, classificar > aplicar
● Botão direto do mouse na camada ‘Produtividade_Milho_MT_2016_S680.shp’ >
propriedades > simbologia > alterar de símbolo simples para graduado com valor
Moisture (umidade do grão) > modo: quartil, classificar > aplicar
● Observar os padrões de comportamento dos pontos, principalmente nas
“cabeceiras” onde tem manobras, com dados de produtividade menores nas
bordaduras, esses dados contêm muitos outliers, são dados brutos e necessitam de
um processamento, como a filtragem, visando eliminar dados discrepantes.