MODELAGEM VARIOGRÁFICA. Estatística espacial Estatística clássica variáveis independentes sem...
32
MODELAGEM VARIOGRÁFICA
Transcript of MODELAGEM VARIOGRÁFICA. Estatística espacial Estatística clássica variáveis independentes sem...
MODELAGEM VARIOGRÁFICA
Estatística espacial
Estatística clássica variáveis independentes sem continuidade espacial
Estatística espacial valores associados à localização no espaço e/ou no tempo distribuição contínua dos valores
Geoestatística variáveis regionalizadas: fenômeno natural aspecto estrutural (determinístico) aspecto errático (casual) correlação espacial
Procedimentos em geoestatística análise exploratória dos dados calculo do variograma experimental modelagem krigagem: estimativa e interpolação simulação
Geoestatística
Metodologia para determinar se existe uma
autocorrelação espacial entre dois pontos amostrados
quantificar o efeito da localização espacial sobre a variabilidade amostral
Baseia-se no pressuposto que pontos mais próximos normalmente estão mais relacionados entre si do que pontos muito afastados.
A hipótse intrínsica
A hipótese intrínsica estabelece que a distribuição das diferenças entre pares de valores obtidos em locais de amostragem é a mesma em toda a área, dependendo apenas da distância e da orientação entre os pontos.
Em outras palavras, as diferenças devem ser consistentes, NÃO constantes por toda a área.
Variograma e krigagemGeoestatística: variograma e krigagem.
Variograma: descrição matemática da relação entre a variância de pares de observação (pontos de amostragem) e a distância que separa essas observações
Krigagem: a autocorrelação espacial pode ser usada para estimar valores em pontos não amostrados.
para a utilização do semivariograma as seguintes suposições básicas são requeridas:
a) as diferenças entre pares de valores de amostras são determinadas apenas pela orientação espacial relativa dessas amostras;
b) o interesse é enfocado apenas na média e na variância das diferenças, significando que esses dois parâmetros dependem unicamente da orientação (hipótese intrínseca);
c) por conveniência assume-se que os valores da área de interesse não apresentam tendência que possa afetar os resultados e assim a preocupação será apenas com a variância das diferenças entre valores das amostras.
Semivariograma mostra, pela análise estrutural, o comportamento espacial da variável regionalizada ou de seus resíduos, quando na presença de tendência:
tamanho da zona de influência em torno de uma amostra; toda amostra cuja distância ao ponto a ser estimado for menor ou igual ao alcance, fornece informações sobre o ponto
anisotropia, quando os semivariogramas se mostram diferentes para diferentes direções de linhas de amostragem;
continuidade, pela forma do variograma quando para h=0 (h) já apresenta algum valor (efeito pepita/nugget); pode ser atribuído à erros de medição ou ao fato de que os dados não foram coletados a intervalos suficientemente pequenos para mostrar o comportamento espacial subjacente do fenômeno em estudo.
Variograma
semivariograma experimental mínimo de 30 pares remoção de valores anômalos maior h, a metade da maior distância existente
entre os pontos. grau de casualidade dos dados, E = Co/C
E<0,15: componente aleatória pequena 0,15 < E < 0,30: componente aleatória significante E > 0,30: componente aleatória muito significativa
extremo do grau de casualidade é o modelo de pepita pura, onde não ocorre covariância entre os valores e a análise semivariográfica não se aplica
iniciar com semivariograma omnidirecional semivariograma teórico
modelagem: processo que envolve várias tentativas e no qual a experiência pesa muito
pode-se optar por um ajuste manual, mais sujeito à erros, ou com o auxílio de algoritmos
Ajuste do variograma experimental a um modelo variográfico teórico Comparação visual Técnicas de ajuste automático:
Método dos mínimos quadradosCritério AIC (Akaike Information
Criterion)Critério CressieCritério Variowin, etc.
Validação cruzada
Modelos variográficos
As funções matemáticas dos modelos devem permitir que a matriz de covariâncias, nele baseada, possa ser invertida, como ocorre na krigagem. Somente certos modelos podem ser usados.
Dados
1ª. Variável: cádmio (example.dat/GEOEAS)2ª. Variable: valor de cádmio+100As coordenadas X e Y são as mesmas para ambas as variáveis ID X Y Cd Cd1001 288 311 11.5 111.52 285.6 288 8.5 108.53 273.6 269 7 1074 280.8 249 10.7 110.75 273.6 231 11.2 111.26 276 206 11.6 111.67 285.6 182 7.2 107.28 288 164 5.7 105.79 292.8 137 5.2 105.210 278.4 119 7.2 107.211 360 315 3.9 103.912 355.2 291 9.5 109.513 367.2 272 8.9 108.914 367.2 250 11.5 111.515 352.8 226 10.7 110.716 350.4 203 8.3 108.3…60 ...
Distribuição dos pontos tem as mesmas coordenadas.Valores não.
11. 5
8. 5
7. 0
10. 7
11. 2
11. 6
7. 2
5. 7
5. 2
7. 2
3. 9
9. 5
8. 9
11. 5
10. 7
8. 3
6. 1
6. 7
6. 2
0. 0
5. 5
4. 0
7. 0
5. 3
11. 6
9. 0
14. 5
12. 1
0. 9
0. 0
3. 2
1. 2
1. 7 1. 2
7. 6
11. 6
8. 7
5. 8
3. 8
10. 4
10. 0
7. 1
4. 4
10. 4
1. 6
15. 0
3. 4
6. 8
10. 814. 9 9. 9 11. 6 6. 5 10. 1 11. 8 11. 0 16. 7 11. 6 6. 99. 9
300 350 400 450
150
200
250
300
111. 5
108. 5
107. 0
110. 7
111. 2
111. 6
107. 2
105. 7
105. 2
107. 2
103. 9
109. 5
108. 9
111. 5
110. 7
108. 3
106. 1
106. 7
106. 2
100. 0
105. 5
104. 0
107. 0
105. 3
111. 6
109. 0
114. 5
112. 1
100. 9
100. 0
103. 2
101. 2
101. 7 101. 2
107. 6
111. 6
108. 7
105. 8
103. 8
110. 4
110. 0
107. 1
104. 4
110. 4
101. 6
115. 0
103. 4
106. 8
110. 8114. 9 109. 9 111. 6 106. 5 110. 1 111. 8 111. 0 116. 7 111. 6 106. 9109. 9
300 350 400 450
150
200
250
300
Cd
Cd100
Determinação de h
Maior “h”: metade da maior distância entre pontos
E-W: 492-254= 238 N-S: 315-118 = 197Maior distância = (2382 + 1972 )1/2 = 309Maior “h” = 150
Extensão do menor “h”? Extensão do intervalo x
número de pares no intervalo
h = 2 (extensão = 75)
h = 75 (extensão = 2)
h = 15 (extensão = 10)
h = 10 (extensão = 15)
Nuvem variográficadireção NE-SW
h = 5 (extensão = 30)
Krigagem para a estimativa de um ponto (Xo)
1
)xx(
:
)xx(
)xx(
:
01...11
1)xx(...)xx()xx(
:::::
1)xx(...)xx()xx(
1)xx(...)xx()xx(
0n
02
01
n
2
1
nn2n1n
n22212
n12111
)xx()zzVar( 0i
n
1ii0
*0
2K
Esta equação pode ser expressa também em termos de covariância (h) = C(0) - C(h)
Cálculo da variância(desvio padrão) associada(o) ao valor obtido por estimativa
Modelagem da variável Cd pelo Variowin
Modelagem da variável Cd100 pelo Variowin
Krigagem da variável Cd e respectivo mapa de desvios padrão dos valores estimados pela krigagem
Modelo adotado: exponencial, de acordo com a indicação da qualidade do ajuste calculado pelo Variowin
3 0 0 3 5 0 4 0 0 4 5 0
1 5 0
2 0 0
2 5 0
3 0 0
- 0 . 5 0
0 . 5 0
1 . 5 0
2 . 5 0
3 . 5 0
4 . 5 0
5 . 5 0
6 . 5 0
7 . 5 0
8 . 5 0
9 . 5 0
1 0 . 5 0
1 1 . 5 0
1 2 . 5 0
1 3 . 5 0
1 4 . 5 0
3 0 0 3 5 0 4 0 0 4 5 0
1 5 0
2 0 0
2 5 0
3 0 0
1 . 6 5
1 . 7 5
1 . 8 5
1 . 9 5
2 . 0 5
2 . 1 5
2 . 2 5
2 . 3 5
2 . 4 5
2 . 5 5
2 . 6 5
2 . 7 5
2 . 8 5
2 . 9 5
3 . 0 5
C D
Krigagem da variável Cd100 e respectivo mapa de desvios padrão dos valores estimados pela krigagem
Modelo adotado: exponencial, de acordo com a indicação da qualidade do ajuste calculado pelo Variowin
300 350 400 450
150
200
250
300
99. 5
100. 5
101. 5
102. 5
103. 5
104. 5
105. 5
106. 5
107. 5
108. 5
109. 5
110. 5
111. 5
112. 5
113. 5
114. 5
300 350 400 450
150
200
250
300
1. 65
1. 75
1. 85
1. 95
2. 05
2. 15
2. 25
2. 35
2. 45
2. 55
2. 65
2. 75
2. 85
2. 95
3. 05
Cd100
Escolha do modelo linear (default do SURFER) para a variável Cd
0 20 40 60 80 100 120 140 160
Lag D istance
0
2
4
6
8
10
12
14
16
18
20
Var
iogr
am
D irection: 0.0 Tolerance: 90.0C olum n D : C d
Escolha do modelo linear (default do SURFER) para a variável Cd100
0 20 40 60 80 100 120 140 160
Lag D istance
0
2
4
6
8
10
12
14
16
18
20
Var
iogr
am
D irection: 0.0 Tolerance: 90.0C olum n E : C d100
Krigagem da variável Cd e respectivo mapa de desvios padrão dos valores estimados pela krigagemModelo adotado: linear (SURFER)
3 0 0 3 5 0 4 0 0 4 5 0
1 5 0
2 0 0
2 5 0
3 0 0
01234567891 01 11 21 31 41 51 61 7
3 0 0 3 5 0 4 0 0 4 5 0
1 5 0
2 0 0
2 5 0
3 0 0
0
0 . 5
1
1 . 5
2
2 . 5
3
3 . 5
4
4 . 5
5
5 . 5
6
6 . 5
7
C d : m o d e l o l i n e a r
Krigagem da variável Cd100 e respectivo mapa de desvios padrão dos valores estimados pela krigagemModelo adotado: linear (SURFER)
300 350 400 450
150
200
250
300
0
0. 5
1
1. 5
2
2. 5
3
3. 5
4
4. 5
5
5. 5
6
6. 5
7
300 350 400 450
150
200
250
300
99. 0100. 0101. 0102. 0103. 0104. 0105. 0106. 0107. 0108. 0109. 0110. 0111. 0112. 0113. 0114. 0115. 0116. 0117. 0
Cd100: modelo linear
CONCLUSÕES
O aspecto geral dos dois mapas para a variável Cd, resultantes da krigagem, é praticamente o mesmo.
Idem em relação à variável Cd100
No primeiro caso, tendo sido ajustado o melhor modelo (exponencial), o resultado, para ambas as variáveis, com relação ao intervalo de valores de desvios padrão é o mesmo (1,65-3,05), independentemente da escala dos valores em Cd e Cd100.
Quando, porem, é adotado um modelo que não reflete o comportamento espacial dos dados (modelo linear), os desvios padrão se apresentam com valores maiores (0-7), para ambas as situações.
Os menores valores de desvio-padrão estão, sempre, associados aos locais com pontos de amostragem.
Isso porque: O semivariograma mostra a medida do
grau de dependência espacial entre valores e é uma medida da variabilidade em relação à distância. A krigagem usa essas informações para encontrar os pesos ótimos a serem associados às amostras que irão estimar um ponto.
A variância da krigagem é independente dos valores dos pontos usados para obter os estimadores Zi* e mede somente a configuração espacial dos dados. Sendo a krigagem baseada apenas no variograma, que é global, os valores da variância independe dos valores locais dos pontos de amostragem.
