Source code for bmtk.simulator.filternet.modules.record_rates

import os
import csv
import pandas as pd
import h5py
import numpy as np

from .base import SimModule


[docs]class RecordRates(SimModule): def __init__(self, csv_file=None, h5_file=None, tmp_dir='output'): csv_file = csv_file if csv_file is None or os.path.isabs(csv_file) else os.path.join(tmp_dir, csv_file) self._save_to_csv = csv_file is not None self._tmp_csv_file = csv_file if self._save_to_csv else os.path.join(tmp_dir, '__tmp_rates.csv') self._tmp_csv_fhandle = open(self._tmp_csv_file, 'w') self._tmp_csv_writer = csv.writer(self._tmp_csv_fhandle, delimiter=' ') self._tmp_csv_writer.writerow(['node_id', 'population', 'timestamps', 'firing_rates']) h5_file = h5_file if h5_file is None or os.path.isabs(h5_file) else os.path.join(tmp_dir, h5_file) self._save_to_h5 = h5_file is not None self._h5_file = h5_file self._timestamps = None
[docs] def save(self, sim, cell, times, rates): if self._timestamps is None: self._timestamps = times for t, r in zip(times, rates): self._tmp_csv_writer.writerow([cell.gid, cell.population, t, r]) self._tmp_csv_fhandle.flush()
[docs] def finalize(self, sim): self._tmp_csv_fhandle.flush() self._tmp_csv_fhandle.close() if self._save_to_h5: try: rates_df = pd.read_csv(self._tmp_csv_file, sep=' ') rates_h5 = h5py.File(self._h5_file, 'w') rates_grp = rates_h5.create_group('/firing_rates') for pop, pop_table in rates_df.groupby('population'): pop_grp = rates_grp.create_group(pop) node_ids = pop_table['node_id'].unique().astype(np.uint) n_nodes = len(node_ids) n_timestamps = len(self._timestamps) pop_grp.create_dataset('mapping/node_ids', data=node_ids) pop_grp.create_dataset('mapping/timestamps', data=self._timestamps) pop_grp.create_dataset( 'data', data=np.reshape(pop_table['firing_rates'].values.astype(np.float), (n_nodes, n_timestamps)).T ) except Exception as e: print(e) print('Unable to save rates to hdf5') if not self._save_to_csv: os.remove(self._tmp_csv_file)