blimpy/guppi.py
def plot_spectrum(self, filename=None, plot_db=True):
        """ Do a (slow) numpy FFT and take power of data """
        header, data = self.read_next_data_block()

        print("Computing FFT...")
        d_xx_fft = np.abs(np.fft.fft(data[..., 0]))
        d_xx_fft = d_xx_fft.flatten()

        # Rebin to max number of points
        dec_fac_x = 1
        if d_xx_fft.shape[0] > MAX_PLT_POINTS:
            dec_fac_x = d_xx_fft.shape[0] / MAX_PLT_POINTS

        d_xx_fft = rebin(d_xx_fft, dec_fac_x)

        print("Plotting...")
        if plot_db:
            plt.plot(10 * np.log10(d_xx_fft))
            plt.ylabel("Power [dB]")
        else:
            plt.plot(d_xx_fft)
            plt.ylabel("Power")
            plt.xlabel("Channel")
        plt.title(self.filename)
        if filename:
            plt.savefig(filename)
        plt.show()
def plot_spectrum(self, filename=None, plot_db=True):
        """ Do a (slow) numpy FFT and take power of data """
        header, data = self.read_next_data_block()

        print("Computing FFT...")
        d_xx_fft = np.abs(np.fft.fft(data[..., 0]))
        d_xx_fft = d_xx_fft.flatten()

        # Rebin to max number of points
        dec_fac_x = 1
        if d_xx_fft.shape[0] > MAX_PLT_POINTS:
            dec_fac_x = d_xx_fft.shape[0] / MAX_PLT_POINTS

        d_xx_fft = rebin(d_xx_fft, dec_fac_x)

        print("Plotting...")
        if plot_db:
            plt.plot(10 * np.log10(d_xx_fft))
            plt.ylabel("Power [dB]")
        else:
            plt.plot(d_xx_fft)
            plt.ylabel("Power")
            plt.xlabel("Channel")
        plt.title(self.filename)
        if filename:
            plt.savefig(filename)
        plt.show()
def plot_spectrum ( self , filename = None , plot_db = True ) : header , data = self . read_next_data_block ( ) print ( "Computing FFT..." ) d_xx_fft = np . abs ( np . fft . fft ( data [ ... , 0 ] ) ) d_xx_fft = d_xx_fft . flatten ( ) # Rebin to max number of points dec_fac_x = 1 if d_xx_fft . shape [ 0 ] > MAX_PLT_POINTS : dec_fac_x = d_xx_fft . shape [ 0 ] / MAX_PLT_POINTS d_xx_fft = rebin ( d_xx_fft , dec_fac_x ) print ( "Plotting..." ) if plot_db : plt . plot ( 10 * np . log10 ( d_xx_fft ) ) plt . ylabel ( "Power [dB]" ) else : plt . plot ( d_xx_fft ) plt . ylabel ( "Power" ) plt . xlabel ( "Channel" ) plt . title ( self . filename ) if filename : plt . savefig ( filename ) plt . show ( )
Do a (slow) numpy FFT and take power of data
Do a ( slow ) numpy FFT and take power of data