-
Notifications
You must be signed in to change notification settings - Fork 0
/
Copy pathscrape_data.py
127 lines (93 loc) · 3.46 KB
/
scrape_data.py
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
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
# -*- coding: UTF-8 -*-
#!/usr/bin/env
from __future__ import (unicode_literals, print_function)
# ---------------------------------------------------------------------------
# Name: Parser Geocoder from Excel file
# Purpose: todo
#
# Author: Julien Moura (https://twitter.com/geojulien)
#
# Python: 2.7.x
# Created: 14/05/2015
# Updated: 15/05/2015
#
# Licence: GPL 3
# ----------------------------------------------------------------------------
# ############################################################################
# ########## Libraries #############
# ##################################
# Standard library
# Python 3 backported
# 3rd party libraries
from geojson import dump, Feature, FeatureCollection, Point
import overpy
# set path
import os
print(os.getcwd())
# Custom modules
# ############################################################################
# ########## Main program ##########
# ##################################
# liste pour stocker les différents points
li_openBeerMap = []
# création de l'instace de l'API overpass
api = overpy.Overpass()
# zone dans laquelle rechercher
# Paris
# bbox = (48.658291, 2.086790, 49.046940, 2.637910)
# Melbourne
bbox = (-37.8140000 - 0.195, 144.96332 - 0.195, -37.8140000 + 0.195, 144.96332 + 0.195) # Melbourne
# stockages des résultats pour trier plus facilement ensuite
bars = api.query(str("node{0}[amenity=bar];out;".format(bbox)))
pubs = api.query(str("node{0}[amenity=pub];out;".format(bbox)))
brasseries_art = api.query(str("node{0}[craft=brewery];out;".format(bbox)))
# quelques petites stats
print(len(bars.nodes))
print(len(pubs.nodes))
print(len(brasseries_art.nodes))
x = 1
# liste des bars où la bière est précisée
for node in bars.nodes:
x += 1
# création de l'objet géographique
point = Point((round(node.lon, 4), round(node.lat, 4)))
obj = Feature(geometry=point,
id=x,
properties={"NAME": node.tags.get('name'),
"TYPE": node.tags.get('amenity'),
"OSM_ID": node.id
})
li_openBeerMap.append(obj)
# liste des pubs où la bière est précisée
for node in pubs.nodes:
x += 1
# création de l'objet géographique
point = Point((round(node.lon, 4), round(node.lat, 4)))
obj = Feature(geometry=point,
id=x,
properties={"NAME": node.tags.get('name'),
"TYPE": node.tags.get('amenity'),
"OSM_ID": node.id
})
li_openBeerMap.append(obj)
# liste des pubs où la bière est précisée
for node in brasseries_art.nodes:
x += 1
# création de l'objet géographique
point = Point((round(node.lon, 4), round(node.lat, 4)))
obj = Feature(geometry=point,
id=x,
properties={"NAME": node.tags.get('name'),
"TYPE": 'Brasserie',
"OSM_ID": node.id
})
li_openBeerMap.append(obj)
# sérialisation
featColl = FeatureCollection(li_openBeerMap)
with open("./melbourne_bar_viz.geojson", "w") as outfile:
dump(featColl, outfile, sort_keys=True)
# ############################################################################
# ##### Stand alone program ########
# ##################################
# if __name__ == '__main__':
# """ standalone execution """