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:

library(sen2r)

Wczytamy też regionem dla którego pozyskamy dane z użyciem funkcji pakietu sf:

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:

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:

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:

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.