gwpy
gwpy copied to clipboard
gwpy
Package for analysing and characterising gravitational wave data in python
Using the ligo-py37 environment (in Jupyter on CIT): ```python from gwpy.timeseries import TimeSeries tmp3f = TimeSeries.find('H1:IMC-VCO_FREQUENCY.rms,m-trend',\ '2019-04-29 05:00:00',\ '2019-04-29 22:00:00',\ frametype='H1_M', nproc=6) print(tmp3f.duration, (22-5)*3600) ``` gives 64800.0 s 61200 It's...
In a discussion with @emaros, he has recommended not using `IFrameFStream.GetNumberOfFrames()` and https://github.com/gwpy/gwpy/blob/001c27fc4237b1c5c445c81ebeacffd9fe3bb76f/gwpy/timeseries/io/gwf/framecpp.py#L195-L208 to get the epoch of each frame. Instead he suggested using `FrTOC.GetGTime{S,N}` which will return arrays of...
I want to plot something from a frame builder behind a router, so I have opened up an SSH tunnel on a custom port. The `gwpy-plot` utility seems to lack...
To whom will concern it, I would like to ask that whether there is any constraint on the usage of the Q-transform from gwpy or pycbc? In the attached docx...
Searching the spectrogram for the minimum nonzero value exposed a functional style indexing problem with 2D GWpy objects. This snippet reproduces it. ```from gwpy.timeseries import TimeSeries import numpy as np...
It would be nice to have the ability to specify that missing data be padded in the CLI. Right now I don't see a way to do this.
```python In [1]: import gwpy In [2]: from gwpy.timeseries import TimeSeries In [3]: ts = TimeSeries.fetch('L1:GRD-IFO_OK', 1226370000, 1226370001) In [4]: ts == 1 Out[4]: False In [5]: ts == 2.0...
I would like to consider moving gwpy to a BSD-style license, away from GPL, mainly after reading [this post](http://nipy.sourceforge.net/software/license/johns_bsd_pitch.html) from John Hunter (creator of Matplotlib). cc @gwpy/gwpy-developers
I would like to deprecate the lalframe GWF API in `gwpy.timseries.io.gwf.lalframe` and double-down on using frameCPP as the API. Maintaining two different APIs is just more work. I don't want...
The code snippet below demonstrates the issue which traces to gwpy/signal/filter_design.py:245, 247 ``` zeros = numpy.array(zeros, dtype=float, copy=False) zeros = zeros[numpy.isfinite(zeros)] poles = numpy.array(poles, dtype=float, copy=False) ``` It produces the...