2017-01-17 74 views
0

我试图让从KNMI德precipation雷达数据在R.传单的数据位于这里加载检索WMS图像:http://adaguc.knmi.nl/contents/webservices/WebServices_RADNL_OPER_R___25PCPRR_L3.html如何从特定时间/日期

到目前为止,我已经设法让当前的雷达图像与下面的代码瓣叶层:

require("leaflet") 
require("dplyr") 
leaflet() %>% addTiles() %>% 
    addWMSTiles(
    "http://geoservices.knmi.nl/cgi-bin/RADNL_OPER_R___25PCPRR_L3.cgi", 
    layers = "RADNL_OPER_R___25PCPRR_L3_KNMI", 
    options = WMSTileOptions(format = "image/png", transparent = TRUE), 
    attribution = "KNMI" 
) 

到目前为止好...

但其实我想,是从给定获取雷达图像时间戳。该系列可追溯至2009年(请参阅KNMI查看器中的示例:http://geoservices.knmi.nl/viewer2.0/?srs=EPSG%3A3857&bbox=-725.7203842048766,6500000,1220725.7203842048,7200000&service=http%3A%2F%2Fgeoservices.knmi.nl%2Fcgi-bin%2FRADNL_OPER_R___25PCPRR_L3.cgi%3F&layer=RADNL_OPER_R___25PCPRR_L3_COLOR%24image%2Fpng%24true%24default%241%240&selected=0&dims=time $ current & baselayers = streetmap $ world_line)。

在网络上,我发现了几个支持日期/时间的小册子的示例(如:http://apps.socib.es/Leaflet.TimeDimension/examples/example14.html),但不幸的是我缺乏将此示例转换为R版小册子的技能。

可以这样做吗?例如,我可以从2017-01-07T04:00检索降水雷达图像吗?如果是的话:我会在哪里开始?

或者如果不能用传单完成;任何其他建议?

回答

0

这是更容易比预期;-)

我不得不在WMSTileOptions添加随时间变化的。

这工作:

require("leaflet") 
require("dplyr") 
leaflet() %>% addTiles() %>% 
    addWMSTiles(
    "http://geoservices.knmi.nl/cgi-bin/RADNL_OPER_R___25PCPRR_L3.cgi", 
    layers = "RADNL_OPER_R___25PCPRR_L3_COLOR", 
    options = WMSTileOptions(format = "image/png", transparent = TRUE, TIME = "2017-01-07T04:00:00"), 
    attribution = "KNMI neerslagradar" 
)