Exemplos de noiseASDF em Python, exemplos de noisedbase.noiseASDF em Python

Exemplo n.º 1

0

Exibir arquivo

Arquivo: xcorr_Alaska_Love_find_bad.py Projeto: cugyyh/NoisePy

import noisedbase
import numpy as np
import timeit

#-----------------------------------------
# New xcorr database
#-----------------------------------------

dset = noisedbase.noiseASDF(
    '/work1/leon/ALASKA_work/ASDF_data/xcorr_Alaska_RTZ_lov_20190314_TA_AK.h5')
# dset.copy_stations(inasdffname = '/work1/leon/ALASKA_work/ASDF_data/xcorr_Alaska_20190218.h5', channel='LH*')

# dset.compute_xcorr(datadir = '/work2/leon/temp_working_2001_2005',
#         startdate='20010101', enddate='20051231', nprocess=10, fastfft=True)

# dset.xcorr_stack(datadir='/work2/leon/COR_ALASKA_dir', startyear=2001, startmonth=1, endyear=2019, endmonth=2, \
#                  inchannels=['LHN', 'LHE', 'LHZ'], fnametype=1)

# dset.xcorr_rotation()
# # #
# # # dset.xcorr_prephp(outdir='/work1/leon/ALASKA_work/xcorr_working_dir/pre_disp', mapfile='./MAPS_ALASKA/smpkolya_phv')
# # # # #
# dset.xcorr_aftan(prephdir='/work1/leon/ALASKA_work/xcorr_working_dir/pre_disp_L', channel='TT')
# # # # # #
# dset.interp_disp(channel='TT')
# # dset.xcorr_raytomoinput(outdir='/work1/leon/ALASKA_work/xcorr_working_dir/raytomo_input_20190131_selected_three_lambda', staxml='/home/leon/code/DataRequest/ALASKA.xml',\
# #         netcodelst=['AK', 'TA', 'PO', 'XR', 'AV', 'XN', 'XY', 'CN', 'US'])
# dset.xcorr_raytomoinput(outdir='/work1/leon/ALASKA_work/xcorr_working_dir_Love/raytomo_input_20190131_all_three_lambda', snr_thresh=10., channel='TT')
# exclude_stalst='ex_sta.lst', \
dset.xcorr_raytomoinput_debug(outdir='/work1/leon/ALASKA_work/xcorr_debug', exclude_stalst='ex_sta.lst', \
                              snr_thresh=10., channel='TT', pers=np.array([50.]))

Exemplo n.º 2

0

Exibir arquivo

Arquivo: xcorr_Alaska_Love_find_bad.py Projeto: NoisyLeon/NoisePy

import noisedbase
import numpy as np
import timeit

#-----------------------------------------
# New xcorr database
#-----------------------------------------

dset = noisedbase.noiseASDF('/work1/leon/ALASKA_work/ASDF_data/xcorr_Alaska_RTZ_lov_20190314_TA_AK.h5')
# dset.copy_stations(inasdffname = '/work1/leon/ALASKA_work/ASDF_data/xcorr_Alaska_20190218.h5', channel='LH*')

# dset.compute_xcorr(datadir = '/work2/leon/temp_working_2001_2005',
#         startdate='20010101', enddate='20051231', nprocess=10, fastfft=True)

# dset.xcorr_stack(datadir='/work2/leon/COR_ALASKA_dir', startyear=2001, startmonth=1, endyear=2019, endmonth=2, \
#                  inchannels=['LHN', 'LHE', 'LHZ'], fnametype=1)

# dset.xcorr_rotation()
# # # 
# # # dset.xcorr_prephp(outdir='/work1/leon/ALASKA_work/xcorr_working_dir/pre_disp', mapfile='./MAPS_ALASKA/smpkolya_phv')
# # # # # 
# dset.xcorr_aftan(prephdir='/work1/leon/ALASKA_work/xcorr_working_dir/pre_disp_L', channel='TT')
# # # # # # 
# dset.interp_disp(channel='TT')
# # dset.xcorr_raytomoinput(outdir='/work1/leon/ALASKA_work/xcorr_working_dir/raytomo_input_20190131_selected_three_lambda', staxml='/home/leon/code/DataRequest/ALASKA.xml',\
# #         netcodelst=['AK', 'TA', 'PO', 'XR', 'AV', 'XN', 'XY', 'CN', 'US'])
# dset.xcorr_raytomoinput(outdir='/work1/leon/ALASKA_work/xcorr_working_dir_Love/raytomo_input_20190131_all_three_lambda', snr_thresh=10., channel='TT')
# exclude_stalst='ex_sta.lst', \
dset.xcorr_raytomoinput_debug(outdir='/work1/leon/ALASKA_work/xcorr_debug', exclude_stalst='ex_sta.lst', \
                              snr_thresh=10., channel='TT', pers=np.array([50.]))
