General Utilities¶

gatherDebugging.py

Description:	Script to gather information about the Pytho/numpy/LSL install to help with troubleshooting.
Usage:	gatherDebugging.py
Options:	None

makeWisdom.py

Description:	Build LSL-specific FFTW wisdom and save it to a file within the LSL distribution.
Usage:	makeWisdom.py
Options:	None

updateLSLSSMIF.py

Description:

Update the internal LWA1 SSMIF used by LSL.

Usage:

updateLSLSSMIF.py [options]

Options:

`-h, --help`	Display this help information
`-u, --update`	Update the default LWA1 SSMIF
`-r, --revert`	Revert the default LWA1 SSMIF to an older version

astroevents.py

Description:	Application to display rise, transit, and set times for various astronomical sources from LWA-1 for the current date.
Usage:	astroevents.py
Options:	None

astrostatus.py

Description:

Application to calculate real-time ephemeris for a LWA site.

Usage:

astrostatus.py [options]

Options:

`-h, --help`	show this help message and exit
`-s SITE, --site=SITE`
	site name (default LWA-1)

driftcurve.py

Description:

Generate a drift curve for a dipole at LWA-1 observing at a given frequency in MHz.

Usage:

driftcurve.py [OPTIONS]

Options:

`-h, --help`	Display this help information
`-f, --freq`	Frequency of the observations in MHz (default = 74 MHz)
`-p, --polarization`
	Polarization of the observations (NS or EW; default = EW)
`-l, --lf-map`	Use LF map instead of GSM
`-t, --time-step`
	Time step of simulations in minutes (default = 10)
`-x, --do-plot`	Plot the driftcurve data
`-v, --verbose`	Run driftcurve in vebose mode

getIonosphericRM.py

Description:

Estimate the ionospheric contribution to the RM for an observation using the IGS final product and the IGRF.

Usage:

getIonosphericRM.py [options] RA Dec Start Stop

RA: J2000 right ascension in HH:MM:SS[.SSS]

Dec: J2000 declination in sDD:MM:SS[.SSS]

Start: YYYY/MM/DD HH:MM:SS start time in UTC

Stop: YYYY/MM/DD HH:MM:SS stop time in UTC

Options:

`-h, --help`	Display this help information
`-n, --n-samples`
	Number of samples to take between the start and stop times (default = 11)

lwa_cat_view.py

Description:

Simple LWDA astronomical source catalogue display application.

Usage:

lwa_cat_view.py [options]

Options:

`-h, --help`	show this help message and exit
`-s SITE, --site=SITE`
	site name (default LWDA)
`-p PERIOD, --period=PERIOD`
	update period in seconds (default 5)

inspectTarball.py

Description:	Given a MCS metadata tarball, print out details of the associated observations.
Usage:	inspectTarball.py metaData
Options:	None

plotAntenna.py

Description:

Plot the relative dipole response for both polarizations of an isolated LWA-1 antenna at a particular frequency.

Usage:

plotAntenna.py [OPTIONS]

Options:

`-h, --help`	Display this help information
`-f, --freq`	Frequency of the observations in MHz (default = 74 MHz)
`-v, --verbose`	Run plotAntenna in vebose mode

plotStands.py

Description:

Plot the x, y, and z locations of stands at LWA-1. Also, mark and label particular stands, if requested.

Usage:

plotStands.py [OPTIONS] [stand1 [stand2 [...]]]

Options:

`-h, --help`	Display this help information
`-m, --metadata`	Name of SSMIF or metadata tarball file to use for mappings
`-l, --label`	Label the stands with their ID numbers (default = No)
`-v, --verbose`	Run plotStands in vebose mode
`-o, --output`	Filename to save the plot to (default = do not save)

plotUVCoverage.py

Description:

Randomly select 20 antennae from LWA-1 and plot the uv-plane coverage for a zenith snapshot and the expected beam. Alternatively, select some FRACTION of the stands with installed FEEs to use or use the specified list of stands.

