sacfile_decimateコマンド マニュアル

(The documentation of sacfile_decimate command)
Last Update: 2024/8/7

◆機能・用途(Purpose)

SAC時系列データのサンプル間隔を間引く。
Resample the samples of a SAC time series data.

◆ソースコード(Source code)

$YMAEDA_OPENTOOL_DIR/sac_data/src/sacfile_decimate.c


◆使用方法(Usage)

コマンドライン引数でパラメータを指定する。 パラメータの一覧を下表に示す。
Specify parameters by command-line arguments. The table below shows a list of parameters.

●「-」から始まらない引数 (Arguments not beginning with “-”)

引数
Argument 		与える値
Quantity to be given
第1引数
1st argument 		入力SACデータファイル名。
The input SAC data file name.
第2引数
2nd argument 		出力する間引き後のSACデータファイル名。
Name of the output SAC data file into which the resampled time series data will be written.


●1つの「-」から始まる引数 (Arguments beginning with a single “-”)

このコマンドでは1つの「-」から始まる引数は存在しない。
This command does not have arguments beginning with a single “-”.

●「--パラメータ名=パラメータ値」の形式の引数 (Arguments of a form “--Parameter name=Parameter Value”)

「--パラメータ名=パラメータ値」の形式の引数は自由な順番で指定できる。 「-」から始まらない引数の間に挿入しても良い。 相反する指定がなされた場合には後の指定が優先される。 デフォルト値を持つパラメータは省略できる。
Arguments of a form “--Parameter name=Parameter Value” can be placed in an arbitrary order. They can even be inserted between arguments not beginning with “-”. In case of conflicting options being specified, the latter option has a higher priority. Parameters that have default values can be omitted.
パラメータ名
Parameter name 		意味
Meaning 		可能なパラメータ値
Allowed parameter values 		デフォルト値
Default value
new_dt 間引き後の新しいサンプリング間隔\(\Delta t_{new}\)(s)。
The new sampling interval \(\Delta t_{new}\) (s) after the resampling. 		入力SACデータにおけるサンプリング間隔の整数倍となる正の実数。
A positive real number that is an integer multiple of the sampling interval of the input SAC data. 		省略不可
Cannot be omitted
mode リサンプル後の時系列データに使用する値。
Values used in the resampled time series data.
  • normal
    残すサンプル時刻での元々の時系列データの値。
    The value of the original time series data at each remained sample time.

  • ave
    区間\([t-\Delta t_{new}/2,t+\Delta t_{new}/2]\) での平均値を各サンプル時刻\(t\)での値とする。
    Use the average in section \([t-\Delta t_{new}/2,t+\Delta t_{new}/2]\) as the value at each sample time \(t\).

  • abs_ave
    区間\([t-\Delta t_{new}/2,t+\Delta t_{new}/2]\) での絶対値平均値を各サンプル時刻\(t\)での値とする。
    Use the absolute value average in section \([t-\Delta t_{new}/2,t+\Delta t_{new}/2]\) as the value at each sample time \(t\).

  • max
    区間\([t-\Delta t_{new}/2,t+\Delta t_{new}/2]\) での最大値を各サンプル時刻\(t\)での値とする。
    Use the maximum value in section \([t-\Delta t_{new}/2,t+\Delta t_{new}/2]\) as the value at each sample time \(t\).

  • min
    区間\([t-\Delta t_{new}/2,t+\Delta t_{new}/2]\) での最小値を各サンプル時刻\(t\)での値とする。
    Use the minimum value in section \([t-\Delta t_{new}/2,t+\Delta t_{new}/2]\) as the value at each sample time \(t\).

  • abs_max
    区間\([t-\Delta t_{new}/2,t+\Delta t_{new}/2]\) での絶対値最大値を各サンプル時刻\(t\)での値とする。
    Use the absolute maximum value in section \([t-\Delta t_{new}/2,t+\Delta t_{new}/2]\) as the value at each sample time \(t\).

 normal


◆動作(Behaviour)

第1引数で指定したSACファイルを読み込み、 サンプリング間隔がパラメータnew_dtで指定した値になるように間引いた上で 第2引数で指定したSACファイルに出力する。 先頭時刻は変えず、フィルタリングや補間も行わずに 単純な間引きのみを行う。 但し、パラメータmodeが指定された場合はその指定内容に従う。
Read the SAC file specified by the 1st argument, resample the data to a sampling interval specified by parameter new_dt, and output the result into the SAC file specified by the 2nd argument. The beginning time of the data is kept unchanged, and the data samples are simply reduced without filtering nor interpolation. If parameter mode is specified, the output data is based on that parameter.

◆使用例(Example)

sacfile_decimate NU.NKY.U.sac NU.NKY.U.10Hz.sac --new_dt=0.1


◆使用上の注意(Note)

このコマンドではアンチエイリアスフィルターは使用せずに 単純なサンプルの間引きのみを行う。 必要であればこのコマンドの前にフィルタリングを行うこと。
An anti-alias filter is not applied by this command; data samples are simply reduced. If necessary, apply the filter separately before executing this command.

◆補足(Additional information)

引数modeを用意した目的は火山爆発に伴う空振パルスのような短時間の強いシグナルを 長時間の連続波形データから目視で見つけ出すことを容易にすることである。 長時間の連続波形データをプロットする際にはサンプルを間引かないと処理に時間がかかる。 しかし単純に間引くと瞬間的な強いシグナルが消えてしまう。 そこで、このコマンドに--mode=maxあるいは--mode=abs_maxのオプションを付けて実行することで 瞬間的な強いシグナルは残したままサンプルを間引くことができる。
The purpose of argument mode is to facilitate manual finding of short-duration intense signals, such as an infrasonic pulse caused by volcanic explosions, in long continuous waveform data. When plotting the long continuous waveforms, downsampling the data significantly shortens the processing time. However, simply downsampling the data results in erasing the short but intense signals. Using --mode=max or --mode=absmax options in this command, the data can be downsampled remaining the short but intense signals.