# dset.xcorr_get_field( channel='TT', snr_thresh=10., netcodelst=['AK', 'TA'])

Exemplo n.º 3

0

Exibir arquivo

Arquivo: xcorr_Alaska.py Projeto: weijias-fork/NoisePy

import noisedbase
import numpy as np
import timeit

#-----------------------------------------
# New xcorr database
#-----------------------------------------
# dset = noisedbase.noiseASDF('/work1/leon/ALASKA_work/ASDF_data/xcorr_Alaska_20180221.h5')
# dset = noisedbase.noiseASDF('/work1/leon/ALASKA_work/ASDF_data/xcorr_Alaska_20190218.h5')
dset = noisedbase.noiseASDF(
    '/work1/leon/ALASKA_work/ASDF_data/xcorr_Alaska_RTZ_ray_20190314_all_snr_10.h5'
)
#
#
# # dset.xcorr_stack(datadir='/scratch/summit/life9360/ALASKA_work/COR_work_dir', startyear=1991, startmonth=1, endyear=2018, endmonth=1)
#
# dset.xcorr_prephp(outdir='/work1/leon/ALASKA_work/xcorr_working_dir/pre_disp', mapfile='./MAPS_ALASKA/smpkolya_phv')
# # # # #
# dset.xcorr_aftan(prephdir='/work1/leon/ALASKA_work/xcorr_working_dir/pre_disp_R')
# # # # # #
# dset.interp_disp()
# # # dset.xcorr_raytomoinput(outdir='/work1/leon/ALASKA_work/xcorr_working_dir/raytomo_input_20190131_selected_three_lambda', staxml='/home/leon/code/DataRequest/ALASKA.xml',\
# # #         netcodelst=['AK', 'TA', 'PO', 'XR', 'AV', 'XN', 'XY', 'CN', 'US'])
# dset.xcorr_raytomoinput(outdir='/work1/leon/ALASKA_work/xcorr_working_dir_Rayleigh/raytomo_input_20190131_all_three_lambda_snr_10', snr_thresh=10.)
# dset.xcorr_get_field(snr_thresh=10.)
# dset.write_stationtxt('sta_alaska.lst')

#-----------------------------------------
# Appending xcorr database
#-----------------------------------------
# dset = noisedbase.noiseASDF('/work1/leon/ALASKA_work/ASDF_data/xcorr_Alaska_20190312.h5')

Exemplo n.º 4

0

Exibir arquivo

import noisedbase
import numpy as np
import timeit

#-----------------------------------------
# New xcorr database
#-----------------------------------------

dset = noisedbase.noiseASDF('/work1/leon/US_work/ASDF_data/temp_2004_2008.h5')

dset.compute_xcorr(datadir='/work2/leon/temp_working_2004_2008',
                   startdate='20060801',
                   enddate='20081231',
                   nprocess=15,
                   fastfft=True)

Exemplo n.º 5

0

Exibir arquivo

import noisedbase
import numpy as np
import timeit

dset = noisedbase.noiseASDF('xcorr_Alaska.h5')

Exemplo n.º 6

0

Exibir arquivo

Arquivo: xcorr_USCON.py Projeto: NoisyLeon/NoisePy

import noisedbase
import numpy as np
import timeit

#-----------------------------------------
# New xcorr database
#-----------------------------------------

dset = noisedbase.noiseASDF('/work1/leon/US_work/ASDF_data/xcorr_TA.h5')
dset.copy_stations(inasdffname = '/work1/leon/US_work/ASDF_data/USCON_TA.h5_cat_inv', channel='LH*')

