eseis (version 0.5.0)

spatial_amplitude: Locate the source of a seismic event by modelling amplutide attenuation

Description

The function fits a model of signal amplitude attenuation for all grid cells of the distance data sets and returns the residual sum as measure of the most likely source location of an event.

Usage

spatial_amplitude(data, coupling, d_map, aoi, v, q, f, a_0, normalise = TRUE)

Arguments

data

Numeric matrix or eseis object, seismic signals to work with. Since the function will calculate the maxima of the data it is usually the envolopes of the data that should be used here.

coupling

Numeric vector, coupling efficiency factors for each seismic station. The best coupled station (or the one with the highest amplification) must receive 1, the others must be scaled relatively to this one.

d_map

List object, distance maps for each station (i.e., SpatialGridDataFrame objects). Output of distance_map.

aoi

SpatialGridDataFrame or raster object that defines which pixels are used to locate the source. If omitted, the entire distance map extent is used. aoi and d_map objects must have the same extents, projections and pixel sizes. The aoi map must be of logical values.

v

Numeric value, mean velocity of seismic waves (m/s).

q

Numeric value, quality factor of the ground.

f

Numeric value, frequency for which to model the attenuation.

a_0

Logical value, start parameter of the source amplitude, if not provided, a best guess is made as 100 times the maximum amplitude value of the data set.

normalise

Logical value, option to normalise sum of residuals between 0 and 1. Default is TRUE.

Value

A raster object with the sums of squared model residuals for each grid cell.

Examples

# NOT RUN {
# }
# NOT RUN {
L <- spatial_amplitude(data = s_f, 
coupling = c(1, 1, 1, 1),
d_map = D_50$maps, 
v = 900, 
q = 100, 
f = 7.5)

# }
# NOT RUN {
# }