Seis.jl

Seis.jl is an open, fast and flexible framework for analysing seismic data in Julia.

For detailed instructions on using Seis, see the documentation via the links below.

This readme file gives a brief introduction to installing and using Seis.

Installation

At present, Seis is unregistered and it and its dependencies must be installed manually. This will change once Seis is registered in the General registry.

To install Seis, do:

julia> ] # As soon as you press ']', you enter Pkg mode...

(v1.3) pkg> add https://github.com/anowacki/Geodesics.jl https://github.com/anowacki/Seis.jl

Basic introduction

Seis is based around single traces of evenly-sampled continuous data where time is the independent variable. (In future we may support unevenly-sampled data and traces with gaps, and we plan to soon support frequency-domain traces.) Each trace holds its sampling interval (the property .delta) and a start time (.b). Traces have the type Trace.

Traces also optionally hold information about the event (.evt) and station (.sta) associated with this recording. Events and Stations are the corresponding types, which also contain useful properties.

For more information, see the docstrings for Trace, Event and Station. This is easily done in the Julia REPL like so:

julia> using Seis

julia> ? # As soon as you press '?', the prompt changes to 'help?>'

help?> Trace
search: Trace trace backtrace AbstractSet AbstractVector AbstractVecOrMat

  Trace

  Evenly-sampled time series recorded at a single seismic station. The start
  time of the trace, in s, is in the b field, whilst the sampling interval, in
  s, is delta. The trace itself is accessed using the trace method, like
  trace(t).

  [...]

All three types above also hold a .meta property, which contains any extra metadata you want to associate with the trace, event or station.

Common conventions

f vs f!

Traces are mutable structs, and therefore for all functions which potentially modify a trace, there are two versions. Firstly, an in-place function (e.g., bandpass!) which as per Julia convention ends with an exclamation mark, and modifies the trace. Secondly, for convenience, there is always a copying version without the exclamation mark (e.g., bandpass) which returns a modified copy of the input trace.

Arrays of traces

There is no special type for holding multiple traces. Instead, we operate on arrays of traces. For instance, reading multiple traces from the same event, we can access all the station names like so:

julia> t = sample_data(:array);

julia> t.sta.sta
60-element Array{String,1}:
 "ABA"
 "APA"
 "AWI"
 "BBH"
 "BBO"
 "BDL"
 "BTA"
 "BWH"
 "CRA"
 "CSF"
 "EAB"
 "EAU"
 "EBH"
 ⋮    
 "PMS"
 "TSA"
 "WAL"
 "WCB"
 "WME"
 "WPM"
 "XAL"
 "XDE"
 "YEL"
 "YLL"
 "YRC"
 "YRE"

Or to get the full channel code from available header information:

julia> channel_code.(t)
60-element Array{String,1}:
 ".ABA..SHZ"
 ".APA..SHZ"
 ".AWI..SHZ"
 ".BBH..SHZ"
 ".BBO..SHZ"
 ".BDL..SHZ"
 ".BTA..SHZ"
 ".BWH..SHZ"
 ".CRA..SHZ"
 ".CSF..SHZ"
 ".EAB..SHZ"
 ".EAU..SHZ"
 ".EBH..SHZ"
 ⋮          
 ".PMS..SHZ"
 ".TSA..SHZ"
 ".WAL..SHZ"
 ".WCB..SHZ"
 ".WME..SHZ"
 ".WPM..SHZ"
 ".XAL..SHZ"
 ".XDE..SHZ"
 ".YEL..SHZ"
 ".YLL..SHZ"
 ".YRC..SHZ"
 ".YRE..SHZ"

Note the use of the broadcasted . operation (channel_code.(t)) which applied the 'scalar' function (channel_code) to each trace in the array t.

IO

Currently, SAC and miniSEED data are read and written. SAC files may be either bigendian (SAC/BRIS convention) or little-endian (usual IRIS SAC convention)).

Use the read_sac and write_sac functions for SAC-formatted IO, and read_mseed and write_mseed for miniSEED files.

Future work will add support for reading many more formats and format auto-detection.

Plotting

Makie.jl

Seis.jl enables the use of the Makie plotting package for plotting seismic data. To enable the functionaliy, load a Makie backend package, such as CairoMakie or GLMakie. Then plotting functions such as plot, plot_section and plot_hodogram can be used.

As an example:

julia> using Seis, CairoMakie

julia> t = sample_data(:array);

julia> plot_section(t)

produces:

See the online documentation for more information.

Plots.jl

Plotting using Plots.jl is supported as a legacy option. You can enable this by doing import Pkg; Pkg.add("Plots"). You then need to do using Plots when you want to start using Seis.jl's plotting routines.

Basic processing

Basic time series processing of traces is possible using functions such as integrate, bandpass, remove_trend. See the online documentation for a full description of the functions available and how to use them.

Maintainers

At present, Andy Nowacki (@anowacki) is the primary maintainer of Seis.jl.

Contributing

Contributions to Seis are most welcome. Please open issues in this repo with bug reports or feature requests. Pull requests for the same are also very welcome.