Usage:

plotUVCoverage.py [FRACTION | STAND LIST]

Options:

`-h, --help`	Display this help information
`-f, --frequency`
	Frequency in MHz to compute the uv coverage (default 50 MHz)
`-m, --metadata`	Name of SSMIF or metadata tarball file to use for mappings
`-o, --output`	Filename to save the plot to (default = do not save)

Data Reading and Writing¶

splitTBN.py

Description:

Split a TBN file containing multiple seconds into several files

Usage:

splitTBN.py [options] file

Options:

`-h, --help`	Display this help information
`-c, --count`	Number of seconds to keep
`-o, --offset`	Number of seconds to skip before splitting
`-d, --date`	Label the split files with a date rather than a squence number

Note

This script does not use a lsl.reader.buffer buffer to try to re-order or verify all packets and simply splits files based on size.

splitDRX.py

Description:

Split a DRX file containing multiple seconds into several files

Usage:

splitDRX.py [options] file

Options:

`-h, --help`	Display this help information
`-c, --count`	Number of seconds to keep
`-o, --offset`	Number of seconds to skip before splitting
`-d, --date`	Label the split files with a date rather than a squence number

readTBW.py

Description:	Example script for reading in TBW data and writing it to a TSFITS file.
Usage:	readTBW.py file
Options:	None

readTBN.py

Description:	Example script for reading in TBN data and writing it to a TSFITS file.
Usage:	readTBN.py file
Options:	None

readDRX.py

Description:	Example script for reading in DRX data and writing it to a SD-FITS file.
Usage:	readDRX.py file
Options:	None

splitSession.py

Description:	Given a MCS metadata tarball and a session DRX recording, split the session recording into the individual observations.
Usage:	splitSession.py metaData data
Options:	None

plotMapper.py

Description:	Read and plot the NOSTA_MAPPER table in a FITS IDI file writen by `lsl.writer.fitsidi` if it exists.
Usage:	plotMapper.py file
Options:	None

Data Analysis¶

tbwSpectra.py

Description:

Given a TBW file, plot the time averaged spectra for each digitizer input.

Usage:

tbwSpectra.py [OPTIONS] file

Options:

`-h, --help`	Display this help information
`-m, --metadata`	Name of SSMIF or metadata tarball file to use for mappings
`-t, --bartlett`	Apply a Bartlett window to the data
`-b, --blackman`	Apply a Blackman window to the data
`-n, --hanning`	Apply a Hanning window to the data
`-q, --quiet`	Run tbwSpectra in silent mode
`-l, --fft-length`
	Set FFT length (default = 4096)
`-g, --gain-correct`
	Correct signals for the cable losses
`-s, --stack`	Stack spectra in groups of 6 (if ‘-g’ is enabled only)
`-d, --disable-chunks`
	Display plotting chunks in addition to the global average
`-o, --output`	Output file name for spectra imag

Warning

tbwSpectra.py currently assumed that the system it is running on has enough memory to read in a full TBW capture. Due to data representation and processing overheads this amounts to about 16 GB.

tbnSpectra.py

Description:

Given a TBN file, plot the time averaged spectra for each digitizer input.

Usage:

tbnSpectra.py [OPTIONS] file

Options:

`-h, --help`	Display this help information
`-m, --metadata`	Name of SSMIF or metadata tarball file to use for mappings
`-t, --bartlett`	Apply a Bartlett window to the data
`-b, --blackman`	Apply a Blackman window to the data
`-n, --hanning`	Apply a Hanning window to the data
`-s, --skip`	Skip the specified number of seconds at the beginning of the file (default = 0)
`-a, --average`	Number of seconds of data to average for spectra (default = 10)
`-q, --quiet`	Run tbwSpectra in silent mode
`-l, --fft-length`
	Set FFT length (default = 4096)
`-d, --disable-chunks`
	Display plotting chunks in addition to the global average
`-o, --output`	Output file name for spectra image

drxSpectra.py

