예제 #1
0
파일: fsua.py 프로젝트: amerand/PRIMA
def findFringes(filename):
    """
    Find fringes in the scanning file 'filename'. Plots the results and returns
    the value of the offset.
    """
    f = pyfits.open(filename)
    fringes = (f['IMAGING_DATA_FSUA'].data['DATA1']-f['IMAGING_DATA_FSUA'].data['DATA3'])[:,0]
    snr = f['IMAGING_DATA_FSUA'].data['PDSNR']
    opd = np.interp(f['IMAGING_DATA_FSUA'].data['TIME'],
                    f['OPDC'].data['TIME'], f['OPDC'].data['RTOFFSET'])
    s = opd.argsort()
    fringes = fringes[s]
    snr = snr[s]
    opd = opd[s]
    
    plt.figure(0)
    plt.clf()
    ax=plt.subplot(211)
    plt.plot(opd, fringes, label='A-C')
    plt.legend()
    plt.subplot(212, sharex=ax)
    rms = slidop.slidingStd(opd,fringes, 10e-6)
    rms/= rms.mean()
    plt.plot(opd, rms, 'r', label='A-C RMS')
    plt.plot(opd, snr, 'g', label='PDSNR')
    plt.legend()
    f.close()
    return opd[rms.argmax()]
예제 #2
0
파일: fsua.py 프로젝트: amerand/PRIMA
def ftkPerfo(filename, minState=5):
    """
    estimate OPDC performance as OPD residuals
    """
    f = pyfits.open(filename)
    x = f['OPDC'].data.field('TIME')*1e-6
    y = f['OPDC'].data.field('UWPHASE')
    s = f['OPDC'].data.field('STATE')
    f.close()
     
    x = x[np.where(s>=minState)]
    y = y[np.where(s>=minState)]
    
    y = y[np.where(1-np.isnan(x))]
    x = x[np.where(1-np.isnan(x))]
    
    x = x[np.where(1-np.isnan(y))]
    y = y[np.where(1-np.isnan(y))]
    
    x = x[::3]
    y = y[::3]
   
    plt.figure(0)
    plt.clf()
    #plt.plot(x,y,'+', label='raw')
    all_dx = np.logspace(-1,0.5,10)
    jitt = []
    for dx in all_dx:
        try:
            yp = slidop.slidingStd(1.0*x,1.0*y,dx)*2.2/(2*np.pi)
        except:
            print np.diff(x)
            y_,yp,x_ = slidop.sliding_avg_rms(y*2.2/(2*np.pi),x,dx)
            print yp
            yp = np.interp(x,x_,yp)
        jitt.append(yp.mean())
    plt.plot(all_dx,jitt, 'ok-', label='residual jitter')
    plt.ylabel('opd residual RMS (um)')
    plt.xlabel('time window')
    plt.xscale('log')
    plt.vlines(21*0.040, plt.ylim()[0], plt.ylim()[1], label='MIDI frame')
    plt.legend()
    #plt.yscale('log')
    return