def ewavelet(wavelet,custom,par):
for i in range(2):
Plot(wavelet+'-'+str(i+1),wavelet,
'window n2=1 f2=%d | transp | window |'%i +
waveplot('%d %s'% (i,custom) ,par))
pplot.p1x2(wavelet,wavelet+'-1',wavelet+'-2',0.5,0.5,-11.5)
def eframe(frame,
movie,
index,
custom,
axis,
par,
xscale=0.75,
yscale=0.75,
shift=-8.25):
Flow([movie + '-plt', movie + '-bar'], movie,
'byte bar=${TARGETS[1]} gainpanel=a pclip=100 %s' % custom)
for i in range(2):
Plot(frame + '-' + str(i), movie + '-plt',
'window n3=1 f3=%d n4=1 f4=%d |' % (i, index) + cgrey('', par))
# Result(frame+'-'+str(i),movie+'-plt',
# 'window n3=1 f3=%d n4=1 f4=%d |' % (i,index)
# + cgrey('',par))
# if(axis==1):
# pplot.p2x1(frame,frame+'-1',frame+'-0',yscale,xscale,shift)
# else:
# pplot.p1x2(frame,frame+'-0',frame+'-1',yscale,xscale,shift)
pplot.p1x2(frame, frame + '-0', frame + '-1', yscale, xscale, shift)
def cliptogetherK(plot,file1,file2,title1,title2,axis,custom,par):
Flow(plot+'-all',[file1,file2],
'''
cat axis=3 space=n ${SOURCES[1]} |
byte gainpanel=a pclip=100|
put label1=%(lkz)s label2=%(lkx)s
unit1=%(ukz)s unit2=%(lkx)s
o1=%(okz)f d1=%(dkz)f o2=%(okx)f d2=%(dkx)f
''' %par)
if(axis==1):
if(not par.has_key('ys')): par['ys']=0.75
if(not par.has_key('xs')): par['xs']=0.75
if(not par.has_key('xc')): par['xc']=-8.25
Plot(file1,plot+'-all','window n3=1 f3=0 |' + fdmod.cgrey(custom+' title=%s'%title1,par))
Plot(file2,plot+'-all','window n3=1 f3=1 |' + fdmod.cgrey(custom+' title=%s label1= unit1= '%title2,par))
pplot.p1x2(plot,file1,file2,par['ys'],par['xs'],par['xc'])
else:
if(not par.has_key('ys')): par['ys']=.75
if(not par.has_key('xs')): par['xs']=.75
if(not par.has_key('xc')): par['xc']=-10
Plot(file1,plot+'-all','window n3=1 f3=1 |' + fdmod.cgrey(custom+' title=%s wantaxis2=n label2= unit2='%title1,par))
Plot(file2,plot+'-all','window n3=1 f3=0 |' + fdmod.cgrey(custom+' title=%s '%title2,par))
p2x1(plot,file1,file2,par['ys'],par['xs'],par['xc'])
def waveplotE2d(wav, custom, par):
Plot(wav + '_V', wav,
'window n2=1 f2=0 | transp | window |' + waveplot(custom, par))
Plot(wav + '_H', wav,
'window n2=1 f2=1 | transp | window |' + waveplot(custom, par))
pplot.p1x2(wav, wav + '_V', wav + '_H', 0.5, 0.5, -11.5)
Result(wav, wav, 'Overlay')
def waveplotE2d(wav,custom,par):
Plot(wav+'_V',wav,
'window n2=1 f2=0 | transp | window |'+
waveplot(custom,par))
Plot(wav+'_H',wav,
'window n2=1 f2=1 | transp | window |'+
waveplot(custom,par))
pplot.p1x2(wav,wav+'_V',wav+'_H',0.5,0.5,-11.5)
Result(wav,wav,'Overlay')
def edata(plot,data,custom,par):
Flow([plot+'_plt',plot+'_bar'],data,
'scale axis=123 | byte bar=${TARGETS[1]} gainpanel=a pclip=100 %s' % custom)
for i in range(2):
ctag = "-%d"%i
Plot( plot+ctag,[plot+'_plt',plot+'_bar'],
'window n2=1 f2=%d bar=${SOURCES[1]} | transp |' % i
+ dgrey('screenratio=1.75 screenht=20 %s' % custom,par))
pplot.p1x2(plot,plot+'-0',plot+'-1',0.5,0.5,-11)
def oneirST_box(inp, zdel, xdel, m, n, box, custom, par):
par['nzspk'] = 2 * m + 1
par['nxspk'] = 2 * n + 1
par['fzspk'] = (par['nzspk'] - 1) / 2 + 1
par['fxspk'] = (par['nxspk'] - 1) / 2 + 1
Flow(
inp, None, '''
spike nsp=1 mag=1
n1=%(nzspk)d o1=%(oz)g d1=%(dz)g k1=%(fzspk)d
n2=%(nxspk)d o2=%(ox)g d2=%(dx)g k2=%(fxspk)d
''' % par)
for j in (zdel, xdel):
ddapply(inp + '-' + j, inp, j, 'y', par)
Flow(
inp + '-all', [inp + '-' + zdel, inp + '-' + xdel], '''
cat axis=3 space=n ${SOURCES[1]} |
scale axis=123 |
byte gainpanel=a pclip=99.9|
put o1=0 d1=1 o2=0 d2=1 label1=Sample# label2=Sample# unit1= unit2=
%s
''' % custom)
Plot(
inp + '-' + zdel + '-', inp + '-all', 'window n3=1 f3=0 |' + opplot(
"\F5 L\_z ",
' wantaxis=y o2num=0 d2num=10 n2tic=6 o1num=0 d1num=10 n1tic=6 '
+ custom, par))
Plot(
inp + '-' + xdel + '-', inp + '-all', 'window n3=1 f3=1 |' + opplot(
"\F5 L\_x ",
' wantaxis=y o2num=0 d2num=10 n2tic=6 o1num=0 d1num=10 n1tic=6 unit1= label1= wantaxis1= '
+ custom, par))
# Plot(inp+'-'+zdel,inp+'-all','window n3=1 f3=0 |' + opplot('',''))
# Plot(inp+'-'+xdel,inp+'-all','window n3=1 f3=1 |' + opplot('',''))
Plot(inp + '-' + zdel, [inp + '-' + zdel + '-', box], 'Overlay')
Plot(inp + '-' + xdel, [inp + '-' + xdel + '-', box], 'Overlay')
if ('ys' not in par): par['ys'] = 0.75
if ('xs' not in par): par['xs'] = 0.75
if ('xc' not in par): par['xc'] = -8.25
pplot.p1x2(inp, inp + '-' + zdel, inp + '-' + xdel, par['ys'], par['xs'],
par['xc'])
def eframe(frame,movie,index,custom,axis,par,xscale=0.75,yscale=0.75,shift=-8.25):
Flow([movie+'-plt',movie+'-bar'],movie,
'byte bar=${TARGETS[1]} gainpanel=a pclip=100 %s' % custom)
for i in range(2):
Plot(frame+'-'+str(i),movie+'-plt',
'window n3=1 f3=%d n4=1 f4=%d |' % (i,index)
+ cgrey('',par))
# Result(frame+'-'+str(i),movie+'-plt',
# 'window n3=1 f3=%d n4=1 f4=%d |' % (i,index)
# + cgrey('',par))
if(axis==1):
pplot.p2x1(frame,frame+'-1',frame+'-0',yscale,xscale,shift)
else:
pplot.p1x2(frame,frame+'-0',frame+'-1',yscale,xscale,shift)
def oneirST_box(inp,zdel,xdel,m,n,box,custom,par):
par['nzspk']=2*m+1
par['nxspk']=2*n+1
par['fzspk']=(par['nzspk']-1)/2+1
par['fxspk']=(par['nxspk']-1)/2+1
Flow(inp,None,
'''
spike nsp=1 mag=1
n1=%(nzspk)d o1=%(oz)g d1=%(dz)g k1=%(fzspk)d
n2=%(nxspk)d o2=%(ox)g d2=%(dx)g k2=%(fxspk)d
''' % par)
for j in (zdel,xdel):
ddapply(inp+'-'+j,inp,j,'y',par)
Flow(inp+'-all',[inp+'-'+zdel,inp+'-'+xdel],
'''
cat axis=3 space=n ${SOURCES[1]} |
scale axis=123 |
byte gainpanel=a pclip=99.9|
put o1=0 d1=1 o2=0 d2=1 label1=Sample# label2=Sample# unit1= unit2=
%s
''' % custom)
Plot(inp+'-'+zdel+'-',inp+'-all','window n3=1 f3=0 |'
+ opplot("\F5 L\_z ",' wantaxis=y o2num=0 d2num=10 n2tic=6 o1num=0 d1num=10 n1tic=6 '+custom,par))
Plot(inp+'-'+xdel+'-',inp+'-all','window n3=1 f3=1 |'
+ opplot("\F5 L\_x ",' wantaxis=y o2num=0 d2num=10 n2tic=6 o1num=0 d1num=10 n1tic=6 unit1= label1= wantaxis1= '+custom,par))
# Plot(inp+'-'+zdel,inp+'-all','window n3=1 f3=0 |' + opplot('',''))
# Plot(inp+'-'+xdel,inp+'-all','window n3=1 f3=1 |' + opplot('',''))
Plot(inp+'-'+zdel,[inp+'-'+zdel+'-',box],'Overlay')
Plot(inp+'-'+xdel,[inp+'-'+xdel+'-',box],'Overlay')
if(not par.has_key('ys')): par['ys']=0.75
if(not par.has_key('xs')): par['xs']=0.75
if(not par.has_key('xc')): par['xc']=-8.25
pplot.p1x2(inp,inp+'-'+zdel,inp+'-'+xdel,par['ys'],par['xs'],par['xc'])
def oneirST_S(inp, zdel, xdel, m, n, custom, par):
par['nzspk'] = 2 * m + 1
par['nxspk'] = 2 * n + 1
par['fzspk'] = (par['nzspk'] - 1) / 2 + 1
par['fxspk'] = (par['nxspk'] - 1) / 2 + 1
Flow(
inp, None, '''
spike nsp=1 mag=1
n1=%(nzspk)d o1=%(oz)g d1=%(dz)g k1=%(fzspk)d
n2=%(nxspk)d o2=%(ox)g d2=%(dx)g k2=%(fxspk)d
''' % par)
for j in (zdel, xdel):
## ddapply(inp+'-'+j,inp,j,'y',par)
Flow(inp + '-' + j, inp, 'math output="input"')
Flow(
inp + '-all', [inp + '-' + zdel, inp + '-' + xdel], '''
cat axis=3 space=n ${SOURCES[1]} |
scale axis=123 |
byte gainpanel=a pclip=98
%s
''' % custom)
Plot(
inp + '-' + zdel, inp + '-all', 'window n3=1 f3=0 |' +
opplot("\F15 \\v148 \F3 /\F15 \\v148 \F3 z", custom, par))
Plot(
inp + '-' + xdel, inp + '-all', 'window n3=1 f3=1 |' +
opplot("\F15 \\v148 \F3 /\F15 \\v148 \F3 x", custom, par))
# Plot(inp+'-'+zdel,inp+'-all','window n3=1 f3=0 |' + opplot('',''))
# Plot(inp+'-'+xdel,inp+'-all','window n3=1 f3=1 |' + opplot('',''))
if ('ys' not in par): par['ys'] = 0.75
if ('xs' not in par): par['xs'] = 0.75
if ('xc' not in par): par['xc'] = -8.25
pplot.p1x2(inp, inp + '-' + zdel, inp + '-' + xdel, par['ys'], par['xs'],
par['xc'])
def cliptogether(plot,file1,file2,title1,title2,axis,custom,par):
Flow(plot+'-all',[file1,file2],
'''
cat axis=3 space=n ${SOURCES[1]} |
byte gainpanel=a %s
''' %custom)
if(axis==1):
if(not par.has_key('ys')): par['ys']=0.75
if(not par.has_key('xs')): par['xs']=0.75
if(not par.has_key('xc')): par['xc']=-8.25
Plot(file1,plot+'-all','window n3=1 f3=0 |' + fdmod.cgrey(custom+' title=%s'%title1,par))
Plot(file2,plot+'-all','window n3=1 f3=1 |' + fdmod.cgrey(custom+' title=%s wantaxis1=n label1= unit1= '%title2,par))
pplot.p1x2(plot,file1,file2,par['ys'],par['xs'],par['xc'])
else:
if(not par.has_key('ys')): par['ys']=.75
if(not par.has_key('xs')): par['xs']=.75
if(not par.has_key('yc')): par['yc']=-5
Plot(file1,plot+'-all','window n3=1 f3=1 |' + fdmod.cgrey(custom+' title=%s wantaxis2=n label2= unit2='%title1,par))
Plot(file2,plot+'-all','window n3=1 f3=0 |' + fdmod.cgrey(custom+' title=%s '%title2,par))
p2x1(plot,file1,file2,par['ys'],par['xs'],par['yc'])
def oneirST(inp,zdel,xdel,m,n,custom,par):
par['nzspk']=2*m+1
par['nxspk']=2*n+1
par['fzspk']=(par['nzspk']-1)/2+1
par['fxspk']=(par['nxspk']-1)/2+1
Flow(inp,None,
'''
spike nsp=1 mag=1
n1=%(nzspk)d o1=%(oz)g d1=%(dz)g k1=%(fzspk)d
n2=%(nxspk)d o2=%(ox)g d2=%(dx)g k2=%(fxspk)d
''' % par)
for j in (zdel,xdel):
ddapply(inp+'-'+j,inp,j,'y',par)
Flow(inp+'-all',[inp+'-'+zdel,inp+'-'+xdel],
'''
cat axis=3 space=n ${SOURCES[1]} |
scale axis=123 |
byte gainpanel=a pclip=98
%s
''' % custom)
Plot(inp+'-'+zdel,inp+'-all','window n3=1 f3=0 |'
+ opplot("\F5 L\_z",custom,par))
Plot(inp+'-'+xdel,inp+'-all','window n3=1 f3=1 |'
+ opplot("\F5 L\_x",custom,par))
# Plot(inp+'-'+zdel,inp+'-all','window n3=1 f3=0 |' + opplot('',''))
# Plot(inp+'-'+xdel,inp+'-all','window n3=1 f3=1 |' + opplot('',''))
if(not par.has_key('ys')): par['ys']=0.75
if(not par.has_key('xs')): par['xs']=0.75
if(not par.has_key('xc')): par['xc']=-8.25
pplot.p1x2(inp,inp+'-'+zdel,inp+'-'+xdel,par['ys'],par['xs'],par['xc'])
def dgreyE2d(data,dbyt,custom,par,xscale=0.5,yscale=0.5,shift=-11):
Plot(data+'_V',dbyt,'window n2=1 f2=0 | transp |'+ dgrey2d('',par))
Plot(data+'_H',dbyt,'window n2=1 f2=1 | transp |'+ dgrey2d('',par))
pplot.p1x2(data,data+'_V',data+'_H',xscale,yscale,shift)
def iovlE2d(out,inp,par,xscale=0.5,yscale=0.5,shift=-11):
Plot(out+'_V',inp,'Overlay')
Plot(out+'_H',inp,'Overlay')
pplot.p1x2(out,out+'_V',out+'_H',xscale,yscale,shift)
def ifrmE2d(wfrm,wbyt,index,custom,par,xscale=0.5,yscale=0.5,shift=-11):
Plot(wfrm+'_V',wbyt,'window n3=1 f3=0 |'+ ifrm2d(index,'',par))
Plot(wfrm+'_H',wbyt,'window n3=1 f3=1 |'+ ifrm2d(index,'',par))
pplot.p1x2(wfrm,wfrm+'_V',wfrm+'_H',xscale,yscale,shift)
def dgreyE2d(data, dbyt, custom, par, xscale=0.5, yscale=0.5, shift=-11):
Plot(data + '_V', dbyt, 'window n2=1 f2=0 | transp |' + dgrey2d('', par))
Plot(data + '_H', dbyt, 'window n2=1 f2=1 | transp |' + dgrey2d('', par))
pplot.p1x2(data, data + '_V', data + '_H', xscale, yscale, shift)
def iovlE2d(out, inp, par, xscale=0.5, yscale=0.5, shift=-11):
Plot(out + '_V', inp, 'Overlay')
Plot(out + '_H', inp, 'Overlay')
pplot.p1x2(out, out + '_V', out + '_H', xscale, yscale, shift)
def ifrmE2d(wfrm, wbyt, index, custom, par, xscale=0.5, yscale=0.5, shift=-11):
Plot(wfrm + '_V', wbyt, 'window n3=1 f3=0 |' + ifrm2d(index, '', par))
Plot(wfrm + '_H', wbyt, 'window n3=1 f3=1 |' + ifrm2d(index, '', par))
pplot.p1x2(wfrm, wfrm + '_V', wfrm + '_H', xscale, yscale, shift)