import obspy
from obspy.core.trace import Trace
import numpy as np

# Creating trace in obspy
npts = 10000
trace_data = np.random.rand(npts)

or_time = obspy.UTCDateTime(2019,10,15,10,0,0) # Origin Time
dt = 1/20.

Trace = Trace(trace_data, header={"starttime": or_time, 'delta': dt})

start = obspy.UTCDateTime(2019,10,15,10,1,0)
end = obspy.UTCDateTime(2019,10,15,10,1,40)

Trace.trim(starttime=start, endtime=end, pad=True, nearest_sample=True, fill_value=0.0)

a=1
Beispiel #2
0
        statsN.sampling_rate = 1.0 / sampling_rate_x
        statsN.delta = sampling_rate_x
        statsN.starttime = starttime
        statsN.npts = len(traceN.data)
        statsN.network = net
        statsN.station = stationname
        statsN.location = ''
        statsN.channel = channelnameN
        traceN.stats = statsN
        traceN.stats.sac = obspy.core.AttribDict()
        traceN.stats.sac.back_azimuth = plottheta  # use this as azimuth of station

        #---applying filters---#
        traceN.filter('bandpass', freqmin=freqmin, freqmax=freqmax)
        tN = traceN.stats.starttime
        traceN.trim(starttime=tN, endtime=tN + trim_end_time)
        traceN.taper(0.05, side='right')

        #----------------------#
        stemp += Stream(traceN)
        # for EW components
        statsE = Stats()
        statsE.sampling_rate = 1.0 / sampling_rate_z
        statsE.delta = sampling_rate_z
        statsE.starttime = starttime
        statsE.npts = len(traceN.data)
        statsE.network = net
        statsE.station = stationname
        statsE.location = ''
        statsE.channel = channelnameE
        traceE.stats = statsE