Commit 3f014db4 authored by Sebastian Heimann's avatar Sebastian Heimann
Browse files

prepare things for tutorial

parent 64e7c25a
......@@ -7,7 +7,7 @@
First, install [Pyrocko](http://emolch.github.io/pyrocko/v0.3/install.html),
then install Lassie:
```
```bash
git clone https://gitext.gfz-potsdam.de/heimann/lassie.git
cd lassie
sudo python setup.py install
......@@ -24,7 +24,7 @@ stations from the
experiment and from 5 permanent stations from regional networks [REFERENCE] are
included. A region of roughly 100 x 100 km is covered by the available data.
### Dataset setup
#### Dataset preparation
```bash
# download data and unpack
......@@ -32,7 +32,7 @@ wget http://kinherd.org/lassie-example-alentejo.tar
tar -xf lassie-example-alentejo.tar
cd lassie-example-alentejo
ls
# data/ - raw data in mseed format
# data/ - raw waveforms in mseed format
# stations.txt - text file with station coordinates
# confirmed-events.txt - catalog with some known events
```
......@@ -45,4 +45,30 @@ application which is part of [Pyrocko](http://emolch.github.io/pyrocko/).
snuffler --stations=stations.txt --events=confirmed-events.txt data/
```
Hints, when preparing your own dataset:
* NET.STA.LOC codes in waveform files must match exactly what is given in the
stations file.
* Waveform files and directory layout do not have to follow any specific
convention, but with large files, data access can become less efficient. As a
rule of thumb, use files smaller than 1MB each. To cut a dataset into a new
chunksize, use the Jackseis program which is part of
[Pyrocko](http://emolch.github.io/pyrocko/):
```bash
# cut whole dataset to hour files:
jackseis INPUT_DIR --tinc=3600 --output-dir=OUTPUT_DIR
```
#### Lassie configuration
Lassie reads its configuration from a text file in the
[YAML](https://en.wikipedia.org/wiki/YAML) format. To create an initial
configuration, run
```bash
lassie init > config.lassie
# or filling in already the names of the stations file and the waveform directory:
lassie init --stations=stations.txt data > config.lassie
```
......@@ -153,24 +153,22 @@ data_paths:
%(data_path)s
## Processing time interval (default: use time interval of available data)
# tmin: '2016-02-05 16:00:00'
# tmax: '2016-02-05 17:40:00'
## Search grid, if not given here, a default grid will be chosen
# grid:
# !coarse.Carthesian3DGrid
# lat: 10.
# lon: 12.
# xmin: 2360.0
# xmax: 2450.0
# ymin: 7270.0
# ymax: 7360.0
# zmin: -415.0
# zmax: -400.0
# dx: 2.0
# dy: 2.0
# dz: 2.0
# tmin: '2012-02-06 04:20:00'
# tmax: '2012-02-06 04:30:00'
## Search grid; if not given here (if commented), a default grid will be chosen
# grid: !lassie.Carthesian3DGrid
# lat: 38.7
# lon: -7.9
# xmin: -70e3 # in [m]
# xmax: 70e3
# ymin: -70e3
# ymax: 70e3
# zmin: 0.0
# zmax: 0.0
# dx: 2.5e3
# dy: 2.5e3
# dz: 2.5e3
## Size factor to use when automatically choosing a grid size
autogrid_radius_factor: 1.5
......@@ -184,10 +182,10 @@ ifcs:
name: 'P'
weight: 1.0
fmin: 1.0
fmax: 50.0
fmax: 15.0
fsmooth_factor: 0.1
shifter: !lassie.VelocityShifter
velocity: 3600.
velocity: 3500.
- !lassie.WavePacketIFC
name: 'S'
weight: 1.0
......@@ -201,7 +199,7 @@ ifcs:
sharpness_normalization: true
## Threshold on detector function
detector_threshold: 200.0
detector_threshold: 100.0
## Output filename for detections
detections_path: 'detections.txt'
......
......@@ -163,7 +163,7 @@ def scan(
ydata_window = tr_stackmax.chop(wmin, wmax, inplace=False).get_ydata()
logger.info('min %g, max %g, median %g' % (
logger.info('CF stats: min %g, max %g, median %g' % (
num.min(ydata_window),
num.max(ydata_window),
num.median(ydata_window)))
......
......@@ -89,8 +89,6 @@ class WavePacketIFC(IFC):
# sumtr.ydata -= 1.0
sumtr.chop(wmin - tpad_new, wmax + tpad_new)
print sumtr.max()
dataset.append((nsl, sumtr))
return dataset
......
Markdown is supported
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment