Dynamic Spectrum and Imaging: Difference between revisions

From EOVSA Wiki
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. You can select a subset of data by specifying a [https://casa.nrao.edu/Release3.3.0/docs/UserMan/UserMansu112.html time range], [https://casaguides.nrao.edu/index.php/Selecting_Spectral_Windows_and_Channels spectral windows/channels],  [https://casaguides.nrao.edu/index.php/Antenna/Baseline_Selection_Syntax_with_or_without_Autocorrelations antenna baseline].


<pre>
<pre>

Revision as of 04:35, 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. You can select a subset of data by specifying a time range, spectral windows/channels, antenna baseline.

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')

Imaging with SunCASA