From 8ecbcc01286d0af9518d4f7322774358ff5444f9 Mon Sep 17 00:00:00 2001
From: Carl <cboettig@gmail.com>
Date: Thu, 7 Dec 2023 22:40:04 +0000
Subject: [PATCH] cleaner

---
 drafts/sst.qmd | 22 ++++++++++++++--------
 1 file changed, 14 insertions(+), 8 deletions(-)

diff --git a/drafts/sst.qmd b/drafts/sst.qmd
index da9b49a..7aa8bcc 100644
--- a/drafts/sst.qmd
+++ b/drafts/sst.qmd
@@ -27,6 +27,7 @@ tm_shape(World) + tm_borders() +
 
 
 ```{r}
+bench::bench_time({
 start =  min(turtles$date) # "2018-01-01" #
 end =  max(turtles$date) # "2018-12-31"   #
 items <- stac("https://cmr.earthdata.nasa.gov/stac/POCLOUD") |>
@@ -35,28 +36,33 @@ items <- stac("https://cmr.earthdata.nasa.gov/stac/POCLOUD") |>
               datetime = paste(start,end, sep = "/")) |>
   post_request() |>
   items_fetch()
+})
 ```
 
 
 
 ```{r}
 stac_dates <- rstac::items_datetime(items) |> as.Date()
-matched <- items[ stac_dates %in% dates ]
+matched <- items$features[ stac_dates %in% dates ]
 
 ```
 
 
 
 ```{r}
-library(stars)
+earthdatalogin::edl_set_token()
 
-feature <- matched$
+urls <- map_chr(matched, list("assets", "data", "href"))
 
-for(feature in matched$features) {
-  url <- feature$assets$data$href
-  ex <- read_stars(paste0("/vsicurl/", url), "analysed_sst")
-  turtle_temp <- st_extract(ex, turtles)
-}
+# smaller test set:
+urls <- urls[1:100]
+ex <- read_stars(paste0("/vsicurl/", urls), "analysed_sst", quiet=TRUE)
+st_crs(ex) <- 4326
+turtle_temp <- st_extract(ex, turtle)
+```
+
+```{r}
+image_collection(urls)
 ```