avgDemeanObspy = tsum / niter #print("Obspy average time for demean %e (s)"%(avgDemeanObspy)) tsum = 0 for i in range(niter): traceObspy.data = a + 1 start = time.time() traceObspy.detrend('linear') tsum = tsum + time.time() - start y = traceObspy.data avgDetrendObspy = tsum / niter tsum = 0 for i in range(niter): traceObspy.data = a + 1 start = time.time() traceObspy.filter("lowpass", corners=2, freq=0.2, zerophase=False) tsum = tsum + time.time() - start avgLPSosFilterObspy = tsum / niter tsum = 0 for i in range(niter): traceObspy.data = a + 1 start = time.time() traceObspy.filter("bandpass", corners=4, freqmin=0.1, freqmax=0.4, zerophase=True) tsum = tsum + time.time() - start avgBPSosFilterObspy = tsum / niter
# for NS components statsN = Stats() 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 = ''