**R – Fronkonstin**, and kindly contributed to R-bloggers)

Hi Hillary, It’s Donald, would you like to have a beer with me in La Cabra Brewing, in Berwyn, Pensilvania? (Hypothetical utilization of The Meeting Point Locator)

Finding a place to have a drink with someone may become a difficult task. It is quite common that one of them does not want to move to the other’s territory. I am sure you have faced to this situation many times. With *The Meeting Point Locator* this will be no longer an issue, since it will give you a list of equidistant bars and coffees to any two given locations. Let’s see an example.

I do not know if Hillary Clinton and Donald Trump have met each other after the recent elections in United States, but the will probably do. Let’s suppose Hillary doesn’t want to go to The White House and that Donald prefers another place instead Hillary’s home. No problem at all. According to this, Hillary lives in **Chappaqua, New York **and Donald will live in **The White House, Washington **(although he supposedly won’t do full time as he announced recently). These two locations are the only input that *The Meeting Point Locator* needs to purpose equidistant places where having a drink. This is how it works:

- Generates a number of coordinates on the great circle which passes through the midpoint of the original locations and is orthogonal to the rhumb defined by them; the number of points depends on the distance between the original locations.
- Arranges these coordinates according to the distance to the original locations, from the nearest to the most distant.
- Depending also on the distance of the original locations, defines a radius to search around each point generated on the great circle (once calculated, this radius is constant for all searches).
- Starting from the nearest point, looks for a number of places (20 by default) to have a drink using the radius calculated previously. To do this, it calls to the Google Places API. Once the number of locations is reached, the proccess stops.

This map shows the places purposed for Hillary and Donald (blue points) as well as the original locations (red ones). You can make *zoom in* for details:

These are the 20 closest places to both of them:

Some other examples of the *utility* of *The Meeting Point Locator*:

**Pau Gasol**(who lives in San Antonio, Texas) and**Marc Gasol**(in Memphis, Tennessee) can meet into have a beer while watching a NBA match. It is 537 kilometers far from both of them.*The Draft Sports Bar,*in Leesville (Louisiana)**Bob Dylan**(who lives in Malibu, California) and**The Swedish Academy**(placed in Stockholm, Sweden) can*smooth things over*drinking a caipirinha in,*Bar São João*, in Tremedal (Brasil)*only*9.810 kilometers far from both of them.**Spiderman**(placed in New York City) and**Doraemon**(in Tokio, Japan) can meet in**Andreyevskaya, in Stroitel (Russia)**to have a have a hot drink. Since they are*superheroes*, they will cover the 9.810 kilometers of separation in no time at all.

I faced with two challenges to do this experiment: how to generate the orthogonal great circle from two given locations and how to define radius and number of points over this circle to do searchings. I will try to explain in depth both things in the future in another post.

You will find the code below. To make it work, do not forget to get your own key for Google Places API Web Service here. I hope this tool will be helpful for someone; if yes, do not hesitate to tell it to me.

library(httr) library(jsonlite) library(dplyr) library(ggmap) library(geosphere) library(DT) library(leaflet) # Write both addresses here (input) place1="Chappaqua, New York, United States of America" place2="The White House, Washington DC, United States of America" # Call to Google Maps API to obtain coordinates of previous addresses p1=geocode(place1, output = "latlon") p2=geocode(place2, output = "latlon") # To do searchings I need a radius radius=ifelse(distGeo(p1, p2)>1000000, 10000, ifelse(distGeo(p1, p2)>100000, 2500, 1000)) # And a number of points npoints=ifelse(distGeo(p1, p2)>1000000, 2002, ifelse(distGeo(p1, p2)>100000, 7991, 19744)) # Place here the Google Places API Key key="PLACE_YOUR_OWN_KEY_HERE" # Build the url to look for bars and cafes with the previous key url1="https://maps.googleapis.com/maps/api/place/nearbysearch/json?location=lat,lon&radius=" url2="&types=cafe|bar&key=" url=paste0(url1,radius,url2,key) # This is to obtain the great circle orthogonal to the rhumb defined by input locations # and which passes over the midpoint. I will explain this step in the future mid=midPoint(p1, p2) dist=distGeo(p1, p2) x=p1 y=p2 while(dist>1000000) { x=midPoint(mid, x) y=midPoint(mid, y) dist=distGeo(x, y) } bea=bearingRhumb(x, y) points=greatCircle(destPoint(p=mid, b=bea+90, d=1), mid, n=npoints) # Arrange the points dependning on the distance to the input locations data.frame(dist2p1=apply(points, 1, function (x) distGeo(p1, x)), dist2p2=apply(points, 1, function (x) distGeo(p2, x))) %>% mutate(order=apply(., 1, function(x) {max(x)})) %>% cbind(points) %>% arrange(order) -> points # Start searchings nlocs=0 # locations counter (by default stops when 20 is reached) niter=1 # iterations counter (if greater than number of points on the great circle, stops) results=data.frame() while(!(nlocs>=20 | niter>npoints)) { print(niter) url %>% gsub("lat", points[niter, 'lat'], .) %>% gsub("lon", points[niter, 'lon'], .) %>% GET %>% content("text") %>% fromJSON -> retrieve df=data.frame(lat=retrieve$results$geometry$location$lat, lng=retrieve$results$geometry$location$lng, name=retrieve$results$name, address=retrieve$results$vicinity) results %>% rbind(df)->results nlocs=nlocs+nrow(df) niter=niter+1 } # I prepare results to do a Data Table data.frame(dist2p1=apply(results, 1, function (x) round(distGeo(p1, c(as.numeric(x[2]), as.numeric(x[1])))/1000, digits=1)), dist2p2=apply(results, 1, function (x) round(distGeo(p2, c(as.numeric(x[2]), as.numeric(x[1])))/1000, digits=1))) %>% mutate(mx=apply(., 1, function(x) {max(x)})) %>% cbind(results) %>% arrange(mx) %>% mutate(rank=row_number()) %>% select(-mx)-> resultsDT # This is the Data table datatable(resultsDT, class = 'cell-border stripe', rownames = FALSE, options = list(pageLength = 5), colnames = c('Distance to A (Km)', 'Distance to B (Km)', 'Latitude', 'Longitude', 'Name', 'Address', 'Rank')) # Map with the locations using leaflet resultsDT %>% leaflet() %>% addTiles() %>% addCircleMarkers( lng=resultsDT$lng, lat=resultsDT$lat, radius = 8, color = "blue", stroke = FALSE, fillOpacity = 0.5, popup=paste(paste0("", resultsDT$name, ""), resultsDT$address, sep=" ") ) %>% addCircleMarkers( lng=p1$lon, lat=p1$lat, radius = 10, color = "red", stroke = FALSE, fillOpacity = 0.5, popup=paste("Place 1", place1, sep=" ") )%>% addCircleMarkers( lng=p2$lon, lat=p2$lat, radius = 10, color = "red", stroke = FALSE, fillOpacity = 0.5, popup=paste("Place 2", place2, sep=" ") )

**leave a comment**for the author, please follow the link and comment on their blog:

**R – Fronkonstin**.

R-bloggers.com offers

**daily e-mail updates**about R news and tutorials on topics such as: Data science, Big Data, R jobs, visualization (ggplot2, Boxplots, maps, animation), programming (RStudio, Sweave, LaTeX, SQL, Eclipse, git, hadoop, Web Scraping) statistics (regression, PCA, time series, trading) and more...