R语言实现根据经纬度得出地理位置

有没有大神能用R语言根据经纬度，得出每个经纬度对应的地理位置？（数据集有上万条观测）
跪求分享代码，要求能切实可运行的。

1. AK <- "你自己的AK"
   
2. lng <- "116.32212329"
   
3. lat <- "39.997986202"
   
4. 
   
5. location <- paste(lat,",",lng,sep="")
   
6. url <- paste("http://api.map.baidu.com/geocoder/v2/?callback=renderReverse&location=",location,"&output=json&pois=1&ak=",AK, sep = "")
   
7. url_string <- URLencode(url)
   
8. connect <- url(url_string)
   
9. temp_geo <- fromJSON(substr(paste(readLines(connect,encoding = "UTF-8",warn = F),collapse = ""),
   
10.                             30,nchar(paste(readLines(connect,encoding = "UTF-8",warn = F),collapse = ""))-1))
    
11. temp_geo

复制代码

这样就行了

http://www.ppvke.com/Blog/archives/37944

这段代码是不可行的，请问你还有更好的方法吗？或者修改一下这个代码？
很感激你的回复！

1. AK <- "你自己的AK"
   
2. lng <- "116.32212329"
   
3. lat <- "39.997986202"
   
4. 
   
5. location <- paste(lat,",",lng,sep="")
   
6. url <- paste("http://api.map.baidu.com/geocoder/v2/?callback=renderReverse&location=",location,"&output=json&pois=1&ak=",AK, sep = "")
   
7. url_string <- URLencode(url)
   
8. connect <- url(url_string)
   
9. temp_geo <- fromJSON(substr(paste(readLines(connect,encoding = "UTF-8",warn = F),collapse = ""),
   
10.                             30,nchar(paste(readLines(connect,encoding = "UTF-8",warn = F),collapse = ""))-1))
    
11. temp_geo

复制代码

radar000 发表于 2017-3-23 17:00
这样就行了

谢谢，感激不尽！

radar000 发表于 2017-3-23 17:00
这样就行了

for (location in lat_lng) {
#生成规则的url地址(具体参数可参考Geocoding API文档)
  url <- paste("http://api.map.baidu.com/geocoder/v2/?callback=renderReverse&location=",location,"&output=json&pois=1&ak=",AK, sep = "")
  url_string <- URLencode(url)
  connect <- url(url_string)
  temp_geo <- fromJSON(substr(paste(readLines(connect,encoding = "UTF-8",warn = F),collapse = ""),
                            30,nchar(paste(readLines(connect,encoding = "UTF-8",warn = F),collapse = ""))-1))
  address<-temp_geo$result$formatted_address
  baidu_address <- c(baidu_address,address)
}
result <- data.frame(lat_lng, baidu_address)

result
                             lat_lng                baidu_address
1  37.8051643371582,112.564300537109 山西省太原市小店区产业路21号
2   37.899299621582,112.557022094727 山西省太原市小店区产业路21号
3  37.8539428710938,112.593185424805 山西省太原市小店区产业路21号
4  37.8113784790039,112.587127685547 山西省太原市小店区产业路21号
5  37.8410720825195,112.540687561035 山西省太原市小店区产业路21号
6  37.8662376403809,112.583618164062 山西省太原市小店区产业路21号
7   37.870849609375,112.631248474121 山西省太原市小店区产业路21号
8  37.8324890136719,112.560577392578 山西省太原市小店区产业路21号
9   37.757740020752,112.574089050293 山西省太原市小店区产业路21号
10 37.7993812561035,112.509735107422 山西省太原市小店区产业路21号

请问我用这样for循环想一下子导出多个经纬度的地址，为什么最后只有一个地址？

radar000 发表于 2017-3-23 17:00
这样就行了

for (location in lat_lng) {
#生成规则的url地址(具体参数可参考Geocoding API文档)
  url <- paste("http://api.map.baidu.com/geocoder/v2/?callback=renderReverse&location=",location,"&output=json&pois=1&ak=",AK, sep = "")
  url_string <- URLencode(url)
  connect <- url(url_string)
  temp_geo <- fromJSON(substr(paste(readLines(connect,encoding = "UTF-8",warn = F),collapse = ""),
                            30,nchar(paste(readLines(connect,encoding = "UTF-8",warn = F),collapse = ""))-1))
  address<-temp_geo$result$formatted_address
  baidu_address <- c(baidu_address,address)
}
result <- data.frame(lat_lng, baidu_address)


请问，我想用for 循环导出多个经纬度的具体地址，但是为什么最后只有一个地址？

> result
                             lat_lng                baidu_address
1  37.8051643371582,112.564300537109 山西省太原市小店区产业路21号
2   37.899299621582,112.557022094727 山西省太原市小店区产业路21号
3  37.8539428710938,112.593185424805 山西省太原市小店区产业路21号
4  37.8113784790039,112.587127685547 山西省太原市小店区产业路21号
5  37.8410720825195,112.540687561035 山西省太原市小店区产业路21号
6  37.8662376403809,112.583618164062 山西省太原市小店区产业路21号
7   37.870849609375,112.631248474121 山西省太原市小店区产业路21号
8  37.8324890136719,112.560577392578 山西省太原市小店区产业路21号
9   37.757740020752,112.574089050293 山西省太原市小店区产业路21号
10 37.7993812561035,112.509735107422 山西省太原市小店区产业路21号

新技能get

radar000 发表于 2017-3-22 21:01
这样就行了

可以加你好友吗？求指导。拜托拜托。