overpass_python_crawler/osm_parser.py at master · ftheeten/overpass_python_crawler · GitHub

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
#code heavily based on the QuickOSM pluging of Quantum GIS
import requests

save_folder="C:\\WORK_CORONA\\osm_congo\\"

root_url="http://overpass-api.de/api/interpreter?data="
size_degree=4

x_min=12
x_max=32
y_min=-14
y_max=6


cpt=1
cpt_min=1
for x in range(x_min,x_max, size_degree):
    print(x)
    for y in range(y_min,y_max, size_degree):
        if cpt>=cpt_min:
            cur_x_min=str(x)
            cur_x_max=x+size_degree
            if cur_x_max > x_max:
                cur_x_max=x_max
            cur_x_max=str(cur_x_max)

            cur_y_min=str(y)
            cur_y_max=y+size_degree
            if cur_y_max > y_max:
                cur_y_max=y_max
            cur_y_max=str(cur_y_max)
            cur_bbox=[cur_x_min, cur_y_min, cur_x_max, cur_y_max]
            print(cur_bbox)
            xml_query='<osm-script output="xml" timeout="25"><union><query type="node"><has-kv k="waterway"/><bbox-query s="'+cur_y_min+'" w="'+cur_x_min+'" n="'+cur_y_max+'" e="'+cur_x_max+'"/></query><query type="way"><has-kv k="waterway"/><bbox-query s="'+cur_y_min+'" w="'+cur_x_min+'" n="'+cur_y_max+'" e="'+cur_x_max+'"/></query><query type="relation"><has-kv k="waterway"/><bbox-query s="'+cur_y_min+'" w="'+cur_x_min+'" n="'+cur_y_max+'" e="'+cur_x_max+'"/></query></union><union><item/><recurse type="down"/></union><print mode="body"/></osm-script>'
            #print(xml_query)
            full_query=root_url+xml_query
            print(full_query)
            resp = requests.get(full_query)
            #print(resp.text)
            file=open(save_folder+'output_osm_'+str(cpt)+".osm", 'w', encoding="utf-8")
            file.write(resp.text)
            file.close()
        cpt=cpt+1