obspy.core.stream.Stream.print_gaps

Stream.print_gaps(min_gap=None, max_gap=None)[source]

Print gap/overlap list summary information of the Stream object.

Parameters:
  • min_gap – All gaps smaller than this value will be omitted. The value is assumed to be in seconds. Defaults to None.

  • max_gap – All gaps larger than this value will be omitted. The value is assumed to be in seconds. Defaults to None.

Example

Our example stream has no gaps:

>>> from obspy import read, UTCDateTime
>>> st = read()
>>> st.get_gaps()
[]
>>> st.print_gaps()  
Source            Last Sample                 Next Sample ...
Total: 0 gap(s) and 0 overlap(s)

So let’s make a copy of the first trace and cut both so that we end up with a gappy stream:

>>> tr = st[0].copy()
>>> t = UTCDateTime("2009-08-24T00:20:13.0")
>>> st[0].trim(endtime=t)  
<...Trace object at 0x...>
>>> tr.trim(starttime=t+1)  
<...Trace object at 0x...>
>>> st.append(tr)  
<...Stream object at 0x...>
>>> st.get_gaps()  
[[..., UTCDateTime(2009, 8, 24, 0, 20, 13), ...
>>> st.print_gaps()  
Source            Last Sample                 ...
BW.RJOB..EHZ      2009-08-24T00:20:13.000000Z ...
Total: 1 gap(s) and 0 overlap(s)

And finally let us create some overlapping traces:

>>> st = read()
>>> tr = st[0].copy()
>>> t = UTCDateTime("2009-08-24T00:20:13.0")
>>> st[0].trim(endtime=t)  
<...Trace object at 0x...>
>>> tr.trim(starttime=t-1)  
<...Trace object at 0x...>
>>> st.append(tr)  
<...Stream object at 0x...>
>>> st.get_gaps()  
[[...'EHZ', UTCDateTime(2009, 8, 24, 0, 20, 13), ...
>>> st.print_gaps()  
Source            Last Sample                 ...
BW.RJOB..EHZ      2009-08-24T00:20:13.000000Z ...
Total: 0 gap(s) and 1 overlap(s)