convert calcium trace to spikes
Deconvolve calcium imaging data to spikes. A wrapper around Chenkov and McColgan so it’s easier to use their model and weights. Assuming you have a imageJ roi measurement csv file, and a known sample_rate of your two photon recording:
import numpy as np from ca2spike import deconvolve calcium_data = np.genfromtxt(file_path, delimiter=',')[1:, 1::4].T spike_data, timestamps, new_sample_rate = deconvolve(data, sample_rate)
Here the new sample rate is just a rounded int, coz otherwise your analysis experience will suck.