{ "cells": [ { "cell_type": "code", "execution_count": 1, "metadata": {}, "outputs": [], "source": [ "import pandas as pd\n", "import xarray as xr\n", "import numpy as np\n", "import geopandas as gpd\n", "\n", "import os\n", "import time" ] }, { "cell_type": "code", "execution_count": 2, "metadata": {}, "outputs": [], "source": [ "fp = '/Users/alinawalch/Documents/OneDrive - epfl.ch/Master Projects/Romain_Sibuet/inputs/'" ] }, { "cell_type": "code", "execution_count": 3, "metadata": {}, "outputs": [ { "name": "stdout", "output_type": "stream", "text": [ "hyd_mod_return_temperature.csv\n", "ashp_hourly_operation.csv\n", "ashp_hourly_dhw_cop.csv\n", "gshp_hourly_dhw_cop.csv\n", "eh_hourly_operation.csv\n", "xhp_hourly_maximum_load.csv\n", "ashp_hourly_sh_cop.csv\n", "hyd_mod_supply_temperature.csv\n", "community.csv\n", "gshp_hourly_sh_cop.csv\n", "pv_hourly_potential.csv\n", "elec_hourly_demand.csv\n", "gshp_hourly_operation.csv\n" ] } ], "source": [ "data = {}\n", "for file in os.listdir(fp):\n", " name = os.path.splitext(file)[0]\n", " \n", " if name == 'heat_hourly_demand':\n", " new_data = pd.read_csv(os.path.join(fp, file), index_col = 0, skiprows = 1)\n", " data['dhw_demand_W'] = new_data.iloc[:,:int(new_data.shape[1]/2)]\n", " data['sh_demand_W'] = new_data.iloc[:,int(new_data.shape[1]/2):]\n", " continue\n", " \n", " data[name] = pd.read_csv(os.path.join(fp, file), index_col = 0)\n", " \n", " print(file)" ] }, { "cell_type": "code", "execution_count": 4, "metadata": { "scrolled": false }, "outputs": [], "source": [ "data_xr = []\n", "for name in data.keys():\n", " if name == 'community':\n", " data[ name ].index = data[ name ].index.astype('str')\n", " data_xr.append( data[ name ].to_xarray() )\n", " continue\n", " \n", " data[ name ].index.name = 'timestamp_ID'\n", " for col in data[name].columns:\n", " data[name] = data[name].rename({ col : col.split('.')[0] }, axis = 1)\n", " \n", " data_xr.append( data[ name ].reset_index()\n", " .melt(var_name = 'EGID', value_name = name, id_vars = 'timestamp_ID')\n", " .set_index(['EGID', 'timestamp_ID']).to_xarray() )\n", "data_xr = xr.merge(data_xr)" ] }, { "cell_type": "code", "execution_count": 5, "metadata": { "scrolled": false }, "outputs": [ { "data": { "text/html": [ "
\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "
<xarray.Dataset>\n",
       "Dimensions:                     (EGID: 518, timestamp_ID: 288)\n",
       "Coordinates:\n",
       "  * EGID                        (EGID) object '1028823' '1028824' ... '9082379'\n",
       "  * timestamp_ID                (timestamp_ID) int64 0 1 2 3 ... 284 285 286 287\n",
       "Data variables:\n",
       "    hyd_mod_return_temperature  (EGID, timestamp_ID) float64 27.51 ... 29.86\n",
       "    ashp_hourly_operation       (EGID, timestamp_ID) float64 nan nan ... 362.3\n",
       "    ashp_hourly_dhw_cop         (EGID, timestamp_ID) float64 nan nan ... 2.633\n",
       "    gshp_hourly_dhw_cop         (EGID, timestamp_ID) float64 2.615 2.615 ... nan\n",
       "    eh_hourly_operation         (EGID, timestamp_ID) float64 0.0 0.0 ... nan nan\n",
       "    xhp_hourly_maximum_load     (EGID, timestamp_ID) float64 3.0 3.0 ... nan nan\n",
       "    ashp_hourly_sh_cop          (EGID, timestamp_ID) float64 nan nan ... 3.273\n",
       "    dhw_demand_W                (EGID, timestamp_ID) float64 8.644 ... 424.1\n",
       "    sh_demand_W                 (EGID, timestamp_ID) float64 797.3 ... 2.185e+03\n",
       "    hyd_mod_supply_temperature  (EGID, timestamp_ID) float64 24.85 ... 26.26\n",
       "    sector                      (EGID) object 'RES' 'RES' 'RES' ... 'RES' 'SER'\n",
       "    BLD_ID                      (EGID) float64 1.5e+08 1.5e+08 ... 1.5e+08\n",
       "    footprint                   (EGID) float64 97.0 97.0 96.0 ... 64.0 167.0\n",
       "    N_floors                    (EGID) float64 1.0 1.0 1.0 1.0 ... 3.0 2.0 1.0\n",
       "    floor_area                  (EGID) float64 97.0 97.0 96.0 ... 128.0 167.0\n",
       "    geometry                    (EGID) object 'POINT (2496755.5 1117207.75)' ...\n",
       "    dhw_demand_kWh              (EGID) float64 2.016e+03 2.16e+03 ... 1.102e+04\n",
       "    sh_demand_kWh               (EGID) float64 8.577e+03 9.189e+03 ... 1.104e+04\n",
       "    Phi                         (EGID) float64 3.0 3.0 4.0 7.0 ... 2.0 2.0 4.0\n",
       "    K_rad                       (EGID) float64 100.0 100.0 133.3 ... 66.67 133.3\n",
       "    geo_potential_kWh           (EGID) float64 3.035e+04 1.482e+04 ... 1.554e+04\n",
       "    elec_demand_kWh             (EGID) float64 3.699e+03 3.699e+03 ... 0.0\n",
       "    pv_potential_kWh            (EGID) float64 1.239e+04 1.106e+04 ... 7.233e+03\n",
       "    radiation_kWh_m2            (EGID) float64 2.238e+03 2.041e+03 ... 1.251e+03\n",
       "    available_area              (EGID) float64 84.46 82.14 85.59 ... 6.343 42.85\n",
       "    N_roofs                     (EGID) int64 2 2 2 2 2 2 2 2 ... 0 1 0 3 1 0 0 1\n",
       "    roof_surface                (EGID) float64 144.4 146.2 140.0 ... 72.34 166.1\n",
       "    tes_capacity_L              (EGID) float64 400.0 400.0 ... 800.0 2.4e+03\n",
       "    gshp_hourly_sh_cop          (EGID, timestamp_ID) float64 4.472 4.428 ... nan\n",
       "    pv_hourly_potential         (EGID, timestamp_ID) float64 0.0 0.0 ... 0.0 0.0\n",
       "    elec_hourly_demand          (EGID, timestamp_ID) float64 298.2 ... 1.472e+03\n",
       "    gshp_hourly_operation       (EGID, timestamp_ID) float64 181.6 189.9 ... nan
" ], "text/plain": [ "\n", "Dimensions: (EGID: 518, timestamp_ID: 288)\n", "Coordinates:\n", " * EGID (EGID) object '1028823' '1028824' ... '9082379'\n", " * timestamp_ID (timestamp_ID) int64 0 1 2 3 ... 284 285 286 287\n", "Data variables:\n", " hyd_mod_return_temperature (EGID, timestamp_ID) float64 27.51 ... 29.86\n", " ashp_hourly_operation (EGID, timestamp_ID) float64 nan nan ... 362.3\n", " ashp_hourly_dhw_cop (EGID, timestamp_ID) float64 nan nan ... 2.633\n", " gshp_hourly_dhw_cop (EGID, timestamp_ID) float64 2.615 2.615 ... nan\n", " eh_hourly_operation (EGID, timestamp_ID) float64 0.0 0.0 ... nan nan\n", " xhp_hourly_maximum_load (EGID, timestamp_ID) float64 3.0 3.0 ... nan nan\n", " ashp_hourly_sh_cop (EGID, timestamp_ID) float64 nan nan ... 3.273\n", " dhw_demand_W (EGID, timestamp_ID) float64 8.644 ... 424.1\n", " sh_demand_W (EGID, timestamp_ID) float64 797.3 ... 2.185e+03\n", " hyd_mod_supply_temperature (EGID, timestamp_ID) float64 24.85 ... 26.26\n", " sector (EGID) object 'RES' 'RES' 'RES' ... 'RES' 'SER'\n", " BLD_ID (EGID) float64 1.5e+08 1.5e+08 ... 1.5e+08\n", " footprint (EGID) float64 97.0 97.0 96.0 ... 64.0 167.0\n", " N_floors (EGID) float64 1.0 1.0 1.0 1.0 ... 3.0 2.0 1.0\n", " floor_area (EGID) float64 97.0 97.0 96.0 ... 128.0 167.0\n", " geometry (EGID) object 'POINT (2496755.5 1117207.75)' ...\n", " dhw_demand_kWh (EGID) float64 2.016e+03 2.16e+03 ... 1.102e+04\n", " sh_demand_kWh (EGID) float64 8.577e+03 9.189e+03 ... 1.104e+04\n", " Phi (EGID) float64 3.0 3.0 4.0 7.0 ... 2.0 2.0 4.0\n", " K_rad (EGID) float64 100.0 100.0 133.3 ... 66.67 133.3\n", " geo_potential_kWh (EGID) float64 3.035e+04 1.482e+04 ... 1.554e+04\n", " elec_demand_kWh (EGID) float64 3.699e+03 3.699e+03 ... 0.0\n", " pv_potential_kWh (EGID) float64 1.239e+04 1.106e+04 ... 7.233e+03\n", " radiation_kWh_m2 (EGID) float64 2.238e+03 2.041e+03 ... 1.251e+03\n", " available_area (EGID) float64 84.46 82.14 85.59 ... 6.343 42.85\n", " N_roofs (EGID) int64 2 2 2 2 2 2 2 2 ... 0 1 0 3 1 0 0 1\n", " roof_surface (EGID) float64 144.4 146.2 140.0 ... 72.34 166.1\n", " tes_capacity_L (EGID) float64 400.0 400.0 ... 800.0 2.4e+03\n", " gshp_hourly_sh_cop (EGID, timestamp_ID) float64 4.472 4.428 ... nan\n", " pv_hourly_potential (EGID, timestamp_ID) float64 0.0 0.0 ... 0.0 0.0\n", " elec_hourly_demand (EGID, timestamp_ID) float64 298.2 ... 1.472e+03\n", " gshp_hourly_operation (EGID, timestamp_ID) float64 181.6 189.9 ... nan" ] }, "execution_count": 5, "metadata": {}, "output_type": "execute_result" } ], "source": [ "data_xr" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "# Select EGID" ] }, { "cell_type": "code", "execution_count": null, "metadata": {}, "outputs": [], "source": [] }, { "cell_type": "code", "execution_count": null, "metadata": {}, "outputs": [], "source": [] }, { "cell_type": "code", "execution_count": null, "metadata": {}, "outputs": [], "source": [] } ], "metadata": { "kernelspec": { "display_name": "Python 3", "language": "python", "name": "python3" }, "language_info": { "codemirror_mode": { "name": "ipython", "version": 3 }, "file_extension": ".py", "mimetype": "text/x-python", "name": "python", "nbconvert_exporter": "python", "pygments_lexer": "ipython3", "version": "3.7.6" } }, "nbformat": 4, "nbformat_minor": 4 }