Ćwiczenie 1: Geostatystyczna analiza danych
Zadanie samodzielne
1 Zadanie samodzielne
Wczytaj dane z pliku
cw_dane.csv. Stwórz obiektpunkty_nizawierający wyłącznie lokalizacje zmiennej Ni_los.Przeprowadź podstawową eksploracyjną analizę danych.
- Wykonaj mapę lokalizacyjną pokazującą przestrzenny rozkład stężenia niklu w lokalizacjach próbki losowej (Ni_los).
- Wykonaj histogram dla zmiennej Ni_los. Czy w danych występują wartości globalnie odstające lub wartości anomalne?
- Analiza i modelowanie struktury przestrzennej.
Wykonaj semiwariogram empiryczny wykorzystując wszystkie 270 pomiarów zmiennej Ni_los.
Usuń wartości anomalne (Lp => 434,973,974,974,976,977) oraz wartości odstające i ponownie wykonaj semiwariogram empiryczny. Sprawdź różne ustawienia szerokości odstępu.
Dopasuj model do semiwariogramu empirycznego.
- Porównaj model dopasowany “manualnie” oraz “półautomatycznie” z użyciem funkcji fit.variogram(). Który model wybierzesz do dalszych analiz?
- Określ parametry modelu.
- Wykonaj kriging prosty oraz kriging zwykły w siatce o parametrach xmin=616550, xmax=643850, ymin=5785850, ymax = 5822150, oraz rozdzielczości 100m. Zwizualizuj wyniki.