예제 #1
0
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)
예제 #2
0
파일: fdm.py 프로젝트: jcapriot/src
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)
예제 #3
0
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'])
예제 #4
0
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')
예제 #5
0
파일: wplot.py 프로젝트: nicklinyi/src
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')
예제 #6
0
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)
예제 #7
0
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'])
예제 #8
0
파일: fdm.py 프로젝트: 1014511134/src
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)
예제 #9
0
파일: fdd.py 프로젝트: 1014511134/src
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'])
예제 #10
0
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'])
예제 #11
0
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']) 
예제 #12
0
파일: fdd.py 프로젝트: 1014511134/src
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'])
예제 #13
0
파일: wplot.py 프로젝트: nicklinyi/src
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)
예제 #14
0
파일: wplot.py 프로젝트: nicklinyi/src
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)
예제 #15
0
파일: wplot.py 프로젝트: nicklinyi/src
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)
예제 #16
0
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)
예제 #17
0
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)
예제 #18
0
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)