Compare commits
No commits in common. "5aafac2025802b898f537e61de75ad81852ea5e5" and "cda2195419d7e28b5dc6336a11b05239040cfb5a" have entirely different histories.
5aafac2025
...
cda2195419
2
.gitignore
vendored
2
.gitignore
vendored
@ -1,2 +0,0 @@
|
|||||||
geoapi.ini
|
|
||||||
.idea/
|
|
||||||
@ -1,9 +0,0 @@
|
|||||||
<?xml version="1.0" encoding="UTF-8"?>
|
|
||||||
<module type="PYTHON_MODULE" version="4">
|
|
||||||
<component name="NewModuleRootManager" inherit-compiler-output="true">
|
|
||||||
<exclude-output />
|
|
||||||
<content url="file://$MODULE_DIR$" />
|
|
||||||
<orderEntry type="inheritedJdk" />
|
|
||||||
<orderEntry type="sourceFolder" forTests="false" />
|
|
||||||
</component>
|
|
||||||
</module>
|
|
||||||
@ -1,12 +0,0 @@
|
|||||||
[google]
|
|
||||||
apikey=<<API KEY HERE>>
|
|
||||||
|
|
||||||
[mapquest]
|
|
||||||
apikey=<<CONSUMER KEY HERE>>
|
|
||||||
secret=<<CONSUMER SECRET HERE>>
|
|
||||||
|
|
||||||
[geocodio]
|
|
||||||
apikey=<<API KEY HERE>>
|
|
||||||
|
|
||||||
[locationiq]
|
|
||||||
token=<<TOKEN HERE>>
|
|
||||||
77
src/dcalc.py
77
src/dcalc.py
@ -1,77 +0,0 @@
|
|||||||
# DCalc.py: basic test of distance calculation from the Boulder office
|
|
||||||
|
|
||||||
import sys
|
|
||||||
import argparse
|
|
||||||
import json
|
|
||||||
import math
|
|
||||||
import configparser
|
|
||||||
import urllib.parse
|
|
||||||
import urllib.request
|
|
||||||
|
|
||||||
|
|
||||||
class GeocodingError(RuntimeError):
|
|
||||||
pass
|
|
||||||
|
|
||||||
|
|
||||||
def geocode_google(config, address):
|
|
||||||
apikey = config['google']['apikey']
|
|
||||||
if not apikey:
|
|
||||||
raise GeocodingError("Google API key not specified")
|
|
||||||
query = { 'key': apikey, 'address': address, 'region': 'us'}
|
|
||||||
url = 'https://maps.googleapis.com/maps/api/geocode/json?' + \
|
|
||||||
urllib.parse.urlencode(query, quote_via=urllib.parse.quote)
|
|
||||||
with urllib.request.urlopen(url) as response:
|
|
||||||
if response.status == 200:
|
|
||||||
apireturn = json.loads(response.read())
|
|
||||||
stat = apireturn['status']
|
|
||||||
if stat == 'OK':
|
|
||||||
results = apireturn['results']
|
|
||||||
if len(results) > 1:
|
|
||||||
raise GeocodingError(f"Google API returned ambiguous results (total count {len(results)})")
|
|
||||||
coords = results[0]['geometry']['location']
|
|
||||||
return coords['lat'], coords['lng']
|
|
||||||
elif stat == 'ZERO_RESULTS':
|
|
||||||
return None
|
|
||||||
else:
|
|
||||||
raise GeocodingError(f"Google API returns status of {stat}")
|
|
||||||
else:
|
|
||||||
raise GeocodingError(f"Google API returns {response.status} HTTP status code")
|
|
||||||
|
|
||||||
|
|
||||||
OFFICE_LOCATION = (40.0187905, -105.2764775) # office location - 1433 Pearl Street, Boulder, CO
|
|
||||||
RADIUS = 6371.0 * 1000.0 # Earth radius in meters
|
|
||||||
METERS_PER_MILE = 1852.0 # number of meters per mile
|
|
||||||
|
|
||||||
|
|
||||||
def distance_miles(point1, point2):
|
|
||||||
Φ1 = math.radians(point1[0])
|
|
||||||
Φ2 = math.radians(point2[0])
|
|
||||||
ΔΦ = math.radians(point2[0] - point1[0])
|
|
||||||
Δλ = math.radians(point2[1] - point1[1])
|
|
||||||
a = math.sin(ΔΦ / 2.0) * math.sin(ΔΦ / 2.0) + math.cos(Φ1) * math.cos(Φ2) * math.sin(Δλ / 2.0) * math.sin(Δλ / 2.0)
|
|
||||||
c = 2 * math.atan2(math.sqrt(a), math.sqrt(1.0 - a))
|
|
||||||
return RADIUS * c / METERS_PER_MILE
|
|
||||||
|
|
||||||
|
|
||||||
cmdline_parser = argparse.ArgumentParser()
|
|
||||||
cmdline_parser.add_argument('address', nargs='+', help='The address to be calculated')
|
|
||||||
cmdline_parser.add_argument('-C', '--config', default='geoapi.ini', help='The geocoding API configuration file')
|
|
||||||
|
|
||||||
|
|
||||||
def main(args):
|
|
||||||
opts = cmdline_parser.parse_args(args)
|
|
||||||
config = configparser.ConfigParser()
|
|
||||||
config.read(opts.config)
|
|
||||||
my_address = ' '.join(opts.address)
|
|
||||||
print(f"Address: '{my_address}'")
|
|
||||||
coords = geocode_google(config, my_address)
|
|
||||||
if coords:
|
|
||||||
print(f"Coordinates: Latitude {coords[0]}, longitude {coords[1]}")
|
|
||||||
dist = distance_miles(OFFICE_LOCATION, coords)
|
|
||||||
print(f"Distance from Boulder office: {dist} miles")
|
|
||||||
else:
|
|
||||||
print("Location was not found.")
|
|
||||||
|
|
||||||
|
|
||||||
if __name__ == '__main__':
|
|
||||||
sys.exit(main(sys.argv[1:]))
|
|
||||||
124
src/gcoder.py
124
src/gcoder.py
@ -1,124 +0,0 @@
|
|||||||
# GCoder.py: basic geocoding test
|
|
||||||
|
|
||||||
import sys
|
|
||||||
import argparse
|
|
||||||
import json
|
|
||||||
import configparser
|
|
||||||
import urllib.parse
|
|
||||||
import urllib.request
|
|
||||||
|
|
||||||
|
|
||||||
class GeocodingError(RuntimeError):
|
|
||||||
pass
|
|
||||||
|
|
||||||
|
|
||||||
def geocode_google(config, address):
|
|
||||||
apikey = config['google']['apikey']
|
|
||||||
if not apikey:
|
|
||||||
raise GeocodingError("Google API key not specified")
|
|
||||||
query = { 'key': apikey, 'address': address, 'region': 'us'}
|
|
||||||
url = 'https://maps.googleapis.com/maps/api/geocode/json?' + \
|
|
||||||
urllib.parse.urlencode(query, quote_via=urllib.parse.quote)
|
|
||||||
with urllib.request.urlopen(url) as response:
|
|
||||||
if response.status == 200:
|
|
||||||
apireturn = json.loads(response.read())
|
|
||||||
stat = apireturn['status']
|
|
||||||
if stat == 'OK':
|
|
||||||
results = apireturn['results']
|
|
||||||
if len(results) > 1:
|
|
||||||
raise GeocodingError(f"Google API returned ambiguous results (total count {len(results)})")
|
|
||||||
coords = results[0]['geometry']['location']
|
|
||||||
return coords['lat'], coords['lng']
|
|
||||||
elif stat == 'ZERO_RESULTS':
|
|
||||||
return None
|
|
||||||
else:
|
|
||||||
raise GeocodingError(f"Google API returns status of {stat}")
|
|
||||||
else:
|
|
||||||
raise GeocodingError(f"Google API returns {response.status} HTTP status code")
|
|
||||||
|
|
||||||
|
|
||||||
def geocode_mapquest(config, address):
|
|
||||||
apikey = config['mapquest']['apikey']
|
|
||||||
if not apikey:
|
|
||||||
raise GeocodingError("MapQuest API key not specified")
|
|
||||||
query = { 'key': apikey, 'inFormat': 'kvp', 'outFormat': 'json', 'location': address,
|
|
||||||
'thumbMaps': 'false', 'maxResults': 1}
|
|
||||||
url = 'https://www.mapquestapi.com/geocoding/v1/address?' + urllib.parse.urlencode(query)
|
|
||||||
with urllib.request.urlopen(url) as response:
|
|
||||||
if response.status == 200:
|
|
||||||
apireturn = json.loads(response.read())
|
|
||||||
stat = apireturn['info']['statuscode']
|
|
||||||
if stat == 0:
|
|
||||||
results = apireturn['results']
|
|
||||||
if len(results) > 1:
|
|
||||||
raise GeocodingError(f"MapQuest API returned ambiguous results (total count {len(results)})")
|
|
||||||
locations = results[0]['locations']
|
|
||||||
if len(locations) > 1:
|
|
||||||
raise GeocodingError(f"MapQuest API returned ambiguous locations (total count {len(locations)})")
|
|
||||||
coords = locations[0]['latLng']
|
|
||||||
return coords['lat'], coords['lng']
|
|
||||||
else:
|
|
||||||
raise GeocodingError(f"MapQuest API returns status of {stat}")
|
|
||||||
else:
|
|
||||||
raise GeocodingError(f"MapQuest API returns {response.status} HTTP status code")
|
|
||||||
|
|
||||||
|
|
||||||
def geocode_geocodio(config, address):
|
|
||||||
apikey = config['geocodio']['apikey']
|
|
||||||
if not apikey:
|
|
||||||
raise GeocodingError("Geocodio API key not specified")
|
|
||||||
query = {'q': address, 'api_key': apikey}
|
|
||||||
url = 'https://api.geocod.io/v1.6/geocode?' + urllib.parse.urlencode(query)
|
|
||||||
with urllib.request.urlopen(url) as response:
|
|
||||||
if response.status == 200:
|
|
||||||
apireturn = json.loads(response.read())
|
|
||||||
coords = apireturn['results'][0]['location']
|
|
||||||
return coords['lat'], coords['lng']
|
|
||||||
else:
|
|
||||||
raise GeocodingError(f"Geocodio API returns {response.status} HTTP status code")
|
|
||||||
|
|
||||||
|
|
||||||
def geocode_locationiq(config, address):
|
|
||||||
apikey = config['locationiq']['token']
|
|
||||||
if not apikey:
|
|
||||||
raise GeocodingError("LocationIQ API key not specified")
|
|
||||||
query = {'q': address, 'key': apikey, 'format': 'json', 'limit': 1}
|
|
||||||
url = 'https://us1.locationiq.com/v1/search.php?' + urllib.parse.urlencode(query, quote_via=urllib.parse.quote)
|
|
||||||
with urllib.request.urlopen(url) as response:
|
|
||||||
if response.status == 200:
|
|
||||||
apireturn = json.loads(response.read())
|
|
||||||
place = apireturn[0]
|
|
||||||
return place['lat'], place['lon']
|
|
||||||
else:
|
|
||||||
raise GeocodingError(f"LocationIQ API returns {response.status} HTTP status code")
|
|
||||||
|
|
||||||
|
|
||||||
geocoding_procs = {
|
|
||||||
'google': geocode_google,
|
|
||||||
'mapquest': geocode_mapquest,
|
|
||||||
'geocodio': geocode_geocodio,
|
|
||||||
'locationiq': geocode_locationiq
|
|
||||||
}
|
|
||||||
|
|
||||||
cmdline_parser = argparse.ArgumentParser()
|
|
||||||
cmdline_parser.add_argument('address', nargs='+', help='The address to be geocoded')
|
|
||||||
cmdline_parser.add_argument('-C', '--config', default='geoapi.ini', help='The geocoding API configuration file')
|
|
||||||
cmdline_parser.add_argument('-g', '--geocoder', default='google', choices=geocoding_procs.keys(),
|
|
||||||
help='Geocoding processor to use to get coordinates (default: google)')
|
|
||||||
|
|
||||||
|
|
||||||
def main(args):
|
|
||||||
opts = cmdline_parser.parse_args(args)
|
|
||||||
config = configparser.ConfigParser()
|
|
||||||
config.read(opts.config)
|
|
||||||
my_address = ' '.join(opts.address)
|
|
||||||
print(f"Address: '{my_address}'")
|
|
||||||
coords = geocoding_procs[opts.geocoder](config, my_address)
|
|
||||||
if coords:
|
|
||||||
print(f"Coordinates: Latitude {coords[0]}, longitude {coords[1]}")
|
|
||||||
else:
|
|
||||||
print("Location was not found.")
|
|
||||||
|
|
||||||
|
|
||||||
if __name__ == '__main__':
|
|
||||||
sys.exit(main(sys.argv[1:]))
|
|
||||||
Loading…
x
Reference in New Issue
Block a user