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)
评论列表
文章目录