Skip to content
GitLab
Menu
Projects
Groups
Snippets
Help
Help
Support
Community forum
Keyboard shortcuts
?
Submit feedback
Contribute to GitLab
Sign in
Toggle navigation
Menu
Open sidebar
Sebastian Heimann
lassie
Commits
2f88f082
Commit
2f88f082
authored
Jun 10, 2016
by
Sebastian Heimann
Browse files
rename config entry ifcs to image_function_contributions
parent
f53bb807
Changes
4
Hide whitespace changes
Inline
Side-by-side
README.md
View file @
2f88f082
...
...
@@ -129,8 +129,8 @@ autogrid_radius_factor: 1.5
## Grid density factor used when automatically choosing a grid
autogrid_density_factor
:
10.0
##
Image function contributions
i
fc
s
:
##
Composition of image function
i
mage_function_contribution
s
:
-
!lassie.WavePacketIFC
name
:
'
P'
weight
:
1.0
...
...
@@ -170,10 +170,10 @@ coda wave packets travelling with characteristic speeds, packet durations and
frequency content through the network. To detect such events automatically, we
can compose an image function with two contributions, one for each dominant
phase. Such image function contributions (IFCs) can be defined under the
`image_function_contributions`
key in Lassie's configuration
[TODO: to be
renamed]. Different types
of IFCs are available and can be combined. Here, we
use the
`WavePacketIFC`
which applies a normalized envelope based
pre-processing to the waveforms
before travel-time compensated stacking.
`image_function_contributions`
key in Lassie's configuration
. Different types
of IFCs are available and can be combined. Here, we
use the
`WavePacketIFC`
which applies a normalized envelope based
pre-processing to the waveforms
before travel-time compensated stacking.
#### Running the detector
...
...
@@ -181,8 +181,7 @@ To test the configuration, select a short processing time span around one of
the known events in the dataset. The time span can be restricted with the
`tmin`
and
`tmax`
entries in Lassie's configuration or with corresponding
command line options
`--tmin='YYYY-MM-DD HH:MM:SS'`
and
`--tmax='YYYY-MM-DD
HH:MM:SS.XXX'`
[TODO: to be implemented], which override configuration
settings.
HH:MM:SS.XXX'`
, which override configuration settings.
Now run the detector with
...
...
apps/lassie
View file @
2f88f082
...
...
@@ -193,8 +193,8 @@ autogrid_radius_factor: 1.5
## Grid density factor used when automatically choosing a grid
autogrid_density_factor: 10.0
## Image function co
ntribu
tion
s
i
fc
s:
## Image function co
mposi
tion
i
mage_function_contribution
s:
- !lassie.WavePacketIFC
name: 'P'
weight: 1.0
...
...
src/config.py
View file @
2f88f082
...
...
@@ -54,7 +54,7 @@ class Config(Object):
help
=
'grid density factor used when automatically choosing a grid '
'spacing'
)
i
fc
s
=
List
.
T
(
i
mage_function_contribution
s
=
List
.
T
(
ifc
.
IFC
.
T
(),
help
=
'image function contributions'
)
...
...
@@ -95,9 +95,11 @@ class Config(Object):
receivers
=
self
.
get_receivers
()
fsmooth_max
=
max
(
ifc
.
get_fsmooth
()
for
ifc
in
self
.
ifcs
)
ifc
.
get_fsmooth
()
for
ifc
in
self
.
image_function_contributions
)
vmin
=
min
(
ifc
.
shifter
.
get_vmin
()
for
ifc
in
self
.
ifcs
)
vmin
=
min
(
ifc
.
shifter
.
get_vmin
()
for
ifc
in
self
.
image_function_contributions
)
spacing
=
vmin
/
fsmooth_max
/
self
.
autogrid_density_factor
...
...
src/core.py
View file @
2f88f082
...
...
@@ -37,9 +37,11 @@ def scan(
norm_map
=
gridmod
.
geometrical_normalization
(
grid
,
receivers
)
ifcs
=
config
.
image_function_contributions
shift_tables
=
[]
tshift_maxs
=
[]
for
ifc
in
config
.
ifcs
:
for
ifc
in
ifcs
:
shift_tables
.
append
(
ifc
.
shifter
.
get_table
(
grid
,
receivers
))
tshift_maxs
.
append
(
shift_tables
[
-
1
].
max
())
...
...
@@ -47,9 +49,9 @@ def scan(
tinc
=
tshift_max
*
10
tpad
=
max
(
ifc
.
get_tpad
()
for
ifc
in
config
.
ifcs
)
+
tshift_max
tpad
=
max
(
ifc
.
get_tpad
()
for
ifc
in
ifcs
)
+
tshift_max
fsmooth_min
=
min
(
ifc
.
get_fsmooth
()
for
ifc
in
config
.
ifcs
)
fsmooth_min
=
min
(
ifc
.
get_fsmooth
()
for
ifc
in
ifcs
)
blacklist
=
set
(
tuple
(
s
.
split
(
'.'
))
for
s
in
config
.
blacklist
)
...
...
@@ -93,7 +95,7 @@ def scan(
pdata
=
[]
trs_debug
=
[]
shift_maxs
=
[]
for
iifc
,
ifc
in
enumerate
(
config
.
ifcs
):
for
iifc
,
ifc
in
enumerate
(
ifcs
):
dataset
=
ifc
.
preprocess
(
trs
,
wmin
,
wmax
,
tshift_max
)
if
not
dataset
:
continue
...
...
@@ -207,8 +209,8 @@ def scan(
f
.
close
()
if
show_detections
:
fmin
=
min
(
ifc
.
fmin
for
ifc
in
config
.
ifcs
)
fmax
=
min
(
ifc
.
fmax
for
ifc
in
config
.
ifcs
)
fmin
=
min
(
ifc
.
fmin
for
ifc
in
ifcs
)
fmax
=
min
(
ifc
.
fmax
for
ifc
in
ifcs
)
plot
.
plot_detection
(
grid
,
receivers
,
frames
,
tmin_frames
,
deltat_cf
,
imax
,
iframe
,
fsmooth_min
,
xpeak
,
ypeak
,
zpeak
,
...
...
Write
Preview
Supports
Markdown
0%
Try again
or
attach a new file
.
Attach a file
Cancel
You are about to add
0
people
to the discussion. Proceed with caution.
Finish editing this message first!
Cancel
Please
register
or
sign in
to comment