Exemplo n.º 1
0
import rf

# First get data by running this command in the console
# obspyDMT --identity 'TA.Z30A.*.BHZ' --min_mag '5.8' --min_date '2011-01-01' --max_date '2011-01-10' --event_catalog 'IRIS' --arc 'N'

# Set data path of obspyDMT here or in conf.py
rf.set_paths('~/obspyDMT-data/2011-01-01_2011-01-10_5.8_9.9')

# Convert pickled events to catalog file events.xml
rf.convert_dmteventfile()

# Select events for RF
rf.create_rfeventsfile(filters=[])

# Calculate RFs. They will be in the RF directory of the corresponding event
rf.rf('dmt', downsample=10)
Exemplo n.º 2
0
import obspy.iris
import rf

# Create getwaveform function which has to be passed to rf function
client = obspy.iris.Client()
def getwaveform(station, t1, t2):
    return client.getWaveform('TA', station, '', 'BH?', t1, t2)

# Set output path here or in config file
rf.set_paths('~/obspyDMT-data/client_test')

# Create event file from given events
rf.create_rfeventsfile('./events.xml')

#Calculate receiver functions. Station coordinates are given in stations.txt
rf.rf('client', getwaveform, './stations.txt', deconvolve='freq')