野辺山偏波計（NoRP）解析マニュアル

ver. 1.2
野辺山宇宙電波観測所

2015 年 4 月 1 日

1

目次
1

はじめに

32

解析ソフトを使うには

33

実際の解析

3.1
3.2

データベース . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .

334

3.3

3.2.1 イベントリスト . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
3.2.2 一日分データの平均値 . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
データのありか . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .

444

キャリブレーション済データの XDR(IDL save) ファイル . . . . . . . . . . . . . . .

45

解析の開始

3.3.1
3.3.2. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .

生データ . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .

3.43.53.6

生データの読み込み — 野辺山偏波計 . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .

3.7スペクトル解析 . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .

3.83.9

生データの読み込み — 旧豊川偏波計 . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
グラフ表示. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .

3.7.1

光学的に薄い非熱ジャイロシンクロトロン放射 （試験公開） . . . . . . . . . . . . .

3.7.2

光学的に薄い熱制動放射 （試験公開） . . . . . . . . . . . . . . . . . . . . . . . . .

具体的なデータ解析の例 . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .

3.8.1

イベント選択 . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .

3.8.23.8.3

光度曲線 . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
スペクトル解析 . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .

その他 . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .

3.9.1

生データを直接よみこむ . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .

A NoRP データアーカイブ匿名 FTP サイト

55667888899992はじめに1

このマニュアルは、国立天文台野辺山宇宙電波観測所にて運用している偏波計（NoRP; Nakajima et al.

1985）と旧豊川偏波計（TYKW; Torii et al. 1979; Shibasaki et al. 1979）とのデータ解析ためのマニュアルである。
このパッケージに関する質問や要望は以下に電子メールを送ってほしい。
norp-help[--at--]solar.nro.nao.ac.jp[--at--] -> @解析ソフトを使うには

2

ここで紹介する偏波計解析ソフトを使うには以下の設定が必要である。

(1) SolarSoftware のインストール「SolarSoftware（SSW）」をサブパッケージ「Radio/NoRP」を含めてインストールする。これ以外にも Radio/NoRH と Yohkoh/SXT とを含めることを強く推奨する。これらのパッケージがインストールされていない場合はコンピュータ管理者に相談してほしい。
「SolarSoftware」の配布元は、http://www.lmsal.com/solarsoft/である。このマニュアルでは SSW をインストールしたトップディレクトリを${SSW} と定義する。また NoRP
サブパッケージのトップディレクトリを${NORP}（通常は${SSW}/radio/norp）と定義する。

(2) NoRP データベースのインストール
国立天文台 NoRP データアーカイブ ftp://solar-pub.nao.ac.jp/pub/nsro/norp/dataからデータベースをコピーしてインストールする。${NORP}/data にインストールするのがデフォルト。

(3) 個人環境の設定
環境設定ファイル~/.cshrc 中でsetenv SSW SSW のディレクトリ1
setenv SSW INSTR ‘‘norp’’source ${SSW}/gen/setup/setup.ssw
を適当な箇所に書き加える。ただし環境変数${SSW} はサイトによって異なるので注意。また他の観測機器も同時解析する場合はsetenv SSW INSTR ‘‘norp norh sxt’’などとする。

3

実際の解析

3.1

解析の開始

この節で紹介する NoRH 用 IDL プロシジャはすべて SolarSoftware(SSW) に依存している。前節の設定を行った後実際に解析する際は SSW/IDL を立ち上げる。

unix% sswidl <CR>
SDAS で SSW/IDL を立ち上げる場合は、以下のコマンドを利用する。
unix% SSWidl <CR>

1 SDAS では /ssw

3

3.2

データベース

3.2.1

イベントリスト

1990 年から 1992 年 5 月までのイベントについてはデータベースがある (’evx’ データ) 2 。たとえば 1991
年 6 月 9 日 UT のイベントのリストを表示するには
IDL> norp pr evx,’1991-06-09’ <CR>
1991 年 6 月 5 日 UT から 6 月 8 日 UT までであればIDL> st day=’1991-06-05’ <CR>
IDL> ed day=’1991-06-08’ <CR>
IDL> norp pr evx,st day,ed day <CR>
17GHz で太陽放射束密度 1000SFU 以上のイベントのみの場合はIDL> norp pr evx,st day,ed day,criterion=1000 <CR>

