三、基礎地圖

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"
? )

待續未完。

總結

以上是生活随笔為你收集整理的R语言在线地图神器：Leaflet for R包（二）基础底图的全部內容，希望文章能夠幫你解決所遇到的問題。

如果覺得生活随笔網站內容還不錯，歡迎將生活随笔推薦給好友。