from sctrace.rawtrace import Record def display_record(t, trace): plot(t, trace) xlabel('time (ms)') ylabel('Current (pA)') title('Cluster example') grid(True) show() if __name__ == "__main__": filename = "./sctrace/samples/cluster.abf" record = Record(filename, filter_f=3000) end = len(record.trace) * record.dt t = arange(0.0, end, record.dt) print('trace: ', record.trace) print('time: ', t) display_record(t, record.trace) new_segement = record.slice(0.1, 1.15, dtype='time') new_cluster = new_segement.find_cluster() popen = new_cluster.Popen() print(new_cluster) print(popen) print(new_cluster.open_level) #plot(new_cluster.display_trace())
def setUp(self): self.filename = "./samples/cluster.abf" self.record = Record(self.filename)