|
八、栅格数据
栅格数据是当前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"),[email protected]$Pop_2009)
leaflet(poly) %>% addTiles() %>%
addPolygons(color=~pal([email protected]$Pop_2009),fillOpacity = 0.8,weight=1)%>%
addLegend(pal = pal, values = [email protected]$Pop_2009,position="bottomright",
title = "2009年人口数量（万人）")
pal <- colorNumeric("Greens",[email protected]$Pop_2009)
leaflet(poly) %>% addTiles() %>%
addPolygons(color=~pal([email protected]$Pop_2009),fillOpacity = 0.8,weight=1)%>%
addLegend(pal = pal, values = [email protected]$Pop_2009,position="bottomright",
title = "2009年人口数量（万人）")
pal <- colorBin(c("darkgreen", "yellow", "orangered"),[email protected]$Pop_2009,10)
leaflet(poly) %>% addTiles() %>%
addPolygons(color=~pal([email protected]$Pop_2009),fillOpacity = 0.8,weight=1)%>%
addLegend(pal = pal, values = [email protected]$Pop_2009,position="bottomright",
title = "2009年人口数量分级")
pal <- colorBin("Greens",[email protected]$Pop_2009,10)
leaflet(poly) %>% addTiles() %>%
addPolygons(color=~pal([email protected]$Pop_2009),fillOpacity = 0.8,weight=1)%>%
addLegend(pal = pal, values = [email protected]$Pop_2009,position="bottomright",
title = "2009年人口数量分级")
|
十、图例：
图例是对地图可视化之后相关颜色或者符号的映射说明——当然，在leaflet里面，主要是针对颜色创建图例。
主要是方法如下：
addLegend（地图，属性，颜色列表，标记值……）
具体的属性方法，大家可以查阅帮助文档，下面给出一个示例：
|
pal <- colorNumeric("Greens",[email protected]$Pop_2009)
leaflet(poly) %>% addTiles() %>%
addPolygons(color=~pal([email protected]$Pop_2009),fillOpacity = 0.8,weight=1)%>%
addLegend(pal = pal, values = [email protected]$Pop_2009,
title = "2009年人口数量（万人）</br>默认，右上角")%>%
addLegend(pal = pal, values = [email protected]$Pop_2009,bins=5,position="bottomleft",
title = "2009年人口数量（万人）</br>图例分级五级，左下角")%>%
addLegend(pal = pal, values = [email protected]$Pop_2009,position="topleft",
labFormat = labelFormat(suffix = " 万"),
title = "2009年人口数量（万人）</br>加单位，左上角")%>%
addLegend(pal = pal, values = [email protected]$Pop_2009,position="bottomright",
opacity = 1,
title = "2009年人口数量（万人）</br>色带条不透明，右下角")
|
好了，到此为止，leafletfor R 的所有方法都介绍完了，正常进行地图可视化，那是绝对够用了……当然如果有哪位同学要在上面玩一些黑科技（想做出一些特别的效果来），也欢迎和我交流。
上面所有的脚本和数据，我将在整理完成之后，放出下载。
|
客服电话
APP下载
官方微信
请发表评论