2017-01-17 75 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" 
)