我正在尝试自动获取波浪高度预测的数据并构建类似于此的图.
使用R,我可以下载数据并使用以下方法绘制它:
library(rgdal)
library(fields)
ftp.string <- "ftp://polar.ncep.noaa.gov/pub/waves//20130205.t00z/nww3.HTSGW.grb"
#this link may become broken with time, as folders are removed after some time. just edit the date to reflect the most recent day at the time you run these lines
download.file(ftp.string, "foo.grb", mode="wb")
grib <- readGDAL("foo.grb")
is.na(grib$band1) <- grib$band1 > 100
image(grib, col=(tim.colors(15)), attr=1)
但是,如果你仔细看看我上面发布的链接,你会发现一个微妙的差异:链接的图表经度超过360度.
这对我正在做的事情非常重要,因为它允许我轻松检查同一地块内所有海洋上的膨胀 - 如果一次只显示360度,这就更难了,因为这会强制导致其中一个海洋切.
尽管我付出了最大的努力,但我找不到绘制超过360度的方法,因为GRIB格式"太聪明"不允许(这不仅仅是偏移数据,而是重复其中的一部分).
任何见解将不胜感激.干杯