def importBandpass(filename, name):
'''Given the name of a bandpass read it from file into Bandpass class
assumes filters are differentiated by #NAME filter_name filter comments'''
filter = Bandpass()
data = []
with open(filename, 'r') as input:
for line in input:
line = line.strip().split()
if (line[0].startswith("#NAME") and (line[1] == name)):
for line in input:
if line.startswith("#NAME"):
break
else:
data.append(map(float, line.split()))
data = np.asarray(data).T
# convert to nanometers and set limits for wavelength
filter.setWavelenLimits(data[0].min() / 10., data[0].max() / 10., 1.)
try:
filter.setBandpass(data[0] / 10., data[1])
except IndexError:
print "Unable to find data values for ", name
exit()
#resample to
filter.resampleBandpass(wavelen_min=300.0,
wavelen_max=1150.0,
wavelen_step=0.1)
return filter
