Skip to content

Commit

Permalink
Browse files Browse the repository at this point in the history
update
  • Loading branch information
trichter committed Feb 14, 2013
1 parent a6439cd commit 257529f
Show file tree
Hide file tree
Showing 17 changed files with 432 additions and 642 deletions.
2 changes: 1 addition & 1 deletion bin/console/plotday.py
Expand Up @@ -56,7 +56,7 @@
kwargs['scale'] = 1.
else:
station = args.file_station
if station.startswith('PB') or station == 'LVC':
if station.startswith('PB') or station == 'LVC' or station.endswith('CX'):
from sito.data import IPOC
data = IPOC(xcorr_append=args.xcorr_append)
elif station == 'PKD':
Expand Down
2 changes: 2 additions & 0 deletions bin/dayplots_overview.py
Expand Up @@ -14,12 +14,14 @@
# 13 stations, 7 days, 2007, 2008
# days = '2006-12-01 2007-02-01 2007-04-01 2007-06-01 2007-08-01 2007-10-01 2007-12-01'
days = '2008-01-01 2008-02-01 2008-04-01 2008-06-01 2008-08-01 2008-10-01 2008-12-01'
days = '%d-01-01 %d-02-01 %d-04-01 %d-06-01 %d-08-01 %d-10-01 %d-12-01' % ((2012,) * 7)
stations = 'PB01 PB02 PB03 PB04 PB05 PB06 PB07 PB08 HMBCX MNMCX PATCX PSGCX LVC'
# 20 stations, 10 days, 2008-2012
# days = '2008-01-01 2008-07-01 2009-01-01 2009-07-01 2010-01-01 2010-07-01 2011-01-01 2011-07-01 2012-01-01 2012-07-01'
# stations = 'PB01 PB02 PB03 PB04 PB05 PB06 PB07 PB08 PB09 PB10 PB11 PB12 PB13 PB14 PB15 HMBCX MNMCX PATCX PSGCX LVC'
pic_dir = '/home/richter/Data/IPOC/raw_pictures/'
output = '/home/richter/Results/IPOC/raw_pics/test.pdf'
output = '/home/richter/Results/IPOC/raw_pics/13stations_7days_2012.pdf'
scale_range = 2000
scale_LVC = 20000

Expand Down
31 changes: 22 additions & 9 deletions bin/maps_ipoc.py
Expand Up @@ -19,7 +19,7 @@ def plot_some_events():
2008-02-15T16:54:04 5.5 -23.00 -70.20 32
2008-02-04T17:01:30 6.6 -20.20 -70.00 36
2007-12-16T08:09:16 7.1 -22.80 -70.00 14
2007-11-14T15:40:51 7.8 -22.34 -70.06 37""" #GEOFON:-22.30 -69.80
2007-11-14T15:40:51 7.8 -22.34 -70.06 37""" #GEOFON:-22.30 -69.80
events = []
for eq in eqs.split('\n'):
time, mag, lat, lon, depth = eq.split()
Expand Down Expand Up @@ -55,9 +55,9 @@ def plot_results1():
ipoc_no_data = ipoc.pick('PB09 PB10 PB11 PB12 PB13 PB14 PB15', replace=False)
#1-3Hz
ipoc_bad_data = ipoc.pick('LVC MNMCX ', replace=False)
ipoc_no_signal = ipoc.pick('PB06 PB08 PSGCX', replace=False) #<0.2%
ipoc_small_signal = ipoc.pick('HMBCX PATCX PB01 PB02 PB03 PB07', replace=False) #0.2-0.5%
ipoc_big_signal = ipoc.pick('PB04 PB05', replace=False) #>0.5%
ipoc_no_signal = ipoc.pick('PB06 PB08 PSGCX', replace=False) #<0.2%
ipoc_small_signal = ipoc.pick('HMBCX PATCX PB01 PB02 PB03 PB07', replace=False) #0.2-0.5%
ipoc_big_signal = ipoc.pick('PB04 PB05', replace=False) #>0.5%

