´ Uvod Vizualizace dat Variogramy Ordinary Kriging Geostatistika v R projektu 18. studentsk´ a konference o matematice Martin Dzurov Krist´ yna Kitzbergerov´ a Lucie ˇ Sindel´ aˇ rov´ a ˇ Cesk´ e vysok´ e uˇ cen´ ı technick´ e v Praze Fakulta stavebn´ ı Katedra matematiky 16. – 18. ˇ cerven 2010 Martin Dzurov, Krist´ yna Kitzbergerov´ a, Lucie ˇ Sindel´ aˇ rov´ a Geostatistika v R projektu
This document is posted to help you gain knowledge. Please leave a comment to let me know what you think about it! Share it to your friends and learn new things together.
Transcript
Uvod Vizualizace dat Variogramy Ordinary Kriging
Geostatistika v R projektu18. studentska konference o matematice
Martin Dzurov Kristyna Kitzbergerova Lucie Sindelarova
Ceske vysoke ucenı technicke v PrazeFakulta stavebnı
Katedra matematiky
16. – 18. cerven 2010
Martin Dzurov, Kristyna Kitzbergerova, Lucie Sindelarova Geostatistika v R projektu
4 Ordinary KrigingTeoreticky uvod do problematikyVlastnı zpracovanı dat
Martin Dzurov, Kristyna Kitzbergerova, Lucie Sindelarova Geostatistika v R projektu
Uvod Vizualizace dat Variogramy Ordinary Kriging Cıl projektu R geoR
Uvod
Cıl projektu
Cılem projektu je vytvorenı ucebnıho textu, ktery ma priblızitproblematiku geostatistiky, predevsım vyznam interpolacnı metodykriging.
modelova data: nadmorske vysky pro CR v merıtku1 : 500 000
zpracovano v prostredı R
pouzit balıcek geoR
Martin Dzurov, Kristyna Kitzbergerova, Lucie Sindelarova Geostatistika v R projektu
Uvod Vizualizace dat Variogramy Ordinary Kriging Cıl projektu R geoR
projekt R
opensource software dostupny proplatformy UNIX, Windows aMacOS
efektivnı nastroj pro manipulaci azpracovanı dat
poskytuje nastroje pro grafickezobrazenı dat
obsahuje radu naimplementovanychfunkcı
moznost rozsırenı pomocı balıcku
Martin Dzurov, Kristyna Kitzbergerova, Lucie Sindelarova Geostatistika v R projektu
Uvod Vizualizace dat Variogramy Ordinary Kriging Cıl projektu R geoR
Balıcek geoR
poskytuje funkce pro geostatistickou analyzu dat v prostredı R
data musı byt ulozena v objektu typu geodata:
objekt musı obsahovat alespon jednu souradnici a hodnotuatributuobjekt se vytvorı pomocı funkce read.geodata(), ktera ztextoveho souboru nacte potrebna data
Parametry odhadnute metodami WLS a MLpro exponencialnı model variogramu
Martin Dzurov, Kristyna Kitzbergerova, Lucie Sindelarova Geostatistika v R projektu
Uvod Vizualizace dat Variogramy Ordinary Kriging teorie zpracovanı dat
Interpolacnı metoda Ordinary Kriging
zakladnı geostatisticka metoda pro urcovanı lokalnıho odhadu
BLUP – nejlepsı linearnı nestranny prediktor
h0 =∑n
i=1 wihi + w0 je linearnı prediktor
h0 je nestranny prediktor, pokud platı E{h0 − h0
}= 0
h0 je nejlepsı, pokud platı E{‖ h0 − h0 ‖2
}je minimalnı
Nastınenı problematiky:
chceme odhadnout hodnotu nadmorske vysky z v mıste x0,kde nemame provedeno merenıpouzijeme linearnı kombinaci hodnot ze sousednıch bodu avah, ktere jsou temto bodum prideleny
z(x0) =
n∑i=1
wiz(xi),
n∑i=1
wi = 1 (7)
Martin Dzurov, Kristyna Kitzbergerova, Lucie Sindelarova Geostatistika v R projektu
Uvod Vizualizace dat Variogramy Ordinary Kriging teorie zpracovanı dat
odhad variance:
E(‖ z0 − z0 ‖2
)= σ2 +
n∑i=1
n∑j=1
wiwjCij − 2n∑
i=1
wiCi0 (8)
kde Cij = cov(zi, zj) = cov(‖ xi, xj ‖) (9)
minimalizacı odhadu variance za podmınky∑n
i=1 wi − 1 = 0
zıskame system rovnic pro Ordinary KrigingC11 · · · C1n 1
... . . ....
...Cn1 · · · Cnn 11 · · · 1 0
︸ ︷︷ ︸
Cn
w1
...wn
λ
︸ ︷︷ ︸
wn
=
C10
...Cn0
1
︸ ︷︷ ︸
dn
(10)
resenım rovnic zıskame neznamou
wn = C−1n dn (11)
odhad variance interpolacnı metody Ordinary Kriging
V ar(z0 − z0) = C00 −n∑
i=1
wiCi0 − λ (12)
Martin Dzurov, Kristyna Kitzbergerova, Lucie Sindelarova Geostatistika v R projektu
Uvod Vizualizace dat Variogramy Ordinary Kriging teorie zpracovanı dat
Vlastnosti
presny interpolator – pokud je hodnota pro vypocet odhadux0 identicka s mıstem, kde bylo provedeno merenı, budevyinterpolovana hodnota totozna s hodnotou namerenou
z(x0) = z(xα), x0 = xα (13)
pouzıva se pro interpolaci dat na pravidelne mrızce pomocınerovnomerne rozmıstenych dat merenı
balıcek geoR obsahuje naimplementovanou funkcikrige.conv()
Martin Dzurov, Kristyna Kitzbergerova, Lucie Sindelarova Geostatistika v R projektu
Uvod Vizualizace dat Variogramy Ordinary Kriging teorie zpracovanı dat
Vlastnı zpracovanı dat
pravidelna mrızka, na ktere bude proveden vypocet
Martin Dzurov, Kristyna Kitzbergerova, Lucie Sindelarova Geostatistika v R projektu
Uvod Vizualizace dat Variogramy Ordinary Kriging teorie zpracovanı dat
Vysledky
predikce nadmorskych vysek metodou OK
Martin Dzurov, Kristyna Kitzbergerova, Lucie Sindelarova Geostatistika v R projektu
Uvod Vizualizace dat Variogramy Ordinary Kriging teorie zpracovanı dat
Vysledky
smerodatna odchylka odhadnutych nadmorskych vysek
Martin Dzurov, Kristyna Kitzbergerova, Lucie Sindelarova Geostatistika v R projektu
Uvod Vizualizace dat Variogramy Ordinary Kriging teorie zpracovanı dat
Dekujeme za pozornost!
Tato prace byla podporena grantem Studentske grantove souteze CVUT
c.OHK1-065/10.
Martin Dzurov, Kristyna Kitzbergerova, Lucie Sindelarova Geostatistika v R projektu