VaST [Variability Search Toolkit]A software package for detection of variable objects
| ||||||||||
VaST is a program aimed to find objects of variable brightness on a series of astronomical images. The images (CCD frames or photographic plate scans) must be taken with the same instrument using the same filter and saved in FITS format. The images may be shifted and/or rotated with respect to each other, but they have to have the same scale (arcsec/pix). VaST performs object detection and aperture photometry using SExtractor on each image, cross-matches lists of detected stars, performs magnitude calibration with respect to the first (reference) image and constructs a lightcurve for each object. The simple sigma-magnitude plot may be used to visually identify variable star candidates. If data permit, period-search techniques ANOVA and BLS may also be employed to identify periodic variable stars among the candidates. All lightcurves are being constructed in an arbitrary magnitudes. If there are reference stars with known magnitudes in the field of view, all magnitudes can be later converted to the standard system. Unlike software based on the image subtraction method, VaST can be used in case of unstable PSF (e.g., bad guiding, digitized wide-field photographic images). VaST is not tied up to any catalog and WCS, so it can be successfully used on images taken with a telephoto lens as well as with a 2.6-m class telescope. VaST is written in C (and partly in BASH scripting language) for GNU/Linux operating system. It is in a continual, albeit slow, state of development. VaST is free software: you can redistribute it and/or modify it under the terms of the GNU General Public License. To compile and use VaST you'll need:
You can find at this page the latest version of VaST
(download via
ftp
or
http, ~1.6M)
along with some sample data
(ftp
or http, 24M)
Special note for Ubuntu users! The following installation procedure is recommended:
ScreenshotsDifferent subroutines in VaST may be started automatically from the main program, or they may be started manually by user from a terminal. Names of the corresponding executable files are given under each screenshot. Click on images to see them in high resolution:
Video (screencast) background comments in Russian only, sorry... I hope to prepare an English version soon...
Compiling VaSTFirst, make sure that you have PGPLOT, CFITSIO and GSL libraries installed in your system. Then unpack the program
tar -xvjf vast-latest.tar.bz2
And try to compile it by running
cd vast-1.0rc48
make
If compilation fails, read carefully compilers output. Most probably it
can't find some necessary libraries. Explore the Makefile with your favorite
text editor and change the libraries location to match your environment.
Using VaSTTo try VaST, download the sample data set, unpack it
wget -c ftp://saistud.sai.msu.ru/pub/vast/sample_data.tar.bz2
change to the VaST installation directory and run the program
tar -xvjf sample_data.tar.bz2
cd vast-1.0rc48
or
./vast ../sample_data/*fit
cd vast-1.0rc48
After a brief computation the variability search window will be opened. Click
on any star on the sigma-magnitude plot to inspect its lightcurve. Click on
any point on the lightcurve plot to see an image from which this point
comes. The star on the image will be marked with a red cross, a red circle
around the star corresponds to the aperture used to measure this image.
./vast -p ../sample_data/*fit If you quit the program without deleting any data, you can restart it without doing all the computations again by running:
./find_candidates aa
To calibrate the magnitude scale using comparison stars in the field of view use this script:
util/magnitude_calibration.sh
It will ask you to specify one or more comparison stars and their
magnitudes. After the comparison stars have been specified, close the image
window with the right mouse click and inspect the relation between
instrumental and catalog magnitudes. This relation may be fitted by a
linear function with the slope fixed to 1, linear function with a free slope
or second degree polynome. Use (P) key on your keyboard to change the fitting
function. For most CCD data sets, linear function with the fixed slope
provides a good fit. The data points may be weighted according to their
estimated errors or weighted equally. Use (W) key to change weighting.
Close the dialog with the right mouse click to apply the calibration to all
stars when you are satisfied with the fit.In order to save the current work, use the script
util/save.sh my_favorite_field_name
If the field name was not provided, the script will ask for it. To restore previous
work, use
util/load.sh my_favorite_field_name
or just copy all saved files from directory [my_favorite_field_name] back to the VaST directory.To delete all files from the previous work, you can use:
util/clean_data all
New feature! There is a new tool to fully automatically identify variable stars detected with VaST. It requires working Internet connection at runtime. The tool may or may not work for your field depending on luck. If the field of view is greater than, say, 30' and the image quality is fine - the odds that the automatic identification will work for your field are pretty good. Here is how to use it... To identify the star number 00190 in the current dataset use the following command:
util/identify.sh out00190.dat 25
where "out00190.dat" is the lightcurve file of the star and "25" is the
estimated field of view size in arcminutes. It is wise to set this value a little bit smaller
than your best-guess field of view. If no field of value is given, the
default value of 40' will be assumed. Note, that the field of view
identification may take really long time (up to ~30 minutes), especially for small
fields. Please, do not interrupt the process - just wait for a result. If
the field identification failed, try to restart it with smaller estimated
field of view size. On successful identification, the script will print the
equatorial coordinates (J2000) of the star and will attempt to query popular
databases to check if this star is already known as a variable. Note, that
the coordinates returned by the identification script may be up to 10" off.
We strongly advise to use the Aladin
interactive sky atlas to verify the identification and obtain more
precise star coordinates from some astrometric catalog (e.g., USNO-B1.0).
The automatic identification function relies on tools provided by Astrometry.net. For more details visit
http://astrometry.net/. Note, that
according to the
use conditions of the index files required to run the Astrometry.net
code, you need to be ready to share the images solved using this software
with any other interested party including the Astrometry.net team. According
to the same use
conditions: "If the indexes are used in any scholarly work,
research grant proposals or other scientific or academic
publications, the Astrometry.net project must be properly cited
and acknowledged. The paper to cite is the Lang et al paper on our
bibliography page http://astrometry.net/biblio.html.
The acknowledgment text is "This research made use of tools provided
by Astrometry.net." The VaST software developers assume, that this statement
concerns the work which uses "util/identify.sh" script from the VaST package.
Using VaST with scanned photographic platesIf you have to deal with digitized photographic materials, you'll want first to convert TIFF images (which are typically produced by scanner software) to the FITS format. It can be easily done with the tiff2fits converter. Do not forget to change your images to positive (white stars on black sky). For example:
./tiff2fits -i input.tiff output.fits
Also, you'll have to enter the information about the date and time of
observation in fits header. This can be done, among other ways, using the
modhead program from the CFITSIO
examples page, MissFITS or edhead from the WCSTools package. You may
enter the observation date and time information in the image header by generating
the usual "DATE-OBS" and "EXPTIME" header keys, or by simply putting the
middle-of-exposure Julian Date into the "JD" key.
You will most likely need to modify the default SExtractor settings stored
in the "default.sex" file located in the VaST directory. Please use
"default.sex.beta_Cas_photoplates" and "default.sex.PHOTO" as working examples
(these files are also found in the VaST directory).Say, you have digitized photographic images in "photo_data" directory. To run VaST on them, change to the VaST directory and run the program
cd vast-1.0rc54
The "-o" parameter tells VaST to use formulas (1) and (3) from
Bacher et al. (2005, MNRAS, 362,
542) for magnitude calibration. This technique provide slightly better
results for the aperture photometry on photographic data compared to the
second-order polynome used in VaST by default for magnitude calibration.
However, the standard calibration technique should also provide reasonably
good results (see Figure 2 in this paper).
./vast -o ../photo_data/*fit VaST mailing listIf you would like to receive updates about new versions of VaST, enter your e-mail address in the form below or just drop a line to kirx[at]scan.sai.msu.ruSome switches which you can use with VaST"-9" — use DS9 instead of VaST's own viewer pgfv to view FITS files."-n" — do NOT perform magnitude calibration. Useless for variable star search but good for doing differential photometry with util/differential_photometry.sh script. "-u" — always assume UTC time system, don't perform conversion to TT. "-e" — failsafe mode: only stars detected on the reference frame will be processed. "-p" — tell VaST that it deals with linear detector. This flag is useful for well calibrated CCD images, but for some CCD images it is better not to use it. Test what works best for your particular dataset! "-l" — perform sigma-filtering of lightcurves (can be extremely CPU-intensive) "-r" — assume that images are not rotated with respect to one another by more than 3 degrees. To get the full list of the available command line arguments use:
./vast -h
List of publications which make use of VaSTIf you are aware of a publication which is not listed here, please, drop us a line to kirx[at]scan.sai.msu.ru so we can include it in the list!
A few other (more or less) useful toolsWeb-based World Coordinate System (WCS) image calibration - a tool based on the Astrometry.net engine that allows one to plate-solve a sky image using minimal a priori information about it. Please use the username vastuser and password vastpassword to access this page.Web-based lightcurve analysis tool - allows one to perform a period search using Lafler & Kinman (1965, ApJS, 11, 216) method. Plase use the username vastuser and password vastpassword to access this page. Times of minima calculator - a simple web-based tool to calculate times of minima from light elements. TIFF2FITS - a simple C program (based on CFITSIO and libtiff) which converts TIFF images to 16-bit FITS format. May be useful for processing scanned photographic plates. OMC2ASCII - a simple C program (based on CFITSIO) which converts INTEGRAL/OMC lightcurves from FITS table to simple ASCII table. May be useful if you don't want to install the Off-line Science Analysis (OSA) software just to get to INTEGRAL/OMC data. If you have any questions fell free to contact authors via e-mail kirx[at]scan.sai.msu.ru
Good luck in using VaST! VaST developers team. | ||||||||||
|
VaST developers team, 2005-2012
|