Geoestatística no R – Lição 22: Fenômeno Isotrópico

Neste tutorial, vamos ilustrar como calcular o variograma experimental assumindo que o fenômeno de interesse é isotrópico.

Um fenômeno é isotrópico quando a estrutura de correlação espacial é a mesma em todas as direções. Neste caso, a função variograma depende unicamente da magnitude do vetor, , e não da direção.

Usaremos o mesmo banco de dados, como na lição anterior, u, jura.pred, e a função variogram(). Mas serão aumentados os argumentos da função variogram() e estes serão justificado ao longo da lição.

O que é fenômeno isotrópico?

Quando a função variogram() não se altera com a direção, diz-se que o fenômeno é isotrópico e é representado em 2D geometricamente por um circulo, e em 3D, por uma esfera.

Variograma experimental (omnidirecional)

Na ausência de anisotropia na área, o variograma experimental é calculado em uma unica direção e é denominado variograma omnidirecional. A função variogram() calcula o variograma experimental da seguinte forma:

 

Usando as funções observadas em lições anteriores do pacote ggplot2, personalizamos o variograma omnidirecional.

 

Usando a função plot(), você pode adicionar propriedades como plot.number = TRUE para adicionar ao gráfico os pares de pontos.

 

Uma segunda maneira de explorar a correlação espacial é desenhar a nuvem variográfica, aumentando o comando cloud =TRUE na função variogram(), como segue:

 

 

 

 

Interessado por geoestatística? Este treinamento vai te ajudar a dominar as técnicas:

 

Série de tutoriais: Geoestatística no R – veja a lista

Lição 1 – Como instalar o R

Lição 2 – Objetos e estruturas

Lição 3 – Leitura e gravação de dados

Lição 4 – Estrutura gráfica em R

Lição 5 – Uso dos gráficos no R

Lição 6 – Paleta de cores

Lição 7 – Análises descritivas dos dados

Lição 8 – Distribuições discretas

Lição 9 – Distribuições contínuas

Lição 10 – Regressão linear

Lição 11 – ggplot2 (histograma)

Lição 12 – ggplot2 (boxplot)

Lição 13 – ggplot2 (scatterplot)

Lição 14 – ggplot2 (qq-plot)

Lição 15 – ggplot2 (série temporal)

Lição 16 – ggplot2 (diagrama de área)

Lição 17 – ggplot2 (diagrama circular)

Lição 18 – Pacote sp

Lição 19 – Pacote raster

Lição 20 – Amostragens

Lição 21 – Variograma

•Lição 22 – Fenômeno isotrópico

Lição  23 – Modelos teóricos de variogramas

Lição 24 – Fenômeno isotrópico ou anisotrópico?

Lição 25 – Ajuste do variograma experimental com anisotropia

Lição 26 – Ajuste do variograma experimental omnidirecional

Lição 27 – Pacote GeoR

Lição 28 – Ajuste do variograma experimental com GeoR

Lição 29 – eyefit e vario4 do pacote GeoR

Lição 30 – Design dos dados

Lição 31 – Função gridded

Lição 32 – Validação cruzada

Lição 33 – Krigagem ordinária