第四章通过REmap包完成基于map分布图示例，前面提到REmap基于Echart2.0, 一方面在移动终端适应效果差，另一方面REmap提供的热力图仅支持全国及省市大版块map，基于上面的原因，参考 http://echarts.baidu.com/examples.html#chart-type-heatmap ，利用第四章清洗后的数据结合Echart3热力图控件开发完整可缩放地图的热力分布。

目录：

  • 效果图
  • 在线构建Echart3脚本包
  • R热力图页面模板
  • R脚本

效果图：

在线构建Echart3脚本包:

  • 在浏览器打开 http://echarts.baidu.com/builder.html， 选择 "下载" -->"在线构建", 在打开的页面选择要打包的图表，点击 “下载”
  • 简单起见，你可以选择所有控件后打包， build 完后保存为 “echarts3.min.js”

R热力图页面模板：

  • 参考Echart 热力图文档，在R脚本定义模板及Class, 保存命名：DemoTemp.R
  • 示例代码如下： 
    setClass("my.map",
    
         representation(
    
           id = "character",
    
           option = "character",
    
           content = "character"
    
         ))

my.mapH = function(data) {

  if (.Platform$OS.type == "windows") {
    
    locate = Sys.getlocale("LC_CTYPE")
    
    Sys.setlocale("LC_CTYPE", "eng")
    
  }

  if (!is.data.frame(data)) {
    
    stop("Map data should be a data frame.")
    
  }

  if (ncol(data) < 3 | nrow(data) == 0) {
    
    stop("Data should have at least 3 columns and 1 row")
    
  }

  if (!is.numeric(data[1, 1]) |
    
      !is.numeric(data[1, 2]) |
    
      !is.numeric(data[1, 3])) {
    
    data[, 1] = as.numeric(data[, 1])
    
    data[, 2] = as.numeric(data[, 2])
    
    data[, 3] = as.numeric(data[, 3])
    
  }

  heatdata = apply(data, 1, function(x) {
    
    out = sprintf('[%s,%s,%s]', x[1], x[2], x[3])
    
    return(out)
    
  })

  heatdata = paste(heatdata, collapse = ",\n")

  ### write remap object
    
  output = new("my.map")
    
  output@id = paste('ID', format(Sys.time(), "%Y%m%d%H%M%S"),
    
                    round(proc.time()[3] * 100), sep = "_")

  output@option = html.data.H$option
    
  head = html.data.H$head
    
  foot = html.data.H$foot

  if (.Platform$OS.type == "windows") {
    
    Sys.setlocale("LC_CTYPE", "chs")
    
  }

  output@option = sub("forChange", "һ", output@option)

  output@option = sub("heatMapData",
    
                      heatdata, output@option)

  ## optionNameData
    
  output@option = sub("optionNameData",
    
                      paste0("option", output@id), output@option)
    
  outputFoot = sub("optionNameData",
    
                   paste0("option", output@id), foot)

  output@option = strsplit(output@option, "kkkmmm")[[1]][2]
    
  output@content = paste(head, output@option, outputFoot, sep = "\n")

  if (.Platform$OS.type == "windows") {
    
    Sys.setlocale("LC_CTYPE", locate)
    
  }
    
  return(output)
    
}