3.2.2

一日分データの平均値

一日分データのうち静穏時のみを選びだした平均値がデータベース化されている（1987 年 11 月以後）。
プリバーストの太陽放射を得たり、長周期変動の研究に役立ててほしい。
旧豊川偏波計（1・2・3.75・9.4GHz）：1987 年 11 月から 1994 年 2 月まで
旧豊川偏波計（3.75GHz のみ）：1994 年 3 月から 1994 年 4 月まで
旧野辺山偏波計（17・35・80GHz）：1990 年 3 月から 1994 年 2 月まで
野辺山偏波計（1・2・3.75・9.4・17・35・80GHz）：1994 年 5 月から現在も稼働中
ただし、35GHz（2400 SFU）と 80GHz（9000 SFU）とは装置の設計上、静穏太陽放射束密度は固定された値になっているので注意。
読み込みは、たとえば 1999 年 1 月 1 日 UT から 6 月 30 日 UT までであればIDL> st day=’1999-01-01’ <CR>
IDL> ed day=’1999-06-30’ <CR>
IDL> norp rd avg,st day,ed day,timavg,fiavg,fvavg <CR>
表示は、1GHz の場合
IDL> utplot,timavg,fiavg(0,*) <CR>
2GHz の場合
IDL> utplot,timavg,fiavg(1,*) <CR>

3.3データのありか

3.3.1

キャリブレーション済データの XDR(IDL save) ファイル

NoRP データアーカイブ (第 A 節参照) にキャリブレーションが済んで XDR ( IDL save ) フォーマットになったファイルがいくらか置いてある。読み込むには IDL コマンド restore を使う。まずはこれらを転送して解析するのをすすめる。定常観測については時間分解能 1 sec のデータが、重要なイベントについては時間分解能 0.1 sec のデータがおいてある。

IDL> file=getenv(’NORP XDR’)+’/1999/08/norp19990828 0056.xdr’ <CR>
IDL> restore,file <CR>
IDL> help <CR>

2 1992 年 6 月以後のイベントについては電波ヘリオグラフのイベントリストを参照のこと

4

DAYFI
FIAVG
FREQ
FV
FVAVG
MVD
TIM

STRING
FLOAT
FLOAT
FLOAT
FLOAT
FLOAT
BYTESTRUCT= ’1999-08-28’= Array[7, 7590]= Array[7]= Array[7]= Array[7, 7590]= Array[7]= Array[7, 7590]= -> ANYTIM2INTS Array[7590]

freq に観測周波数が GHz 単位で（ここでは 7 周波数 1・2・3.75・9.4・17・35・80GHz）、tim が観測時
刻、fi・fv に放射束密度 I(R+L) 成分・V(R–L) 成分（単位 SFU）が、fiavg・fvavg は一日分の静穏時の平均値がはいっている。また mvd は、各データ点が有効・無効をあらわす配列で有効なとき 1、無効なときゼロ。

3.3.2

生データ

NoRP データアーカイブ (第 A 節参照) に生データの一部が公開されている。

3.4

生データの読み込み — 野辺山偏波計

野辺山偏波計データは、豊川偏波計との統合前の旧フォーマットと統合後の新フォーマットとがあるが解
析プログラムは共通である。旧フォーマットデータには 1、2、3.75、9.4 GHz のデータは含まれていないのでこれらに関するオプションは無視される。
実際に読む込むには、読み込みたい日付 (JST) を指定してIDL> day=’2000-4-23’
IDL> norp rd dat,day,mvd,tim,fi,fv,freq <CR>
または、直接ファイル名を指定してIDL> file=’./pl921102’
IDL> norp rd dat,file,mvd,tim,fi,fv,freq <CR>とする。freq に観測周波数が GHz 単位で（ここでは 7 周波数 1・2・3.75・9.4・17・35・80GHz）、tim が観測時刻、fi・fv に放射束密度 I 成分・V 成分（単位 SFU）がはいっている。また mvd は、各データ点が有効・無効をあらわす配列で有効なとき 1、無効なときゼロ。
ある時間区間だけを選びだして読み込むには

