AIPS を用いた
VERA データの解析方法
2011 年 9 月 28
28～～29 日
水沢 VLBI 観測所 ユユーザーズミーティング
ザ ズミ ティング
＠国立天文台三鷹
倉山 智春 (鹿児島大学鹿児島大学))
2011/09/28--29
2011/09/28

水沢VLBI水沢VLBI観測所
観測所 UM1

Target(narrow-band data)
FITLD
(1)

論文に…
論文に…

Re

INDXR [+ FITLD CL1](3)

ACCOR [+ FITLD SN1]

(4)CLCAL [+ FITLD CL2]

(6)

ANTAB [+ FITLD GC1]

(7)APCAL [+ FITLD SN2](8)CLCAL [+ FITLD CL3]

TBIN

(22)

FRING

(23)

[+ AVER SN5] SNEDT(24)[+ AVER SN4]
docalib = 2; gainuse = 2; doband = 0
p = 20, 0
flux = 0; nmaps = 1; refant = 1
solint = 4; solsub = 10; aparm = 0, 0, 0, 0, 1, 2, 3, 0
dparm = 1, 10, 30, 0; snver = 4; invers = 1

yver = 0; gcver = 1

editting large delay data
(25)

TBIN

frequency conversion[+ FITLD SN5]

frequency conversion

(26)

TBIN

[+ FITLD SN6]

docalib = 1; solint = 1; bpassprm = 1, 0;
outvers = 1; smooth = 1
1, 0;
timerang = (8 min. at the highest elevation)

(27)

TBIN

[+ FITLD SN7](28)

CLCAL [+ FITLD CL5, CL6, CL7]
CLCALp
CLCAL snver = 5; gainver = 4; gainuse = 5snver = 6; gainver = 5; gainuse = 6snver = 7; gainver = 6; gainuse = 7(29)

TABED [rewriting FITLD AN]

(30)

SETJY

[+ FITLD SN3]
reading the recalculation file

SNEDT [+ FITLD SN4]

(12)

CLCAL [+ FITLD CL4]

editting large delay values

kkeyp
snver = 4; gainver = 3; gainuse = 4[+ AVSPC whole file]
gainuse = 4
doband = 2; bpver = 1; smooth = 0

AVSPC

(13)

avopt li[+ AVER whole file] UVAVG
1; doacor = 1
1 yinc1;
i =44; zinci =1(14)FRING(15)[+ AVER SN1]
1
doband = 1; refant = 1; solint = 2; solsub = 0aparm = 0, 0, 0, 0, 1, 2, 4, 0; dparm = 0; snver = 1[+ AVER SN2] SNEDT

CLCAL

(17)

[+ AVER SN3]

FRING

(18)

CLCAL

(19)

SPLIT

(20)

ystringt i[rewriting FITLD SU]

(31)

CVEL

(32)

INDXR

(33)

IMAGR [+ image file][+ CVEL whole file]doband = 2; bpver = 1; gainuse = 0docalib = 2; gainuse = 7; doband = 1;
bchan, echan = (three channels centered at peak)cellsize = 1e 4, 1e 4; imsize = 512, 512
rashift, decshift = (values in Table 2)niter = 10000(16)

psnver = 2; gainver = 1; gainuse = 2

TBOUT

sysvel = 57.3; restfreq = 2.223E10, 5080000
veltypparm = 256, 0

editing large delay data[+ AVER CL2]

TBOUT

inext = SN ; invers = 2[+ FITLD BP1]

(11)

( )
(34)

JMFIT
blc, trc = (range of maser radiation, one frequency channel)
ngauss = 1; domax = 1; dopos = 1; dowidth = 1
niter = 40

docalib = 2; gainuse = 2
doband = 1; refant = 1; solint = 4; solsub = 0
aparm = 0, 0, 0, 0, 1, 2, 3, 0
dparm = 1, 10, 30, 0; snver = 3[+ AVER CL3]psnver = 3; gainver = 2; gainuse = 3