ipoc_no_data.plot(m, mfc='w', ms=8, zorder=20)
ipoc_bad_data.plot(m, mfc='#00FFFF', ms=8, zorder=20)
Expand All @@ -69,9 +69,9 @@ def plot_results1():
for st in ipoc:
ipoc[st].longitude += 0.15
ipoc_bad_data = ipoc.pick('PB04', replace=False)
ipoc_no_signal = ipoc.pick('HMBCX MNMCX PB01 PB06 PSGCX LVC', replace=False) #<0.2%
ipoc_small_signal = ipoc.pick('PB02 PB03 PB05 PB07 PB08', replace=False) #0.2-0.5%
ipoc_big_signal = ipoc.pick('PATCX', replace=False) #>0.5%
ipoc_no_signal = ipoc.pick('HMBCX MNMCX PB01 PB06 PSGCX LVC', replace=False) #<0.2%
ipoc_small_signal = ipoc.pick('PB02 PB03 PB05 PB07 PB08', replace=False) #0.2-0.5%
ipoc_big_signal = ipoc.pick('PATCX', replace=False) #>0.5%

ipoc_bad_data.plot(m, marker='s', mfc='#00FFFF', ms=8, zorder=20, annotate=False)
ipoc_no_signal.plot(m, marker='s', mfc='y', ms=8, zorder=20, annotate=False)
Expand All @@ -85,5 +85,18 @@ def plot_results1():
labels, numpoints=1, loc='lower right')
plt.gcf().savefig('/home/richter/Documents/pics/maps/ipoc/ipoc_map_results.png')

plot_results1()
#plt.show()
def plot_small_map_for_F():
from sito.stations import IPOCStations
from matplotlib.lines import Line2D

ipoc = IPOCStations()
ipoc.pick('PATCX')
m = map.createIPOCMap(show=False, ll=(-24, -70.9), ur=(-18, -68.5),
earthquake='Tocopilla_position', stations=ipoc,
cities='Tocopilla')



#plot_results1()
plot_small_map_for_F()
plt.show()
43 changes: 0 additions & 43 deletions bin/noise/convert2sac.py

This file was deleted.

29 changes: 15 additions & 14 deletions bin/noise/noise_s_final_autocorr2.py
Expand Up @@ -13,13 +13,14 @@

def main():
stations = 'PB01 PB02 PB03 PB04 PB05 PB06 PB07 PB08 HMBCX MNMCX PATCX PSGCX LVC'
stations = 'PB09 PB10 PB11 PB12 PB13 PB14 PB15 PB16'
stations2 = None


components = 'Z'
# TOcopilla earthquake: 2007-11-14 15:14
t1 = UTC('2006-02-01')
t2 = UTC('2008-12-31')
t2 = UTC('2012-10-01')

shift = 100
correlations = get_correlations(stations, components, stations2, only_auto=True)
Expand All @@ -29,16 +30,16 @@ def main():
data = IPOC(xcorr_append='/' + method, use_local_LVC=False)
data.setXLogger('_' + method)

# pool = Pool()
# prepare(data, stations.split(), t1, t2, component=components,
# filter=(1, 3, 2, True), downsample=20,
# eventremoval='waterlevel_env2', param_removal=(10, 0),
# whitening=False,
# normalize='1bit', param_norm=None,
# pool=pool)
# noisexcorrf(data, correlations, t1, t2, shift, pool=pool)
# pool.close()
# pool.join()
pool = Pool()
prepare(data, stations.split(), t1, t2, component=components,
filter=(1, 3, 2, True), downsample=20,
eventremoval='waterlevel_env2', param_removal=(10, 0),
whitening=False,
normalize='1bit', param_norm=None,
pool=pool)
noisexcorrf(data, correlations, t1, t2, shift, pool=pool)
pool.close()
pool.join()

# plotXcorrs(data, correlations, t1, t2, start=None, end=None, plot_overview=True, plot_years=False, use_dlognorm=False,
# plot_stack=True, plot_psd=False, add_to_title='', downsample=None)
Expand All @@ -49,9 +50,9 @@ def main():
# plotXcorrs(data, correlations, t1=None, t2=None, start=None, end=None, plot_overview=True, plot_years=False, use_dlognorm=False,
# plot_stack=True, plot_psd=False, add_to_title='', downsample=None,
# stack=('10days', '2days'))
plotXcorrs(data, correlations, t1=None, t2=None, start=0, end=20, plot_overview=True, plot_years=False, use_dlognorm=False,
plot_stack=True, plot_psd=False, add_to_title='', downsample=None,
stack=('10days', '2days'), ext='_hg.png', vmax=0.1)
# plotXcorrs(data, correlations, t1=None, t2=None, start=0, end=20, plot_overview=True, plot_years=False, use_dlognorm=False,
# plot_stack=True, plot_psd=False, add_to_title='', downsample=None,
# stack=('10days', '2days'), ext='_hg.png', vmax=0.1)
# util.checkDir(data.getPlotX(('', ''), t1))
#for correlation in correlations:
# stations = correlation[0][:-1], correlation[1][:-1]
Expand Down
41 changes: 29 additions & 12 deletions bin/noise/noise_s_final_autocorr3.py
Expand Up @@ -13,45 +13,62 @@