IDL> timerange=[’1992-11-2 2:00’,’1992-11-2 4:00’] <CR>
IDL> norp rd dat,day,mvd,tim,fi,fv,freq,timerange=timerange <CR>
この場合キャリブレーションは必要な区間だけに行われるので読み込みが速くなる。
ある周波数だけ読み込むには以下のようにする。たとえば 17GHz だけであればIDL> rdfreq=[0,0,0,0,1,0,0] <CR>
IDL> norp rd dat,day,mvd,tim,fi,fv,freq,rdfreq=rdfreq <CR>
1GHz と 9.4GHz とだけであればIDL> rdfreq=[1,0,0,1,0,0,0] <CR>
IDL> norp rd dat,day,mvd,tim,fi,fv,freq,rdfreq=rdfreq <CR>
つまり rdfreq の 7 つの各要素が順に 1・2・3.75・9.4・17・35・80GHz に対応しており値が 1 のときに読み込む。

3.5

生データの読み込み — 旧豊川偏波計

旧豊川偏波計データは、 1、2、3.75、9.4 GHz のデータが含まれている。

5

実際に読む込むには、読み込みたい日付 (JST) を指定してIDL> day=’1992-11-2’
IDL> tykw rd dat,day,mvd,tim,fi,fv,freq <CR>とする。freq に観測周波数が GHz 単位で（ここでは 4 周波数 1・2・3.75・9.4GHz）、tim が観測時刻、fi・fv に放射束密度 I 成分・V 成分（単位 SFU）がはいっている。また mvd は、各データ点が有効・無効をあらわす配列で有効なとき 1、無効なときゼロ。
直接ファイル名を指定してIDL> file=’./ty921102.01i’ <CR> ; 1GHz I 成分
IDL> file0pa=’./ty921102.0pa’ <CR> ; キャリブレーションデータ
IDL> filestt=’./ty921102.stt’ <CR> ; ステータスデータ
IDL> tykw rd dat,file,file0pa,filestt,mvd,tim,data <CR>とすることもできる。この場合は、各周波数の I・V 成分のどちらかしか一度に読めない。また生データのほかにキャリブレーションデータ・ステータスデータが含まれたファイル名も同時に指定しなければならない。
ある時間区間だけを選びだして読み込むには

IDL> timerange=[’1992-11-2 2:00’,’1992-11-2 4:00’] <CR>
IDL> tykw rd dat,day,mvd,tim,fi,fv,freq,timerange=timerange <CR>
この場合キャリブレーションは必要な区間だけに行われるので読み込みが速くなる。
ある周波数だけ読み込むには以下のようにする。たとえば 2GHz だけであればIDL> rdfreq=[0,1,0,0] <CR>
IDL> tykw rd dat,file,mvd,tim,fi,fv,freq,rdfreq=rdfreq <CR>
1GHz と 9GHz とだけであればIDL> rdfreq=[1,0,0,1] <CR>
IDL> tykw rd dat,file,mvd,tim,fi,fv,freq,rdfreq=rdfreq <CR>
つまり rdfreq の 4 つの各要素が順に 1、2、3.75、9.4 GHz に対応しており値が 1 のときに読み込む。

3.6

グラフ表示

グラフ表示するには、

IDL> mfreq=0 <CR>
IDL> norp plot,mfreq,mvd,tim,fi <CR>
ここで mfreq は、fi（複数の周波数分を含む）にはいっているデータのうち表示したいデータの周波数の番号である。
別の周波数のデータを重ねるには

IDL> mfreq1=1 <CR>
IDL> norp plot,mfreq1,mvd,tim,fi,/over <CR>とする。直接指定してIDL> utplot,tim(where(mvd(mfreq,*))),fi(mfreq,where(mvd(mfreq,*))) <CR>とすることもできる。

3.7

スペクトル解析複数の周波数データを組み合わせることでスペクトルを求めることができる。求まる量は α で定義はFν ∝ ν α6である。実際は α は全周波数帯にわたって一価ではなく、マイクロ波帯のスペクトルは次のような関数でだいたいフィッテングできる。

