virusw_special_reduction.py 文件源码

python
阅读 24 收藏 0 点赞 0 评论 0

项目:Panacea 作者: grzeimann 项目源码 文件源码
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)
评论列表
文章目录


问题


面经


文章

微信
公众号

扫码关注公众号