查看原文
其他

R实现地理位置与经纬度相互转换

2017-01-10 华青莲 R语言中文社区

本实例要实现目标通过输入城市名或者地名,然后找出其经度纬度值,以及通过可视化展现其线路流向以及周边地图展示

address_list数据:

山西省太原市小店区亲贤北街77号

贵州省贵阳市云岩区书香门第B栋3单元

北京市通州区神树商业街168号

贵州省贵阳市南明区兴关路51

北京市东城区长巷二条乙5号

山西省太原市杏花岭区北肖墙12号

北京

北京市通州区

北京市的的的的的

天津市武清区


1、生成可视化地图

library(devtools)
#install_github('lchiffon/REmap')
library(REmap)
city_vec = c("北京","Shanghai","广州")
get_city_coord("Shanghai")
get_geo_position (city_vec)
set.seed(125)
origin = rep("杭州",5)
destination = c('台州','丽水','金华','温州','乡宁')
dat = data.frame(origin,destination)
out = remap(dat,title = "REmap实例数据",subtitle = "theme:Dark")
plot(out)
mapNames("杭州")

结果:



2、根据地名,生成其经纬度

#############使用百度API######################
############地理位置→经纬度###################
##############################################
library(REmap)
#library(RCurl)
#library(rjson)
#get_city_coord('杭州')
#导入地址列表
address <- read.table('address_list.txt',header = F, col.names=c('address'),as.is = c(1))
address <- address$address   #转化为向量格式,备for循环使用
head(address)
#建立备用向量,包括空向量及百度地图api秘钥
baidu_lng <- c()
baidu_lat <- c()
ak <- 'V5dWHCeG7jdMhcElzW自己的APIAK'     #百度地图api的秘钥,需自己申请
#加载包
library(rjson)
library(RCurl)
location <-address[1]
#循环解析过程
for (location in address) {
  #生成规则的url地址(具体参数可参考Geocoding API文档)
  url <- paste('http://api.map.baidu.com/geocoder/v2/?ak=',ak,'&callback=renderOption&output=json&address=',location,sep='')
  #利用URLencode()转换为可解析的URL地址
  url_string <- URLencode(url)
  #通过readLines读取URL地址,并解析JSON格式的结果
  json<- readLines(url_string, warn=F)
  geo <- fromJSON(substr(json,regexpr('\\(',json)+1,nchar(json)-1))
  #在解析结果中提取经纬度
  lng<-geo$result$location$lng
  lat<-geo$result$location$lat
  #存储到已经建好的字段中
  baidu_lng <- c(baidu_lng,lng)
  baidu_lat <- c(baidu_lat,lat)
}
#整理结果
result <- data.frame(address=address,longitude=baidu_lng,latitude=baidu_lat)
result

结果(对应address_list文件位置经纬度):



3、获取输入名字的百度地图展现

library(devtools)
#install_github('badbye/baidumap')
#install_github('lchiffon/REmap')
library(baidumap)
library(ggmap)
## Loading required package: ggplot2
#获取浙江中医药大学的地图信息
####http://www.gpsspg.com/maps.htm
q <- getBaiduMap(c(120.2230040000,30.2154520000), width=600, height=600, zoom=18, scale = 2, messaging=FALSE)
ggmap(q) #绘制地图
q <- getBaiduMap('浙江中医药大学', width=600, height=600, zoom=18, scale = 2, messaging=FALSE)
ggmap(q) #绘制地图

结果:


华青莲日常点滴,方便自己,成长他人!!!



微信回复关键字即可学习

回复 R              R语言快速入门免费视频 
回复 统计          统计方法及其在R中的实现
回复 用户画像   民生银行客户画像搭建与应用 
回复 大数据      大数据系列免费视频教程
回复 可视化      利用R语言做数据可视化
回复 数据挖掘   数据挖掘算法原理解释与应用
回复 机器学习   R&Python机器学习入门 

您可能也对以下帖子感兴趣

文章有问题?点此查看未经处理的缓存