cνFν = F( ν )αtk {νb[ ( )]} { cαFν (ν/bν ) tk
ν αtn −αtk
1 − exp −≈αcν (ν/bνbFν ) tnfor ν  νbfor ν  νb(1)cν はターンオーバー放射束密度・αtk は
ここでフィッテングのパラメータは、νb はターンオーバー周波数・F
光学的に厚い側（低周波数側）のべき指数・αtn は光学的に薄い側（高周波数側）のべき指数での 4 変数である。実際のフィッテングは以下のようにおこなう。

IDL> day=’2000-4-8’
IDL> norp rd dat,day,mvd,tim,fi,fv,freq <CR>
IDL> norp rd avg,day,timavg,fiavg,fvavg <CR>
IDL> for m=0,6 do fi(m,*)=fi(m,*)-fiavg(m) <CR>
IDL> norp alpha,freq,fi,mvd,mvdfit,alpha tk,alpha tn,freqpk,fluxpk <CR>
フレア前静穏状態の放射束密度（一日の平均値）を読み込んだデータから差し引いてフレア成分だけをとりだし、それを norp alpha にわたす。mvdfit は得られたフィッティングが有効な場合 1 が無効な場合（たとえば信号が弱いなどの理由で）についてはゼロがはいっている。
注意： 通常はこのような単純な手順ではフィッティングできない。実際は (1) 積分時間を長めにとってS/N を良くする、(2) 上記の関数から明らかにはずれた点をフィッティングからはずす（理由：低周波数は放射メカニズムが異なる、高周波数は S/N が良くない）などの手続が必要である。手順は節 3.8.3 を参照のこと。

3.7.1

光学的に薄い非熱ジャイロシンクロトロン放射 （試験公開）

注意：改良についてのコメント歓迎。
放射領域物理量から放射を求めることができる。

Dulk (1985) の方法と Ramaty (1969) の方法とが用意されている。ここでは電子の分布関数が以下のようなものであると仮定している。

dN V (E)=KdE(EE0)−δ

E は電子のエネルギー（keV）
・N V (E) はエネルギー E 以上の電子数（particles）。この分布はパラメータδ と E0 と K とによって記述されるが、K のかわりに N V (E0 ) = K/(δ − 1)/E0δ−1 を与える方が物理的な意味が明確になるのでそちらを実際は用いる。以下のプロシジャでは Dulk (1985) にならい、E0 = 10 keVで固定している。

(a) 物理量から放射を求める
放射領域の、電子分布べき指数 δ・磁場強度（G）
・視線と磁場とのなす角度（度）、さらに（エネルギー10keV 以上の）非熱電子数（個）を入力すると、放射束密度（SFU）・円偏波率が以下のように得られる。
Dulk (1985) の方法ではIDL> dulk gysy,delta,bb,theta,nv,freq,fi,rc <CR>
さらに放射源の大きさ（立体角であたえる。単位は sterad）をあたえると光学的厚さが求まる。

IDL> dulk gysy,delta,bb,theta,nv,freq,fi,rc,omega,tau <CR>
同様に Ramaty (1969) の方法ではIDL> ramaty gysy,delta,bb,theta,nv,freq,fi,rc <CR>
IDL> ramaty gysy,delta,bb,theta,nv,freq,fi,rc,omega,tau <CR>7

(b) 放射から物理量を求める
Dulk (1985) の近似モデルにもとづいて、放射から放射領域物理量を求めることができる。
α をもとめたあと（節 3.7 参照）、非熱電子の分布関数のべき指数がただちに求められる。
IDL> norp alpha,freq,fi,mvd,mvdfit,alpha tk,alpha tn,freqpk,fluxpk <CR>
IDL> norp alpha2delta,alpha tn,delta <CR>

3.7.2

光学的に薄い熱制動放射 （試験公開）

注意：改良についてのコメント歓迎。

(a) 物理量から放射を求める
Dulk (1985) の近似モデルにもとづいて、放射領域物理量から放射を求めることができる。
放射領域の電子温度（K）
・視線方向磁場強度（G）
・ボリュームエミッションメジャ（cm−3 ）を入力すると、放射束密度（SFU）・円偏波率が以下のように得られる。

IDL> dulk frfr,te,b loc,vem,freq,fi,rc <CR>
さらに放射源のサイズ（sterad）をあたえると光学的厚みも得られる。

IDL> dulk frfr,te,b loc,vem,freq,fi,rc,omega,tau <CR>

3.8

具体的なデータ解析の例

この節では解析の例を記す。

3.8.1

イベント選択重要なイベントについては、偏波計のウェブサイトに紹介されている。まずはそこをみて概要をつかむのがのぞましい。URL は

http://solar.nro.nao.ac.jp/norp/

3.8.2

光度曲線

イベントを選びだしたら、光度曲線をみてみよう。重要なイベントについてはキャリブレーションの済んだデータが XDR（IDL save）フォーマットで保存してあるのでそれを読む。

IDL> file=getenv(’NORP XDR’)+’/1999/08/norp19990828 0056.xdr’ <CR>
IDL> restore,file <CR>
グラフ表示するには、
IDL> mfreq=0 <CR>
IDL> norp plot,mfreq,mvd,tim,fi <CR>
ここで mfreq は、表示したいデータの周波数の freq の中での順番である。
別の周波数のデータを重ねるには

IDL> mfreq1=1 <CR>
IDL> norp plot,mfreq1,mvd,tim,fi,/over <CR>とする。

8

3.8.3スペクトル解析

スペクトルを計算する際は積分時間を長くとった方がよい。300 点（30 秒間）分のデータを積分するには

IDL> mint=300 <CR>
IDL> norp mkint,mint,mvd,tim,fi,fv,mvdav,timav,fiav,fvav <CR>
次にフレア前静穏電波放射束を差し引く。
IDL> for m=0,6 do fiav(m,*)=fiav(m,*)-fiavg(m) <CR>
えられた結果をスペクトルフィッティングする。
IDL> norp alpha,freq,fiav,mvdav,mvdfit,alpha tk,alpha tn,freqpk,fluxpk <CR>
放射がジャイロシンクロトロン放射であると仮定すると非熱電子のべき δ がただちにもとまる。
IDL> norp alpha2delta,alpha tn,delta <CR>

また、フィッティング関数（方程式 1）から明らかにはずれた観測点があるとフィッティングがうまくいかないことがある。理由は、低周波数では放射メカニズムが異なる（ジャイロシンクロトロンではなくプラズマ放射）のと、高周波数は S/N が良くない（特に 80GHz はよくない）ためである。そのようなときはその観測点をフィッティングから以下のようにしてはずす。

IDL> mvfreq=[2,3,4,5,6] <CR>
IDL> norp alpha,freq(mvfreq),fiav(mvfreq,*),mvdav(mvfreq,*) $ <CR>
IDL> ,mvdfit,alpha tk,alpha tn,freqpk,fluxpk <CR>

フィッティングした結果を確かめるために、フィッティング関数と観測点とを重ね合わせてみる。たとえば n ステップ目のデータを見るには

IDL> norp funcp,alpha tk(n),alpha tn(n),freqpk(n) ,fluxpk(n),freqfit,fifit <CR>
IDL> plot,freq,fiav(*,n),/xlog,/ylog,psym=5 <CR>
IDL> oplot,freqfit,fifit <CR>

3.9

その他

3.9.1

生データを直接よみこむ

生データを直接よみこむ方法は、野辺山データはIDL> norp rd rdt,file,rdata,date,version <CR>
豊川データはIDL> tykw rd rdt,file,rdata <CR>
データに対応する時刻（JST）を得るには

IDL> timej=norp rdt2timej(rdata) <CR>A

NoRP データアーカイブ匿名 FTP サイト

観測データベースファイルは日々蓄積されていくので適宜更新する必要がある。国立天文台天文データセンターでは生データを含むデータベースを匿名 FTP サーバで公開している。NoRP データの URL は

ftp://solar-pub.nao.ac.jp/pub/nsro/norp

9謝辞
以下の方々に感謝します。

SolarSoftWare は、太陽観測データ解析用の IDL パッケージで、NASA のさまざまなプロジェクトのもとで支援を受けています。IDL は Research Systems Inc. が製造・販売しているデータ解析ソフトウェアです。

参考文献
Dulk, G. A., 1985, ARAA, 23, 169
Nakajima et al. ”The Radiometer and Polarimeters at 80, 35, and 17 GHz for Solar Observations atNobeyama”, PASJ, 37, 163 (1985)
Ramaty, 1969, ApJ, 158, 753
Shibasaki et al. ”Solar Radio Data Acquisition and Communication System (SORDACS) of Toyokawa
Observatory”, Proc. of the Res. Inst. of Atmospherics, Nagoya Univ., 26, 117 (1979)
Torii et al. ”Full-Automatic Radiopolarimeters for Solar Patrol at Microwave Frequencies”, Proc. of theRes. Inst. of Atmospherics, Nagoya Univ., 26, 129 (1979)10