-
Notifications
You must be signed in to change notification settings - Fork 0
/
Copy pathget_data.py
61 lines (49 loc) · 1.57 KB
/
get_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
from datetime import datetime
import requests
import matplotlib.pyplot as plt
getObservationPayload = {
"request": "GetObservation",
"service": "SOS",
"version": "2.0.0",
"procedure": "Penlee",
"offering": "PenleeMetData",
"observedProperty": "http://mmisw.org/ont/cf/parameter/air_temperature",
"temporalFilter": {
"during": {
"ref": "om:phenomenonTime",
"value": [
"2015-01-02T13:50:00.000Z",
"2016-01-02T13:50:00.000Z"
]
}
}
}
getDummyObservationPayload = {
"request": "GetObservation",
"service": "SOS",
"version": "2.0.0",
"procedure": "http://www.52north.org/test/procedure/9",
"offering": "http://www.52north.org/test/offering/9",
"observedProperty": "http://www.52north.org/test/observableProperty/9_3",
"featureOfInterest": "http://www.52north.org/test/featureOfInterest/9"
}
def getObservation(server):
resp = requests.post(server, json=getObservationPayload)
#plotData(resp.json())
print resp.json()
def plotData(data):
times = []
values = []
locations = []
for item in data['observations']:
print "adding result"
print times
times.append(datetime.strptime(item['phenomenonTime'], "%Y-%m-%dT%H:%M:%S.%fZ"))
values.append(item['result']['value'])
locations.append(item['featureOfInterest']['geometry']['coordinates'])
#print times
#print values
#plt.plot(times,values)
#plt.show()
if __name__ == "__main__":
getObservation('https://rsg.pml.ac.uk/sensorweb/service')