Homec4science
Diffusion solar_potential (master) Browse

solarPOT/Hourly_PV_public/Technical_potential/Panel_performance/efficiency_mapping.py
e573bf139de5master

Last Change
Raw File

efficiency_mapping.py
View Options

import numpy as np
import xarray as xr
import pandas as pd
from pvlib import pvsystem
import os
import time
import util
## === SPECIFY FILE PATHS === ##
tilted_irrad_fp = # TILTED RADIATION FILE PATH
performance_fp = # PANEL PERFORMANCE FILE PATH
T_max_mmd_fp = # MAX TEMPERATURE MMD FILE PATH
ROOFS_fp = # ROOFTOP CHARACTERISTICS FILE PATH
OUTFILE = # FILE PATH FOR OUTPUT FILE
## === LOAD DATA === ##
tilted_irrad = xr.open_mfdataset( fp, chunks = {'DF_UID' : 100000} )
performance = xr.open_dataset( fp )
T_max_mmd = xr.open_dataset( fp )
roofs = pd.read_csv(ROOFS_fp, usecols = ['DF_UID', 'XCOORD', 'YCOORD'])
OUTPATH = os.path.split( OUTFILE )[0]
try:
os.mkdir( os.path.join(OUTPATH, 'tmp') )
except:
print('tmp directory already exists')
# get pixel coordinates of roofs, by rounding to the coordinate of the closest pixel
roofs['chx'] = util.round_to_interval(roofs['XCOORD'] - 500, 1000) + 500
roofs['chy'] = util.round_to_interval(roofs['YCOORD'] - 500, 1000) + 500
## ==== MAP PANEL PERFORMANCE TO ROOFS ==== ##
# treat every timestamp independently
for i, t in enumerate(tilted_irrad.timestamp):
tt = time.time()
# get month for maximum temperature
month = int(t.dt.month)
# get tilted irradiance for given timestamp
Gt = tilted_irrad.sel(timestamp = t).tilted_irradiance.to_dataframe().reset_index()
Gt['G_t'] = np.round(Gt.tilted_irradiance)
# get maximum temperature
Tmax = T_max_mmd.TmaxD.sel(month = month).to_dataframe().dropna().reset_index().drop(columns = ['lon', 'lat', 'month'])
Tmax['T_amb'] = np.round(Tmax['TmaxD'])
# get efficiency dataframe
eff = performance.to_dataframe().reset_index().drop(columns = ['T_cell', 'losses', 'inverter_eff', 'PF'])
# create dataset with all merged information
tmp = roofs.merge(Tmax, on = ['chx', 'chy'], how = 'left')
tmp = tmp.merge(Gt, on = 'DF_UID', how = 'left')
roofs_w_eff = tmp.merge(eff, on = ['T_amb', 'G_t'], how = 'left')
# create xarray dataset from the relevant information and save intermediate file
eff_xr = roofs_w_eff.set_index(['DF_UID', 'timestamp'])[['TmaxD', 'module_eff', 'PF_mean']].to_xarray()
eff_xr.to_netcdf( os.path.join(OUTPATH, 'tmp', 'efficiency_mapped_%d.nc' %i) )
print('Finished iteration %d in %.3fs' %(i, time.time()-tt))
## ==== MERGE INTERMEDIATE FILES ==== ##
eff = []
for i in range(10):
eff.append( xr.open_mfdataset( os.path.join(OUTPATH, 'tmp', 'efficiency_mapped_%d*.nc' %i),
chunks = {'DF_UID': 100000}) )
eff_xr = xr.concat(eff, dim = 'timestamp').sortby('timestamp')
eff_xr.to_netcdf( OUTFILE )
c4science · Help