r - NASA 瓷砖与 R 中的传单

Question

我想就传单包寻求一些帮助。当您绘制交互式 map 时，您可以执行以下操作。

library(leaflet)
library(magrittr)

m <- leaflet() %>% 
     setView(lng = -71.0589, lat = 42.3601, zoom = 8) %>%
     addTiles() 

m

如果您想添加第三方磁贴，您也可以这样做。以下链接提供了第三方磁贴的选项 ( http://leaflet-extras.github.io/leaflet-providers/preview/index.html ) 以下图像是使用 OpenWeatherMap.Precipitation 创建的.

### They work
m %>% addProviderTiles("MtbMap")
m %>% addProviderTiles("HikeBike.HikeBike")
m %>% addProviderTiles("OpenWeatherMap.Precipitation")

链接中的一些图块选项包括 NASA 的图块。我想使用其中之一。所以我尝试了以下代码。不幸的是，他们都没有工作。

### The default map appears, then a black layer appears on top of the default layer.
m %>% addProviderTiles("NASAGIBS.ModisTerraTRUEColorCR")
m %>% addProviderTiles("NASAGIBS.ModisTerraBands367CR")

唯一有效的选项如下。

m %>% addProviderTiles("NASAGIBS.ViirsEarthAtNight2012")

我的下一次尝试是使用带有 addTiles() 的自定义 URL 模板。 . URL 来自上面的链接。但是，这也没有成功。没有出现错误消息，但磁贴没有变化。

m %>%addTiles(urlTemplate = "http://map1.vis.earthdata.nasa.gov/wmts-webmerc/MODIS_Terra_CorrectedReflectance_Bands367/default/{time}/{tilematrixset}{maxZoom}/{z}/{y}/{x}.{format}",
              tileOptions(minZoom = 1, maxZoom = 8))

我的最后一次尝试如下。这显示了默认 map ，但也没有出现额外的图块。

leaflet() %>%
addTiles() %>%
setView(lng = -71.0589, lat = 42.3601, zoom = 8) %>%
addTiles(urlTemplate = "http://map1.vis.earthdata.nasa.gov/wmts-webmerc/MODIS_Terra_CorrectedReflectance_Bands367/default/{time}/{tilematrixset}{maxZoom}/{z}/{y}/{x}.{format}",
         tileOptions(minZoom = 1, maxZoom = 8))

我的问题是这是否是 NASA 磁贴的潜在错误。或者，我需要在这些脚本中进行哪些修订？提前感谢您的帮助。

更新

我找到了一个 website它使用相同的 NASA 瓷砖。我指定了 NASAGIBS.ModisTerraTRUEColorCR并得到以下图像。该图像显示了邮件是如何从美国传送到瑞典的。如您所见，美国和欧洲都没有图像。我想这可能是我看到黑色瓷砖的原因。我想知道是否有人知道NASA瓷砖的一些细节。我选择了可​​以看到 NASA 图像的区域。但是，我没有运气。

### I expected to see Japan area this time.
foo <- leaflet() %>% 
       setView(lng = 137.37, lat = 35.93, zoom = 5) %>%
       addTiles() 

foo %>% addProviderTiles("NASAGIBS.ModisTerraTRUEColorCR")

更新 2

今天，我又打了一针。此时，我设法获得了以下图像。我拍的时候缩小了一点。在 UPDATE ，我提供了一张你看不到美国和欧洲的 map 。在新图像中，您会看到美国西海岸是黑色的。鉴于所有观察结果，在我看来，人们可能无法始终获得 NASA 某个位置的图像。根据您询问 NASA 磁贴的时间，您可能/可能没有您想要的图像。

m <- leaflet() %>% 
     setView(lng = -71.0589, lat = 42.3601, zoom = 8) %>%
     addTiles()

m %>% addProviderTiles("NASAGIBS.ModisTerraBands367CR")

Answer 1

您的最终结论是正确的:根据您请求图像的位置和请求时间，卫星可能已获取图像，也可能尚未获取图像。所以你可能会得到一个实际的图像或只是一个空的。 (这也在 GIBS API documentation 中说明。)

但是，您可以通过 addProviderTiles() 的“时间”选项指定请求图像的日期。通过指定最近的日期，如果您愿意，您可以获得所有位置的非空图像。

这是语法:

> library(leaflet)
> library(magrittr)

> m <- leaflet() %>% 
       setView(lng = 4.5, lat = 51, zoom = 1) %>%
       addTiles() %>% 
       addProviderTiles("NASAGIBS.ModisTerraTrueColorCR",
                        options = providerTileOptions(time = "2015-08-31", opacity = 0.5))

> m

在撰写本文时 (2015-08-31)，我得到以下结果:



英格兰的阴天，谁能想到？

大多数数据恰好已经存在，但还没有阿拉斯加的图像。另一方面，如果我指定昨天的日期

options = providerTileOptions(time = "2015-08-30", opacity = 0.5)

我们得到完整的图像:



最后说下原因

m %>% addProviderTiles("NASAGIBS.ModisTerraTRUEColorCR")

没有工作可能是因为打字错误。应该是

m %>% addProviderTiles("NASAGIBS.ModisTerraTrueColorCR")