東西・南北成分の波形をradialとtransverse成分に回転する
(Last Modified: December 2, 2014)

■南北成分のSACファイルに必要なヘッダー情報をいれる

SAC> r 南北成分のSAC形式ファイル名
”データの読み込み”

SAC> ch stla 32.8795
SAC> ch stlo 131.0763
SAC> ch stel 0
”SACファイルのヘッダ部分に、本堂の観測点の緯度、経度、高度情報を入力” 

SAC> ch evla -15.49
SAC> ch evlo -172.10
SAC> ch evdp 18
”SACファイルのヘッダ部分に、該当する地震の緯度、経度、深さ情報を入力” 
上の値はあくまでも例です。適当な値に変更してください。

SAC> ch cmpaz 0
SAC> ch cmpinc 90
”SACファイルのヘッダ部分に、 データ成分の方位(cmaz)と鉛直からの角度(cmpinc)を入力” 

SAC> w tmp.sac.ns
"出来た南北成分の波形をtmp.sac.nsに保存する”

■東西成分のSACファイルに必要なヘッダー情報をいれる

SAC> r 東西成分のSAC形式ファイル名
”データの読み込み”

SAC> ch stla 32.8795
SAC> ch stlo 131.0763
SAC> ch stel 0
”SACファイルのヘッダ部分に、本堂の観測点の緯度、経度、高度情報を入力” 

SAC> ch evla -15.49
SAC> ch evlo -172.10
SAC> ch evdp 18
”SACファイルのヘッダ部分に、該当する地震の緯度、経度、深さ情報を入力” 
上の値はあくまでも例です。適当な値に変更してください。

SAC> ch cmpaz 90
SAC> ch cmpinc 90
”SACファイルのヘッダ部分に、 データ成分の方位(cmaz)と鉛直からの角度(cmpinc)を入力” 

SAC> w tmp.sac.ew
"出来た東西成分の波形をtmp.sac.ewに保存する”

■東西・南北成分を、radial成分(大円方向の水平成分)と transverse成分(大円方向に直交する水平成分)に座標回転した 波形をつくる

SAC> r tmp.sac.ns tmp.sac.ew
”南北・東西成分の地データの読み込み”

SAC> rotate
”座標を回転する”

SAC> w tmp.sac.r tmp.sac.t
”出来たradialとtransverse成分のデータを保存”
tmp.sac.r tmp.sac.tは、自分の好きなファイル名でよい。

SAC> p1
”出来たradialとtransverse成分のデータをみる”
上がradial成分、 下がtransverse成分。