Skip to content

Commit

Permalink
Merge pull request #142 from alurban/use-outseg
Browse files Browse the repository at this point in the history
Adjust the size of Q-transform spectrograms
  • Loading branch information
Alex L. Urban committed Sep 26, 2018
2 parents 5e0a43f + d334834 commit 93f8415
Showing 1 changed file with 8 additions and 6 deletions.
14 changes: 8 additions & 6 deletions bin/gwdetchar-omega
Original file line number Diff line number Diff line change
Expand Up @@ -33,6 +33,7 @@ use('agg') # nopep8
from gwpy.utils import gprint
from gwpy.time import tconvert
from gwpy.table import EventTable
from gwpy.segments import Segment
from gwpy.timeseries import TimeSeriesDict
from gwpy.detector import (Channel, ChannelList)
from gwpy.signal.qtransform import QTiling
Expand Down Expand Up @@ -355,15 +356,16 @@ for block in blocks[:]:
rtable = eventgram(gps, hpseries, frange=table.frange, qrange=(Q, Q),
snrthresh=c.snrthresh, mismatch=c.mismatch)

# compute Q-transforms
# compute Q-transform spectrograms
outseg = Segment(gps - max(c.pranges)/2., gps + max(c.pranges)/2.)
tres = min(c.pranges) / 500
fres = c.frange[0] / 5
qscan = wseries.q_transform(qrange=(Q, Q), frange=c.frange,
tres=tres, fres=fres, gps=gps,
search=0.25, whiten=False)
fres = c.frange[0] / 20
qscan = wseries.q_transform(qrange=(Q, Q), frange=c.frange, tres=tres,
fres=fres, gps=gps, search=0.25,
whiten=False, outseg=outseg)
rqscan = hpseries.q_transform(qrange=(Q, Q), frange=c.frange,
tres=tres, fres=fres, gps=gps,
search=0.25, whiten=False)
search=0.25, whiten=False, outseg=outseg)

# prepare plots
if args.verbose:
Expand Down

0 comments on commit 93f8415

Please sign in to comment.