set.seed(1345)
st_bbox(L7_ETMs) |>
st_as_sfc() |>
st_sample(10) -> pts
(x <- st_cells(L7_ETMs, pts))
# get the pixel values (first band only):
st_as_stars(L7_ETMs)[[1]][x]
# get pixel values for all bands:
st_as_stars(L7_ETMs) |> split() |> sapply(`[`, x)
# compare with st_extract():
st_as_stars(L7_ETMs) |> split() |> st_extract(pts)
Run the code above in your browser using DataLab