# dset.compute_xcorr(datadir = '/work2/leon/temp_working_2001_2005',
#         startdate='20010101', enddate='20051231', nprocess=10, fastfft=True)

# dset.xcorr_stack(datadir='/work2/leon/COR_ALASKA_dir', startyear=2001, startmonth=1, endyear=2019, endmonth=2, \
#                  inchannels=['LHN', 'LHE', 'LHZ'], fnametype=1)

# dset.xcorr_rotation()
# # # 
# # # dset.xcorr_prephp(outdir='/work1/leon/ALASKA_work/xcorr_working_dir/pre_disp', mapfile='./MAPS_ALASKA/smpkolya_phv')
# # # # # 
# dset.xcorr_aftan(prephdir='/work1/leon/ALASKA_work/xcorr_working_dir/pre_disp_L', channel='TT')
# # # # # # 
# dset.interp_disp(channel='TT')
# # dset.xcorr_raytomoinput(outdir='/work1/leon/ALASKA_work/xcorr_working_dir/raytomo_input_20190131_selected_three_lambda', staxml='/home/leon/code/DataRequest/ALASKA.xml',\
# #         netcodelst=['AK', 'TA', 'PO', 'XR', 'AV', 'XN', 'XY', 'CN', 'US'])
# dset.xcorr_raytomoinput(outdir='/work1/leon/ALASKA_work/xcorr_working_dir_Love/raytomo_input_20190131_all_three_lambda', snr_thresh=10., channel='TT')
# dset.xcorr_get_field( channel='TT', snr_thresh=10., netcodelst=['AK', 'TA'])
# dset.write_stationtxt('sta_alaska.lst')

#-----------------------------------------

Exemplo n.º 7

0

Exibir arquivo

Arquivo: xcorr_USCON.py Projeto: cugyyh/NoisePy

import noisedbase
import numpy as np
import timeit

#-----------------------------------------
# New xcorr database
#-----------------------------------------

dset = noisedbase.noiseASDF('/work1/leon/US_work/ASDF_data/xcorr_TA.h5')
dset.copy_stations(
    inasdffname='/work1/leon/US_work/ASDF_data/USCON_TA.h5_cat_inv',
    channel='LH*')

# dset.compute_xcorr(datadir = '/work2/leon/temp_working_2001_2005',
#         startdate='20010101', enddate='20051231', nprocess=10, fastfft=True)

# dset.xcorr_stack(datadir='/work2/leon/COR_ALASKA_dir', startyear=2001, startmonth=1, endyear=2019, endmonth=2, \
#                  inchannels=['LHN', 'LHE', 'LHZ'], fnametype=1)

# dset.xcorr_rotation()
# # #
# # # dset.xcorr_prephp(outdir='/work1/leon/ALASKA_work/xcorr_working_dir/pre_disp', mapfile='./MAPS_ALASKA/smpkolya_phv')
# # # # #
# dset.xcorr_aftan(prephdir='/work1/leon/ALASKA_work/xcorr_working_dir/pre_disp_L', channel='TT')
# # # # # #
# dset.interp_disp(channel='TT')
# # dset.xcorr_raytomoinput(outdir='/work1/leon/ALASKA_work/xcorr_working_dir/raytomo_input_20190131_selected_three_lambda', staxml='/home/leon/code/DataRequest/ALASKA.xml',\
# #         netcodelst=['AK', 'TA', 'PO', 'XR', 'AV', 'XN', 'XY', 'CN', 'US'])
# dset.xcorr_raytomoinput(outdir='/work1/leon/ALASKA_work/xcorr_working_dir_Love/raytomo_input_20190131_all_three_lambda', snr_thresh=10., channel='TT')
# dset.xcorr_get_field( channel='TT', snr_thresh=10., netcodelst=['AK', 'TA'])
# dset.write_stationtxt('sta_alaska.lst')

Exemplo n.º 8

0

Exibir arquivo

Arquivo: temp_2004_2008.py Projeto: NoisyLeon/NoisePy

import noisedbase
import numpy as np
import timeit

#-----------------------------------------
# New xcorr database
#-----------------------------------------

dset = noisedbase.noiseASDF('/work1/leon/US_work/ASDF_data/temp_2004_2008.h5')

