Ćwiczenie 1: Geostatystyczna analiza danych

Zadanie samodzielne

Author

Anna Dmowska

1 Zadanie samodzielne

  1. Wczytaj dane z pliku cw_dane.csv. Stwórz obiekt punkty_ni zawierający wyłącznie lokalizacje zmiennej Ni_los.

  2. 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?
  1. 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.
  1. 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.