def facility_vs_building_set(s):
print s
conn = uo.connect('all')
if s == 'AI':
ids = gbs.get_cat_set(['A', 'I'], conn)
elif s == 'covered':
ids = gbs.get_covered_set()
df_f = pd.read_csv(os.getcwd() + '/input/FY/covered/Covered_Facilities_All Energy mmBTUs_FY14_EISA07Sec432_input.csv')
facility_eisa = set(df_f['Facility_Number'].tolist())
facility_eisa = [x[:8] for x in facility_eisa if type(x) != float]
f_ids = [x for x in ids if '0000' in x]
# for x in sorted(f_ids):
# print x
b_ids = [x for x in ids if not '0000' in x]
print 'total {0}, facility {1}, building {2}'.format(len(ids), len(f_ids), len(b_ids))
bf_ids = ['{0}0000{1}'.format(x[:2], x[-2:]) for x in b_ids]
print len(common)
common = (set(bf_ids).intersection(f_ids))
for y in common:
print y
print [x for x in b_ids if '{0}0000{1}'.format(x[:2], x[-2:]) == y]
if s == 'covered':
print 'eisa facility', len(facility_eisa)
print 'common ids from eisa', len(set(f_ids).intersection(facility_eisa))
print 'different ids from eisa', (set(f_ids).difference(facility_eisa))
common = (set(f_ids).difference(facility_eisa))
for y in common:
print y
print [x for x in b_ids if '{0}0000{1}'.format(x[:2], x[-2:]) == y]
def facility_vs_building_set(s):
print s
conn = uo.connect('all')
if s == 'AI':
ids = gbs.get_cat_set(['A', 'I'], conn)
elif s == 'covered':
ids = gbs.get_covered_set()
df_f = pd.read_csv(
os.getcwd() +
'/input/FY/covered/Covered_Facilities_All Energy mmBTUs_FY14_EISA07Sec432_input.csv'
)
facility_eisa = set(df_f['Facility_Number'].tolist())
facility_eisa = [x[:8] for x in facility_eisa if type(x) != float]
f_ids = [x for x in ids if '0000' in x]
# for x in sorted(f_ids):
# print x
b_ids = [x for x in ids if not '0000' in x]
print 'total {0}, facility {1}, building {2}'.format(
len(ids), len(f_ids), len(b_ids))
bf_ids = ['{0}0000{1}'.format(x[:2], x[-2:]) for x in b_ids]
print len(common)
common = (set(bf_ids).intersection(f_ids))
for y in common:
print y
print[x for x in b_ids if '{0}0000{1}'.format(x[:2], x[-2:]) == y]
if s == 'covered':
print 'eisa facility', len(facility_eisa)
print 'common ids from eisa', len(
set(f_ids).intersection(facility_eisa))
print 'different ids from eisa', (set(f_ids).difference(facility_eisa))
common = (set(f_ids).difference(facility_eisa))
for y in common:
print y
print[x for x in b_ids if '{0}0000{1}'.format(x[:2], x[-2:]) == y]
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 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 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 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)
