def rms(tr, params):
d = tr.data()
time, endtime, samprate, nsamp = tr.getv("time", "endtime", "samprate",
"nsamp")
dt = 1.0 / samprate
nsmps = int(params["twin"] * samprate)
inds = []
for i in range(int(nsamp / nsmps)):
istart = i * nsmps
iend = istart + nsmps
if iend > len(d):
break
rms = cs.rms(d[istart:iend])
if abs(rms) < params["thresh"]:
inds.append((istart, iend))
if len(inds) == 0:
return None
inds = _flatten_index_tuples(inds)
ret = []
for i in inds:
ret.append({
"meastype": params["meastype"],
"tmeas": time + dt * i[0],
"twin": dt * (i[1] - i[0]),
"auth": "auto_qc",
})
return ret
def rms(tr, params):
d = tr.data()
time, endtime, samprate, nsamp = tr.getv("time", "endtime", "samprate", "nsamp")
dt = 1.0 / samprate
nsmps = int(params["twin"] * samprate)
inds = []
for i in range(int(nsamp / nsmps)):
istart = i * nsmps
iend = istart + nsmps
if iend > len(d):
break
rms = cs.rms(d[istart:iend])
if abs(rms) < params["thresh"]:
inds.append((istart, iend))
if len(inds) == 0:
return None
inds = _flatten_index_tuples(inds)
ret = []
for i in inds:
ret.append(
{"meastype": params["meastype"], "tmeas": time + dt * i[0], "twin": dt * (i[1] - i[0]), "auth": "auto_qc"}
)
return ret
