Luna包getNASA函数无法获取NASA MCD19A3D v061产品的求助

• 修正产品名称参数
  你要下载的是带D后缀的逐日产品MCD19A3D，原代码使用的MCD19A3是不同时间分辨率的产品，这会导致查询无结果。直接修改product参数为完整产品名，同时补全未定义的path变量：

  library(luna)
  library(terra)
  library(geodata)
  
  # 定义数据存储路径
  path <- "./nasa_mcd_data"
  dir.create(path, showWarnings = FALSE)
  
  senegal <- geodata::gadm("Senegal", level=1, path=".")
  
  available <- getNASA(product = "MCD19A3D", version = "061", 
                       start_date = "2024-01-01", 
                       end_date = "2024-03-01",
                       aoi = senegal, 
                       download = FALSE,
                       path = path)
  

• 验证函数支持性
  如果修改产品名后仍无结果，说明getNASA函数可能未收录MCD19A3D产品。此时可以用LP DAAC官方的STAC接口下载，步骤如下（需先注册NASA Earthdata账号）：

  library(earthdata)
  library(rstac)
  library(purrr)
  
  # 登录Earthdata账号（首次运行会提示输入用户名密码）
  edl_login()
  
  # 初始化STAC客户端
  client <- stac("https://cmr.earthdata.nasa.gov/stac/LPCLOUD/")
  
  # 搜索目标数据
  items <- client %>%
    stac_search(collections = "MCD19A3D.061",
                bbox = terra::ext(senegal) %>% as.numeric(),
                datetime = "2024-01-01/2024-03-01") %>%
    get_request()
  
  # 下载数据到指定路径
  if (length(items$features) > 0) {
    urls <- map_chr(items$features, ~.x$assets$`001`$href)
    walk(urls, ~download.file(.x, file.path(path, basename(.x)), mode = "wb"))
  }
  

• 检查数据发布状态
  部分NASA产品存在发布延迟，若2024年初的数据未上线，可尝试缩小时间范围（比如改用2023年的日期）测试是否能查询到结果。

内容的提问来源于stack exchange，提问作者cha116