Segui tutoriais on-line para krigagem espacial com ambos geoR
e gstat
(e também automap
). Posso executar krigagem espacial e compreendo os principais conceitos por trás disso. Eu sei como construir um semivariograma, como ajustar um modelo a ele e como executar krigagem comum.
O que não entendo é como são determinados os pesos dos valores medidos circundantes. Eu sei que eles derivam do semivariograma e dependem da distância do local da previsão e do arranjo espacial dos pontos medidos. Mas como?
Alguém poderia, por favor, criar um modelo comum de krigagem (não bayesiana) com 3 pontos aleatórios medidos e 1 local de previsão? Seria esclarecedor.
prediction
spatial
weights
kriging
Pigna
fonte
fonte
Respostas:
Primeiro, descreverei a krigagem comum com três pontos matematicamente. Suponha que temos um campo aleatório intrinsecamente estacionário.
Kriging comum
Estamos tentando prever o valor usando os valores conhecidos Z = ( Z ( x 1 ) , Z ( x 2 ) , Z ( x 3 ) ) A previsão queremos é da forma Z ( x 0 ) = λ T Z onde λ = ( λ 1 , λ 2 , λ 3 )Z( x0 0) Z= ( Z( x1) , Z( x2) , Z( x3) ))
O comportamento preciso dos pesos é difícil de ver apenas a partir da equação, mas pode -se dizer muito grosso modo :
Você pode brincar usando a
clickppp
função spatstat :Aqui estão alguns exemplos
Pontos próximos entre si compartilharão os pesos
Ponto próximo "roubando" os pesos
É possível obter pesos negativos
Espero que isso lhe dê uma idéia de como os pesos funcionam.
fonte