def test_decovnolution(path_to_config):
yass.set_config('tests/config_nnet.yaml')
clear_scores, spike_index_clear, spike_index_collision = preprocess.run()
(spike_train_clear, templates,
spike_index_collision) = process.run(clear_scores, spike_index_clear,
spike_index_collision)
deconvolute.run(spike_train_clear, templates, spike_index_collision)
clean_tmp()
def test_process(path_to_config):
yass.set_config(path_to_config)
clear_scores, spike_index_clear, spike_index_collision = preprocess.run()
(spike_train_clear, templates,
spike_index_collision) = process.run(clear_scores, spike_index_clear,
spike_index_collision)
def test_new_process(path_to_config):
cfg = yass.Config.from_yaml(path_to_config)
pp = Preprocessor(cfg)
score, clr_idx, spt = pp.process()
set_config(path_to_config)
(spike_train_clear, templates,
spike_index_collision) = process.run(score, clr_idx, spt)
def _run_pipeline(config,
output_file,
logger_level='INFO',
clean=True,
output_dir='tmp/',
complete=False):
"""
Run the entire pipeline given a path to a config file
and output path
"""
# load yass configuration parameters
set_config(config)
CONFIG = read_config()
ROOT_FOLDER = CONFIG.data.root_folder
TMP_FOLDER = path.join(ROOT_FOLDER, output_dir)
# remove tmp folder if needed
if os.path.exists(TMP_FOLDER) and clean:
shutil.rmtree(TMP_FOLDER)
# create TMP_FOLDER if needed
if not os.path.exists(TMP_FOLDER):
os.makedirs(TMP_FOLDER)
# load logging config file
logging_config = load_logging_config_file()
logging_config['handlers']['file']['filename'] = path.join(
TMP_FOLDER, 'yass.log')
logging_config['root']['level'] = logger_level
# configure logging
logging.config.dictConfig(logging_config)
# instantiate logger
logger = logging.getLogger(__name__)
# run preprocessor
(score, spike_index_clear,
spike_index_collision) = preprocess.run(output_directory=output_dir)
# run processor
(spike_train_clear, templates,
spike_index_collision) = process.run(score,
spike_index_clear,
spike_index_collision,
output_directory=output_dir)
# run deconvolution
spike_train = deconvolute.run(spike_train_clear,
templates,
spike_index_collision,
output_directory=output_dir)
# save metadata in tmp
path_to_metadata = path.join(TMP_FOLDER, 'metadata.yaml')
logging.info('Saving metadata in {}'.format(path_to_metadata))
save_metadata(path_to_metadata)
# save config.yaml copy in tmp/
path_to_config_copy = path.join(TMP_FOLDER, 'config.yaml')
shutil.copy2(config, path_to_config_copy)
logging.info('Saving copy of config: {} in {}'.format(
config, path_to_config_copy))
# save templates
path_to_templates = path.join(TMP_FOLDER, 'templates.npy')
logging.info('Saving templates in {}'.format(path_to_templates))
np.save(path_to_templates, templates)
path_to_spike_train = path.join(TMP_FOLDER, output_file)
np.save(path_to_spike_train, spike_train)
logger.info('Spike train saved in: {}'.format(path_to_spike_train))
# this part loads waveforms for all spikes in the spike train and scores
# them, this data is needed to later generate phy files
if complete:
STANDARIZED_PATH = path.join(TMP_FOLDER, 'standarized.bin')
PARAMS = load_yaml(path.join(TMP_FOLDER, 'standarized.yaml'))
# load waveforms for all spikes in the spike train
logger.info('Loading waveforms from all spikes in the spike train...')
explorer = RecordingExplorer(STANDARIZED_PATH,
spike_size=CONFIG.spikeSize,
dtype=PARAMS['dtype'],
n_channels=PARAMS['n_channels'],
data_format=PARAMS['data_format'])
waveforms = explorer.read_waveforms(spike_train[:, 0])
path_to_waveforms = path.join(TMP_FOLDER, 'spike_train_waveforms.npy')
np.save(path_to_waveforms, waveforms)
logger.info('Saved all waveforms from the spike train in {}...'.format(
path_to_waveforms))
# score all waveforms
logger.info('Scoring waveforms from all spikes in the spike train...')
path_to_rotation = path.join(TMP_FOLDER, 'rotation.npy')
rotation = np.load(path_to_rotation)
main_channels = explorer.main_channel_for_waveforms(waveforms)
path_to_main_channels = path.join(TMP_FOLDER,
'waveforms_main_channel.npy')
np.save(path_to_main_channels, main_channels)
logger.info('Saved all waveforms main channels in {}...'.format(
path_to_waveforms))
waveforms_score = dim_red.score(waveforms, rotation, main_channels,
CONFIG.neighChannels, CONFIG.geom)
path_to_waveforms_score = path.join(TMP_FOLDER, 'waveforms_score.npy')
np.save(path_to_waveforms_score, waveforms_score)
logger.info('Saved all scores in {}...'.format(path_to_waveforms))
# score templates
# TODO: templates should be returned in the right shape to avoid .T
templates_ = templates.T
main_channels_tmpls = explorer.main_channel_for_waveforms(templates_)
path_to_templates_main_c = path.join(TMP_FOLDER,
'templates_main_channel.npy')
np.save(path_to_templates_main_c, main_channels_tmpls)
logger.info('Saved all templates main channels in {}...'.format(
path_to_templates_main_c))
templates_score = dim_red.score(templates_, rotation,
main_channels_tmpls,
CONFIG.neighChannels, CONFIG.geom)
path_to_templates_score = path.join(TMP_FOLDER, 'templates_score.npy')
np.save(path_to_templates_score, templates_score)
logger.info(
'Saved all templates scores in {}...'.format(path_to_waveforms))
def test_new_process_shows_error_if_empty_config():
with pytest.raises(ValueError):
process.run(None, None, None)
import logging
import yass
from yass import preprocess
from yass import process
from yass import deconvolute
# configure logging module to get useful information
logging.basicConfig(level=logging.INFO)
# set yass configuration parameters
yass.set_config('tests/config_nnet.yaml')
# run preprocessor
score, clr_idx, spt = preprocess.run()
# run processor
spike_train_clear, templates, spike_index_collision = process.run(
score, clr_idx, spt)
# run deconvolution
spike_train = deconvolute.run(spike_train_clear, templates,
spike_index_collision)
import logging
import yass
from yass import preprocess
from yass import process
# configure logging module to get useful information
logging.basicConfig(level=logging.DEBUG)
# set yass configuration parameters
yass.set_config('config_sample.yaml')
# run preprocessor
score, clr_idx, spt = preprocess.run()
# run processor
(spike_train_clear, templates,
spike_index_collision) = process.run(score, clr_idx, spt)
import logging
import yass
from yass import preprocess
from yass import process
from yass import deconvolute
# configure logging module to get useful information
logging.basicConfig(level=logging.INFO)
# set yass configuration parameters
yass.set_config('tests/config_nnet.yaml')
# run preprocessor
score, clr_idx, spt = preprocess.run()
# run processor
spike_train, spikes_left, templates = process.run(score, clr_idx, spt)
# run deconvolution
spikes = deconvolute.run(spike_train, spikes_left, templates)
def test_decovnolute_new_pipeline(path_to_config):
yass.set_config('tests/config_nnet.yaml')
score, clr_idx, spt = preprocess.run()
spike_train, spikes_left, templates = process.run(score, clr_idx, spt)
deconvolute.run(spike_train, spikes_left, templates)