Aula 5 - Algoritimos de Interpolação
Aula 5 - Algoritimos de Interpolação
Aula 5 - Algoritimos de Interpolação
• Curvatura mínima
• Regressão polinomial
• Krigagem
A diferença entre os diversos métodos está em como, a partir de valores conhecidos, podem ser
estimados os nós da rede regular. Além disso, a krigagem se destaca por levar em consideração a
variabilidade espacial dos dados.
Estimativa dos valores da grelha
Modelos
Seja uma situação hipotética (Fig. abaixo) em que 5 pontos foram amostrados para
teores de U3O8 e se pretende a partir deles estimar o valor no ponto A (Clark
1979).Aplicando, inicialmente, o método do inverso da distância, onde os pesos são
corrigidos para obter uma soma igual a um, é obtido o resultado :
Inverso do quadrado da distancia
Neste caso, para estimar um valor desconhecido a partir de valores conhecidos importa
apenas as distâncias euclidianas para a atribuição dos pesos e isso é válido para todos os
outros algoritmos, com exceção da krigagem. Na krigagem há necessidade de encontrar,
preliminarmente, a variabilidade espacial dos dados o que é conseguido por uma análise
variográfica.
b)b) Cálculo médio de uma variável regionalizada para um volume maior que o suporte
geométrico como, por exemplo, no cálculo do teor médio de uma jazida a partir de
informações obtidas de testemunhas de sondagens;
Em todas essas situações o método fornece, além dos valores estimados, o “erro”
associado a tal estimativa, o que o distingue dos demais algoritmos à disposição.
Krigagem
Em todas essas situações o método fornece, além dos valores estimados, o “erro”
associado a tal estimativa, o que o distingue dos demais algoritmos à disposição.
Krigagem
Na krigagem para a estimativa de A são necessários não apenas os valores das distâncias de cada
um dos 5 pontos em relação a A, mas também as distâncias entre os pontos. Não são as distâncias
euclidianas, porém, que são utilizadas, mas sim os valores resultantes da análise espacial pelo
variograma modelado. Isso pode ser obtido graficamente colocando os valores no semivariograma ou
analiticamente utilizando a fórmula de um modelo, no caso, esférico.
A Figura abaixo mostra as relações entre amostras (em preto) e entre as amostras e o ponto a ser
estimado (em vermelho). Assim a distância euclidiana entre o ponto 1, o mais próximo, e A é de
21,54 pés, a qual corresponde, pelo variograma modelado, ao valor 322,7 (Fig. abaixo). Esse valor
também pode ser encontrado segundo:
V* = p1 x1 + p2x2 + p3x3 + … + pnxn, onde “os pi são os pesos atributos a cada amostra i.”
É evidente que existe associado a esse estimador um erro ε=V-V* e que se, teoricamente, diversas estimativas forem feitas a
média de erros é zero. Se os erros, portanto, apresentarem valores próximos a zero, o estimador é de confiança e isso pode ser
verificado pela distribuição desses valores. A maneira mais simples de medir estatisticamente tal distribuição é via o desvio
padrão ou a variância. No caso em questão, porém, a variância não pode ser obtida porque não se conhece o valor real que se
esta estimando e, portanto, também não se sabe qual o erro associado.
Variância dos erros = = desvios ao quadrado em relação ao erro médio = média de (V-V*)2.
Para encontro da variância pode-se, porém, utilizar o semivariograma, em que são medidas as diferenças ao quadrado. Num
semivariograma, previamente calculado, dada uma distância h entre os pontos, pode-se estimar a variância simplesmente lendo
o valor no eixo dos γ’s e multiplicando-o por 2.
Estimativa de um ponto utilizando a Krigagem
Se a soma dos pesos for igual a 1 e não ocorrer tendência local dos valores, esse estimador é o melhor e não
tendencioso, pois a partir dos pesos atribuídos a cada amostra, minimiza a estimativa da variância.
Isso é obtido construindo-se um sistema de n equações com n incógnitas (λ1, λ2, λ3, …λn) e havendo a
restrição de que Σλi=1, passa-se a n+1 equações. Como se tem apenas n incógnitas desconhecidas, introduz-
se uma outra, também desconhecida, para balancear o sistema, ou seja, o chamado multiplicador de
Lagrange, μ.
Estimativa de um ponto utilizando a Krigagem
O objetivo da krigagem é procurar pelo conjunto ótimo de ponderadores de modo que a variância
do erro de estimativa seja a menor possível. Para isso, é organizado um sistema de equações com
n+1 incógnitas, para a estimativa de um ponto (So):
Estimativa de um ponto utilizando a Krigagem
Em notação matricial:
A matriz [Si,Si] contém os valores obtidos no semivariograma referentes às distâncias entre as amostras estimadoras; o vetor
[Si,S0] contém os valores obtidos no semivariograma referentes às distâncias entre cada amostra e o ponto (So) a ser estimado
e o vetor [λi] contém os ponderadores a serem calculados. Resolvido o sistema de equações, obtém-se os pesos λi e o
multiplicador de Lagrange, μ, segundo:
Para o ponto S0 a ser estimado, obtém-se uma combinação linear dos valores dos pontos vizinhos e respectivos pesos: So =
ΣλiSi
Para o cálculo da variância (σ²) associada ao valor So obtido por estimativa usa-se a expressão: