Mapa de Parâmetros de solo - Krigagem com Smart Map

O que é interpolação no Smart Map?

Interpolação é o processo de transformar dados pontuais em uma superfície contínua.

Exemplo clássico no agro:

Você coleta:

• pH do solo

• Fósforo (P)

• Potássio (K)

• Produtividade

• NDVI em pontos amostrais

Só que esses valores existem apenas onde você coletou.

A interpolação usa matemática + estatística espacial pra estimar os valores nas áreas onde você não coletou.

Como gerar a krigagem via Smart map.

1. Importação dos pontos amostrais

Inicialmente, foi carregada no QGIS a camada contendo os pontos de amostragem do talhão, oriundos da planilha de análise de solo.

No plugin Smart-Map: Decision Support System for Precision Agriculture, foram definidos:

• Input Layer: planilha de campo contendo os pontos georreferenciados

• Variável Z: Fósforo (P) em mg/dm³

• Sistema de Coordenadas: EPSG:32722 (WGS 84 / UTM Zona 22S)

• Opção Eliminate Outliers ativada para remoção de valores discrepantes

Após a importação, os dados de coordenadas X, Y e os teores de fósforo foram visualizados para conferência.

2. Definição da grade de interpolação

Na aba Grid, foi criada a malha de interpolação com:

• Pixel Size X = 5,0 m

• Pixel Size Y = 5,0 m

Em seguida, foi inserido o polígono do talhão (talhao_oficial) como limite da interpolação (Outline Polygon).

O Smart-Map apresentou:

• Distribuição espacial dos pontos

• Índice de Moran = 1,000

• p-valor = 0,003

Esse resultado indica forte dependência espacial dos dados de fósforo, validando o uso da geoestatística (krigagem).

3. Ajuste do Semivariograma

Na aba Interpolation → Ordinary Kriging, foi ajustado o modelo geoestatístico:

• Modelo: Exponential

• Maximum Distance: 450,000

• Lag (h): 99,999

  
  

Parâmetros do modelo:

Parâmetro	Valor
Nugget (Co)	0
Sill (Co+C)	1023,873
Range (A)	424,104
R²	0,924

O alto coeficiente de determinação (R²) indica bom ajuste do modelo aos dados.

4. Validação Cruzada

Foi realizada a validação cruzada para verificar a confiabilidade da interpolação:

• RMSE = 21,800

• R² = 0,409

A dispersão observada entre valores observados e preditos indica variabilidade espacial considerável no teor de fósforo do talhão.

5. Geração do mapa interpolado

Após o ajuste do semivariograma, foram definidos:

• Número de vizinhos: 16

• Raio de busca: 894,427 m

Em seguida, foi executada a interpolação (Interpolate), gerando o mapa contínuo de fósforo (mg/dm³).

O resultado apresentou:

• Zonas com baixos teores de P (vermelho)

• Zonas com teores médios (amarelo)

• Zonas com altos teores (verde)

Esse mapa evidencia a variabilidade espacial da fertilidade dentro do talhão.

O resultado obtido pode ser reclassificado a fim de classificarmos que área está com quantidade boa de fósforo e que área necessita de atenção de acordo com a porcentagem de argila (%) no solo conforme documentação da Embrapa.