def main():
stations = 'PB01 PB02 PB03 PB04 PB05 PB06 PB07 PB08 HMBCX MNMCX PATCX PSGCX LVC'
stations = 'PB09 PB10 PB11 PB12 PB13 PB14 PB15 PB16'
stations = 'PATCX'
stations2 = None


components = 'Z'
# TOcopilla earthquake: 2007-11-14 15:14
t1 = UTC('2006-02-01')
t2 = UTC('2008-12-31')
#t2 = UTC('2008-12-31')
#t2 = UTC('2012-10-01')
t2 = UTC('2011-12-31')
# t1 = UTC('2009-05-01')
# t2 = UTC('2009-05-03')

shift = 100
correlations = get_correlations(stations, components, stations2, only_auto=True)
shift = 60
#correlations = get_correlations(stations, components, stations2, only_auto=True)
correlations = get_correlations(stations, components, stations2)
print correlations

method = 'FINAL_filter4-6_1bit_auto'
method = 'FINAL_filter4-6_1bit_auto_3C'
method = 'FINAL_filter3-5'

data = IPOC(xcorr_append='/' + method, use_local_LVC=False)
data.setXLogger('_' + method)

# pool = Pool()
pool = Pool()
# prepare(data, stations.split(), t1, t2, component=components,
# filter=(4, 6, 2, True), downsample=20,
# filter=(2, 4, 2, True), downsample=50,
# eventremoval='waterlevel_env2', param_removal=(10, 0),
# whitening=False,
# normalize='1bit', param_norm=None,
# pool=pool)
# noisexcorrf(data, correlations, t1, t2, shift, pool=pool)
# pool.close()
# pool.join()
# noisexcorrf(data, correlations, t1, t2, shift, period=24 * 3600, pool=pool)

# noisexcorrf(data, correlations, t1, t2, shift, period=5 * 60, pool=pool,
# max_preload=1000)
pool.close()
pool.join()

# plotXcorrs(data, correlations, t1, t2, start=None, end=None, plot_overview=True, plot_years=False, use_dlognorm=False,
# plot_stack=True, plot_psd=False, add_to_title='', downsample=None)
plotXcorrs(data, correlations, t1, t2, start=0, end=20, plot_overview=True, plot_years=False, use_dlognorm=False,

plotXcorrs(data, correlations, t1, t2, start= -20, end=20, plot_overview=True, plot_years=False, use_dlognorm=False,
plot_stack=True, plot_psd=False, add_to_title='', downsample=None, ext='_hg.png', vmax=0.1)

# stack(data, correlations, dt= -1)
# stack(data, correlations, dt=10 * 24 * 3600, shift=2 * 24 * 3600)
#stack(data, correlations, dt=60 * 60, period=5 * 60)
# stack(data, correlations, dt=24 * 60 * 60, period=5 * 60)
# plotXcorrs(data, correlations, t1=None, t2=None, start=None, end=None, plot_overview=True, plot_years=False, use_dlognorm=False,
# plot_stack=True, plot_psd=False, add_to_title='', downsample=None,
# stack=('10days', '2days'))
plotXcorrs(data, correlations, t1, t2, start=0, end=20, plot_overview=True, plot_years=False, use_dlognorm=False,
plot_stack=True, plot_psd=False, add_to_title='', downsample=None,
stack=('10days', '2days'), ext='_hg.png', vmax=0.1)
# plotXcorrs(data, correlations, t1, t2, start=0, end=20, plot_overview=True, plot_years=False, use_dlognorm=False,
# plot_stack=True, plot_psd=False, add_to_title='', downsample=None,
# period=60 * 5, stack=(60 * 60, None), ext='_hg.png', vmax=0.1)

# util.checkDir(data.getPlotX(('', ''), t1))
#for correlation in correlations:
Expand Down

0 comments on commit 257529f

Please sign in to comment.