使用RStudio和wgrib2重现天气研究,下载和提取文件

时间:2015-08-25 12:41:28

标签: r rstudio grib

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的信息在文章中提供,并且可用hereWindows 7的来源。

我希望你们中的一些人觉得这很有趣!

0 个答案:

没有答案