Description:

Given a DRX file, plot the time averaged spectra for each beam output.

Usage:

drxSpectra.py [OPTIONS] file

Options:

`-h, --help`	Display this help information
`-t, --bartlett`	Apply a Bartlett window to the data
`-b, --blackman`	Apply a Blackman window to the data
`-n, --hanning`	Apply a Hanning window to the data
`-s, --skip`	Skip the specified number of seconds at the beginning of the file (default = 0)
`-a, --average`	Number of seconds of data to average for spectra (default = 10)
`-q, --quiet`	Run tbwSpectra in silent mode
`-l, --fft-length`
	Set FFT length (default = 4096)
`-d, --disable-chunks`
	Display plotting chunks in addition to the global average
`-o, --output`	Output file name for spectra image

drSpecSpectra.py

Description:

Given a DR spectrometer file, plot the time averaged spectra for each beam output.

Usage:

drSpecSpectra.py [OPTIONS] file

Options:

`-h, --help`	Display this help information
`-s, --skip`	Skip the specified number of seconds at the beginning of the file (default = 0)
`-a, --average`	Number of seconds of data to average for spectra (default = 10)
`-q, --quiet`	Run drSpecSpectra in silent mode
`-d, --disable-chunks`
	Display plotting chunks in addition to the global average
`-o, --output`	Output file name for spectra image

correlateTBW.py

Description:

Cross-correlate data in a TBW file

Usage:

correlateTBW.py [OPTIONS] file

Options:

`-h, --help`	Display this help information
`-m, --metadata`	Name of SSMIF or metadata tarball file to use for mappings
`-l, --fft-length`
	Set FFT length (default = 2048)
`-q, --quiet`	Run correlateTBW in silent mode
`-x, --xx`	Compute only the XX polarization product (default)
`-y, --yy`	Compute only the YY polarization product
`-2, --two-products`
	Compute both the XX and YY polarization products

correlateTBN.py

Description:

Example script that reads in TBN data and runs a cross-correlation on it. The results are saved in the Miriad UV format.

Usage:

correlateTBN.py [OPTIONS] file

Options:

`-h, --help`	Display this help information
`-m, --metadata`	Name of SSMIF or metadata tarball file to use for mappings
`-f, --fft-length`
	Set FFT length (default = 256)
`-t, --avg-time`	Window to average visibilities in time (seconds; default = 6 s)
`-s, --samples`	Number of average visibilities to generate (default = 10)
`-o, --offset`	Seconds to skip from the beginning of the file
`-q, --quiet`	Run correlateTBN in silent mode
`-x, --xx`	Compute only the XX polarization product (default)
`-y, --yy`	Compute only the YY polarization product
`-2, --two-products`
	Compute both the XX and YY polarization products
`-4, --four-products`
	Compute all for polariation products: XX, YY, XY, and YX.

possm.py

Description:	Script that takes a FITS IDI file and mimics the AIPS task POSSM by plotting average cross-power spectra for all baselines in the FITS IDI file.
Usage:	possm.py file
Options:	None

imageIDI.py

Description:

Script that takes a FITS IDI file and images the data.

Usage:

imageIDI.py file

Options:

`-h, --help`	Display this help information
`-1, --freq-start`
	First frequency to image in MHz (Default = 10 MHz)
`-2, --freq-stop`
	Last frequency to image in MHz (Default = 88 MHz)
`-s, --dataset`	Data set to image (Default = All)
`-m, --uv-min`	Minimun baseline uvw length to include (Default = 0 lambda at midpoint frequency)
`-n, --no-labels`
	Disable source and grid labels
`-g, --no-grid`	Disable the RA/Dec grid

General Utilities¶

Data Reading and Writing¶

Data Analysis¶

Table Of Contents

Previous topic

Next topic

This Page

Navigation

General Utilities¶

Data Reading and Writing¶

Data Analysis¶

Table Of Contents

Previous topic

Next topic

This Page

Quick search

Navigation