得到OPenDAP ncml需要身份验证，使用R汇总数据

Question

我试图从使用R的raster包NASA OPenDAP服务器获取数据TRMM起初我有关于认证一定的难度，但问题是resolved。

NASA OPenDAP服务器发布TRMM 3B42_daily数据作为单独的文件，每天一个和（使用ncml）聚合的年度数据。所以，我现在的问题是，，使用R raster包和认证文件.dodsrc和.netrc我可以下载单独的NetCDF文件，但我不能下载的汇总数据。

所以，这个工作：

library(raster)
single_date_opendap <- 'https://disc2.gesdisc.eosdis.nasa.gov:443/opendap/TRMM_L3/TRMM_3B42_Daily.7/2002/04/3B42_Daily.20020405.7.nc4'
test <- stack(single_date_opendap, varname = 'precipitation')

这不：

library(raster)
url_opendap_no_brkt <- 'https://disc2.gesdisc.eosdis.nasa.gov:443/opendap/ncml/aggregation/TRMM_3B42_Daily.7/TRMM_3B42_daily.7_Aggregation_2001.ncml'
test <- stack(url_opendap_no_brkt, varname = 'precipitation')

而且给我的错误信息：

Error in .local(.Object, ...) :
An error occurred while creating a virtual connection to the DAP server:
Error while reading the URL: https://disc2.gesdisc.eosdis.nasa.gov:443/openda
p/ncml/aggregation/TRMM_3B42_Daily.7/TRMM_3B42_daily.7_Aggregation_2001.ncml.
ver.
The OPeNDAP server returned the following message:
Unauthorized: Contact the server administrator.
Error in .rasterObjectFromFile(x, band = band, objecttype = "RasterLayer",
Cannot create a RasterLayer object from this file. (file does not exist)

是否有可能得到一个发布聚合数据OPenDAP服务器的数据？

Answer 1

与美国航天局的支持，并与安东尼奥的小技巧交流后，发现是R raster包将不会与聚合数据集工作。但是ncdf4::nc_open能够处理它。奇怪，因为，从我的理解，raster包在后台调用nc_open。

不管怎么说，这个工作：

library(ncdf4)
url_opendap <- 'https://disc2.gesdisc.eosdis.nasa.gov:443/opendap/ncml/aggregation/TRMM_3B42_Daily.7/TRMM_3B42_daily.7_Aggregation_2001.ncml'
trmm <- nc_open(url_opendap)

这不

library(raster)
url_opendap <- 'https://disc2.gesdisc.eosdis.nasa.gov:443/opendap/ncml/aggregation/TRMM_3B42_Daily.7/TRMM_3B42_daily.7_Aggregation_2001.ncml'
trmm <- stack(url_opendap, varname = "precipitation")