示例#1
0
from pylab import *

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())
示例#2
0
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())
示例#3
0
 def setUp(self):
     self.filename = "./samples/cluster.abf"
     self.record = Record(self.filename)