我正在尝试使用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]
}
并且
geocode
的ggmap
函数在所有这些代码(“ 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地图,它就可以很好地工作:
有人对这个看似奇怪的现象有任何线索吗?
编辑
遵循以下注释中的一个建议,我尝试在版本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/