EIS WORKSHOPSHINSUKE IMADA

コピーペースト用テキストファイル cmndlist.txt 参照

1. EIS データ解析前に必要な情報
自分が見たい EIS データを探すのに有用な WEB サイト、及びそのデータを簡単に見る（QL）方法を説明する。
まず、どのような観測を探したいか、以下のサイトから Study を探す。
http://msslxr.mssl.ucl.ac.uk:8080/SolarB/EisStudyList.jsp
続いて、以下のサイトからその Study が走っている日時を検索する。
http://solar-b.nrl.navy.mil/cgi-bin/timeline4.cgi
MSSL のサイトでも同じように検索できるが、こちらの方が見やすいのでおすすめ。

DARTS などで、見たい観測データを手に入れたら、IDL>xﬁlesとすると QL が立ち上がり簡単にどのような結果が得られているか見る事ができる。
(Viggo’s presentation を参照)
この他、MSSL のサイトから直接強度、速度、ライン幅の QL を見る事もできる。

EIS user がキャリブレーション情報などを交換または議論しているサイト
http://msslxr.mssl.ucl.ac.uk:8080/eiswiki/Wiki.jsp?page=DataProAnalysis

Date: 2009.07.13.
1

2

SHINSUKE IMADA

2. EIS データのキャリブレーション
EIS データのキャリブレーション方法について（中身については Viggo’s presentation 参照）

dark current, ﬂat ﬁeld 等の基本的なキャリブレーションを eis prep を用いて行う。
IDL>ﬁlename=’eis l0 20061206 191520.ﬁts.gz’
IDL> eis prep, ﬁlename, /default, /save（時間がかかるので今回は省略：もう少し小さいデータで後ほど）とすれば（数分かかる、ものによってはかなりの時間がかかる）、dark current, ﬂat ﬁeld 等のキャリブレーション済みの強度データ、及び誤差のデータが得られる。様々なオプションが存在する。doc library を参照してください。

3. EIS データの読み込み
EIS データの読み込みには様々な方法があるが、ここでは eis getwindata を紹介する。
読み込み方法はきわめて簡単で、読み込みたい波長（大体）を調べ
IDL>ﬁlename2=’eis l1 20061206 191520.ﬁts’
IDL> wd=eis getwindata(ﬁlename2,195.12)
IDL>help, wd, /strとすれば structure 形式で Fe XII 195.12A のデータを読み込む。

IDL> plot image,total(wd.int[*,*,*],1)とすれば強度のイメージを作る事が可能（Fig 1）。

4. EIS データのフィッテング
EIS データを自動でフィッテングし、強度、速度、ライン幅を求める。
（かなり時間がかかる）
IDL>eis auto ﬁt, wd, ﬁt, refwvl=195.12
IDL>help, ﬁt, /str
（時間がかかるので今回は省略）
フィティングした結果を表示するには

IDL>eis ﬁt viewer, wd, ﬁt

EIS WORKSHOP

Figure 1. Intensity Map:

3

4SHINSUKE IMADAとすれば良い。左上から強度、速度、ライン幅のデータが表示される。速度のイメージに縦の線が入っているが、これはフィッティングの計算が間違っている訳ではなく、ひので衛生の軌道周回変動の影響である。また実際には EIS のスリットは CCD に対し多少傾いているため、Artiﬁcial な速度を補正する必要がある。
注）1” slit の tilt は非常に小さいが、2” slit の tilt はやや大きい。
上であげた２つの Artiﬁcial な速度を補正する方法を説明する。
まず、スリットの tilt を補正するには、

IDL>outﬁt=eis tilt correction(wd, ﬁt)を走らせれば補正される。続いて、軌道周回変動を補正するには

IDL>outﬁt2=eis orbit spline(outﬁt, /setzero, /plotvel, spline=spline)を走らせれば良い。これで波長の補正も完了した。この補正方法は基本的にスリット方向にドップラーシフトが平均してゼロになるように補正する。つまり同じスリット位置でほとんどの領域がドップラーシフトしているような領域（例えばコロナホール）では間違った補正がかかる可能性が高い。軌道周回変動の補正の仕方はいくつか別の方法もあるが、現在、衛星のハウスキーピングデータを用いて補正をかけるような EIS のサイエンスデータによらないものも考案中である。
補正したデータを見るには以前と同様に

IDL>eis ﬁt viewer, wd, outﬁt2を用いれば良い。作成した outﬁt2 から自分なりのプロットを作るには、例えば

IDL>window, 0, xs=512, ys=1024
IDL>!P.multi=[0,1,3]
IDL>loadct, 3
IDL>plot image, outﬁt2.int
IDL>load vel
IDL>plot image, outﬁt2.vel, max=-50, min=50
IDL>loadct, 33
IDL>plot image, outﬁt2.wid, max=0.05, min=0.0223としてあげれば良い（Fig2）。
以上が基本的な EIS データの取り扱い方である。

EIS WORKSHOP

Figure 2. Intensity, Velocity, and Width Map:

56SHINSUKE IMADA

5. おまけ5.1. Map の作成
ヘッダーは wd.hdr にある。しかし、map 形式にしようとして、

IDL>index2map,wd.hdr,outﬁt2.int,map intensityとすると失敗する。なので少し面倒ではあるが、IDL>obj=obj new(’eis data’, ﬁlename2)
IDL>index long=ﬁtshead2struct((obj->gethdr())->getbte hdr())
ヘッダーに入っている EIS の pointing 情報は CCDLongWavelength の情報です。少なくとも FeXII のような ShortWavelength のデータは２つの CCD 間のズレを直す必要があります。
Y（スリット）方向に-16”、X（スキャン）方向に-2”。

IDL>index short=index long
IDL>index short.xcen=index short.xcen-2.
IDL>index short.ycen=index short.ycen-16.
map を作成する
IDL>index2map,index shrot,outﬁt2.int,map intensity
IDL>index2map,index short,outﬁt2.vel,map velocity
IDL>map velocity.data=-1.*map velocity.data
IDL>index2map,index short,outﬁt2.wid,map width
描画する

IDL>plot map,map intensity
IDL>plot map,map velocity, dmax=50, dmin=-50
IDL>plot map,map width, dmax=0.05, dmin=0.0223
5.2. スペクトルの表示
cmndlist.txt 参照

5.3. ひとつのウィンドウに２つ以上のラインが入っているケース
それなりに難解。以下を参照ください。
http://msslxr.mssl.ucl.ac.uk:8080/eiswiki/attach/DataProAnalysis/eis auto ﬁt gen.pdf

EIS WORKSHOP

7

6. 小さいデータでの練習問題
２００９年２月５日のデータをやってみよう。スキャンしている領域はディスク中心の静穏領域。2”スリットを使い小さい領域を高速でスキャンしているデータ。班１で使用するデー
タと同じ Study
解答は cmndlist.txt の prctice 以下にあるが自分で考えながらやってみよう。（Fig3）注意）軌道周回変動の補正が少し手を加えなければいけない

6.1. 他のデータとのアラインメント
小さい領域を高速スキャンしているデータは位置合わせが非常に難しい。どのブライトポイ
ントが EIT195 のどれに対応するかすらよくわからない事が多い。そういう場合は、少し前の時間を探し、割と大きい領域をはいているものを探す。そこで位置を合わせ EIS の視野の修正分を小さい視野のデータに適応すれば、5 秒各程度には位置が合ってくる。
２００９年２月５日のデータを EIT とあわせるには、X+41, Y-49
実演して簡単に説明する。

8

SHINSUKE IMADA

Figure 3. Intensity, Velocity, and Width Map: