関数semblance_given_slowness_azimuth マニュアル

(The documentation of function semblance_given_slowness_azimuth)
Last Update: 2025/8/7


◆機能・用途(Purpose)

指定されたスローネスと波の到来方向を用いて時系列データのセンブランスを計算する。
Compute the semblance of time-series data using a given slowness and back azimuth.

◆形式(Format)

#include <array/stack.h>
inline struct sequence semblance_given_slowness_azimuth
(const struct array_station_info stations,
 const struct sequence ∗data_before, const struct sequence ∗data,
 const struct sequence ∗data_after,
 const double slowness,const double back_azimuth,
 const double window_length,const double tinc)


◆引数(Arguments)

stations アレイの観測点情報。
The station information of the array.
data_before センブランス計算の対象とするタイムウインドウの直前の時系列データ を各観測点について並べた配列。
An array composed of the time-series data at every station immediately before the time window for the semblance computation.
data センブランス計算の対象とするタイムウインドウの時系列データ を各観測点について並べた配列。
An array composed of the time-series data at every station for the time window for the semblance computation.
data_after センブランス計算の対象とするタイムウインドウの直後の時系列データ を各観測点について並べた配列。
An array composed of the time-series data at every station immediately after the time window for the semblance computation.
slowness センブランス計算の際に仮定する波のスローネス。 時刻の単位は時系列データに、距離の単位は観測点座標に合わせる。
The slowness of the wave assumed in the computation of the semblance. The unit of the time must be same as that of the time-series data, and the unit of the distance must be same as that of the station coordinates.
back_azimuth センブランス計算の際に仮定する波の到来方向(°)。 東から反時計回りに測る。
The back azimuth (°) of the wave assumed in the computation of the semblance, measured counterclockwise from the east.
window_length 一つのセンブランス値の計算に用いるタイムウインドウ長。
Length of the time window used for the computation of each semblance value.
tinc センブランス値のサンプリング間隔。 dataのサンプリング間隔の整数倍かつ dataのウインドウ長の整数分の1でなければならない。
Sampling interval of the semblance values; this must be a multiple of the sampling interval of the data, and the window length of the data must be dividable by this value.


◆戻り値(Return value)

センブランス値を並べた時系列データ。 センブランス値は引数dataの先頭から末尾までの時刻範囲で、 引数tincの時間刻み幅で計算する。 個々のセンブランス値の計算には その時刻を中心とした±window_length/2のタイムウインドウを用いる。 センブランス値の計算式は \[\begin{equation} S=\frac{\int dt s_1(t)^2}{N\int dt s_2(t)} \label{eq.S} \end{equation}\] を用いる。ここで\(N\)は観測点数、 \(s_1(t)\)は指定されたスローネスと波の到来方向を用いてスタックした時系列データ、 \(s_2(t)\)は指定されたスローネスと波の到来方向を用いて2乗スタックした時系列データである。
A time series data composed of the semblance values, covering the time window of argument data with time steps given by argument tinc. Each semblance value is computed using a time window of window_length/2 centered on the time of the data sample of the semblance. Eq. (\ref{eq.S}) is used as the definition of the semblance, where \(N\) is the number of stations, and \(s_1(t)\) and \(s_2(t)\) are the stacked and square-stacked time-series data, respectively, using the given slowness and back azimuth.

◆使用例(Example)

struct sequence semblance=semblance_given_slowness_azimuth
     (stations,data_before,data,data_after,0.001,30.0,5.0,2.0);