中級編-2: 2次元データ解析1 （オーロラ画像、 SuperDARNレーダーデータ等)堀智昭(名古屋大・STE研) Hori, T., 2次元データ解析1, IUGONET解析講習会@NIPR    8/20/2014 1 | T | { H Science Cente
♪　IDLが起動したままになっていたら,一旦終了→再起動させて下さい.
 　メモリーをクリア,カラーテーブル初期化♪ IDLのコマンド検索パスの順番で, SDツールフォルダ、UDASフォ
ルダがTDASフォルダより先、になっていることを確認して下さい。 › この講習で使用するSDツールはSPEDAS-j githubサイトにあります。 ▷  https://github.com/spedas-j/erg_sd_tool/tree/iugws_201408またコマンドリストも以下にありますので、コピペして実行したりできるようにWebブラウザ等で開いておいて下さい。 ›  https://github.com/spedas-j/erg_sd_tool/blob/iugws_201408/iugonet_ws_advanced_2_crib.pro今回の講習後に自習したり辞書的に使ってもらえるように、この講
習資料は網羅的に書いてあります。ただ今日の講習では時間が限られているので、RTIプロット、2D地図プロットのそれぞれ重要な箇所についてのみやります。 >  2
11 R F Science ( ▷ 2次元データの可視化・解析をやってみる。 １では2次元データの例としてSuperDARNデータを用い△る。(2では全天力メラデータ)典型的な可視化形態であるRTIプロット、2Dマッププロットを作成できるようにする。地図上にデータをプロットする方法を学ぶ。使用するデータはERG Science Center (ERG-SC)) (http://ergsc.stelab.nagoya-u.ac.jp)から公開されているCDF形式のCommon time fitacfデータ詳細はHori et al., Adv. Polar Sci., 2013を参照のこと。 http://journal.polar.gov.cn/EN/10.3724/SP.J.1085.2013.00069 SuperDARNデータは各レーダーのPIから提供されている。この講習で使うIDLプログラムはERG-SCとIUGONETとの共同開発C plug-inは、 THEMISチームによって開発されたTHEMIS Inalysis Software suite ( TDAS, http://themis.ssl.berkeley.edu/software.shtml)をベHori, T., 2次元データ解析1, IUGONET解析講習会@NIPR 8/20/2014謝謝辞3 FRG Science Center
iugonet △ 4 Hori, T., 2次元データ解析1, IUGONET解析講習会©NIPR 8/20/2014
IUGONEI HF帯・VHF帯電波のray pathと電離圏電子密》 HF帯電波は電子密度の勾配で屈折する。電離層中のdecameter (~10¹m)スケールのirregularityでBragg散乱を受け、極一部が後方散乱される。散乱されて戻ってきた電波のドップラーシフトを計測することで電離圏フラズマの水平ドリフト速度を知ることができる。この密度縞が無いと工コーが返ってこないので計測できない。 △ 5
Hori, T., 2次元データ解析1, IUGONET解析講習会@NIPR
8/20/2014 FFR G Science Center
1 beamを3秒、扇型を一通りスキャン》ドップラー速度などの計測を、送信ビームに沿って75-110 range gateに分割した各ピクセル毎に行う(最近は24 Beamのレーダーも)。 › ビームを方位角方向に振ることで扇型の領域内をスキャンする。 △ 6 Hori, T., 2次元データ解析1, IUGONET解析講習会@NIPR　 8/20/2014 ERG Science Center From http://vt.superdarn.org/tiki-index.php?page=Radar+Overview Hori, T., 2次元データ解析1, IUGONET解析講習会@NIPR     8/20/2014 ★ 7 ERG Science Center
iugonet ERG Science Center
INGONEI ★ 9 Hori, T., 2次元データ解析1, IUGONET解析講習会@NIPR 8/20/2014 11:11 Science Cen
個々のレーダーによる視線方向Doppler速度、工コー強度、スペクトル幅の観測(~10秒—1分) 1ビームに3秒→ 2次元スキヤン(通常16ビーム)は1
分で得られるデルERG-SC plug-inのSDライブラリ、及びこの講習資料はこっちのプロット形態用。マップポテンシャルプロットを作るにはVirginia Techのサイト(この資料の最後にある補足資料を参照)でやるのがよい。またSDレーダーPIに依頼してもよい。 Hori, T., 2次元データ解析1, IUGONET解析講習会@NIPR    8/20/2014 [Cousins+2010 (the original technique developed by Ryohoniemi+1998)]マップボテンシャルと呼ばれる1   10 FRG Science Center G
特に個々のレーダーのデSC/SI-induced transient flows [Hori+2012] ionospheric disturbance ▷  11 (LSTID) [Hayashi+2019]ori, T., 2次元データ解析1, IUGONET解析講習会©NIPR 8/20/2014 FR F Science Center
[Liu+2011] [Hori+2012]緯度-経度グリッド上での2D地図プ"* [Ebihara**2008]ロットの例Range-Time-Intensity (RTI)プロットある時刻の2D空間プロファイルを見たの例り、他の地上・衛星観測との位置関係変動のタイミング・周期を詳細に見るたを見るために使われる。めに使われる-- Hori, T., 2次元データ解析1, IUGONET解析講習会@NIPR 8/20/2014 A 12日本人Science Center

データ読込
▷ erg_load_sdfit SuperDARNデータをダウンロード・読込データ加工・RTIプロット▷ splitbeam tplot変数をビーム別に分割ムset_coords縦軸の座標系を変更ムloadct_sdカラーデニブルをセットget_fixed_pixel_graph Yつの歴クセルの値を独立したtplot変数にするA get_sd_ave _(*)ある緯度経度範囲の平均値を計算get_sd_lat_profile (*)純度プロファイルの万オグラムを作成》 asciidump_scan
データをアスキーダンプ2D地図プロット環境設定map2d_init (旧sd_init) 2D地図描画関係の変数初期化▷ map2d_time (旧sd_time) 2D地図プロット用の時刻をセット▷ map2d_set (旧sd_map_set)地理座標及びAACGM座標緯度-経度グリッドを定義▷ map2d coord 2D地図プロット用の座標系をセット▷ map2d_grid緯度-経度グリッドを点線で描画
SuperDARNデータの2D地図プロット
▷ plot_map_sdfit SDデータを2D砲図上にプロットするall in oneコマンド▷   overlay_map_sdfit
データを緯度-経度グリッド上にプロットト  overlay_map_sdfov
レーダーの視野をチロット▷ make_fanplot_pictures 2Dプロットをまとめて作成2D地図プロット関連overlay_map_coast
地図(海岸線)をプロットoverlay_map_vec (*)
2D地図上に線を播画(*)今回は資料に入れられなかったので、興味のある方はソースコード先頭のドキュメントをご覧下さい。 Hori, T., 2次元データ解析1, IUGONET解析講習会@NIPR     8/20/2014 > 13 14 Hori, T., 2次元データ解析1, IUGONET   8/20/2014

erq_load_sdfitを使う
get_support_dataキーワードをセットすることで、窓換に必要なパラメータなども一緒に読み込むsiteキーワードでレーダー名を指定する(以下以外もある) hok:北海道-陸別短波レーダー ksr: King Salmonレーダー (アラスカ) sye:  Syowa Eastレーダー (南極昭和基地) sys: Syowa Southレーダー ( pyk: Pykkvibaer (Iceland East)レーダー (アイスランド) han: Hankasalmiレーダー (フィンランド) erg_load_sdfitコマンド1つでfitacfデータをCDFファイルにしたものを自動ダウンロード(ローカルPC上に保存される)各種パラメータをtplot変数としてIDL上に読み込むをやってくれるデータの保存先はWindows: c:¥data Unix系: ~/data     の下THEMIS> timespan, '2013-09-23' THEMIS› erg_load_sdfit, site='pyk', /get_support_data Range-Time-Intensity (RTI)プHori, T., 2次元データ解析1, IUGONET  8/20/2014解析講習会@NIPR 16
THEMIS> tplot, 'sd_pyk_pwr_2' TPLOT(398):   5 sd_pyk_pwr_2 THEWIS> THEMIS> tplot, ['sd_pyk_pwr_2', 'sd_pyk_vlos_2','sd_pyk_spec_width_2'] TPLOT(410):   5 sd_pyk_pwr_2 TPLOT(410):   9 sd_pyk_vlos_2 TPLOT(410):   7 sd_pyk_spec_width_2 THEWIS> Beam0, Beam1, Beam2, …, Beam15の順に1ビーム約3秒ずつ観測していく。右図ではビーム順に時間方向に並べてプロットしている。この時Pykkvibaerレーダーは70レンジゲートモードだった。 Hori, T., 2次元データ解析1, IUGONET解析講習会@NIPR 8/20/2014 1   17
RTIプロットでは通常各beamごとに時間変化を見Beam05, 07, 09を並べてみた。しかしカラーバーのスケールが合ってない... > 18 Hori, T., 2次元データ解析1, IUGONET解析講習会@NIPR 8/20/2014
THEMIS› zlim, 'sd_pyk_vlos_2_azim??', -800, 800 THEMIS> 36): 47 sd_pyk_vlos_2_azim( IPLOT(486): 49 sd_pyk_vlos_2_azimi iH=WiS> zlim, 'sd_pyk_vlos_2_azim??'のようにワイルドカードを使うことで、...azim00 ~ ...azim15全てについて、カラーバーの上下限を設定できる。 * (アスタリスク)もワイルドカードとして使える。カラーバーの上限下限を-800, +800 [m/s]にしたらDoppler速度の変動が浮き出てきた！ Hori, T., 2次元データ解析1, IUGONET解析講習会©NIPR 8/20/2014

THEMIS› tlimit
:ウィンドウに十字スケールが出てくるので、プロットしたい時間帯の:最初の辺に合わせて走クリック１回:最後の辺りに合わせて左クリック１回、 :右の例では20:30, 00:00(右端)辺りをクリックして得られたTHEMIS> tlimit, /last時間で前回プロットした時の時間UREMIS> tlimit. /ful ;時間隔をtimespanで指定してある; THEMIS> tlimit, '2013-09-23/20:30', '2013-09-24/00:00' ;時間幅を自動きして指定することもでtlimit +マウスクリックtlimit, /last tlimit, /full tlimit, [開始時刻,終了時刻]を駆使することで、プロットの時間幅を自由に設定できるよくやるのは、 tlimit+マウスクリック→狙った時刻幅を少し外したのでtlimit,/lastで元に戻す→再度tlimit,
とか.自由自在にプロット時刻幅を変更できる.このように拡大すると, 22 UT直後のドップラー速度が
beam05→07→09と赤色主体から青色になっていることていることがわかる。 →異なるbeam(つまり異なる方向)で速度が反転している！ ★ 20 Hori, T., 2次元データ解析1, IUGONET解析講習会©NIPR 8/20/2014
sd_pyk_vlos_bothscat_2はsd_pyk_vlos_iscat_2とsd_pyk_vlos_qscat_2の２つを指すマルチtplot変数なので、以下のコマンドでiscat, gscat両方をビーム毎に分割することができるTHEMIS› splitbeam, 'sd_pyk_vlos_bothscat_2' ：電離圏エコーのみをプロットTHEMIS> zlim, 'sd_pyk_vlos_*scat*', -800,800 THEMIs› tplot, 'sd_pyk_vlos_iscat_2_azim05' ；電離圏エコーと地上エコー(灰色でマスク)をプロットTHEMIS› loadct_sd, 43     :灰色入りカラーテーブルをセットthemis› tplot, 'sd_pyk_vlos_bothscat_2_azim05'ただしここでの地上エコーの判定は経験式に基づいている.実際にはエコ強度,Doppler速度,スペクトル幅などから総合的に判断する必要がある。 > 21 Hori, T., 2次元データ解析1, IUGONET解析講習会@NIPR 8/20/2014
：縦軸を地理線度に変換THEMIS› set_coords, ['sd_pyk_vlos_2_azim??'], 'glat (縦軸をGeographical Latに変更,という表示) :再プロット(縦軸がrange gateから地理緯度になる) THEMIS> tplot,['sd_pyk_vlos_2_azim05','sd_pyk_vlos_2_azim09'] ;縦軸をAACGM緯度に変換THEMIS> set_coords, ['sd_pyk_vlos_2_azim??'], 'mlat' (縦軸をAACGM Latに変更,という表示) THEMIS> tplot,['sd_pyk_vlos_2_azim05','sd_pyk_vlos_2_azim09'] AACGM緯度~69°を境に、それより高緯度側(beam05)ではレーダーに向かう方向(つまり正の)のDoppler速度、低緯度側(beam09)ではレーダーから離れる向きの(負の)Doppler速度となるようなflow shearが観測されている！ AACGM: Altitude-Adjusted Corrected GeoMagnetic coordinates磁力線マッピングを考慮した地磁気座標の１つ[Baker+, JGR, 1989] ( http://www.aer.com/science-research/space/space- weather/altitude-adjusted-corrected-geomagnetic-coordinates )
ylimコマンドを使えば縦軸の範囲も自由に変更可能！ > 22 Hori, T., 2次元データ解析1, IUGONET解析講習会©NIPR 8/20/2014
THEMIS› set_coords, 'sd_pyk_vlos_2_azim??', 'gate' :再び機軸をrange gatelに変更THEMIS› get_fixed_pixel_graph, 'sd_pyk_vlos_2', $ beam=7, range_gate=18 ; beam07, range gate 18のpixelを選択THEMIS> tplot, ['sd_pyk_vlos_2_azim07'] sd_pyk_vlos_2_bm07rq018 ; get_fixed_pixel_graph ; beam07, range gate 18のpixelの値を取り出して、新しい: tpの繁殖物に整細する速度の値をグラフで確認することができるTHEMIS> get_data, ' sd_cvw_vlos_8_bm17rg032', data=dataで、値を通常のIDL変数(実際には構造体)に代入することも可能Hori, T., 2次元データ解析1, IUGONET解析講習会©NIPR 8/20/2014
/centerを付けるとカラーバーの真ん中辺りを透明にする。→ゼロ付近を目立たないようにできる。他のカラーバーにも適用できる。 THEMIS> loadct_sd, 44 THEMIS> tplot, 'sd_pyk_vlos_2_azim07' THEMIS> loadct_sd, 45, /center THEMIS> tplot, 'sd_pyk_vlos_2_azim07'次項の講習のためにloadct_sd, 43で元に戻しておいて下さい正の値を寒色系、負の値を暖色系で描ゼロのところで寒色系-暖色系を切り替画。SDプロットでは慣習的にこっちをえる→値の正負をわかりやすくする使うことが多い。 Hori, T., 2次元データ解析1, IUGONET解析講習会@NIPR 8/20/2014 ＼ 24
THEMIS> asciidump_scan, 'sd_pyk_vlos_2' The ascii files have been generated in /home/idl/work/IDL8_project/tmp/asci idump Program finished.いを使っと、カレント
リュというディレクトリュヘリにasciidu多数のファイルができるので注意(1スキャン1分→ 1日で1440個!) .今のところLOSVデータ(sd_STN_vlos_?
STNはレーダー名)のみに対応。 Hori, T., 2次元データ解析1, IUGONET解析講習会@NIPR 8/20/2014 2次元地図プロットの作成26 Hori, T., 2次元データ解析1, IUGONET    8/20/2014 FRO Science
▶ある時刻の２次元スキャンのデータを、緯度・経度グリッド(+世界地図)上に描画する観測値の2次元空間分布がわかる他の観測データを重ね描きすれば位置関係を調べることができるある１つの時刻のデータしかプロットすることができない(異なる時刻の複数のプロットを作る必要) Hori, T., 2次元データ解析1, IUGONET解析講習会©NIPR     8/20/2014
THEMIS> map2d_init THEMIS> timespan, '2013-09-23' THEMIS> map2d_time, 2200        ;描画時刻をセットplot_time: 2013-09-23/22:00:00 THEMIS> map2d_set, /erase :デフォルトで地理緯度経度グリッドを定義THEMIS> overlay_map_coast       ;地図を描画THEMIS> plots, [345., 15], [65,78] ; plotsコマンドにより(地理緯度,経度)で; (75°, 345°) ⇒ (78°, 15°)に直線を描くTHEMIS> polyfill，$ [15., 45., 45., 15.] ,  $ [70., 70., 75., 75.],      color = 240 ; (60°,17h), (60°,19h), (70°,19h), (70°,17h)の4点で;囲まれる台形を150番の色(この場合は黄緑)で塗りつぶHori, T., 2次元データ解析1, IUGONET解析講習会@NIPR     8/20/2014 FR G Science Center
iugonet map2d_init, map2d_timeは既に設定済みなので省略THEMIS> map2d_coord, 'aacgm'  ;AACGM座標を指定THEMIS> map2d_set, /erase, /mltlabel :AACGM緯度-LTグリッドを定義THEMIS> overlay_map_coast       ;地図を描画THEMIS> plots,  [0., 9] /24.*360.,   [70., 80.] ; plotsコマンドで(AACGM緯度,LT)で; (70°, 0.0h) → (80°, 9.0h)に直線を描くTHEMIS> polyfill，$ [17., 19., 19., 17.]  /24.*360.,   $ [60., 60., 70., 70.],      color = 150 ; (60°,17h), (60°,19h), (70°,19h), (70°,17h)の4点で;囲まれる台形を150番の色(この場合は黄緑)で塗りつぶす▷  map2d_setの他のオプションはヘッダードキュメント等を参照★ 29 Hori, T., 2次元データ解析1, IUGONET解析講習会@NIPR    8/20/2014 FRG Science Center
map2d_coord,'geo' & map2d_set  (地理座標用) map2d_coord,'aac9m' & map2d_set (AACGM座標用) map2d_setでは↑のような緯度-経度グリッドが定義される。地理座標では、plot, (経度), (緯度)のようにすればそ
のままプロットできる。座標系はmap2d_coordコマンドで指定する。 › AACGM座標の場合、韓度はAACGM緯度、AACGM地方時を0h→0° 9h→90° 12h→180° 18h→270°のように変換
した仮想的な経度にして、plot, (AACGM緯度), (仮想的な経度)とすれば対応する場所に作図できる。 》 overlay_map_????系のコマンドは内部で緯度,経度を上記のように計算して、点・線などを描画している。 ‣つまり同じように緯度,経度を与えるようにすれば、自作のプログラムでも同様なプロットをすることができる。 Hori, T., 2次元データ解析1, IUGONET解析講習会@NIPR     8/20/2014 ▲ 30 ERG Science Center
iugonet ;環境をセットアップするTHEMIS> map2d_init ;プロットする時刻を22:10 UTに指定するTHEMIS> map2d_time, 2210 ;指定時刻のLOS velocityデータを描画するTHEMIS› plot_map_sdfit, 'sd_pyk_vlos_2' ; coastキーワードをセットすると世界地図を重ねて描くTHEMIS> plot_map_sdfit, 'sd_pyk_vlos_2' , /coast ;HANデータをロードして一緒に描画するTHEMIS> erg_load_sdfit, site='han' , /get_support THEMIS> zlim, 'sd_han_vlos_2', -800,800 :スケールを揃えるTHEMIS> plot_map_sdfit, ['sd_pyk_vlos_2', $ 'sd_han_vlos_2'], /coast   : pykとhan両方をプロット；座標を地理座標に変えてプロットTHEMIS> map2d_coord, 'geo'   ;座標系切替コマンドTHEMIS> plot_map_sdfit, , ['sd_pyk_vlos_2', $ 'sd_han_vlos_2'], /coast THEMIS> overlay_map_sdfov, site='pyk han' :レーダーの視野を重ね描きHori, T., 2次元データ解析1, IUGONET解析講習会@NIPR     8/20/2014 FRA Science Center
INGONEI ;clipキーワードをセットするとズームイン。 :レーダー視野を外れることがあるのでcenter_glat ;center_glonキーワードで描画中心の地理雑度経度をJ形:する。 THEMIS> map2d_coord, 'aacgm'  ;AACGMに戻すTHEMIS> plot_map_sdfit, ['sd_pyk_vlos_2', $ 'sd_han_vlos_2'], /coast,  /clip, $ center_glat=75, center_glon=0 ; $(ダラー)を付けると1行を分割できる:拡大するとカラーバーがはみ出ることが多いので、手動で位置を調整するとよい。colorscaleposキーワードで指定する。 THEMIS> plot_map_sdfit, ['sd_pyk_vlos_2', $ 'sd_han_vlos_2'], /coast, /clip, $ center_glat=75, center_glon=0, $ colorscalepos=[0.05, 0.65, 0.08, 0.95] ; MLTラベルを重ね描き(mltlabelキーワード) , THEMIS> plot_map_sdfit, ['sd_pyk_vlos_2', $ 'sd_han_vlos_2'], /coast, /clip, $ center_glat=75, center_glon=0, $ colorscalepos=[0.05, 0.65, 0.08, 0.95], $ /mltlabel緯度線は80°,70°,60°...と10度刻み
Hori, T., 2次元データ解析1, IUGONET解析講習会©NIPR
8/20/2014
THEMIS> map2d_time, 2155 THEMIS> plot_map_sdfit, ['sd_pyk_vlos_2','sd_han_vlos_2'],/clip, /coast, center_glat=75, center_glon=0, position=[0.0,0.5,0.5,1.0] /nocolorscale plot map sdiit r /nocolorscaleまで数行動THEMIS> map2d_time, 2200して１行でタイTHEMIS> plot_map_sdfit, ['sd_pyk_vlos_2','sd_han_vlos_2'],/clip, /coast, center_glat=75, center_glon=0, position=[0.5,0.5,1.0,1.0], /noerase , /nocolorscale THEMIS> map2d_time, 2205 THEMIS> plot_map_sdfit, ['sd_pyk_vlos_2','sd_han_vlos_2'],/clip, /coast, center_glat=75, center_glon=0, position=[0.0,0,0.5,0.5]. /noerase, /nocolorscale THEMIS> map2d_time, 2210 THEMIS> plot_map_sdfit, ['sd_pyk_vlos_2','sd_han_vlos_2'],/clip, /coast, center_glat=75, center_glon=0, position=[0.5,0.0,1.0,0.5], /noerase, colorscalepos=[0.05, 0.65, 0.08, 0.95] ● positionキーワードにnormal座標でのプロットの位置を与える([x0, y0, x1, y1]) ● plot_map_sdfitはデフォルトで描画毎にウィンドウを
し2D地図プロットにすれば、スライド22のRTIプロットで見えていたflow shear構造が見易い。 > 33 Hori, T., 2次元データ解析1, IUGONET解析講習会@NIPR　8/20/2014 [ ] { f  Science Cent
THEMIS> !p.position = [0,0,1,1]    ;画面分割設定を初期化THEMIS> make_fanplot_pictures, $ ['sd_pyk_vlos_2','sd_han_vlos_2'], $ 2155, 2210, /clip, /coast, center_glat=75, $ center_glon=0, prefix='pngdir/sd_pykhan_' 2Dプロットの画像ファイル(png)をいっぺんに作成するコマンド。第2,3引数の2155, 2210の意味→ 21:55-22:10 UTの間のプロットを全て作成。出力先はprefixキーワードで指定する。上の例だと, pngdirというフォルダを作ってその中にsd_pykhan_hhmm.pngというファイル名で出力(hhmmは時刻)。
> 34 Hori, T., 2次元データ解析1, IUGONET解析講習会©NIPR 8/20/2014 RTIプロットや2Dプロットの応用編35 Hori, T., 2次元データ解析1, IUGONET    8/20/2014解析講習会@NIPR FRA Science Center
HIGONEI THEMIS> map2d_init THEMIS> map2d_time, '2013-09-23/22:10' ;編度総度抽画モードをオンにして、グリッドを描画する。 ; aacgm:セットするとAACGM座標で描画; glatc/gloncで描画する際の中心位置を指定。 ; mltlabel  キーワードをセットするとMLTのラベルを描; eraseキーワードをセットすると一度ウィンドウ内を消; scale:プロットのスケールTHEMIS> map2d_coord, 'aacgm' THEMIS> map2d_set, $ glatc=75,glonc=0, scale=35e+6, $ /mltlabel, /erase ;指定時刻のLOS velocityデータを重ね描きするTHEMIS> overlay_map_sdfit, $ ['sd_pyk_vlos_2','sd_han_vlos_2'], $ colorscalepos=[0,05, 0,65, 0,08, 0,95] ;世界地図を重ね描きするTHEMIS› overlay_map_coast実はplot_map_sdfitは内部で•map2d_set •overlay_map_sdfit •overlay_map_coast (/coastの場合)を順に実行している。 ▲ 36 Hori, T., 2次元データ解析1, IUGONET解析講習会@NIPR 8/20/2014 11:11 Science D  
INGONEI ›  map2d_setは内部でIDL標準のmap_setを呼んでいる. › map_setを宣言すると, plot,経度,緯度, …のようにして地図上に描画できるようになる. ♪ map_set上に描画するような他の(自作・他作)プロットルーチンと
組み合わせることで,複数種のデータを地図上に描画できる.
! ！注意！！データの読込でかなりメモリーを消費するので、 erg_load_sdfit, iug_load_asi_niprを実行したところでout of memoryエラーが出る可能性があります。その場合、もっとメモリーに余裕があるPCでやってみて下さい。
23:04 UT
THEMIS> timespan, '2012-03-27/23:00:00',1,/hourのデータの読込および2D地図プロットコマンド。 THEMIS> erg_load_sdfit, site='pyk', /get THEMIS> iug_load_asi_nipr, site='hus'  ≪ THEMIS> map2_time, 2304 & map2d_coord, 'geo' THEMIS> map2d_set, /erase, scale=17e+6, glatc=70, glonc=346 SDのピクセルを0.5倍にしてTHEMIS> loadct, 0  0番のカラーテーブル(白黒)をセットするTHEMIS> overlay_map_asi_nipr, 'nipr_asi_hus_0000', /nocolor, colorr=[20,140] THEMIS> loadct_sd, 44 & overlay_map_sdfit, 'sd_pyk_vlos_2' " pixel=0.5 THEMIS> overlay_map_coast, col=40地図を40番の色(茶色)で描画Hori, T., 2次元データ解析1, IUGONET解析講習会@NIPR ★ 38 8/20/2014 39
ほとんどのレーダーでbeam0が一番西端のbeam (例外有しーダーによってこれらの数
♪主な変数名と中身
pwr:エコー強度vlos: Line-of-sightドップラー速度(LOSV) spec_width:スペクトル幅vnorth: LOSドップラー速度の地理緯度成分(北向き) veast: LOSドップラー速度の地理経度成分(東向き) (vlos|vnorth|veast)_iscat:電離圏エコーのみのデー R (vlos|vnorth|veast)_gscat:地上エコーのみのデータ(vlos|vnorth|veast)_bothscat:電離圏・地上エコー両方のデータelev_angle: elevation angle値echo_flag:電離圏エコーか地上エコーかの判定フラグquality:データのqualityについての情報(0: good, 1以上: poor) quality_flag: quality判定の内訳(詳細は担当者へ) position_tbl:各pixelの四隅の緯度、経度値テーブルpositioncnt_tbl:各pixelの中心の緯度、経度値テーブルcpid: beam毎の観測モードtfreg: beam毎の周波数noise: beam毎のノイズレベル
College of Engineering v Radars Maps/Tables/Links Operating Schedule Radar Finder Ray-Tracing Tool hdw.dat Updater v Software
v Data
Latest Data (10-Plot CP D Inventory Data Requests SAPS Examples 24- hrs Overview Plot Range-Time Plot Time-Series Plot Map Velocity Tool Map Potential Plot Global Plotting Tool ASCII Potential Tool Dynamic FFT Plot Conjugacy Tool GPS/TEC Plot Res Bay 630nm Plot Convection Map Browser V  Diagnostics ACF Plot and Fit Badlags Calculator Badlag Finder Print Fit Records ACF Applet –ダーに関する情報。野プロットetc.データアーカイブ状況。
データリクエスト。各種プロットをオンデマンドで作成できる。 http://vt.superdarn.org | Maps/Tables/Links page最近消滅したJHU/APLのサイトに代わって今ではSuperDARNI口ジェクトのメインポータルサイトになっている.各種プロットをオンラインで作成できる. P b Convection mapプロットやGPS-TECプロットなど各種プロット有り特にConvection図作成はこのサイトが便利v DaViT Hori, T., 2次元データ解析1, IUGONET解析講習会@NIPR 8/20/2014