html.data.H = list(
    
  head = "<html>
    
  <head>
    
  <meta charset=\"utf-8\">
    
  <style type=\"text/css\">
    
  body {
    
  margin: 0;
    
  }
    
  #main {
    
  height: 100%;
    
  }
    
  </style>
    
  </head>
    
  <body>
    
  <div id=\"main\"></div>
    
  <script src=\"./js/echarts.js\"></script>
    
  <script src=\"./js/bmap.js\"></script>
    
  <script src = \"http://api.map.baidu.com/api?v=2.0&ak=密钥\"></script>
    
  <script>
    
  var myChart = echarts.init(document.getElementById(\"main\"));

  ",
    
  option = "forChangekkkmmm
    
  var heatData = [heatMapData];
    
  var optionNameData =
    
  {
    
  animation: false,
    
  bmap: {
    
  center: [103.855096, 36.056805],
    
  zoom: 5,
    
  min: 4,
    
  max:12,
    
  roam: true
    
  },
    
  visualMap: {
    
  show: true,
    
  top: 'top',
    
  min: 1,
    
  max: 50,
    
  seriesIndex: 0,
    
  calculable: true,
    
  inRange: {
    
  color: ['blue', 'blue', 'green', 'yellow', 'red']
    
  }
    
  },
    
  series: [{
    
  type: 'heatmap',
    
  coordinateSystem: 'bmap',
    
  data: heatData,
    
  pointSize: 6,
    
  blurSize: 8
    
  }]
    
  }",
    
  foot = ";
    
  myChart.setOption(optionNameData);
    
  </script>
    
  </body>
    
  </html>"
    
  )

R脚本：

  • R脚本代码，demo.rda 为清洗后保存的数据，示例

    run <- function(...) {
    
  map_name <- my.writeMapH()
    
  out(map_name)
    
}

my.writeMapH <- function() {

  path = c("/var/www/html")
    
  file_name = paste0("3.0/Demo")
    
  full_path = paste0(path, "/", file_name, ".html")
    
  if (file.exists(full_path)) {
    
    return(file_name)
    
  }

  source("/var/FastRWeb/web.R/DemoTemp.R")
    
  tmp <- load("/var/FastRWeb/web.R/demo.rda")
    
  data <- data.frame(pdata$lon,pdata$lat,c(1))
    
  out <- my.mapH(data)

  writeLines(out@content, full_path, useBytes = T)
    
  return(file_name);
    
}
  • 通过FastRWeb框架调用R脚本成功后，返回的是在/var/www/html目录下生成的文件名
  • client再次发起请求，调用html文件
  • 注意： 在/var/www/html 目录下部署脚本引用的 js 文件

R(七): R开发实例-map热力图的更多相关文章

  1. R(四): R开发实例-map分布图

    前几章对R语言的运行原理.基本语法.数据类型.环境部署等基础知识作了简单介绍,本节将结合具体案例进行验证测试. 案例场景:从互联网下载全国三甲医院数据,以地图作为背景,展现各医院在地图上的分布图.全国 ...

  2. 工业级GBDT算法︱微软开源 的LightGBM（R包正在开发....）

    看完一篇介绍文章后,第一个直觉就是这算法已经配得上工业级属性.日前看到微软已经公开了这一算法,而且已经发开python版本,本人觉得等hadoop+Spark这些平台配齐之后,就可以大规模宣传啦~如果 ...

  3. C#软件开发实例.私人订制自己的屏幕截图工具（七）加入放大镜的功能

    上一篇:C#软件开发实例.私人订制自己的屏幕截图工具(六)加入配置管理功能 因为截图时可能须要精确截取某一部分,所以须要放大镜的功能,这样截取的时候才更easy定位截图的位置. 加入PictureBo ...

  4. Android音乐播放器的开发实例

    本文将引导大家做一个音乐播放器,在做这个Android开发实例的过程中,能够帮助大家进一步熟悉和掌握学过的ListView和其他一些组件.为了有更好的学习效果,其中很多功能我们手动实现,例如音乐播放的 ...

  5. Android开发实例之miniTwitter登录界面的实现

    原文: http://www.jizhuomi.com/android/example/134.html 本文要演示的Android开发实例是如何完成一个Android中的miniTwitter登录界 ...

  6. Android开发实例详解之IMF（Android SDK Sample—SoftKeyboard）

    本博前面的文章介绍了Android开发环境的搭建和模拟器的常用操作.本次,将以Android Sample中经典的SoftKeyboard项目为例,详细解析Android上一个小型项目的开发过程和注意 ...

  7. 持久化框架Hibernate 开发实例（一）

    1 Hibernate简介 Hibernate框架是一个非常流行的持久化框架,其中在web开发中占据了非常重要的地位, Hibernate作为Web应用的底层,实现了对数据库操作的封装.HIberna ...

  8. [转贴]JAVA:RESTLET开发实例（二）使用Component、Application的REST服务

    上一篇文章,我们介绍了基于JAX-RS的REST服务,本篇文章我们介绍不基于JAX-RS的模式.JAX-RS其实就是一个简单的 Application服务.和我们接下来介绍的Application基本 ...

  9. Android 应用开发实例之情景模式

    2013-07-01 Android 应用开发实例 1. 情景模式 使用TabHost来实现主界面的布局. 设置一组RadioButton来切换不同的情景模式. 对比普通情景模式,定时情景模式需要加上 ...

随机推荐

  1. xmind的第三天笔记

  2. jQuery显示隐藏动态效果的几种写法

    <script type="text/javascript"> $(document).ready(function() {     /*$("#test1& ...

  3. PHP 将json的stdClass Object转成数组array

    PHP和JS通讯通常都用json,但是PHP要用json的数据,通过json_decode转出来的数组并不是标准的array,所以需要用这个函数进行转换. function object_array( ...

  4. POJ 3114 Countries in War（强连通）（缩点）（最短路）

                                    Countries in War Time Limit: 1000MS   Memory Limit: 65536K Total Sub ...

  5. 英语语法最终珍藏版笔记-6“情态动词＋have+ done”的含义

    “情态动词+have+ done”的含义 1.Must have done的含义.“must have+过去分词”表示对过去的推测,意思是“一定已经,想必已经,准是已经….”,只用于肯定句中.例如: ...

  6. windows 10安装framework 3.5失败的解决方案

    装了两次win 10,全都因为没法安装framework 3.5,用不了老版本的开发环境,又换回了win7. 网上有两种解决方案: a,通过iso安装. 可是拜托,我的系统都是用ghost版本安装的, ...

  7. shell下的作业管理[转]

    作业管理 举例来说,我们在登陆 bash 后, 想要一边复制文件.一边进行数据搜寻.一边进行编译,还可以一边进行 vi 程序撰写! 当然我们可以重复登陆那六个文字介面的终端机环境中,不过,能不能在一个 ...

  8. 类的扩展－－类目和延展－－ios

    person+money.h  这是类目类 #import "Person.h" //这是扩展person类的接口类,独立一个文件 @interface Person (Money ...

  9. android语言适配

    虽然我们现在做的app基本只在国内使用,所以只需要中文,但是我们要把眼光放长远一点,以后我们的APP发展到了国外呢,所以我们就要做语言适配了   很简单 在res目录下,右键  New Resourc ...

  10. PgSQL · 追根究底 · WAL日志空间的意外增长

    问题出现 我们在线上巡检中发现,一个实例的pg_xlog目录,增长到4G,很是疑惑.刚开始怀疑是日志归档过慢,日志堆积在pg_xlog目录下面,未被清除导致.于是检查归档目录下的文件,内容如下.但发现 ...