Darmowe dane wysokościowe SRTM w R – funkcja getData

W tym poście pokażemy jak wykorzystać funkcję getData znanego już nam pakietu raster do pobrania danych SRTM. Funkcja pozwala na pobranie danych z trzech darmowych baz:

• ‘SRTM’ – dane wysokościowe z misji SRTM dostępne również na stronie: http://srtm.csi.cgiar.org/
• ’GADM’ – globalna baza granic administracyjnych państw. Więcej na stronie https://gadm.org/,
• . 'worldclim’ – globalna baza rastrowych danych klimatycznych. Więcej na stronie http://www.worldclim.org/.

Nas interesuje baza z danymi wysokościowymi zdefiniowana jako “SRTM”. Dane SRTM zapisane są w kaflach o wielkości 5 x 5 stopni. Do pobrania danych przy użyciu funkcji konieczne jest również zdefiniowanie współrzędnych geograficznych interesującego nas obiektu, np. szerokość 52 stopni i długość 21 stopni. Nasza funkcja powinna wyglądać następująco:

dem <- getData("SRTM",lon=21,lat=52)

Na nasz komputer zostanie pobrany bezpośrednio ze strony kafelek zawierający wskazane współrzędne:

> dem
class       : RasterLayer
dimensions  : 6001, 6001, 36012001  (nrow, ncol, ncell)
resolution  : 0.0008333333, 0.0008333333  (x, y)
extent      : 19.99958, 25.00042, 49.99958, 55.00042  (xmin, xmax, ymin, ymax)
coord. ref. : +proj=longlat +datum=WGS84 +ellps=WGS84 +towgs84=0,0,0
data source : d:/Program Files/RStudio/srtm_41_02.tif
names       : srtm_41_02
values      : -32768, 32767  (min, max)

Zobaczmy jak wygląda pobrany plik:

plot(dem)

Jak pobrać dane z pozostałych baz dowiecie się z dokumentacji funkcji getData. Powiemy Wam tylko że dla obu pozostałych baz można pobierać dane o różnej szczegółowości. Miłej nauki:)