dset.compute_xcorr(datadir = '/work2/leon/temp_working_2004_2008',
        startdate='20060801', enddate='20081231', nprocess=15, fastfft=True)

Exemplo n.º 9

0

Exibir arquivo

Arquivo: test_xcorr.py Projeto: weijias-fork/NoisePy

import noisedbase
import numpy as np
import time

#-----------------------------------------
# New xcorr database
#-----------------------------------------

dset = noisedbase.noiseASDF(
    '/work1/leon/ALASKA_work/ASDF_data/xcorr_Alaska_20190306_Love.h5')
# dset.copy_stations(inasdffname = '/work1/leon/ALASKA_work/ASDF_data/xcorr_Alaska_20190218.h5', channel='LH*')

# xcorr_lst = dset.compute_xcorr(datadir = '/work2/leon/temp_working_2001_2005', startdate='20010101', enddate='20051231')

dset.compute_xcorr(datadir='/work2/leon/temp_working_2001_2005',
                   startdate='20010101',
                   enddate='20051231',
                   nprocess=10,
                   fastfft=True,
                   subsize=20)

# t1=time.time()
# for i in range(30):
#     xcorr_lst[i*3].convert_amph_to_xcorr(datadir='/work2/leon/temp_working_2001_2005', fastfft=True)
# #
# t2=time.time()
# for i in range(30):
#     xcorr_lst[i*3].convert_amph_to_xcorr(datadir='/work2/leon/temp_working_2001_2005', fastfft=False)
# t3=time.time()
# #
# print t3-t2, t2-t1

Exemplo n.º 10

0

Exibir arquivo

Arquivo: xcorr_Alaska.py Projeto: NoisyLeon/NoisePy

import noisedbase
import numpy as np
import timeit

#-----------------------------------------
# New xcorr database
#-----------------------------------------
# dset = noisedbase.noiseASDF('/work1/leon/ALASKA_work/ASDF_data/xcorr_Alaska_20180221.h5')
# dset = noisedbase.noiseASDF('/work1/leon/ALASKA_work/ASDF_data/xcorr_Alaska_20190218.h5')
dset = noisedbase.noiseASDF('/work1/leon/ALASKA_work/ASDF_data/xcorr_Alaska_RTZ_ray_20190314_all_snr_10.h5')
# 
# 
# # dset.xcorr_stack(datadir='/scratch/summit/life9360/ALASKA_work/COR_work_dir', startyear=1991, startmonth=1, endyear=2018, endmonth=1)
# 
# dset.xcorr_prephp(outdir='/work1/leon/ALASKA_work/xcorr_working_dir/pre_disp', mapfile='./MAPS_ALASKA/smpkolya_phv')
# # # # # 
# dset.xcorr_aftan(prephdir='/work1/leon/ALASKA_work/xcorr_working_dir/pre_disp_R')
# # # # # # 
# dset.interp_disp()
# # # dset.xcorr_raytomoinput(outdir='/work1/leon/ALASKA_work/xcorr_working_dir/raytomo_input_20190131_selected_three_lambda', staxml='/home/leon/code/DataRequest/ALASKA.xml',\
# # #         netcodelst=['AK', 'TA', 'PO', 'XR', 'AV', 'XN', 'XY', 'CN', 'US'])
# dset.xcorr_raytomoinput(outdir='/work1/leon/ALASKA_work/xcorr_working_dir_Rayleigh/raytomo_input_20190131_all_three_lambda_snr_10', snr_thresh=10.)
# dset.xcorr_get_field(snr_thresh=10.)
# dset.write_stationtxt('sta_alaska.lst')

#-----------------------------------------
# Appending xcorr database
#-----------------------------------------
# dset = noisedbase.noiseASDF('/work1/leon/ALASKA_work/ASDF_data/xcorr_Alaska_20190312.h5')
# dset.xcorr_append(inasdffname='/work1/leon/ALASKA_work/ASDF_data/xcorr_Alaska_20180221.h5', inchannels=['LHZ'], verbose=False, \
#                    datadir='/work1/leon/ALASKA_work/COR_work_dir', startyear=2018, startmonth=1, endyear=2019, endmonth=1, fnametype=2)