Prostorová data

Interpretace geodat

Jan Caha

2024-04-15

Balíky

  • sf - Simple Features balík založený na GDAL/OGR pro načítání a zpracování zejména vektorových dat
  • raster a terra (novější iterace) - balíky pro načítání a zpracování rastrových dat
  • stars - balík pro načítání a zpracování zejména rastrů v podobně multidimenzionálních polí (např. i časové řady)
  • tmap - tvorba vizualizací prostorových dat
  • qgisprocess a qgis

sf

  • OGC standard simple features
  • data se načtou jako standatní data.frame či tibble, geometrie je zvláštní sloupec se specificky uloženými daty
  • pro práci s daty je možné použít běžné dplyr funkce (u některých dostáváme “zdarma” prostorové rozšíření - např. kombinace group_by() a summarize() automaticky agreguje i geometrie)
  • zvláštní funkce balíku využívají předponu st_* (spatiotemporal)
  • další funkce v balíku lwgeom

Ukázky sf

raster a terra

  • velice podobné balíky, terra je novější a má některé vylepšení
  • pracuje hlavně s rastrovými daty, ale umí i něco málo s vektorovými
  • rastry se nenačítají do paměti, pokud o to explicitně nežádáme, drží se v principu odkaz na data + metadata
  • data se získavají jako matrix

Ukázky raster a terra

stars

  • postaveno na principu datových kostek (v GIS relativně málo používané)
  • zaměřeno na multidimenzionální data (typicky formát netCDF), často s časovou složkou
  • příklady v vignetách

tmap

  • balík pro tvorbu vizualizací prostorových dat
  • jak statické tak i dynamické vizualizace (s využitím leaflet)
  • odvozeno od ggplot2 podobný zápis, ale nelze kombinovat
  • pracuje s vrstvami
  • ukázky v vignetách

qgisprocess a qgis

  • balíky pro spuštění QGIS procesů z R
  • qgisprocess je obecnější, ale uživatel musí znát QGIS procesy
  • qgis je nadstavba, která vytváří pro každý QGIS nástroj odpovídající funkci v R a to včetně nápovědy

Ukázky qgisprocess a qgis

  • instalace
install.packages(c("qgisprocess", "remotes"))
remotes::install_github("JanCaha/r_package_qgis")
  • maximálně podporují existující balíky, takže lze předávat jako vstupní data např. data.frame se sloupcem geometrie či raster a stars objekty
  • ukázky ve vignetách qgisprocess a qgis

Shrnutí

  • poměrně rozsáhlý ekosystém balíků pro práci s prostorovými daty
  • lze srovnávat s ekosystémem v jazyce Python
  • více zaměřené na zpracování dat
  • lépe řešené vizualizace
  • díky základu jazyka R lepší integrace s ostatními balíky, odpadají konverze dat mezi formáty různých balíků

Dotazy?