La Palma SVSTデータの解析      rev.0.2  Katsukawa (2004/12/20)

(1)ファイルの読み込み
○svt30may98_gband.00100.finを読む場合

IDL> ftemp = '/mwork2/katskwyk/lapalma/30may98/final/svt30may98_gband.######.fin'
IDL> data = anafrd(fns(ftemp,100),h)
IDL> help,data, h
DATA            INT       = Array[2028, 2044]
H               STRING    = '30-May-1998 08:19:49.871    101  2029 2047  16'...

一度に多くのファイルを読み込むとすぐにメモリが一杯になるので注意


(2)ヘッダーに含まれている情報
La Palmaデータのヘッダーにはろくな情報が含まれていない

IDL> print, h
30-May-1998 08:19:24.108    100  2029 2047  16  1100241

rang = 153.366

[時刻]
IDL> time = token(h,1)+' '+token(h,2)
IDL> print, time
30-May-1998 08:19:49.871

[exposur(多分)]
IDL> dur = token(h, 6)
IDL> print, dur
16
(16msだと思う)

[pointing, pixel scale]
自分で調べる必要あり
30-May-98データについては 0.082arcsec/pixel
(ref. Berger et al. 2001)


(3)Mapへ変換
IDL> map = mk_lpmap(data, time=time, dx=0.082, dy=0.082, wave='G-band' $
                    , dur=dur)
IDL> help,map,/st
** Structure <8418f2c>, 12 tags, length=8290544, data length=8290544, refs=1:
   DATA            INT       Array[2028, 2044]
   XC              FLOAT           0.00000
   YC              FLOAT           0.00000
   DX              FLOAT         0.0820000
   DY              FLOAT         0.0820000
   TIME            STRING    '14-Dec-2004 13:51:29.568'
   ID              STRING    'La Palma G-band'
   ROLL_ANGLE      FLOAT           0.00000
   ROLL_CENTER     FLOAT     Array[2]
   DUR             FLOAT           0.00000
   XUNITS          STRING    'arcsecs'
   YUNITS          STRING    'arcsecs'


(4)LCTで速度場を計算 
全視野では決して行なわない、部分像を切り出して行なうこと
○3D cubeを作成
IDL> data = intarr(256,256,20)
IDL> for i=0,19 do data[*,*,i]= (anafrd(fns(ftemp,i+100),h))[768:1023,512:767]

○LCTで速度場を計算
grid間隔が8 pixl, apodization window(ガウス関数)のFWHMは8pixel*2
IDL> tile  = 8
IDL> fwfac = 2
IDL> vel = flowmap(data, tile, fwfac, gx=gx, gy=gy)
IDL> help, vel, gx, gy
VEL             FLOAT     = Array[2, 32, 32, 20]
GX              LONG      = Array[32, 32]
GY              LONG      = Array[32, 32]
gxはgridのX座標を、gyはY座標をそれぞれ表す

○20枚(実際は19枚)の平均
IDL> va = flowaverage(vel)
IDL> help, va
VA              FLOAT     = Array[2, 32, 32]

○速度場をプロットする
IDL> plot_image, data[*,*,0]
IDL> plot_flow,va[0,*,*],va[1,*,*],gx[*,0],gy[0,*],/over