Dane satelitarne Sentinel-2 w R
Zbierane od 2015 roku dane z europejskiego satelity Sentinel stanowią doskonałe bezpłatne źródło satelitarnych danych obrazowych. W niniejszym poście pokażemy Wam jak można pobrać i przetworzyć te dane w R dla określonego rejonu. Wykorzystamy do tego pakiet sen2r. Załadujmy go:
|
1 |
library(sen2r) |
Wczytamy też regionem dla którego pozyskamy dane z użyciem funkcji pakietu sf:
|
1 2 |
library(sf) border <- read_sf("D:/GIS_in_R/warsaw_border.shp") |
Pakiet sen2r umożliwia pobieranie i przetworzenie danych za pomocą GUI. Graficzny Interfejs Użytkownika wywołuje się funkcją sen2r bez atrybutów:
|
1 |
sen2r() |
My jednak pobierzemy i przetworzymy dane Sentinel-2 z użyciem bezpośrednio funkcji tego pakietu. Zapoznanie z GUI pozostawiamy Wam.
Do wylistowania dostępnych dla naszego regionu zobrazowań użyjemy funkcji s2_list. W funkcji tej zdefiniować można wiele atrybutów tj. zakres przestrzenny (spatial_extent), zakres czasu (time_interval), maksymalne pokrycie chmurami (max_cloud) i wiele innych. Nasza funkcja będzie wyglądała następująco:
|
1 2 3 4 5 6 7 |
images_list <- s2_list( spatial_extent = border, time_interval = as.Date(c("2015-05-01", "2020-08-30")), max_cloud = 1 ) images_list<- as.data.frame(images_list) |
Tabelę danych images_list możemy dowolnie filtrować. Z pobranej listy pobierzemy z wykorzystaniem sen2r zobrazowanie z 2020-07-01:
|
1 2 3 4 5 6 7 8 9 |
sen2r(gui = FALSE, timewindow = c(as.Date("2020-07-01")), extent = border, list_prods = c("BOA"), list_indices = c("NDVI"), s2_levels = "l2a", path_l2a = "d:/GIS_in_R/s2_data/l2a/", path_out = "d:/GIS_in_R/s2_data/" ) |
W folderze wyjściowym w BOA i NDVI zapisane zostały zdefiniowane przetworzenia, a w folderze l2a dane źródłowe.
Pakiet sen2r zawiera wiele przydatnych funkcji, które można wykorzystać nie tylko do danych z satelity Sentinel-2. Odkrywanie tego pakietu w szczegółach pozostawiamy Wam.
