我正在尝试使用ggmap(版本2.4)中的以下(基本)代码对意大利的不同IATA机场代码进行地理编码

#list of all IATA codes
geo_apt <- c("AOI", "BGY", "BLQ", "BRI", "CTA", "FCO", "LIN", "MXP", "NAP",
"PMF", "PSA", "PSR", "RMI", "TRN", "VCE", "VRN")

#preparing an empty dataframe to store the geocodes
apt_geo <- data.frame(IATA=rep(NA,16), lon=rep(NA,16), lat=rep(NA,16))

#geocoding the codes
for (i in seq_along(geo_apt)) {
  apt_geo[i,1] <- geo_apt[i]
  apt_geo[i,2] <- (geocode(paste(geo_apt[i],"airport")))[1]
  apt_geo[i,3] <- (geocode(paste(geo_apt[i],"airport")))[2]
}


并且geocodeggmap函数在所有这些代码(“ PSR”除外)上都可以正常工作

   IATA        lon      lat
1   AOI  13.363752 43.61654
2   BGY   9.703631 45.66957
3   BLQ  11.287859 44.53452
4   BRI  16.765202 41.13751
5   CTA  15.065775 37.46730
6   FCO  12.246238 41.79989
7   LIN   9.276308 45.45218
8   MXP   8.725531 45.63006
9   NAP  14.286579 40.88299
10  PMF  10.295935 44.82326
11  PSA  10.397884 43.68908
12  PSR -81.117259 33.94855  #<- doens't work
13  RMI  12.618819 44.02289
14  TRN   7.647867 45.19654
15  VCE  12.339771 45.50506
16  VRN  10.890141 45.40000


我尝试使用revgeocode,这些坐标对应于以下地址:

revgeocode(as.numeric(apt_geo[12,2:3]))
#Information from URL : http://maps.googleapis.com/maps/api/geocode/json?latlng=33.948545,-81.1172588&sensor=false
[1] "Kentucky Avenue, West Columbia, SC 29170, USA"


相反,如果我去Google地图,它就可以很好地工作:

r - 地理编码结果与Google map 不同-LMLPHP

有人对这个看似奇怪的现象有任何线索吗?

编辑

遵循以下注释中的一个建议,我尝试在版本2.4上再次使用geocode(italy PSR airport),而不是抛出更准确的结果甚至相同的结果,这是我得到的警告:

geocode("italy PSR airport")
  lon lat
1  NA  NA
Warning message:
geocode failed with status ZERO_RESULTS, location = "italy PSR airport"


而尝试airport PSR的坐标甚至与PSR airport的坐标不同(至少这次是一个实际的机场,尽管其IATA代码是LEX而不是PSR)。

revgeocode(as.numeric(geocode("airport PSR")))
Information from URL : http://maps.googleapis.com/maps/api/geocode/json?latlng=38.0381454,-84.5970727&sensor=false
[1] "3895 Terminal Drive, Lexington, KY 40510, USA"

最佳答案

整个问题是可能的duplicate

尽管如此,我不知道API和Google地图使用不同数据集的原因...

关于r - 地理编码结果与Google map 不同,我们在Stack Overflow上找到一个类似的问题:https://stackoverflow.com/questions/31720359/

10-12 01:29