def meteostat(lat_degr,long_degr,alt_avrg,start_time,end_time):
#Retrieve nearest weather station
stations = Stations()
stations = stations.nearby(lat_degr, long_degr)
station = stations.fetch(1)
#Use Point to agregate nearby weather stations data for better accuracy
weather_point = Point(lat_degr, long_degr,alt_avrg)
weather_data = Hourly(weather_point, start_time - timedelta(0,7200), end_time + timedelta(0,7200))
weather_data = weather_data.fetch()
#Use weather data from nearest station if Point returns an empty dataset
if weather_data.empty:
weather_data = Hourly(station, start_time - timedelta(0,7200), end_time + timedelta(0,7200))
weather_data = weather_data.fetch()
return weather_data
def fetch(
start=datetime(2020, 1, 1),
end=datetime.now(),
lat=33.5020, # Closest to UAB
lon=-86.8064):
start += timedelta(hours=6) # UTC offset
end += timedelta(hours=6)
end.replace(microsecond=0, second=0)
stations = Stations(lat=lat, lon=lon)
station = stations.fetch(1)
data = Hourly(station, start, end)
data = data.normalize()
data = data.interpolate()
df = data.fetch()
out = {}
last_row = None
for row in df.itertuples():
if last_row:
out.update(interpolate_minutes(last_row, row))
last_row = row
current_time = last_row.time - timedelta(hours=6)
while current_time <= end - timedelta(hours=6):
out[create_datetime_ID(current_time)] = {
"time": current_time,
"temp": last_row.temp
}
current_time += timedelta(minutes=1)
return out
def get_weather(station: str, start: datetime, end: datetime) -> dict:
# Get hourly data
hourly_obs = Hourly(station, start, end)
df = hourly_obs.fetch()
return {
"rain": df["prcp"].fillna(0).sum(),
"temperature": df["temp"].mean(),
"wind_speed": df["wspd"].mean(),
}
from meteostat import Stations, Hourly from datetime import datetime import matplotlib.pyplot as plt # Hourly stations = Stations(lat = 50, lon = 8) station = stations.fetch(1) data = Hourly(station, start = datetime(2010, 1, 1), end = datetime(2020, 1, 1, 23, 59)) data = data.fetch() data.plot(x = 'time', y = ['temp'], kind = 'line') plt.show()
def get_weather(road):
#url='http://pro.openweathermap.org/data/2.5/weather?q=Birmingham,uk&APPID=b68960a402d40a497b69695725ad73af'
# Create Point for Birmingham
#Birmingham = Point(52.49, -1.86, 140)
#Bradford = Point(53.79, -1.75, 168.78)
#Leeds= Point(53.801277, -1.548567, 340)
#print('City : ',City)
month = datetime.today().month
day = datetime.today().day
year = datetime.today().year
hour = datetime.today().hour
#print(year, month, day, hour, City)
start = datetime(year, month, day, hour, 0)
end = datetime(year, month, day, hour, 59)
#print(hour)
if road == 'M606':
#City = Point(53.79, -1.75, 168.78)
data = Hourly(Point(53.79, -1.75, 168.78), start, end)
#City='Bradford'
elif road == 'M621':
#City='Leeds'
City = Point(53.801277, -1.548567, 340)
data = Hourly(Point(53.801277, -1.548567, 340), start, end)
elif road == 'A38(M)':
#City='Birmingham'
City = Point(52.49, -1.86, 140)
data = Hourly(Point(52.49, -1.86, 140), start, end)
#data = Hourly(City, start, end)
weather_df = data.fetch()
print(start)
print(weather_df)
coco = list(weather_df['coco'])
#print(coco[0])
#str1=str(round(coco[0]))
#print(str1)
#print(weather[coco])
weather = {
'1': 'Clear',
'2': 'Fair',
'3': 'Cloudy',
'4': 'Overcast',
'5': 'Fog',
'6': 'Freezing Fog',
'7': 'Light Rain',
'8': 'Rain',
'9': 'Heavy Rain',
'10': 'Freezing Rain',
'11': 'Heavy Freezing Rain',
'12': 'Sleet',
'13': 'Heavy Sleet',
'14': 'Light Snowfall',
'15': 'Snowfall',
'16': 'Heavy Snowfall',
'17': 'Rain Shower',
'18': 'Heavy Rain Shower',
'19': 'Sleet Shower',
'20': 'Heavy Sleet Shower',
'21': 'Snow Shower',
'22': 'Heavy Snow Shower',
'23': 'Lightning',
'24': 'Hail',
'25': 'Thunderstorm',
'26': 'Heavy Thunderstorm',
'27': 'Storm'
}
weather_df = weather_df.reset_index()
#print(weather[str1])
return coco
from meteostat import Stations, Hourly from datetime import datetime # Hourly stations = Stations(lat = 50, lon = 8) station = stations.fetch(1) data = Hourly(station, start = datetime(2020, 1, 1), end = datetime(2020, 1, 1, 23, 59)) print(data.fetch())
def point_hourly():
"""
Return hourly point data in JSON format
"""
# Get query parameters
args = utils.get_parameters(parameters)
# Check if required parameters are set
if args['lat'] and args['lon'] and len(args['start']) == 10 and len(
args['end']) == 10:
try:
# Convert start & end date strings to datetime
start = datetime.strptime(args['start'], '%Y-%m-%d')
end = datetime.strptime(f'{args["end"]} 23:59:59',
'%Y-%m-%d %H:%M:%S')
# Get number of days between start and end date
date_diff = (end - start).days
# Check date range
if date_diff < 0 or date_diff > max_days:
# Bad request
abort(400)
# Caching
now_diff = (datetime.now() - end).days
if now_diff < 3:
cache_time = 60 * 60
elif now_diff < 30:
cache_time = 60 * 60 * 24
else:
cache_time = 60 * 60 * 24 * 3
Hourly.max_age = cache_time
# Create a point
location = Point(args['lat'], args['lon'], args['alt'])
# Get data
data = Hourly(location,
start,
end,
timezone=args['tz'],
model=args['model'])
# Check if any data
if data.count() > 0:
# Normalize data
data = data.normalize()
# Aggregate
if args['freq']:
data = data.aggregate(args['freq'])
# Unit conversion
if args['units'] == 'imperial':
data = data.convert(units.imperial)
elif args['units'] == 'scientific':
data = data.convert(units.scientific)
# Fetch DataFrame
data = data.fetch()
# Convert to integer
data['tsun'] = data['tsun'].astype('Int64')
data['coco'] = data['coco'].astype('Int64')
# DateTime Index to String
data.index = data.index.strftime('%Y-%m-%d %H:%M:%S')
data = data.reset_index().to_json(orient="records")
else:
# No data
data = '[]'
# Inject meta data
meta = {}
meta['generated'] = datetime.now().strftime('%Y-%m-%d %H:%M:%S')
meta['stations'] = location.stations.to_list()
# Generate output string
output = f'''{{"meta":{json.dumps(meta)},"data":{data}}}'''
# Return
return utils.send_response(output, cache_time)
except BaseException:
# Bad request
abort(400)
else:
# Bad request
abort(400)