R语言在线地图神器：Leaflet for R包（二）基础底图

三、基础地图

leaflet包里面，内置了多个基础地图，包括了OpenStreetMap、Esri和CartoDB等等。还可以加载自定义的WMTS，当然，有些地图是需要翻墙的。
设置底图的方法是：

addProviderTiles("地图标号")
默认用的osm。
下面是虾神找到的，可以直接在leaflet包中直接加载的地图标号（不用翻墙啊！再一次热泪盈眶）

###############################
# OpenStreetMap.Mapnik
# OpenStreetMap.BlackAndWhite
# OpenStreetMap.DE
# OpenStreetMap.France
# OpenStreetMap.HOT
# OpenTopoMap
# Thunderforest.OpenCycleMap
# Thunderforest.Transport
# Thunderforest.TransportDark
# Thunderforest.SpinalMap
# Thunderforest.Landscape
# Thunderforest.Outdoors
# Thunderforest.Pioneer
# OpenMapSurfer.Roads
# OpenMapSurfer.Grayscale
# Hydda.Full
# Stamen.Toner
# Stamen.TonerBackground
# Stamen.TonerLite
# Stamen.Watercolor
# Stamen.Terrain
# Stamen.TerrainBackground
# Stamen.TopOSMRelief
# Esri.WorldStreetMap
# Esri.DeLorme
# Esri.WorldTopoMap
# Esri.WorldImagery
# Esri.WorldTerrain
# Esri.WorldShadedRelief
# Esri.WorldPhysical
# Esri.OceanBasemap
# Esri.NatGeoWorldMap
# Esri.WorldGrayCanvas
# MtbMap
# CartoDB.Positron
# CartoDB.PositronNoLabels
# CartoDB.PositronOnlyLabels
# CartoDB.DarkMatter
# CartoDB.DarkMatterNoLabels
# CartoDB.DarkMatterOnlyLabels
# HikeBike.HikeBike
# HikeBike.HillShading
# NASAGIBS.ModisTerraTrueColorCR
# NASAGIBS.ModisTerraBands367CR
# NASAGIBS.ViirsEarthAtNight2012
# NASAGIBS.ModisTerraLSTDay
# NASAGIBS.ModisTerraSnowCover
# NASAGIBS.ModisTerraAOD
# NASAGIBS.ModisTerraChlorophyll
#################################################
挑几个来看看，首先挑Esri提供的四个地图：

leaflet()%>%setView(lng=116.38,lat=39.9,zoom=3)%>%
  addTiles()%>%addProviderTiles("Esri.WorldStreetMap")

leaflet()%>%setView(lng=116.38,lat=39.9,zoom=3)%>%
  addTiles()%>%addProviderTiles("Esri.WorldImagery")

leaflet()%>%setView(lng=116.38,lat=39.9,zoom=3)%>%
  addTiles()%>%addProviderTiles("Esri.WorldTerrain")

leaflet()%>%setView(lng=116.38,lat=39.9,zoom=3)%>%
  addTiles()%>%addProviderTiles("Esri.NatGeoWorldMap")

好了，以后再也不会说缺少底图资源了，更神奇的是，leaflet可以加载大量NASA提供的数据，而且有些数据还可以选择不同时间，如下：

加载NASA提供的2012年全球灯光夜景数据：
leaflet()%>%setView(lng=116.38,lat=39.9,zoom=2)%>%
  addTiles()%>%addProviderTiles("NASAGIBS.ViirsEarthAtNight2012")


加载2015年1月15日、5月15日和8月15日和11月15日的modis数据：
leaflet() %>%
  setView(lng = 110, lat = 30, zoom = 2) %>%
  addTiles() %>%
  addProviderTiles("NASAGIBS.ModisTerraTrueColorCR",
                   options = providerTileOptions(
                     time = "2015-01-15", opacity = 1))

leaflet() %>%
  setView(lng = 110, lat = 30, zoom = 2) %>%
  addTiles() %>%
  addProviderTiles("NASAGIBS.ModisTerraTrueColorCR",
                   options = providerTileOptions(
                     time = "2015-05-15", opacity = 1))


leaflet() %>%
  setView(lng = 110, lat = 30, zoom = 2) %>%
  addTiles() %>%
  addProviderTiles("NASAGIBS.ModisTerraTrueColorCR",
                   options = providerTileOptions(
                     time = "2015-08-15", opacity = 1))

leaflet() %>%
  setView(lng = 110, lat = 30, zoom = 2) %>%
  addTiles() %>%
  addProviderTiles("NASAGIBS.ModisTerraTrueColorCR",
                   options = providerTileOptions(
                     time = "2015-11-15", opacity = 1))


当然，还可以加载自己的WMS服务，不过鉴于国内WMS服务实在不好找，找到的又实在不好用，所以这里用官方的示例：

leaflet() %>% addTiles() %>% setView(-93.65, 42.0285, zoom = 4) %>%
  addWMSTiles(
    "http://mesonet.agron.iastate.edu/cgi-bin/wms/nexrad/n0r.cgi",
    layers = "nexrad-n0r-900913",
    options = WMSTileOptions(format = "image/png", transparent = TRUE),
    attribution = "Weather data ? 2012 IEM Nexrad"
  )

待续未完。

更多的更新的文章，虾神会第一时间在公众号推送，请关注虾神公众号：虾神daxialu