topic Re: Is there a way to Reverse Geocode? Take Lat/Long and convert to an address? in Graphics Programming

Is there a way to Reverse Geocode? Take Lat/Long and convert to an address?

DSouza — Tue, 07 Jul 2015 18:49:31 GMT

Is there a way to take a set of latitude and longitude coordinates and reverse geocode them to get an address?

Re: Is there a way to Reverse Geocode? Take Lat/Long and convert to an address?

Reeza — Tue, 07 Jul 2015 21:05:16 GMT

I don't know if SAS can do it, but there are tools online that can do it:

geocoder.ca: geocoding for North America - USA and Canada

You could probably write an algorithm to pass data to the service via the API - $1 per 400 lookups appears to be the rate.

Probably need to investigate the TOS.

Re: Is there a way to Reverse Geocode? Take Lat/Long and convert to an address?

Darrell_sas — Wed, 08 Jul 2015 14:20:22 GMT

SAS cannot do Reverse Geocoding.

There are some commercial places but there are a few free ones. I can't say how good they are. You might Google "Reverse Geocoding". A few more that show up are:

Free Online Reverse Geocoding

GeoNames reverse address geocoding

Re: Is there a way to Reverse Geocode? Take Lat/Long and convert to an address?

FriedEgg — Wed, 08 Jul 2015 14:38:51 GMT

Here is an example using the google geocoding service. It can be adapted for reverse geocoding.

An interface for using Google Geocode API with PROC GROOVY in SAS

Another implementation of using Google geocoding API in SAS, this time using PROC FCMP

Re: Is there a way to Reverse Geocode? Take Lat/Long and convert to an address?

Darrell_sas — Wed, 08 Jul 2015 15:05:15 GMT

Google prohibits using a service (such as geocoding) without displaying the contents/results on a Google Map. See 10.1.1.h of the Google Maps Terms of Use:

https://developers.google.com/maps/terms

Re: Is there a way to Reverse Geocode? Take Lat/Long and convert to an address?

FriedEgg — Wed, 08 Jul 2015 15:19:38 GMT

A very important note

A direct link:

Google Maps/Google Earth APIs Terms of Service

Re: Is there a way to Reverse Geocode? Take Lat/Long and convert to an address?

bpfields — Mon, 19 Sep 2016 20:08:13 GMT

After much research I landed on a service by TA&M University. You can process batches. I assume this is why you are looking for a SAS solution:

http://geoservices.tamu.edu/Services/Geocode/

Please let us know if you found a SAS solution. Thanks.

Re: Is there a way to Reverse Geocode? Take Lat/Long and convert to an address?

GraphGuy — Tue, 25 Oct 2016 13:19:57 GMT

Do you need to find the closest street address? ... or would the closest city or zipcode be good enough?

(if the closest city or zipcode is good enough, then I could write some SAS code to do that)

Re: Is there a way to Reverse Geocode? Take Lat/Long and convert to an address?

bpfields — Tue, 25 Oct 2016 14:04:57 GMT

Zip Code is fine.

Re: Is there a way to Reverse Geocode? Take Lat/Long and convert to an address?

GraphGuy — Tue, 25 Oct 2016 14:30:42 GMT

Ok - then here's a way to do that, using brute force (and the sashelp.zipcode dataset).

You might refine this to make it a little less code, and more efficient - but I like to take things 1 step at a time, to make it simpler to understand. 🙂

data foo;
input customer_number lat long;
datalines;
1 35.8227817 -78.755733
2 38.8976398 -77.036627
;
run;

/* Create all possible pairs of my data, and zipcode lat/longs */
proc sql;
create table pairs as
select unique foo.*, zipcode.zip, zipcode.city, zipcode.x, zipcode.y
from foo, sashelp.zipcode
where zipcode.x^=. and zipcode.y^=.;
quit; run;

/* calculate distance between each customer and each zipcode centroid */
data pairs; set pairs;
Distance = geodist(lat, long, y, x, 'DM');
run;

/* find the closest zipcode to each customer */
proc sort data=pairs out=pairs;
by customer_number distance;
run;
data pairs; set pairs (keep = customer_number lat long zip city);
by customer_number;
label zip='Closest Zipcode';
if first.customer_number then output;
run;