Dynamic Spectrum and Imaging: Difference between revisions
Jump to navigation
Jump to search
| Line 1: | Line 1: | ||
==Get Dynamic Sepctrum with SunCASA== | ==Get Dynamic Sepctrum with SunCASA== | ||
Within SunCASA, you are using IPython to interact with the system. This does not mean extensive python experience is necessary. Basic Python interactions are straightforward, e.g., assigning parameters, importing modules, running functions. The first module we introduce is ''dspec''. This module allows you to generate a dynamic spectrum from an MS file, and visualize it. | Within SunCASA, you are using IPython to interact with the system. This does not mean extensive python experience is necessary. Basic Python interactions are straightforward, e.g., assigning parameters, importing modules, running functions. The first module we introduce is ''dspec''. This module allows you to generate a dynamic spectrum from an MS file, and visualize it. | ||
<pre> | <pre> | ||
from suncasa.utils import dspec | from suncasa.utils import dspec as ds | ||
bl=' | msfile = 'IDB20170821201020-203020.12s.slfcaled.ms' # the visbility data | ||
spw=' | bl = '' ## | ||
specfile = msfile+'. | spw = '' | ||
timeran=' | specfile = msfile + '.dspec.npz' | ||
timeran = '' | |||
domedian =True | |||
ds.get_dspec(vis=msfile, specfile=specfile, bl=bl, spw=spw, domedian=domedian) | |||
ds.plt_dspec(specfile, pol='XXYY') | |||
</pre> | </pre> | ||
==Imaging with SunCASA == | ==Imaging with SunCASA == | ||
Revision as of 04:21, 15 May 2019
Get Dynamic Sepctrum with SunCASA
Within SunCASA, you are using IPython to interact with the system. This does not mean extensive python experience is necessary. Basic Python interactions are straightforward, e.g., assigning parameters, importing modules, running functions. The first module we introduce is dspec. This module allows you to generate a dynamic spectrum from an MS file, and visualize it.
from suncasa.utils import dspec as ds msfile = 'IDB20170821201020-203020.12s.slfcaled.ms' # the visbility data bl = '' ## spw = '' specfile = msfile + '.dspec.npz' timeran = '' domedian =True ds.get_dspec(vis=msfile, specfile=specfile, bl=bl, spw=spw, domedian=domedian) ds.plt_dspec(specfile, pol='XXYY')