Homec4science
Diffusion solar_potential (aw_solar) Browse

solarPOT/Potential_hourly/scripts/merge_annual_data.py
8cbd68ca1491aw_solar

Last Change
Raw File

merge_annual_data.py
View Options

import numpy as np
import pandas as pd
import xarray as xr
import sys
import os
import time
# # Define file paths
COMM_FP = '/work/hyenergy/output/solar_potential/Sonnendach_Gemeinden_Dachflaechen.csv' # '/Users/alinawalch/Documents/EPFL/data/EMPA_data/Sonnendach_Gemeinden_Dachflaechen.csv'
ROOF_FP = '/work/hyenergy/output/solar_potential/geographic_potential/CH_ROOFS_all_replaced.csv' # '/Volumes/Seagate Backup Plus Drive/DOKUMENTE/EPFL/data/Results/PV_w_uncertainty/CH_ROOFS_all_replaced.csv'
AREA_FP = '/work/hyenergy/output/solar_potential/geographic_potential/available_area/PRED_available_area_unc_filled.csv' # '/Volumes/Seagate Backup Plus Drive/DOKUMENTE/EPFL/data/Results/PV_w_uncertainty/PRED_available_area_unc_filled.csv'
PV_PATH = sys.argv[1]
OUTPATH = sys.argv[2]
## CONSTANTS
GT_IN = 'tilted_irrad_1*.nc'
PV_IN = 'pv_potential_1*.nc'
DF_OUT = '2010_annual_roofs.csv'
HR_OUT = '2010_hourly_comms.nc'
PV_OUT = '2010_hourly_pv_kW.csv'
GT_OUT = '2010_hourly_gt_kW_m2.csv'
FLAT_TILT = 10
MAX_ASPECT = 90
# Print file paths
print('\nCommune mapping file path: %s' %COMM_FP)
print('Available area file path: %s' %AREA_FP)
print('Roof surfaces file path: %s' %ROOF_FP)
print('PV and irradiance input path: %s' %PV_PATH)
print('Output path: %s\n\n' %OUTPATH)
# # Load commune map
tt = time.time()
comm_map = pd.read_csv(COMM_FP, encoding = 'cp1250', usecols = ['DF_UID', 'BFS_NUMMER', 'NAME'] )
# # Load roof data
roofs = pd.read_csv(ROOF_FP, usecols = ['DF_UID', 'NEIGUNG', 'AUSRICHTUNG'])
av_area = pd.read_csv(AREA_FP, usecols = ['DF_UID', 'available_area', 'available_area_std'])
roof_info = roofs.merge(av_area, on = 'DF_UID').merge(comm_map, on = 'DF_UID', how = 'left')
print('Loaded roof info in %.2fs' %(time.time()-tt))
# # Get suitable roofs
tt = time.time()
roof_info['suitable'] = ( ( (roof_info.NEIGUNG <= FLAT_TILT) | (abs(roof_info.AUSRICHTUNG) <= MAX_ASPECT) )
& (roof_info.available_area >= 8) ).astype(int)
# # Load prediction
tilted_irrad = xr.open_mfdataset( os.path.join(PV_PATH, GT_IN), combine = 'by_coords' )
pv_potential = xr.open_mfdataset( os.path.join(PV_PATH, PV_IN), combine = 'by_coords' )
roof_info = roof_info.merge( tilted_irrad.yearly_kWh_m2.to_series().rename('Gt_kWh_m2'), on = 'DF_UID' )
roof_info = roof_info.merge( pv_potential.yearly_PV_kWh.to_series().rename('EPV_kWh'), on = 'DF_UID' )
print('Processed annual PV data in %.2fs\n' %(time.time()-tt))
# ### Save rooftop data
tt = time.time()
roof_save = roof_info.drop(columns = ['NEIGUNG', 'AUSRICHTUNG', 'NAME'])
print(roof_save.head())
roof_save.to_csv( os.path.join(OUTPATH, DF_OUT), index = False )
print('\nSaved annual PV data in %.2fs' %(time.time()-tt))
c4science · Help