Produce a stacked image of all epochs for a VAST field, optionally applying positional and flux corrections. Run with --help
for more details.
Requires SWarp which is also available via conda-forge. Or, on mortimer
:
module use /sharedapps/LS/cgca/modulefiles/
module load swarp
The following will produce a stacked image of the VAST_0012+00A field, assuming that the images are stored following the VAST data release structure.
A spreadsheet of the VAST pilot corrections in Excel format is required to apply the positional and flux corrections to each image before stacking. This can be obtained by download a copy of the corrections Google Sheet in Excel format, see the VAST Wiki for details.
#!/bin/bash
module use /sharedapps/LS/cgca/modulefiles/
module load swarp
conda activate vast-tools
python stack_images.py \
--out output_dir \
--clean \
--type=tiles \
--imagepath="/raid-17/LS/kaplan/VAST/EPOCH*/TILES/STOKESI_IMAGES/" \
--qc "VAST Pilot Corrections.xlsx" \
--beam \
"0012+00A"