Forced photometry task using the LSST stack
Scripts to run forced photometry at user-specified positions. Based on this Twinkles photometry code by Michael Wood-Vasey. See this post for discussion relating the the development of the original code.
The routine is designed to measure forced photometry on a set of images for Stripe 82 data in the Prototype Data Access center.
Example:
get_forcedphot.py /datasets/gapon/data/DC_2013/calexps ./output id1340590_images.json \
'1,25.612807,0.291621' u id1340590_photometry_u.tbl
produces a table in IPAC format id1340590_photometry.tbl at the specified RA=25.612807,
Dec=0.291621, for all the images in id1340590_images.json. The json input file of images
overlapping this point is separately generated by a DAX query.