virusw_special_reduction.py 文件源码

def get_master_sky(wavelength, spectrum, fiber_to_fiber, path, exptime):
    masterwave = []
    masterspec = []
    for wave, spec, ftf in zip(wavelength,spectrum,fiber_to_fiber):
        masterwave.append(wave)
        masterspec.append(np.where(ftf>1e-8, spec/ftf, 0.0))
    masterwave = np.hstack(masterwave)
    ind = np.argsort(masterwave)
    masterwave[:] = masterwave[ind]
    masterspec = np.hstack(masterspec)
    masterspec[:] = masterspec[ind]
    mastersky = medfilt(masterspec, 281)
    wv = np.arange(masterwave.min(),masterwave.max()+0.05,0.05)
    s = np.zeros((len(wv),2))
    s[:,0] = wv
    s[:,1] = np.interp(wv, masterwave, mastersky / exptime)    
    np.savetxt(op.join(path, 'sky_model.txt'), s)