How to GeoCode a simple address using Data Science Toolbox

by

Like this:

library(httr)
library(rjson)

data <- paste0("[",paste(paste0("\"",dff$address,"\""),collapse=","),"]")
url  <- "http://www.datasciencetoolkit.org/street2coordinates"
response <- POST(url,body=data)
json     <- fromJSON(content(response,type="text"))
geocode  <- do.call(rbind,sapply(json,
                                 function(x) c(long=x$longitude,lat=x$latitude)))
geocode
#                                                long      lat
# San Francisco, California, United States -117.88536 35.18713
# Mobile, Alabama, United States            -88.10318 30.70114
# La Jolla, California, United States      -117.87645 33.85751
# Duarte, California, United States        -118.29866 33.78659
# Little Rock, Arkansas, United States      -91.20736 33.60892
# Tucson, Arizona, United States           -110.97087 32.21798
# Redwood City, California, United States  -117.88536 35.18713
# New Haven, Connecticut, United States     -72.92751 41.36571
# Berkeley, California, United States      -122.29673 37.86058
# Hartford, Connecticut, United States      -72.76356 41.78516
# Sacramento, California, United States    -121.55541 38.38046
# Encinitas, California, United States     -116.84605 33.01693
# Birmingham, Alabama, United States        -86.80190 33.45641
# Stanford, California, United States      -122.16750 37.42509
# Orange, California, United States        -117.85311 33.78780
# Los Angeles, California, United States   -117.88536 35.18713

This takes advantage of the POST interface to the street2coordinates API (documented here), which returns all the results in 1 request, rather than using multiple GET requests.

The absence of Phoenix seems to be a bug in the street2coordinates API. If you go the API demo page and try “Phoenix, Arizona, United States”, you get a null response. However, as your example shows, using their “Google-style Geocoder” does give a result for Phoenix. So here’s a solution using repeated GET requests. Note that this runs much slower.

geo.dsk <- function(addr){ # single address geocode with data sciences toolkit
  require(httr)
  require(rjson)
  url      <- "http://www.datasciencetoolkit.org/maps/api/geocode/json"
  response <- GET(url,query=list(sensor="FALSE",address=addr))
  json <- fromJSON(content(response,type="text"))
  loc  <- json['results'][[1]][[1]]$geometry$location
  return(c(address=addr,long=loc$lng, lat= loc$lat))
}
result <- do.call(rbind,lapply(as.character(dff$address),geo.dsk))
result <- data.frame(result)
result
#                                     address         long        lat
# 1        Birmingham, Alabama, United States   -86.801904  33.456412
# 2            Mobile, Alabama, United States   -88.103184  30.701142
# 3           Phoenix, Arizona, United States -112.0733333 33.4483333
# 4            Tucson, Arizona, United States  -110.970869  32.217975
# 5      Little Rock, Arkansas, United States   -91.207356  33.608922
# 6       Berkeley, California, United States   -122.29673  37.860576
# 7         Duarte, California, United States  -118.298662  33.786594
# 8      Encinitas, California, United States  -116.846046  33.016928
# 9       La Jolla, California, United States  -117.876447  33.857515
# 10   Los Angeles, California, United States  -117.885359  35.187133
# 11        Orange, California, United States  -117.853112  33.787795
# 12  Redwood City, California, United States  -117.885359  35.187133
# 13    Sacramento, California, United States  -121.555406  38.380456
# 14 San Francisco, California, United States  -117.885359  35.187133
# 15      Stanford, California, United States    -122.1675   37.42509
# 16     Hartford, Connecticut, United States   -72.763564   41.78516
# 17    New Haven, Connecticut, United States   -72.927507  41.365709

More Related Contents:

  1. rgdal package installation
  2. Plotting bar charts on map using ggplot2?
  3. How to put a geom_sf produced map on top of a ggmap produced raster
  4. R convert zipcode or lat/long to county
  5. ggplot centered names on a map
  6. Fixing maps library data for Pacific centred (0°-360° longitude) display
  7. Choropleth map in ggplot with polygons that have holes
  8. Improve centering county names ggplot & maps
  9. Plot coordinates on map
  10. Labeling center of map polygons in R ggplot
  11. Maps, ggplot2, fill by state is missing certain areas on the map
  12. Convert latitude and longitude coordinates to country name in R
  13. Group by multiple columns and sum other multiple columns
  14. Fitting a density curve to a histogram in R
  15. Conditionally change panel background with facet_grid?
  16. How to delete rows from a dataframe that contain n*NA
  17. How to order data by value within ggplot facets
  18. How to use R with Google Colaboratory?
  19. Replace negative values by zero
  20. run a for loop in parallel in R
  21. double integral in R
  22. R’s read.csv prepending 1st column name with junk text [duplicate]
  23. Delete duplicate rows in two columns simultaneously [duplicate]
  24. Should I use a data.frame or a matrix?
  25. add image in title page of rmarkdown pdf
  26. ggplot: Boxplot of multiple column values [duplicate]
  27. Gantt style time line plot (in base R)
  28. The perils of aligning plots in ggplot
  29. How to create a new variable in a data.frame based on a condition?
  30. What does the %% operator mean in R? [closed]

Leave a Comment