TL;博士
首先,下面提供了重现问题的所有必要代码。我在使用 grib 部分时遇到了一些挑战,但是对于download.file()
还有更基本的东西。问题
关于在R here中使用GFS天气数据的文章不错。这将是进一步分析的一个很好的起点,但我有几个问题。这篇文章包含了重现它所需的所有R代码,但是这里有前几行也是最重要的部分:
#STEP 1 (doesn't work either way because of bad link)
loc=file.path("ftp://ftp.ncep.noaa.gov/pub/data/nccf/com/gfs/prod/gfs.2009121700/gfs.t00z.sfluxgrbf03.grib2")
我的第一个问题是该链接无法正常工作,但这是来自同一来源的相同格式文件的网址,因此这应该有效:
#PART 1 with working link (at least in google chrome):
loc=file.path("ftp://ftp.ncep.noaa.gov/pub/data/nccf/com/gfs/prod/gfs.2015081112/gfs.t12z.sfluxgrbf00.grib2")
download.file(loc,"temp.grb",mode="wb")
#PART 2
shell("wgrib2 -s temp03.grb | grep :LAND: | wgrib2 -i temp00.grb -netcdf LAND.nc",intern=T)
#PART3
library(ncdf)
landFrac <-open.ncdf("LAND.nc")
land <- get.var.ncdf(landFrac,"LAND_surface")
x <- get.var.ncdf(landFrac,"longitude")
y <- get.var.ncdf(landFrac,"latitude")
#PART4
rgb.palette <- colorRampPalette(c("snow1","snow2","snow3","seagreen","orange","firebrick"), space = "rgb")#colors
image.plot(x,y,t2m.mean,col=rgb.palette(200),axes=F,main=as.expression(paste("GFS 24hr Average 2M Temperature",day,"00 UTC",sep="")),axes=F,legend.lab="o C")
contour(x,y,land,add=TRUE,lwd=1,levels=0.99,drawlabels=FALSE,col="grey30")
我的系统
Windows 7,64位 Rstudio 0.99.467
问题1 - download.file()
我遇到了download.file()的问题,即使是在Chrome中运行良好的链接(或ftp路径)也是如此。这是错误消息的一部分:
在download.file中(url = loc,destfile =&#34; temp.grb&#34;,mode =&#34; wb&#34;):
InternetOpenUrl失败:
任何人都知道是什么原因造成的?
问题:2 - wgrib2
我手动下载了一个示例文件以便继续工作,但在第2部分再次陷入困境。以下r命令在命令窗口中会如何直接显示?
shell("wgrib2 -s temp03.grb | grep :LAND: | wgrib2 -i temp00.grb -netcdf LAND.nc",intern=T)
我已经在下面的链接下安装了必要的dll,并且我已经启动并运行了wgrib 2(显然)。有关wgrib2的信息在文章中提供,并且可用here。 Windows 7的来源。
我希望你们中的一些人觉得这很有趣!