We implement an algorithm for detecting and removing artifacts from
astronomical images by means of outlier rejection during stacking. Our method
is capable of addressing both small, highly significant artifacts such as
cosmic rays and, by applying a filtering technique to generate single frame
masks, larger area but lower surface brightness features such as secondary
(ghost) images of bright stars. In contrast to the common method of building a
median stack, the clipped or outlier-filtered mean stacked point-spread
function (PSF) is a linear combination of the single frame PSFs as long as the
latter are moderately homogeneous, a property of great importance for weak
lensing shape measurement or model fitting photometry. In addition, it has
superior noise properties, allowing a significant reduction in exposure time
compared to median stacking. We make publicly available a modified version of
SWarp that implements clipped mean stacking and software to generate single
frame masks from the list of outlier pixels.Comment: PASP accepted; software for download at
http://www.usm.uni-muenchen.de/~dgruen