2011/09/28-29

水沢VLBI観測所 UM
doband =(21)FITLD[+ REFMAP whole file]interpolcalsour
snver = 2; gainver = 2; gainuse = 3(10)(wide-band data)reading self-calibrated images

yver = 1; gcver = 1
snver = 2; aparm = 0; solint = 0BPASS

Re

self-calibration andimaging with DIFMAP

interpolcalsour
snver = 1; gainver = 1; gainuse = 2(9)(continued from lower left)parm = 0, 0, 0.02, 0Same asa target sourcee

PASJ 63 513

Target(narrow-band data)

FITTP

(2)

(5) EXTDEST [ FITLD all GCs]
; inext =inclass =► Kurayama+ 2011

(wide-band data)[+ REF whole file]inclass AVER
gainuse = 3parm = 2, 0(continuing to upper right)2位相補償画像
較r正opple
相 償Do
位相補償償
位位償
画画像合成成画成Self--calibraSelf
c 償 ation
位再計算uvvw
バンドパス較正)天r体(mmaser

目的天体

振振幅較正正振ド ス
参照天体連(続波天体体体波 体)3
水沢VLBI観測所 UM
2011/09/28-29

画 フリ成ンジサーチ

全体的な解析の流れ較正の大分類► 観測量・ビジビリティは
観ビジビ

振幅
周波数

 複素数で(位相と振幅を持つ位相と振幅を持つ))
 時間変化し 周波数にも依存する。
► これに対応して、各較正もれに対応して、各較正も分類できる。

位相周波数

 位相較正 振幅較正 バンドパス較正
2011/09/28-29

時刻

水沢VLBI観測所 UM

時刻

4

振幅較正
► チョッパチョッパーホイール法
ホイ ル法► 各スキャンの開始前に電波吸収体をかぶせる► FITS ファイル内に必要なデータを格納
Scan
Scan
R (1天体組を観測R (1(1天体組を観測))(1天体組を観測天体組を観測))

時刻
2011/09/28-29

水沢VLBI観測所 UM

5

位相較正
► FITS 用 uvw 再計算

 相関器での追尾計算の修正
相関器での追尾計算の修正・高精度化
高精度化► フリンジサーチ► SelfSelflf-calibrationlb► 位相補償 2 天体間の位相の引き算
 機械的位相差の較正
2011/09/28-29

水沢VLBI観測所 UM6

FITS 用 uvw 再計算
局位置について► ソフトウェアのバグ
ソフトウ アのバグ 2011 年の観測に対して、でたらめな局位置で計算を行っていた► 地震の影響 水沢近辺ではここ数年地震が多発►2008.6.14
2008 6 14 岩手・宮城内陸地震
岩手 宮城内陸地震►2011.3.11 東日本大震災►前震・余震前震 余震

 正確な局位置を求めるのは容易ではない
2011/09/28-29

水沢VLBI観測所 UM7

FITS 用 uvw 再計算
EOP 改訂に伴う変更► 地球回転パラメータ (EOP(EOP、、Earth OrientationParameter)) の提供方式が変更された
提供方変更► 座標系の変更► 対応して局位置情報も変更が必要► 2011 年 7 月以降の観測は再計算ができない状態
 現在対応中2011/09/28-29

水沢VLBI観測所 UM8

位相補償► 位相の引き算► 周波数の変換を考慮する必要あり
周波数 変換を考慮する必要あり

 AIPS では基準 (≒先頭先頭)) 周波数の位相 + 傾きで補正
 変換プログラムまたは AIPS task の組み合わせで対応

連続波天体の帯域連続波天体の位相位相メーザー天体の帯域

周波数
2011/09/28-29

水沢VLBI観測所 UM9

メーザー探し
Fringe Rate Map► 細かいところが見える
細か と ろが見える

画像から広い範囲のメーザーを探すのは困難► Fringe Rate Map を作成すると容易に発見できることもある。

2011/09/28-29

水沢VLBI観測所 UM10