3: [1,1,2,2]
}
from gbtpy.makecube import velo_iterator
velo = velo_iterator(datapfits,linefreq=None).next()
interp_polyorder=15
scanrange = [9,54]
off_templates = {}
for percentile in (50,95,99):
for ii,exclude in enumerate(([-58,103,158,193],[-10,70],[])):
off_template,off_template_in,off_poly = make_off_template.make_off(hdul, scanrange=scanrange,
exclude_velo=exclude, interp_vrange=[-150,250],
interp_polyorder=interp_polyorder, sampler=sampler, return_uninterp=True,
return_poly=True,
percentile=percentile, exclude_spectral_ends=10)
off_templates[(percentile,ii)] = off_template,off_template_in,off_poly
oti_med = off_templates[(50,0)][1]
pl.figure(1)
pl.clf()
for ot in off_templates:
off_template,off_template_in,off_poly = off_templates[ot]
ax1 = pl.subplot(4,2,ot[1]+1)
L, = ax1.plot(velo, off_template)
ax1.plot(velo, off_template_in, label=ot, color=L.get_color())
ax1.set_ylim(0.915,1.0)
filename = '/Users/adam/observations/gbt/AGBT12B_221_01/AGBT12B_221_01.raw.acs.fits'
filepyfits = pyfits.open(filename,memmap=True)
datapfits = filepyfits[1].data
dataarr = datapfits.DATA
from gbtpy.makecube import velo_iterator
velo = velo_iterator(datapfits,linefreq=None).next()
sampler = 'A9'
off_templates = {}
for percentile in (50,95,99):
for ii,exclude in enumerate(([-58,103,158,193],[-127,103,158,193],[-127,43,69,103,158,193],[-127,-96,-78,38,62,96,158,193])):
off_template,off_template_in,off_poly = make_off_template.make_off(filepyfits, scanrange=scanrange,
exclude_velo=exclude, interp_vrange=[-350,350],
interp_polyorder=10, sampler=sampler, return_uninterp=True,
return_poly=True,
percentile=percentile, exclude_spectral_ends=10)
off_templates[(percentile,ii)] = off_template,off_template_in,off_poly
oti_med = off_templates[(50,0)][1]
pl.figure(1)
pl.clf()
for ot in off_templates:
off_template,off_template_in,off_poly = off_templates[ot]
ax1 = pl.subplot(4,2,ot[1]+1)
L, = ax1.plot(velo, off_template)
ax1.plot(velo, off_template_in, label=ot, color=L.get_color())
ax1.set_ylim(0.89,0.98)
