狠狠撸

地磁気データ解析讲习
in データ解析講習会
2012/Aug/10

担当：阿部（九大ICSWSE）

2012/08/10 IUGONETデータ解析講習会＠NIPR 1

IUGONETで取り扱っている地磁気?地磁気活動度データ
? 下はIUGONET参加機関の所有する観測機器と範囲を示した図。
地磁気による超高層分野の研究には長い歴史がある。IUGONET
においても複数の参加機関が地磁気観測をおこなっており、地上
全域をカバーしている。



? 現在UDASに組み込まれているIUGONET参加機関の地磁気デー
タとload procedureは以下のとおり（制限のあるものもあり。詳細
は各load procedureやウェブサイトを参考）
? http://www.iugonet.org/software/loadprocedures.html

観測データ提供機関ロードプロシージャ名
AE, Dst, ASY/SYM指数、観測所
地磁気データWDC形式1時間京大WDC iug_load_gmag_wdc
値?1分値

erg_load_gmag_nipr
昭和基地、アイスランド地磁気極地研
(iug_load_gmag_nipr)

erg_load_mag_mm210
210°地磁気観測網名大、九大
(iug_load_gmag_mm210)
MAGDAS地磁気九大 iug_load_mag_serc

鲍顿础厂による厂厂奥イベント（2009年1月24日）时の
IUGONET機関地磁気データプロット

?UDAS使用のおおまかな流れ
?UDASを初期化する
?解析する时间幅を入力する
?load procedureでデータを読む
?迟辫濒辞迟変数を确认する
?そのままプロットする
?データや时间幅を処理してプロットする


解析する时间幅を入力する
解析時間幅は、以下の書式で入力

THEMIS〉timespan,’yyyy-mm-dd/hh:mm:ss’[,n,/<option>]

例1）2012年1月23日から1日分を指定
THEMIS〉timespan,’2012-01-23’ ←（,/1,/dayを続けても可）

例2）2009年1月1日から31日間を指定

THEMIS〉timespan,’2009-01-1’,31,/day

今回は2009年1月24日のSSWイベントを解析するので、
例2の時間幅を使う


load procedureでデータを読む
? データの読み込みは、各機関提供のload procedure
を使う(詳細は各load procedureの中や、crib sheetを
読むこと)
THEMIS〉erg_load_gmag_nipr, site=[‘syo’]
THEMIS〉erg_load_gmag_mm210,site=[’ktb’]
THEMIS〉iug_load_gmag_serc,site=[‘asb’, ‘anc’, ‘dav’]
THEMIS〉kyoto_dst_load ※udasロードプロシージャではありません
? 以上の4操作だけで
– 極地研の昭和観測点
– MAGDASの芦別、アンコン、ダバオ観測点
– 210MMのコトタバン観測点
– 京大WDCのDst指数(リアルタイム値)
のデータを読み込み、極域から赤道に至る
磁場変動と地磁気インデックスを見る準備ができる

load procedureでデータを読む
? データを読み込んだ
際、各機関のrules
of the roadがIDLウィ
ンドウに表示される。
内容を必ず確認して
おくこと！

←データ利用者として
必ず守るべきルール
や謝辞の方法などが
記載されています！

迟辫濒辞迟変数を确认する
THEMIS〉tplot_names THEMIS> tplot_names
% Compiled module: TPLOT_NAMES.
今回の例では。。。 1 nipr_mag_syo_1sec
2 mm210_mag_ktb_1min_hdz
3 mm210_mag_ktb_1h_hdz
4 magdas_mag_anc
5 magdas_mag_asb
6 magdas_mag_dav
7 kyoto_dst

各tplot変数の詳細を知るには
THEMIS〉tplot_names, tplot変数名or 番号,/verbose
THEMIS〉tplot_names, ‘magdas_mag_asb’ ,/verbose

