Tohban EOVSA Imaging Tutorial A-Z
Jump to navigation
Jump to search
Step 1: Downloading raw data (IDB) on pipeline server
On pipeline, CASA
from astropy.time import Time import os trange = Time(['2017-08-21 20:15:00', '2017-08-21 20:25:00']) #### (Optional) change output path, default current directory "./" ##### outpath = './msdata/' if not os.path.exists(outpath): os.makedirs(outpath) ###################################################### msfiles = importeovsa(idbfiles=trange, ncpu=1[, visprefix=outpath])
OR
from suncasa.tasks import task_calibeovsa as calibeovsa from suncasa.tasks import task_importeovsa as timporteovsa from split_cli import split_cli as split import dump_tsys as dt from util import Time import numpy as np import os from glob import glob from eovsapy import util trange = Time(['2020-06-07 21:35:00', '2020-06-07 21:55:00']) idbdir = util.get_idbdir(trange[0]) info = dt.rd_fdb(trange[0]) sidx = np.where( np.logical_and(info['SOURCEID'] == 'Sun', info['PROJECTID'] == 'NormalObserving') & np.logical_and( info['ST_TS'].astype(np.float) >= trange[0].lv, info['ST_TS'].astype(np.float) <= trange[ 1].lv)) filelist = info['FILE'][sidx] outpath = './msdata/' if not os.path.exists(outpath): os.makedirs(outpath) inpath = idbdir + '{}/'.format(trange[0].datetime.strftime("%Y%m%d")) ncpu = 1 msfiles = timporteovsa.importeovsa(idbfiles=[inpath + ll for ll in filelist], ncpu=ncpu, timebin="0s", width=1, visprefix=outpath, nocreatms=False, doconcat=False, modelms="", doscaling=False, keep_nsclms=False, udb_corr=True)
Step 2: Concatenate all the 10 mins data
# This is to set the path/name for the concatenated files concatvis = os.path.basename(msfiles[0])[:11] + '_concat.ms' vis = calibeovsa(msfiles, doconcat=True, concatvis=concatvis[, msoutdir=outpath])
Step 3: Calibration
calibeovsa(vis='IDB20170821202020.ms', caltype=['refpha','phacal'], doimage=True)
Step 4: Self-calibration on inti server
https://github.com/binchensun/casa-eovsa/blob/master/slfcal_example.py