BlackMarbleR Example

blackmarble.R

library(blackmarbler)
library(geodata)
library(sf)

# Set NASA EarthData Token
bearer <- "NASA-EARTHDATA-TOKEN"

# Retrieve GADM polygon of Lebanon
roi <- gadm(country = "LBN", 
                    level=1, 
                    path = tempdir()) |> 
                    st_as_sf()

# Create NTL raster
ntl_r <- bm_raster(roi_sf = roi,
                        product_id = "VNP46A4",
                        date = 2019:2022,
                        bearer = bearer)

# Create NTL dataset from 2012 to 2022 
ntl_df <- bm_extract(roi_sf = roi,
                     product_id = "VNP46A4",
                     date = 2012:2022,
                     bearer = bearer,
                     aggregation_fun = "sum")