如何将 terra R package 与需要身份验证的云优化 geotiffs 一起使用？

Question

我目前正在使用栅格 package 访问 NASA Earthdata 云优化的 geotiffs，但想切换到 terra package。当我尝试使用 terra 时，出现错误。 我认为这是因为我没有传递必要的 GDAL 配置选项。

使用光栅 package 时，我执行以下操作：

library(raster) 

library(rgdal) 

rgdal::setCPLConfigOption(ConfigOption = "GDAL_HTTP_UNSAFESSL", value = "YES") 

rgdal::setCPLConfigOption(ConfigOption = "GDAL_HTTP_COOKIEFILE", value = ".rcookies") 

rgdal::setCPLConfigOption(ConfigOption = "GDAL_HTTP_COOKIEJAR", value = ".rcookies") 

rgdal::setCPLConfigOption(ConfigOption = "GDAL_DISABLE_READDIR_ON_OPEN", value = "EMPTY_DIR") rgdal::setCPLConfigOption(ConfigOption = "CPL_VSIL_CURL_ALLOWED_EXTENSIONS", value = "TIF") 

url <- "/vsicurl/https://data.lpdaac.earthdatacloud.nasa.gov/lp-prod-protected/HLSS30.020/HLS.S30.T10SEJ.2021214T184919.v2.0/HLS.S30.T10SEJ.2021214T184919.v2.0.B8A.tif" 

r <- raster(url) 

一个 NASA Earthdata 帐户，一个带有登录名和密码的 .netrc 文件，以及这些 GDAL 配置选项是这个工作所必需的。

NASA Earthdata 登录/注册： https://urs.earthdata.nasa.gov/

Rscript to create.netrc file: https://git.earthdata.nasa.gov/projects/LPDUR/repos/hls_tutorial_r/browse/Scripts/earthdata.netrc_setup.R

如果我想使用 terra package 做同样的事情，我该如何传递必要的 GDAL 配置选项？ 还是有其他问题？

使用 terra 的示例：

library(terra) 

url <- "/vsicurl/https://data.lpdaac.earthdatacloud.nasa.gov/lp-prod-protected/HLSS30.020/HLS.S30.T10SEJ.2021214T184919.v2.0/HLS.S30.T10SEJ.2021214T184919.v2.0.B8A.tif" 

r <- rast(url) 

Answer 1

你可以做

library(terra) 
setGDALconfig("GDAL_HTTP_UNSAFESSL", "YES")
setGDALconfig("GDAL_HTTP_COOKIEFILE", ".rcookies") 
setGDALconfig("GDAL_DISABLE_READDIR_ON_OPEN", "EMPTY_DIR")
setGDALconfig("CPL_VSIL_CURL_ALLOWED_EXTENSIONS", "TIF") 

url <- "/vsicurl/https://data.lpdaac.earthdatacloud.nasa.gov/lp-prod-protected/HLSS30.020/HLS.S30.T10SEJ.2021214T184919.v2.0/HLS.S30.T10SEJ.2021214T184919.v2.0.B8A.tif" 

r <- rast(url)