Revision 29335c53549d146e7c295b894972e47f42f465c4 authored by Collin Capano on 25 June 2020, 20:26:52 UTC, committed by GitHub on 25 June 2020, 20:26:52 UTC
* update the inference workflow docs * fix some typos
1 parent 2eab5a5
read.py
import pycbc.psd
import pycbc.types
import pylab
filename = 'example_psd.txt'
# The PSD will be interpolated to the requested frequency spacing
delta_f = 1.0 / 4
length = int(1024 / delta_f)
low_frequency_cutoff = 30.0
psd = pycbc.psd.from_txt(filename, length, delta_f,
low_frequency_cutoff, is_asd_file=False)
pylab.loglog(psd.sample_frequencies, psd, label='interpolated')
# The PSD will be read in without modification
psd = pycbc.types.load_frequencyseries('./example_psd.txt')
pylab.loglog(psd.sample_frequencies, psd, label='raw')
pylab.xlim(xmin=30, xmax=1000)
pylab.legend()
pylab.xlabel('Hz')
pylab.show()
# Save a psd to file, several formats are supported (.txt, .hdf, .npy)
psd.save('tmp_psd.txt')
Computing file changes ...
