def PlotSCD( files='*.pk' ):
	import pylab
	if isinstance( files, basestring ):
		if os.path.isdir( files ): files = os.path.join( files, '*.pk' )
		files = glob.glob( files )
	d = DataFiles.load( files )
	r = SigTools.correlate( d['x'], d['y'], axis=0 )
	SigTools.imagesc( r*numpy.abs(r), y=d['channels'], x=SigTools.samples2msec( range( r.shape[1] ), d['fs'] ), aspect='auto', balance=0.0, colorbar=True )
	pylab.title( ', '.join( [ '%d: %d' % ( yi, ( d['y'] == yi ).sum() ) for yi in numpy.unique( d['y'] ) ] ) )
	pylab.draw()
Exemple #2
0
def PlotSCD(files='*.pk'):
    import pylab
    if isinstance(files, basestring):
        if os.path.isdir(files): files = os.path.join(files, '*.pk')
        files = glob.glob(files)
    d = DataFiles.load(files)
    r = SigTools.correlate(d['x'], d['y'], axis=0)
    SigTools.imagesc(r * numpy.abs(r),
                     y=d['channels'],
                     x=SigTools.samples2msec(range(r.shape[1]), d['fs']),
                     aspect='auto',
                     balance=0.0,
                     colorbar=True)
    pylab.title(', '.join([
        '%d: %d' % (yi, (d['y'] == yi).sum()) for yi in numpy.unique(d['y'])
    ]))
    pylab.draw()