Example #1
0
#=========================================================
dicLog = {}#Diccionario que albergará listas de logA por estación
Gamma=[]
heads={}
sfiles=np.loadtxt('sfiles.txt',dtype=str)

for evento in sfiles:
    print evento
    #Lista que albergara las distancias de las estaciones a estudiar
    logr=[]
    r=[]
    logA=[]
    amp_t=[]
    stations=[]
    
    traces,head=cut_traces(evento,cut)#diccionario que alberga trazas cortadas
    dist=distance(evento)#Diccionario que alberga la distancia de cada estacion al evento
    
    for station in traces:
        #print station
        trN=traces[station][0]
        trE=traces[station][1]
        #print trN
    
        seedrespN=respN(station,units)
        #print seedrespN
    
        seedrespE=respE(station,units)
    
        if seedrespE != {} and seedrespN != {}:
#        print station
Example #2
0
cut=1
#file_name=raw_input('06-0503-41L.S201506')

#units=raw_input('Ingrese las unidades de salida de la repuesta instrumental (VEL, DIS, ACC), si NO desea filtrar 0: ')
units='DIS'
sfile = '13-1807-32L.S201405' #Magnitud 3.7

#11-0119-41L.S201305 27-1021-16L.S201403
#31-0253-08L.S201401 13-1807-32L.S201405
#01-2250-03L.S201403 04-2306-20L.S201406
#15-0655-42L.S201403 25-0727-39L.S201406
#17-0430-31L.S201403 20-2047-32L.S201407
#17-1026-05L.S201403 24-0215-50L.S201508
#24-2319-31L.S201403 

traces,head=cut_traces(sfile,cut)#diccionario que alberga trazas cortadas

#Diccionario que alberga la distancia de cada estacion al evento
dist=distance(sfile)

#Lista que albergara las distancias de las estaciones a estudiar
logr=[]
logA=[]
#Diccionario que alberga los logaritmos de las amplitudes para encontrar diferencia entre cut_traces4 y 5

#30 35
#=========================================================
pre_filt = (0.02, 0.03, 20.0, 25.0)
<<<<<<< HEAD
=======
<<<<<<< HEAD
Example #3
0
from cut_traces3 import cut_traces
from resp import respN
from obspy.signal.invsim import seisSim,simulate_seismometer

units = 'DIS'
station = 'GUA'
pre_filt = (0.2, 0.3, 20.0, 25.0)

dic,head=cut_traces('01-2250-03L.S201403',0)
tr=dic[station][0]
tr2 = tr.copy()
tr1 = tr.copy()

seedrespN=respN(station,units)
tr1.simulate(paz_remove=None, pre_filt=pre_filt, seedresp=seedrespN)

seedrespN['date']=tr2.stats.starttime
seedrespN['network']=tr2.stats.network
seedrespN['station']=tr2.stats.station
seedrespN['location']=tr2.stats.location
seedrespN['network']=tr2.stats.network
seedrespN['channel']=tr2.stats.channel
tr2.data=simulate_seismometer(tr2.data, 100, remove_sensitivity=False, simulate_sensitivity=False, 
                 water_level=0, seedresp=seedrespN, pre_filt=None )

tr.plot()
#tr1.plot()
tr2.plot()