def show_covered_exception(b, **kwargs): euas = gbs.get_all_building_set() print b in euas bs = [x for x in euas if '{0}0000{1}'.format(x[:2], x[-2:]) == b] if 'bs' in kwargs: bs = kwargs['bs'] dfs = [] df = uo.view_building(b, 'Electric_(kBtu)') dfs.append(df) if 'year' in kwargs: year = kwargs['year'] df = df[(df['year'] == year) & (df['month'].isin([1, 2, 3]))] print b print df[['Region_No.', 'month', 'Electric_(kBtu)']] for x in bs[:3]: print x df = uo.view_building(x, 'Electric_(kBtu)') if 'year' in kwargs: year = kwargs['year'] df = df[(df['year'] == year) & (df['month'].isin([1, 2, 3]))] dfs.append(df) if 'year' in kwargs: print print df[['Region_No.', 'month', 'Electric_(kBtu)']] df_all = pd.concat(dfs, ignore_index=True) df_all.drop(['year', 'month'], axis=1, inplace=True) df_all.sort(['Fiscal_Year', 'Fiscal_Month', 'Building_Number'], inplace=True) print 'write to {0}.csv'.format(b) df_all.to_csv(homedir + 'question/facility_building/{0}.csv'.format(b), index=False) return
def check_covered_component(): ids = gbs.get_covered_set() f_ids = [x for x in ids if '0000' in x] b_ids = [x for x in ids if not '0000' in x] euas = gbs.get_all_building_set() print len(euas) print 'total {0}, facility {1}, building {2}'.format( len(ids), len(f_ids), len(b_ids)) for f in f_ids: dfs = [] df = uo.view_building(f, 'Electric_(kBtu)') bs = [ x for x in euas if '{0}0000{1}'.format(x[:2], x[-2:]) == f and not '0000' in x ] if len(bs) == 0: print 'no building under {0}'.format(f) continue dfs.append(df) for b in bs: df = uo.view_building(b, 'Electric_(kBtu)') dfs.append(df) df_all = pd.concat(dfs, ignore_index=True) df_all.sort(['year', 'month', 'Building_Number'], inplace=True) print 'write to {0}.csv'.format(f) df_all.to_csv(homedir + 'question/facility_building/{0}.csv'.format(f), index=False)
def eisa_building_in_euas(): df = pd.read_csv(os.getcwd() + '/input/FY/covered/Covered_Facilities_All Energy mmBTUs_FY14_EISA07Sec432_input.csv') buildings = df['Building_Number'].unique() facility = df['Facility_Number'].unique() print 'total building {0}, total facility {1}'.format(len(buildings), len(facility)) euas = gbs.get_all_building_set() print len(euas.intersection(buildings)), len(euas.intersection(facility)) print euas.intersection(buildings)
def euas_covered(): covered = gbs.get_covered_set() euas = gbs.get_all_building_set() good_elec = gbs.get_energy_set('eui_elec') good_gas = gbs.get_energy_set('eui_gas') print len(covered.intersection(euas)) print len(covered.intersection(good_elec)) print len(covered.intersection(good_gas))
def eisa_building_in_euas(): df = pd.read_csv( os.getcwd() + '/input/FY/covered/Covered_Facilities_All Energy mmBTUs_FY14_EISA07Sec432_input.csv' ) buildings = df['Building_Number'].unique() facility = df['Facility_Number'].unique() print 'total building {0}, total facility {1}'.format( len(buildings), len(facility)) euas = gbs.get_all_building_set() print len(euas.intersection(buildings)), len(euas.intersection(facility)) print euas.intersection(buildings)
def check_covered_component(): ids = gbs.get_covered_set() f_ids = [x for x in ids if '0000' in x] b_ids = [x for x in ids if not '0000' in x] euas = gbs.get_all_building_set() print len(euas) print 'total {0}, facility {1}, building {2}'.format(len(ids), len(f_ids), len(b_ids)) for f in f_ids: dfs = [] df = uo.view_building(f, 'Electric_(kBtu)') bs = [x for x in euas if '{0}0000{1}'.format(x[:2], x[-2:]) == f and not '0000' in x] if len(bs) == 0: print 'no building under {0}'.format(f) continue dfs.append(df) for b in bs: df = uo.view_building(b, 'Electric_(kBtu)') dfs.append(df) df_all = pd.concat(dfs, ignore_index=True) df_all.sort(['year', 'month', 'Building_Number'], inplace=True) print 'write to {0}.csv'.format(f) df_all.to_csv(homedir + 'question/facility_building/{0}.csv'.format(f), index=False)