pptx

Report
ひので SOT SP データ解析
データ解析ゼミ
2014.5.15 一本
SOT Spectro-Polarimeter
光球のベクトル磁場(Bx, By, Bz)、視線速度場(v)、温度(T)
の3次元分布(x, y, t)として求める装置。
時間分解能 > 数時間
空間分解能 x 精度 世界最高
SP 参考文献
• Hinode Solar Optical Telescope Data Analysis Guide, 2008.1.23
• ”The Hinode Spectro-Polarimeter”, Lites etal., 2013,
Sol.Phys..283..579L
• “The SP_PREP Data Preparation Package for the Hinode SpectroPolarimeter” , Lites and Ichimoto, 2013, Sol.Phys..283..601L
• 国立天文台ひのでサイエンスセンター
http://hinode.nao.ac.jp/index.shtml
SOT ポラリメータの構成
Collimator
lens unit
(CLU)
Polarization
modulator unit
(PMU)
1 rev./1.6 sec
HDM
CTM-TM
Astigmatism
corrector lens
(ACL)
M2
M1
Mech. Mask wheel
shutter
NFI- Polarization
analyzer
Reimaging lens
Tunable filter
FG/NFI
Non-polarizing
beam splitter
FG-CCD
Blocking filter wheel
SP
Slit scan mirror
Slit
SP- Polarization analyzer
(polarizing beam splitter)
SP-CCD left/right
16 frames/sec
3
Polarization modulation and demodulation
V
U
Q
Waveplate angle [deg.]
PMU segment
SP onboard
demodulation
0
+
+
+
-
1
+
+
-
2
+
-
3
+
+
-
4
+
+
+
+
5
+
+
+
6
+
+
7
+
+
+
8
+
+
+
-
9 10 11 12 13 14 15
+ +
+ + + +
+
- + +
+
+ - + +
- - + + +
+




I’
Q’
U’
V’
4
SOT observables
FG:
- simple image
- Dopplergram (2 wavelengths)
- Stokes IV
- Stokes IQUV
- IVDG (2 wavelengths)
SP:
- normal mode (0.16” step, 4.8s/slit)
- fast mode (0.32” step, 3.2s/slit)
- dynamic mode (0.16” step, 1.6s/slit)
- deep mode (0.16” step, 9.6s/slit)
5
データの探し方
NAOJ ひのでサイエンスセンターQL
http://hinode.nao.ac.jp/QLmovies/
LMSAL operation info.
https://sot.lmsal.com/operations/timeline/
SP 生データ
SP4D20070228_183214.2.fits
I
-Q
-U
level-0
-V
I
Q
U
V
sp_prep.pro
(標準ソフト)
処理内容:
-
ダーク・フラット補正
偏光校正
分散方向の傾き補正
スペクトルの曲がり補正
CCD両サイド(直交スペクトル)の結合
波長反転
周回ドリフト補正
使い方: ex_sp_prep2.pro
HSC にて
@/home/ichimoto/idlpro/solarb/analysis/file_io/sot_filelist.pro
;topdir ='/isas/data/solarb_fits/sot/QuickLook'
topdir ='/isas/data/hinode_fits/sot/level0'
date='20070202' & start_time=date+'_121300' & end_time=date+'_150000‘
outdir='/nfs/scr01/ichimoto/spcal/'+date+'/‘
files=sot_filelist('SP4D',start_time,end_time,topdir=topdir)
sp_prep, files, outdir=outdir, /verbose;
end
9
SP 1次処理データ
SP4D20070228_183214.2C.fits
I
level-1
Q
U
V
Definition of SOT polarization coordinate
This definition is applied to the Stokes vectors obtained after application of
the X-matrix. Raw Stokes products of FPP are not consistent with this definition.
-Q
-U
+U
N
-V
-Q
-U
FPP
+U
W
+Q
E
+Q
+V
-V
S/C +Y
+V
S
S/C +X
View from the top of SOT
View towards the sun
This definition is the same as that used in the analysis of the suntest data of 2004.8 and consistent
with the ASP definition, ie. positive V at blue side of spectral line gives positive magnetic flux. This
is also consistent with the definition of Stokes V: (right circ. – left circ.), where right circular
polarization is positive when electric vector rotates clockwise looking at the source.
簡易マップ作成; simple_iquvmap.pro (一本製)
; smap:[*,*,k]
'continuum', $
'Q', $
'U', $
'V', $
;
'NCP (6302)', $
'Dopp. (CG)', $
'Dopp. (6301 botm)', $
'Dopp. (6301 wing 0.8)', $
'I (6301 core)', $
'total rss(circluar)', $
'total rss(linear)', $
'Dopp. (6302 I+V botm)', $
'Dopp. (6302 I+V wing 0.8)', $
'Dopp. (6302 I-V botm)', $
'Dopp. (6302 I-V wing 0.8)', $
'NCP (6301)', $
'V (6302+0.365)', $
'V (6302-0.276)', $
'V (6301)' $
]
dir=‘level-1 ディレクトリ\'
vname=[
spos=0
0  全体 1 scan とみなす
20  1scan 20 slit pos の繰り返しと見なす
様々な量の2次元マップsmap[*,*,*] をつくり
dir\smap.isv に保存。
例えば Q-map は下の波長領域を積分
l0 -10 ~ l0 -2
l0 +2 ~ l0 +10
l0 =75
k
;0
;1
;2
;3
;4
;5
;6
;7
;8
;9
; 10
; 11
; 12
; 13
; 14
; 15
; 16
; 17
; 18
磁場へ変換 Stokes インバージョン
HSC ME インバージョンコード by Yokoyama 他
sp_meinv.pro
(一本製)
;------ data(*,*,i) -------------------------------------;
0
Field strength
Gauss
;
1
Field inclination
degree
;-----------------------------------------;
2
Field azimuth
degree
date='20071015'
;
3
Doppler shift 630.15nm km/sec
;-----------------------------------------;
4
Doppler shift 630.25nm km/sec
dir='/nfs/scr01/ichimoto/spcal/'+date
5
Doppler width
mA
invdir='/nfs/scr01/ichimoto/spcal/meinv' ;
;
6
Line strength
non-dim
outdir=invdir
;
7
Damping parameter
non-dim
;
8
Source function
non-dim
;
9
Source function gradient non-dim
;
10
Macro-turbulence
km/sec
;
11
Straylight fraction
non-dim
;
12
Straylight shift
mA
;
13
Continuum intensity
DN
;
14
Polarization degree
マップ表示 show_maps.pro
(一本製)

similar documents