1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
|
#!/usr/bin/env python3
import argparse
import urllib.parse
import xdg.BaseDirectory
from OSMPythonTools.cachingStrategy import CachingStrategy, JSON
from OSMPythonTools.api import Api
from shapely.geometry import shape
CachingStrategy.use(JSON, cacheDir=xdg.BaseDirectory.save_cache_path("osm_where"))
import logging
logging.getLogger("OSMPythonTools").setLevel(logging.ERROR)
def parse_name(name, area, lang):
overpass = Overpass()
result = overpass.query(
f'area["ISO3166-1"="{area}"][admin_level=2];relation["name:{lang}"="{name}"](area);out center;'
)
if len(result.elements()) > 0:
bakhmut = result.elements()[0]
return bakhmut.centerLat(), bakhmut.centerLon()
else:
return None
def parse_URL(url):
api = Api()
parsed_url = urllib.parse.urlparse(url)
way = api.query(parsed_url.path)
geom = shape(way.geometry())
return geom.centroid.y, geom.centroid.x
def main():
parser = argparse.ArgumentParser(
prog="osm_where", description="Get geo: URI for given locality from OSM"
)
parser.add_argument("url", nargs="?", help="URL of the location")
parser.add_argument("-n", "--name", help="name of the location")
parser.add_argument("-a", "--area", default="UA", help="ISO 3166 area code")
parser.add_argument("-l", "--lang", default="en", help="ISO 639-1 language code")
args = parser.parse_args()
if args.name:
geo = parse_name(args.name, args.area, args.lang)
else:
geo = parse_URL(args.url)
if geo:
print(f"geo:{geo[0]:.7f},{geo[1]:.7f}")
if __name__ == "__main__":
main()
|
