接着上次的GIS包继续说

不仅仅是geo数据 ggmap还可以绘制任何地方的地图

下面通过绘制近期中国地震地缘分布来练习一下使用
不多说上代码

library(ggmap)
library(mapproj)
## Google啊Google给我China的地图数据吧
map <- get_map(location = 'China', zoom = 4)
ggmap(map)

tables <- readHTMLTable(webpage,stringsAsFactors = FALSE)
raw <- tables[[6]]
data <- raw[-1,c(1,3,4,6)]
names(data) <- c('date','lan','lon','earth_size')
pattern <- '[0-9]+'
gregout <- gregexpr(pattern, data$earth_size)
earth_size <- substr(data$earth_size[1],3,5)
for(i in 2:length(data$earth_size)){
  earth_size <- rbind(earth_size, substr(data$earth_size[i],3,5) )
}
data$earth_size <- as.numeric(earth_size)
data$lan <- as.numeric(data$lan)
data$lon <- as.numeric(data$lon)
data$date <- as.Date(data$date,  "%Y-%m-%d")

#绘制地震点图
mymap <- get_googlemap(center = 'china', zoom=4,maptype='terrain')
ggmap(mymap,extent='device')+
  geom_point(data=data,aes(x=lon,y=lan, size = earth_size),colour = 'red',alpha=0.7, na.rm = T)

图中红点的大小表示地震的级数, 可以看到的是貌似都是小地震, 4级以下偏多, 主要分布在西部貌似, 中部比如武汉基本没啥地震..感觉我们这还是比较安全的= =地震比较多的还是那些盆地比如四川之类的

感觉还是挺有意思的 还有爬虫的东西下次再补吧

发现一个css selector神器啊..都不用自己去找网页源码直接给出节点号了.爬虫无忧啊..发现R真是神器, 比python确实好用些, 就是感觉不好写太大的程序.