example.py

from datetime import date

from weather.weather import WeatherExtractor, WeatherApi


def print_data(weather_data):
    for row in weather_data.iterrows():
        # row is tuple (index, columns)
        measure = row[1]

        print "Measurement of %s at from %s for %s" % (measure['shortName'], measure['validDateTime'], measure['validityDateTime'])
        for lat, lon, val in zip(measure['lats'], measure['lons'], measure['values']):
            print "%f N %f S = %f" % (lat, lon, val)


wa = WeatherApi()

# download forecast data for november 2017
wa.get(from_date=date(2017, 11, 1), to_date=date(2017, 11, 30),
       target='nov2017-forecast.grib', request_type='forecast')

# download actual weather data november 2017
wa.get(from_date=date(2017, 11, 1), to_date=date(2017, 11, 30),
       target='nov2017-actual.grib', request_type='actual')

# query the downloaded data
we = WeatherExtractor()
# load actual and forecasted weather data
we.load(['nov2017-forecast.grib', 'nov2017-actual.grib'])

""" Get forecasted data from 1-11-2017 for 2-11-2017, 3-11-2017 and 4-11-2017 for all grid points. """
weather_data = we.get_forecast(base_date=date(2017, 11, 1), from_date=date(
    2017, 11, 2), to_date=date(2017, 11, 4), aggtime='hour', aggloc='grid')

# print the result
print_data(weather_data)

""" Get forecasted data from 1-11-2017 for 2-11-2017, 3-11-2017 and 4-11-2017 for 
two specific points with latitudes and longitudes: (45.01, 13.00) and (46.00, 12.05) """
points = [{'lat': 45.01, 'lon': 13.0}, {'lat': 46.0, 'lon': 12.05}]
weather_data = we.get_forecast(base_date=date(2017, 11, 1), from_date=date(
    2017, 11, 2), to_date=date(2017, 11, 4), aggtime='hour', aggloc='points', interp_points=points)
# print the result
print_data(weather_data)

""" Get actual weather data for 2-11-2017, 3-11-2017 and 4-11-2017 for 
two specific points with latitudes and longitudes: (45.01, 13.00) and (46.00, 12.05) """
points = [{'lat': 45.01, 'lon': 13.0}, {'lat': 46.0, 'lon': 12.05}]
weather_data = we.get_actual(from_date=date(
    2017, 11, 2), to_date=date(2017, 11, 4), aggtime='hour', aggloc='points', interp_points=points)
# print the result
print_data(weather_data)