提取bz2文件并使用R打开ncdf

时间:2016-07-30 16:48:15

标签: r ftp netcdf bz2

我试图从ftp打开bz2压缩的netcdf文件,但无法让它工作。非常感谢任何帮助。

我尝试使用以下方式直接打开它:

url <- 'ftp://podaac-ftp.jpl.nasa.gov/allData/ghrsst/data/L4/GLOB/NCDC/AVHRR_OI/1982/001/19820101-NCDC-L4LRblend-GLOB-v01-fv02_0-AVHRR_OI.nc.bz2'
nc_open(url)

但这不起作用(我得到:Error in R_nc4_open: No such file or directory)。我想因为压缩,所以我尝试先下载文件并解压缩。

temp <- tempfile()
download.file(url,temp)
nc_open(bzfile(temp, '19820101-NCDC-L4LRblend-GLOB-v01-fv02_0-AVHRR_OI.nc', 'rb'))

但这也不起作用,我得到错误和警告:

Error in bzfile(temp, "19820101-NCDC-L4LRblend-GLOB-v01-fv02_0-AVHRR_OI.nc",  : 
  cannot open the connection
In addition: Warning message:
In bzfile(temp, "19820101-NCDC-L4LRblend-GLOB-v01-fv02_0-AVHRR_OI.nc",  :
  cannot open bzip2-ed file '/var/folders/hs/k9t_8wxs2hn48qq4vp_7xmgm0000gn/T//Rtmpv9UIBr/filef5659606aee', probable reason 'Invalid argument'

任何想法?

由于

1 个答案:

答案 0 :(得分:2)

这似乎按照你想要的方式工作:

A  B  C  D  E  F  G0
A  B  C  D  E  F  G1
A* B* C* D* E* F* G0*
A* B* C* D* E* F* G1*