baseFilterbankFilename, extFilterbankFilename = os.path.splitext(filterbankFilename)
plt.savefig(baseFilterbankFilename + "_dynamic_" + str(dynamic[0]) + "-" + str(dynamic[-1]) + ".png")
sys.exit(0)
elif opts.timeseries:
timeseries = opts.timeseries.split()
channel = int(timeseries[0])
begin = int(timeseries[1])
end = int(timeseries[2])
#print begin,end,channel
if (end >= numberSpectra):
end = numberSpectra-1
print "Selected spectrum exceeds available number of spectra!"
if (channel >= numberChannels):
channel = numberChannels-1
print "Selected channel exceeds available number of channels!"
timeseriesChannel = filterbankFile.getChan(channel)
timeseriesChannel.shape
timeseriesSelected = timeseriesChannel[begin:end]
plt.plot(timeseriesSelected)
plt.xlabel("Spectrum")
plt.ylabel("Intensity (a.u.)")
if not opts.savePlot:
plt.show()
else:
baseFilterbankFilename, extFilterbankFilename = os.path.splitext(filterbankFilename)
plt.savefig(baseFilterbankFilename + "_timeseries_ch" + str(channel) + "_" + str(begin) + "-" + str(end) + ".png")
sys.exit(0)
elif opts.bandpass:
totalBandpass = filterbankFile.bandpass() # calculate bandpass from entire observation
plt.plot(totalBandpass)
plt.xlabel("Channel")
