Permian-Global-Research / chewie

A package to efficiently download GEDI data
Apache License 2.0
5 stars 1 forks source link

Error with grab_gedi #11

Open emilyluisedoyle opened 7 months ago

emilyluisedoyle commented 7 months ago

DAAC1819S_2_search <- find_gedi(DAAC1819S_2,

  • gedi_product = "2A",
  • date_start = "2019-01-01",
  • date_end = "2019-12-31") ✔ Using cached GEDI data Error in [[<-.data.frame(*tmp*, i, value = FALSE) : replacement has 1 row, data has 0

Assuming it is because there are no GEDI files within the AOI, but error code to bypass this also did not work

h-a-graham commented 6 months ago

Hey Em,

I'm struggling to reproduce this one... Does the following work for you:

library(chewie)
#> ✔ NASA Earthdata Credentials already set.
#> ✔ GEDI cache set in the following directory:
#> → "/home/hugh/.chewie"
aoi2 <- sf::st_bbox(c(
  xmin = -67.669, ymin = -9.762,
  xmax = -67.622, ymax = -9.748
)) |>
  sf::st_as_sfc() |>
  sf::st_set_crs(4326)

x2 <- find_gedi(aoi2,
  gedi_product = "2A",
  date_start = "2019-01-01",
  date_end = "2019-12-31",
  cache = FALSE
)

x2
#> 
#> ── chewie.find ───────────────────────────────────────────────────────────────────────────────────────────────────────────────────
#> •  GEDI-2A
#>                        id          time_start            time_end                                                   url
#>                    <char>              <POSc>              <POSc>                                                <char>
#> 1: G2006475899-LPDAAC_ECS 2019-07-10 12:57:16 2019-07-10 14:30:03 https://e4ftl01.cr.usgs.gov//GEDI_L1_L2/GEDI/GEDI0...
#> 2: G2007536358-LPDAAC_ECS 2019-09-11 11:59:57 2019-09-11 13:32:46 https://e4ftl01.cr.usgs.gov//GEDI_L1_L2/GEDI/GEDI0...
#> 3: G2007171172-LPDAAC_ECS 2019-10-02 15:28:58 2019-10-02 17:01:48 https://e4ftl01.cr.usgs.gov//GEDI_L1_L2/GEDI/GEDI0...
#> 2 variables not shown: [geometry <sfc_POLYGON>, cached <lgcl>]
#> 
#> ──────────────────────────────────────────────────────────────────────────────────────────────────────────────────────────────────
plot(x2)

Created on 2024-02-27 with reprex v2.0.2