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")