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

这个例子的目的是输入城市名称或地名，然后找出它的经纬度值，并直观地显示它的路线流向和周围的地图显示

地址列表数据:

山西省太原市小店区秦县北街77号

贵州省贵阳市云岩区书香门第乙座3单元

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

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

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

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

北京

北京市通州区

北京的，北京的

天津武清区

1.生成可视化地图

库(devtools)

# install _ github(' lchiffon/REmap ')

库(REmap)

City_vec = c("北京"、"上海"、"广州")

get_city_coord("上海")

get _ geo _ location(city _ vec)

set.seed( 125)

原点= rep(“杭州”，5)

目的地= c('台州'，'丽水'，'金华'，'温州'，'乡宁')

dat = data.frame(起点，终点)

Out= remap(dat，title = "REmap实例数据"，副标题= "主题:黑暗")

绘制(出)

地图名称(“杭州”)

结果:

2.根据地名，生成其经纬度

# # # # # # # # # # # # # #使用百度API # # # # # # # # # # # # # # #。

# # # # # # # # # # # #地理位置→经纬度# # # # # # # # # # # # # # #。

##############################################

库(REmap)

#库(RCurl)

#library(rjson)

#get_city_coord('杭州')

#导入地址列表

地址<。- read.table( 'address_list.txt '，header = F，col.names=c( 'address ')，as.is = c( 1))

地址<。-地址$address#被转换为矢量格式以供回收

标题(地址)

#创建备用矢量，包括空矢量和百度地图api键

百度_ lng & lt- c()

百度_ lat & lt- c()

ak <。V5dWHCeG7jdMhcElzW自带的apiAK'#百度地图API的秘钥，需要自己申请

#加载包

库(rjson)

库(RCurl)

位置<。-地址[ 1]

#循环解析过程

for(地址中的位置){

#生成规则的网址(具体参数请参考地理编码应用编程接口文档)

url <。-粘贴(' http://API . map . Baidu . com/geocoder/v2/？ak= '，ak '，& amp回调= renderOption & amp输出= json & amp地址= '，位置，sep= ')

#使用URLencode()转换为可解析的网址

url _ string & lt- URLencode(url)

#通过readLines读取网址，并以JSON格式解析结果

json<。- readLines(url_string，warn=F)

geo <。- fromJSON( substr(json，regexpr('('，json)+ 1，nchar(json)- 1))

#从分析结果中提取纬度和经度

lng<。-geo $result$location$lng

lat<。-地理$结果$位置$稍后

#存储在已构建的字段中。

百度_液化天然气<。- c(百度_lng，lng)

百度_ lat & lt- c(百度_lat，lat)

}

#整理结果

结果<。- data.frame(address=address，经度=baidu_lng，纬度=baidu_lat)

结果

结果(对应地址列表文件位置的经纬度):

3.获取输入名称的百度地图显示

库(devtools)

# install _ github(' bad bye/Baidu map ')

# install _ github(' lchiffon/REmap ')

图书馆(baidumap)

库(ggmap)

##加载所需的包: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('浙江中医药大学'，宽= 600，高= 600，缩放= 18，比例= 2，消息传递= false)

Ggmap(q) #绘制地图

结果:

华清莲日常点滴，方便自己，成长他人！！！

——————————————