Dynamic Spectrum and Imaging: Difference between revisions

From EOVSA Wiki
Jump to navigation Jump to search
Line 3: Line 3:
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]. The selection syntax follows the ''CASA'' convention. More information may be found in the CASA guide of  [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] selection pages.
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]. The selection syntax follows the ''CASA'' convention. More information may be found in the CASA guide of  [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] selection pages.


<pre style="background-color: #d8bfd8;">
<source lang="python">
from suncasa.utils import dspec as ds
from suncasa.utils import dspec as ds


Line 32: Line 32:
## this plot the dynamic spectrum.
## this plot the dynamic spectrum.
ds.plt_dspec(specfile, pol='XXYY')
ds.plt_dspec(specfile, pol='XXYY')
</pre>
</source>
 
 
[[file:fig-dspec.png|thumb|center|500px|Figure 1: EOVSA cross power dynamic spectrum at stokes XX and YY]]
[[file:fig-dspec.png|thumb|center|500px|Figure 1: EOVSA cross power dynamic spectrum at stokes XX and YY]]



Revision as of 15:05, 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. The selection syntax follows the CASA convention. More information may be found in the CASA guide of time range, spectral windows/channels, antenna baseline selection pages.

<source lang="python"> from suncasa.utils import dspec as ds

  1. define the visbility data file

msfile = 'IDB20170821201020-203020.12s.slfcaled.ms'

    1. define the output filename of the dynamic spectrum

specfile = msfile + '.dspec.npz'

    1. antenna selection
    2. leave it blank for selecting all baselines

bl =

    1. Spectral Windows and Channels selection

leave it blank for selecting all spectral windows spw =

    1. time range selection
    2. leave it blank for selecting the entire time interval

timeran =

    1. select baselines with a median length (0.2~0.8km)

domedian =True

    1. this step generates a dynamic spectrum and saves it to specfile

ds.get_dspec(vis=msfile, specfile=specfile, bl=bl, spw=spw, domedian=domedian)

    1. this plot the dynamic spectrum.

ds.plt_dspec(specfile, pol='XXYY') </source>


Figure 1: EOVSA cross power dynamic spectrum at stokes XX and YY

Imaging with SunCASA

cd to your working directory where the measurement sets file is located.

from suncasa.utils import qlookplot as ql
msfile = 'IDB20170821201020-203020.12s.slfcaled.ms'
vis = msfile
timerange = '20:21:10~20:21:30'  ## time range selection
spw = '3.4~6.0GHz'  ## Spectral Windows and Channels selection
stokes = 'XXYY'     ## polarizations selection
ql.qlookplot(vis, timerange=timerange, spw=spw, stokes=stokes)