R语言实现莆田医院的可视化

摘要：R语言实现莆田医院的可视化...

摘要：Idea来自于这里，原文作者用python完成了数据爬取,可视化的过程,本文用R实现一遍...

基本的流程如下:

数据来源于凤凰资讯,由于数据隐藏在JS代码中,所以我们使用readLines读取,字符处理的方式获得数据.

raw = readLines("http://news.ifeng.com/mainland/special/ptxyy/", encoding = "UTF-8") rawHospital = raw[123 : 238] rawHospital[1]

[1] " [上海,上海市闵行区中医院<br/>上海华美美莱整形医院<br/>上海万众医院<br/>上海安真医院<br/>上海福华医院<br/>上海玛丽女子医院<br/>上海真爱女子医院<br/>上海心脏病医院<br/>上海远大心胸医院<br/>上海仁爱医院<br/>上海天大医疗美容医院<br/>上海沪申五官科医院<br/>上海博爱医院<br/>上海新虹桥医院<br/>上海九龙男子医院<br/>上海城市女子医院<br/>上海西郊骨科医院<br/>上海真美妇科医院<br/>上海南浦妇科医院<br/>上海虹桥医院<br/>上海健桥医院<br/>上海博爱医院<br/>上海长江医院<br/>解放军411医院<br/>上海阳光中医医院<br/>上海英港泌尿外科医院<br/>上海阿波罗男子医院<br/>上海玫瑰女子医院<br/>上海康新医院<br/>上海圣爱医院<br/>上海同德医院<br/>上海安平医院<br/>江东泌尿外科医院<br/>上海九州泌尿医院<br/>上海青城医院<br/>上海博大医院<br/>上海中亚医院,38],"

接下来,用一个函数提取每个rawHospital中的医院名,并整理成一个字符向量.

extFun = function(x) extFun(rawHospital[1]) hosList = sapply(rawHospital,extFun) hospital = c() for(i in 1:length(hosList)) hospital = append(hospital, hosList[[i]]) ## 获取经纬度 library(baidumap) blackHospital = getCoordinate(hospital, formatted = T) blackHospital = na.omit(blackHospital) plotdata = data.frame(lon = blackHospital[,1], lat = blackHospital[,2], city = rownames(blackHospital)) ## 绘制地图 library(REmap) remapB(markPointData = data.frame(plotdata$city), markPointTheme = markPointControl(symbol = "pin", effect=F, symbolSize = 5, color="red"), geoData = plotdata)

本文作者郎大为为雪晴数据网兼职讲师，博客地址http://www.chiffon.ninja

本文链接：http://www.xueqing.tv/cms/article/199