上篇介绍了使用IDL mir对SMA偏振数据进行数据检查，apply系统温度, bandpass calibration, fix bandpass delay, gain calibration, cross-rx gain calibration的过程。本文介绍使用miriad进行instrumental polarization correction和成图的过程。（这一部分就没有官方manual了，我是根据张其洲老师的script总结的。）

一． 准备

上文中，我们最后得到了pol calibrator和source的UV数据：3C279_lsb，3C279_usb， G28_lsb, G28_usb.

‘

在c-shell中设置变量（b-shell中命令格式略有变化）

set bpcal='3C279' %pol calibrator

set source='G28'

set refant=2 %使用天线2做基准天线

set avetime=2 %每两分钟平均

set bpflux=7.1 %pol calibrator的flux

二． 使用pol calibrator解出pol solution。 （以lsb为例，usb照此办理）

1.     再一次做更精细的gain calibration。（我感觉这一步影响不是很大）

mfcal vis=$bpcal'_lsb' refant=${refant} options=nopassol interval=2

2.     对数据以两分钟为间隔做平均。（这步必须做）

puthd in=$bpcal'_lsb'/evector value=0.7853 %固定步骤，改正偏振角初始角度

puthd in=$bpcal'_lsb'/mount value=4 %固定步骤， 我忘了这步干啥的了

uvaver vis=$bpcal'_lsb' out=$bpcal'_lsb_ave' interval=${avetime} %平均并做某些修正

uvsort vis=$bpcal'_lsb_ave' out=$bpcal'_lsb_ave_sort' %将数据按编号排列好

3.     重新计算视差信息。（这一步不是必需的，如果视差信息已经正确地被存储到了数据里的话）

uvredo vis=$bpcal'_lsb_ave_sort' options=chi out=$bpcal'_lsb_ave_sort_redo'

4.     使用pol calibrator解出instrumental polarization.

blflag vis=$bpcal'_lsb_ave_sort_redo' axis=real,imag options=nobase device=/xs %flag掉坏点

gpcal vis=$bpcal'_lsb_ave_sort_redo' refant=${refant} options=noxy,circular,qusolve interval=${avetime} flux=$bpflux %解出solution

5.     检查pol calibration.

uvplt vis=$bpcal'_lsb_ave_sort_redo' device=2/xs axis=real,imag options= nocal,nopol,noba average=1 stokes=rl,lr %calibrate之前

uvplt vis=$bpcal'_lsb_ave_sort_redo' device=2/xs axis=real,imag options=noba average=1 stokes=rl,lr %calibrate之后，图形的中心与原点重合。

三． 分离源的连续谱和谱线

1.     检查谱线

uvspec vis=G28_lsb stokes=i axis=chan,ampl device=/xs options=avall,nobase interval=10000 line=channel yrange=-1,5 %

2.     分离

uvlin vis=$source'_lsb' out=$source'_lsb_cont' order=0 mode=chan0 options=nowin,nopol chans=200,2250,2400,7850,7950,7990

uvlin vis=$source'_lsb' out=$source'_lsb_line' order=1 mode=line options=nowin,nopol chans=200,2250,2400,7850,7950,7990

四． 将polarization solution 应用到source上

1.     不解释了。固定步骤，照做即可。

puthd in=$source'_lsb_cont'/evector value=0.7853

puthd in=$source'_lsb_cont'/mount value=4

uvsort vis=$source'_lsb_cont' out=$source'_lsb_cont_uvsort'

uvaver vis=$source'_lsb_cont_uvsort' out=$source'_lsb_cont_uvsort_ave' interval=${avetime}

gpcopy vis=$bpcal'_lsb_ave_sort_redo' out=$source'_lsb_cont_uvsort_ave' options=nocal %将pol solution 复制到源中

uvredo vis=$source'_lsb_cont_uvsort_ave' options=chi out=$source'_lsb.pol.out'

五． 成图

set date='merge0203'

1.     将数据从uv空间进行傅里叶变换转化成dirty image

invert vis=$source'_usb.pol.out',$source'_lsb.pol.out' map=$source'_cont'$date.i.mp,$source'_cont'$date.q.mp,$source'_cont'$date.u.mp beam=$source'_cont'$date.bm stokes=i,q,u options=mfs line=chan,1,1 robust=2 cell=0.3 imsize=256

2.     修正第三维

puthd in=$source'_cont'$date.i.mp/cdelt3 value=1.0 type=double

puthd in=$source'_cont'$date.q.mp/cdelt3 value=1.0 type=double

puthd in=$source'_cont'$date.u.mp/cdelt3 value=1.0 type=double

3.     clean

set clcuti=0.005 %cutoff 一般设置为3-5倍rms。可通过对dirty image成图估计rms

clean map=$source'_cont'$date.i.mp, beam=$source'_cont'$date.bm out=$source'_cont'$date.i.clmp niters=1000 gain=0.1 cutoff=$clcuti

4.     restore

restor model=$source'_cont'$date.i.clmp beam=$source'_cont'$date.bm map=$source'_cont'$date.i.mp out=$source'_cont'$date.i.cvmp

六． 生成偏振文件

set qurms =  %Q/U map的rms

set irms=  %I map的rms

impol in=$source'_cont'$date.q.mp,$source'_cont'$date.u.mp,$source'_cont'$date.i.cvmp poli=$source'_cont'$date.poli,$source'_cont'$date.polierr polm=$source'_cont'$date.polm,$source'_cont'$date.polmerr pa=$source'_cont'$date.pa,$source'_cont'$date.paerr sigma=$qurms,$irms sncut=2,5