そのままプロットする
THEMIS〉tplot,[‘tplot変数名’,またはtplot変数番号 …]
window, 0, xsize=1000, ysize=600
tplot_options, 'region', [0.05, 0, 1, 1]
tplot,[‘kyoto_dst‘,'nipr_mag_syo_1sec', 'magdas_mag_asb‘,
'mm210_mag_ktb_1min_hdz', 'magdas_mag_dav,’,
'magdas_mag_anc‘, ]


データや时间幅を処理してプロットする
;MAGDAS地磁気を各成分に分解
;H成分から平均値を差し引く
split_vec,'magdas_mag_asb‘
tsub_average, 'magdas_mag_asb_0'
split_vec,'magdas_mag_anc‘
tsub_average, 'magdas_mag_anc_0'
split_vec,'magdas_mag_dav‘
tsub_average, 'magdas_mag_dav_0' 簡単な処理で、
;コトタバン地磁気を各成分に分解見たい部分だけを
split_vec,'mm210_mag_ktb_1min_hdz' はっきりと表示する。
;昭和地磁気を各成分に分解 H成分のSSWイベント周辺数
;H成分を1分平均値化日だけを切り出してみよう。
split_vec,'nipr_mag_syo_1sec'
avg_data,'nipr_mag_syo_1sec_x',60
; 各観測点のH成分とDst指数をプロット
tplot,['nipr_mag_syo_1sec_x_avg','magdas_mag_as
b_0-d','mm210_mag_ktb_1min_hdz_x','kyoto_dst']
時間軸を変更
tlimit, '2009-01-16/00:00','2012-01-28/00:00'


;正規表現によるtplot。
;コトタバン地磁気各成分（絶対値）
tplot,[‘mm210_mag_ktb_1min_hdz_?’]
変数の指定に
正規表現が可能

;時間軸を直前のもの（1か月間）に戻す
tlimit, /last

;時間微分は deriv_data
deriv_data, ‘magdas_mag_anc_0’
deriv_data, ‘magdas_mag_dav_0’
tplot, ‘magdas_mag_*_0_ddt’


生プロットだけではなく
; 赤道域、低緯度、極域データデータ解析後の表示もできる
をウェーブレット解析して並べるウェーブレット解析をしてみる
wav_data,'magdas_mag_dav_0'
,/kol ,maxpoints=24l*3600*2
wav_data,'magdas_mag_asb_0'
,/kol ,maxpoints=24l*3600*2
wav_data,'nipr_mag_syo_1sec_
x_avg',/kol ,maxpoints=24l*360
0*2
; カラープロットの幅指定
zlim,'*pow', .0001,.1,1
; ウェーブレット表示
tplot,['magdas_mag_dav_0_wv
_pow',
'magdas_mag_asb_0_wv_pow',
'nipr_mag_syo_1sec_x_avg_wv
_pow']
tlimit,/full

; 大量データのウェーブレットは時間がか
かるので、必要なデータを切り出す生プロットだけではなく
newname = 'nipr_mag_syo_1sec_x' + データ解析後の表示もできる
'_clip'
ウェーブレット解析をしてみる
t1=time_double('2009-01-23/00:00')
t2=time_double('2009-01-24/23:59:59')
trange_clip, 'nipr_mag_syo_1sec_x', t1, t2,
newname=newname
;昭和地磁気H成分をウェーブレット解析
wav_data,'nipr_mag_syo_1sec_x_clip',/k
ol ,maxpoints=24l*3600*2
;カラープロットの幅指定
zlim,'*pow', .0001,.1,1
;昭和地磁気H成分と
;ウェーブレット結果の表示
tplot, ['nipr_mag_syo_1sec_x_clip',
'nipr_mag_syo_1sec_x_clip_wv_pow',
'kyoto_dst']
;時間軸を拡大
tlimit, '2009-01-23/00:00','2009-01-
25/00:00'

;昭和地磁気H成分にPi 2帯（40s-150s）で
バンドパスフィルタ
;tplot変数から値を取り出す TDASのプロシージャにない処理でも
get_data, 'nipr_mag_syo_1sec_x_clip', tplot変数から値を取り出して計算できる。
data=x ここではサンプルとしてバンドパスフィルタ
;thm_lsp_filterを逐次実行を手動でかけてみる
flow = 1d/150
fhigh =1d/40
dt=1
db=120.0 ;畳み込みとtplot変数再格納
nyquist = 0.5d/dt x.y = convol(x.y,cofs,/edge_t,/nan)
fhigh = double(fhigh/nyquist) < 1.d store_data, ‘nipr_mag_syo_1sec_x_clip_filt’,
flow = double( flow/nyquist) > 0.d data=x
fmin = min([flow, fhigh])
npts = long(!pi/fmin) > 1 ;SYOのH成分とバンドパスとウェーブレット表示
npts = npts < n_elements(x.y) ylim, ‘nipr_mag_syo_1sec_x_clip_filt', -10, 10
;デジタルフィルタの定義 tplot,['nipr_mag_syo_1sec_x_clip','nipr_mag_sy
cofs = digital_filter(flow,fhigh,db,npts, o_1sec_x_clip_filt’, ’nipr_mag_syo_1sec_x_clip_
/double) wv_pow’]


;フーリエ変換によるパワースペクトル表示
newname2 = 'nipr_mag_syo_1sec_x_clip' + '_psd'
tdpwrspc, 'nipr_mag_syo_1sec_x_clip', newname=newname2, nboxpoint=512
tplot,['nipr_mag_syo_1sec_x_clip','nipr_mag_syo_1sec_x_clip_filt',
'nipr_mag_syo_1sec_x_clip_wv_pow','nipr_mag_syo_1sec_x_clip_psd']
フーリエ変換を用いたパワー
スペクトル計算も可能


狠狠撸

Iugonet 20120810 abe_analysis

Recommended

More Related Content

What's hot (20)

Viewers also liked (6)

Similar to Iugonet 20120810 abe_analysis (20)

More from Iugo Net (20)

Iugonet 20120810 abe_analysis