Esempio n. 1
0
if FIGURE_PRESET.endswith('for poster'):
    FIGURE_PATH = os.path.expanduser('~/Dropbox/MyWork/imperial/PhD/writing'
                                     '/posters/UKERCposter2013/')
else:
    FIGURE_PATH = os.path.expanduser('~/Dropbox/MyWork/imperial/PhD/writing'
                                     '/papers/tetc2013/figures/')

FIGURE_SUFFIX = '.pdf'
DATA_DIR = '/data/mine/vadeec/merged/house1'

if FIGURE_PRESET in ['daily usage histogram', 'daily usage histogram for poster']:
    START_DATE = None # datetime.datetime(year=2013, month=3, day=1)
    END_DATE = None # datetime.datetime(year=2013, month=3, day=1)
    BIN_SIZE = 'T' # D (daily) or H (hourly) or T (minutely)
    TIMESPAN = 'D' # D (daily) or W (weekly)
    spfl.setup(fig_height=8)
    GRID = True
    XTICKS_ON = False
    LATEX_PDF_OUTPUT_FILENAME = os.path.join(FIGURE_PATH,
                                             'daily_usage_histograms'+FIGURE_SUFFIX)
    if FIGURE_PRESET == 'daily usage histogram':
        CHAN_IDS = [2,3,7,17,9,19,25,8,10,11,13,42,14,45,16]
        TITLE_Y = 0.87
    else:
        CHAN_IDS = [2,3,7,25,10,42,14]
        TITLE_Y = 0.8
elif FIGURE_PRESET == 'weekly usage histogram':
    START_DATE = None # datetime.datetime(year=2013, month=3, day=1)
    END_DATE = None # datetime.datetime(year=2013, month=3, day=1)
    BIN_SIZE = 'D' # D (daily) or H (hourly) or T (minutely)
    TIMESPAN = 'W' # D (daily) or W (weekly)
Esempio n. 2
0
import setupPlottingForLaTeX as spfl
import os

ON_DURATION_THRESHOLD = 0.1 # hours

FIGURE_PATH = os.path.expanduser('~/Dropbox/MyWork/imperial/PhD/writing'
                                 '/papers/tetc2013/figures/')

DATA_DIR = DD

# Load metoffice data
print("Opening metoffice data...")
weather = metoffice.open_daily_xls('/data/metoffice/Heathrow_DailyData.xls')

# Setup figure
spfl.setup(columns=2)
fig = plt.figure()

def correlate(chan_id, weather_variable, subplot_index, annotate_y):
    # 25 = lighting circuit # (R^2 = 0.443)
    # 8 = kitchen lights (R^2 = 0.194)
    # 2 = boiler (versus radiation R^2 = 0.052, 
    #             versus mean_temp R^2 = 0.298,
    #             versus max_temp  R^2 = 0.432,
    #             versus min_temp  R^2 = 0.212)
    # 3 = solar (R^2 = 0.798)
    # 12 = fridge vs min_temp R^2 = 0.255 (with on_power_threshold = 20)
    
    print("Opening channel data...")
    channel = Channel(DATA_DIR, chan_id)
Esempio n. 3
0
                                     '/posters/UKERCposter2013/')
else:
    FIGURE_PATH = os.path.expanduser('~/Dropbox/MyWork/imperial/PhD/writing'
                                     '/papers/tetc2013/figures/')

FIGURE_SUFFIX = '.pdf'
DATA_DIR = '/data/mine/vadeec/merged/house1'

if FIGURE_PRESET in [
        'daily usage histogram', 'daily usage histogram for poster'
]:
    START_DATE = None  # datetime.datetime(year=2013, month=3, day=1)
    END_DATE = None  # datetime.datetime(year=2013, month=3, day=1)
    BIN_SIZE = 'T'  # D (daily) or H (hourly) or T (minutely)
    TIMESPAN = 'D'  # D (daily) or W (weekly)
    spfl.setup(fig_height=8)
    GRID = True
    XTICKS_ON = False
    LATEX_PDF_OUTPUT_FILENAME = os.path.join(
        FIGURE_PATH, 'daily_usage_histograms' + FIGURE_SUFFIX)
    if FIGURE_PRESET == 'daily usage histogram':
        CHAN_IDS = [2, 3, 7, 17, 9, 19, 25, 8, 10, 11, 13, 42, 14, 45, 16]
        TITLE_Y = 0.87
    else:
        CHAN_IDS = [2, 3, 7, 25, 10, 42, 14]
        TITLE_Y = 0.8
elif FIGURE_PRESET == 'weekly usage histogram':
    START_DATE = None  # datetime.datetime(year=2013, month=3, day=1)
    END_DATE = None  # datetime.datetime(year=2013, month=3, day=1)
    BIN_SIZE = 'D'  # D (daily) or H (hourly) or T (minutely)
    TIMESPAN = 'W'  # D (daily) or W (weekly)