Я обрабатываю данные объекта в R с дополнительной геоинформацией, такой как название улицы, номер улицы, почтовый индекс и город из Германии. Я хочу использовать функцию для получения информации о долготе и широте для конкретного адреса и обновления моего фрейма данных. С долготой и широтой для каждого объекта из этого фрейма данных я создаю geom_point в qmap.
Любая помощь или информация о возможном решении в R или о том, где можно найти дополнительную информацию, приветствуются. Себастьян
Вы можете использовать Google Maps API для обратного геокодирования, но скорость будет ограничена довольно быстро. Уже есть пакет R, который использует этот API Google, но я написал эту функцию для обратного геокодирования с использованием API openstreetmap. Однако не забудьте указать свою электронную почту в звонке, чтобы они могли связаться с вами, если вы начнете нагружать их сервер.
FindGPS <- function(cityname) {
library("RJSONIO")
cityname <- gsub(" ","\\+",cityname) #add '+' character to URL
#Open Connection
SITE <- paste('http://nominatim.openstreetmap.org/
search?q=',cityname,'&format=json&email=
[email protected]', sep="")
con <- url(SITE)
data.json <- fromJSON(paste(readLines(con,warn=FALSE),
collapse=""))
close(con)
t<-length(data.json)
if(t>0) {
lat <- data.json[[1]]$lat
long <- data.json[[1]]$lon
latlong<-c(lat, long)
return (latlong)
}
}
