Difference between revisions of "Dynamic Spectrum and Imaging"

From EOVSA Wiki
Jump to: navigation, search
(Get Dynamic Sepctrum with SunCASA)
(Get Dynamic Sepctrum with SunCASA)
Line 8: Line 8:
 
msfile = 'IDB20170821201020-203020.12s.slfcaled.ms'  # the visbility data
 
msfile = 'IDB20170821201020-203020.12s.slfcaled.ms'  # the visbility data
 
specfile = msfile + '.dspec.npz'  ## The output filename of the dynamic spectrum  
 
specfile = msfile + '.dspec.npz'  ## The output filename of the dynamic spectrum  
bl = ''                 ##  antenna selection
+
bl = ''   ##  antenna selection
spw = ''               ## Spectral Windows and Channels selection
+
spw = ''   ## Spectral Windows and Channels selection
timeran = ''         ## time range selection
+
timeran = ''   ## time range selection
domedian =True
+
domedian =True   ## select baselines with a median length (0.2~0.8km)
 
ds.get_dspec(vis=msfile, specfile=specfile, bl=bl, spw=spw, domedian=domedian)
 
ds.get_dspec(vis=msfile, specfile=specfile, bl=bl, spw=spw, domedian=domedian)
 
ds.plt_dspec(specfile, pol='XXYY')
 
ds.plt_dspec(specfile, pol='XXYY')

Revision as of 04:49, 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. More information of the selection syntax may be found in the time range, spectral windows/channels, antenna baseline pages.

from suncasa.utils import dspec as ds

msfile = 'IDB20170821201020-203020.12s.slfcaled.ms'  # the visbility data
specfile = msfile + '.dspec.npz'  ## The output filename of the dynamic spectrum 
bl = ''    ##  antenna selection
spw = ''    ## Spectral Windows and Channels selection
timeran = ''    ## time range selection
domedian =True    ## select baselines with a median length (0.2~0.8km)
ds.get_dspec(vis=msfile, specfile=specfile, bl=bl, spw=spw, domedian=domedian)
ds.plt_dspec(specfile, pol='XXYY')

Imaging with SunCASA