将“rgb”图例添加到R传单热图

时间:2017-12-28 17:35:05

标签: r leaflet data-visualization heatmap

我使用传单(特别是 leaflet.extras 包中的 addHeatmap()命令制作了一些交互式热图,并且闪亮即可。创建了所需的地图后,我想添加一个图例。

我感兴趣的是一个“rgb”图例,它基于由纯long / lat坐标中的 addHeatmap()推导出的密度值。 我需要的是这样的地图 - https://www.patrick-wied.at/static/heatmapjs/example-legend-tooltip.html - 遗憾的是我不了解JS并且无法用R语言重写这段代码/为我的问题包含正确的JS代码片段。

我到目前为止尝试的是 addLegend()命令,它没有给出所需的结果,因为在这种情况下我需要指定一个为其准备图例的变量。我还试图从创建的传单对象中提取颜色范围和指定值,但是没有成功。

这是指向运行可重现示例的完整数据的链接: https://drive.google.com/file/d/1h3jL_PU6DGTtdIWBK02Tt37R7IB2ArH9/view

这是前20条记录:

structure(list(latitude = c(30.309522, 30.24429616, 30.30038194, 
                            30.27752338, 30.23294081, 30.23038507, 
                            30.34285933, 30.24962237, 30.26594744, 
                            30.20515821, 30.22363485, 30.2759184, 
                            30.28283226, 30.33816909, 30.26611565, 
                            30.18835401, 30.26704789, 30.27456699, 
                            30.19237135, 30.1925213), 
              longitude = c(-97.73171047, -97.77446858, -97.77885789, 
                            -97.71919076, -97.58937812, -97.76581095, 
                            -97.73598704, -97.72215443, -97.74144275, 
                            -97.8782895, -97.78329845, -97.71321066, 
                            -97.70820152, -97.82413058, -97.7327258, 
                            -97.81606795, -97.68989589, -97.7580592, 
                            -97.7816127, -97.73138523)), 
                .Names =  c("latitude", "longitude"), row.names = 
                                c(NA, 20L), class = "data.frame")      

这是一个示例代码,我想通过上述功能进行扩展:

library(magrittr)
library(leaflet)
library(leaflet.extras)

data <- read.csv('DATA.csv')
leaflet(data) %>%
  addTiles(group="OSM") %>%
  addHeatmap(group="heat", lng = ~longitude, lat = ~latitude, max=.5, blur = 60)

这是该代码的结果(在整个数据集上):

https://i.stack.imgur.com/6VFNC.jpg

总结一下我想做的事情:基于这样的图片,我想提取绘制颜色的范围以及分配给它们的值,并使用该信息绘制图例。

我有什么遗失的吗?这看起来很简单,但过去几个小时我一直在努力寻找解决方案。 在此先感谢您的帮助!

编辑:扩展了可重复的示例。

1 个答案:

答案 0 :(得分:0)

您的样本数据没有任何值可以实际映射密度叠加层。

您可以使用colorBin()指定垃圾箱的数量,然后使用pal函数指定那些垃圾箱。您可以根据需要在data_values发行版中设置不同的垃圾箱。 colorBin()的帮助部分可帮助您确定所需的正确参数。

bins <- c(0,1,2,3,4)
pal <- colorBin("Spectral", domain = data_value, bins = bins, na.color = "transparent")

m <-leaflet() %>% 
      addTiles() %>% 
      addHeatmap(lng= long_cords, lat = lat_cords, intensity = data_value,
                 blur = 20, max = 400, radius = 15, cellSize = 3) %>%     
      addLegend(pal = pal, values = data_value,
                title="Heat map legend")

您必须使用addHeatmap参数来获得正确的透明度和密度设置。