R语言在线地图神器：Leaflet for R包（六,完结篇）栅格、颜色和图例设置

八、栅格数据 栅格数据是当前GIS领域里面最庞大的一类数据（天上N颗卫星24小时无死角的对地球君进行直播……每天生成遥感数据都以PB为单位），leaflet作为地图可视化模块，自然也提供了加载栅格数据的方法。 方法如下：addRasterImage(地图，参数) 要注意的一点是，栅格不像矢量数据，矢量本身就带有空间参考信息，而栅格有些是不带空间参考的，如果你不指定投影参数，那么leaflet会默认把加载的栅格投影成WGS84 Web Mercator (Auxiliary Sphere)（EPSG:3857），也就是谷歌地图用的那种空间参考。 下面通过两个示例来说明，首先是加载一个NC（netCDF）文件：注意，需要安装raster包和ncdf4包

library(raster) r <- raster("data/Leaflet/YSUZ98_KWBN_201312130055-var0-z0-rt0-t0.nc") pal <- colorNumeric(c("transparent",topo.colors(100,alpha = NULL)), values(r),                     na.color = "transparent") leaflet() %>% addTiles() %>%   addRasterImage(r, colors = pal, opacity = 0.8) %>%   addLegend(pal = pal, values = values(r),             title = "USA Snow 2013-12-13")

  

 下面加载一个北京主干道的车流密度分析栅格：

r <- raster("E:/flow.tif") pal <- colorNumeric(c("transparent",topo.colors(20,alpha = NULL)), values(r),                     na.color = "transparent") leaflet() %>% addTiles() %>%   addRasterImage(r, colors = pal, opacity = 0.8) %>%   addLegend(pal = pal, values = values(r),position="bottomleft",             title = "北京交通主干道车流密度分析")

 

  九、颜色设置 空间可视化的最重要部分，就是把空间要素上的附件属性映射到颜色上，通过颜色的对比，来反应变量之间的关系，这些变量可以是数值，也可以是类别。 R语言号称顶级的可视化语言，所以对于颜色设置的方法多到你眼花缭乱，有兴趣的同学可以去参考R语言的相关教程。  leaflet里面的颜色设置一共有个参数，分布是color和fillcolor，从名称就可以看出来，一个是主颜色，一个是填充色，那么对应各种对象来说，color一般就是所谓的边框的颜色了。 在配合leaflet进行数值渲染的R语言调色函数里面，最常用的有以下几个：1、用于连续数值的：colorNumeric，colorBin和colorQuantile;2、用于分类输入，colorFactor。 下面分别简单解释一下：colorNumeric: 将连续数值线性映射为设定的颜色模式的过程，一般来说，会对设定的颜色模式根据数值的变化进行平滑插值。 colorBin:也是将连续数值线性映射为设定的颜色模式的过程，也会对数据进行平滑内插生成颜色，但是与上一种方法不同的是，这个方法会对设定的颜色模式根据数值的变化进行分级设定。 colorQuantile ： 也是线性映射，但是设定的方式是通过百分位数进行分级。 colorFactor： 把分类映射到颜色模式，如果颜色模式设定的数量和分类数量不一样，那么就对颜色模式进行平滑内插。 当然，你也可以使用R语言提供的那些调色板，比如heat.colors、cm.colors、rainbow等等，相关内容请查阅资料。 对于颜色设置的，实际上每一篇博客的demo的时候，都做了颜色渲染，大家有兴趣的时候回头去看，下面给出本章的一个实例： 

poly <- readShapePoly(paste(path,"CNPG_S.shp",sep ="")) pal <- colorNumeric(c("darkgreen", "yellow", "orangered"),poly@data$Pop_2009) leaflet(poly) %>% addTiles() %>%   addPolygons(color=~pal(poly@data$Pop_2009),fillOpacity  = 0.8,weight=1)%>%   addLegend(pal = pal, values = poly@data$Pop_2009,position="bottomright",           title = "2009年人口数量（万人）")   pal <- colorNumeric("Greens",poly@data$Pop_2009) leaflet(poly) %>% addTiles() %>%   addPolygons(color=~pal(poly@data$Pop_2009),fillOpacity  = 0.8,weight=1)%>%   addLegend(pal = pal, values = poly@data$Pop_2009,position="bottomright",             title = "2009年人口数量（万人）")   pal <- colorBin(c("darkgreen", "yellow", "orangered"),poly@data$Pop_2009,10) leaflet(poly) %>% addTiles() %>%   addPolygons(color=~pal(poly@data$Pop_2009),fillOpacity  = 0.8,weight=1)%>%   addLegend(pal = pal, values = poly@data$Pop_2009,position="bottomright",             title = "2009年人口数量分级")   pal <- colorBin("Greens",poly@data$Pop_2009,10) leaflet(poly) %>% addTiles() %>%   addPolygons(color=~pal(poly@data$Pop_2009),fillOpacity  = 0.8,weight=1)%>%   addLegend(pal = pal, values = poly@data$Pop_2009,position="bottomright",             title = "2009年人口数量分级")

 

 十、图例： 图例是对地图可视化之后相关颜色或者符号的映射说明——当然，在leaflet里面，主要是针对颜色创建图例。 主要是方法如下： addLegend（地图，属性，颜色列表，标记值……） 具体的属性方法，大家可以查阅帮助文档，下面给出一个示例： 

pal <- colorNumeric("Greens",poly@data$Pop_2009) leaflet(poly) %>% addTiles() %>%   addPolygons(color=~pal(poly@data$Pop_2009),fillOpacity  = 0.8,weight=1)%>%   addLegend(pal = pal, values = poly@data$Pop_2009,             title = "2009年人口数量（万人）</br>默认，右上角")%>%   addLegend(pal = pal, values = poly@data$Pop_2009,bins=5,position="bottomleft",             title = "2009年人口数量（万人）</br>图例分级五级，左下角")%>%   addLegend(pal = pal, values = poly@data$Pop_2009,position="topleft",           labFormat = labelFormat(suffix = " 万"),           title = "2009年人口数量（万人）</br>加单位，左上角")%>%   addLegend(pal = pal, values = poly@data$Pop_2009,position="bottomright",           opacity = 1,           title = "2009年人口数量（万人）</br>色带条不透明，右下角")

 

 好了，到此为止，leafletfor R 的所有方法都介绍完了，正常进行地图可视化，那是绝对够用了……当然如果有哪位同学要在上面玩一些黑科技（想做出一些特别的效果来），也欢迎和我交流。 上面所有的脚本和数据，我将在整理完成之后，放出下载。