Could use a more sophisticated method for the contamination
searches, perhaps based on source flux ratios and errors
Could have a list of source numbers as input and output a set of
region files based on those sources. The output region files could be
based on the input root plus the source number. The problem with this
is that it would contradict the pipeline philosophy of always
specifying all input and output parameters.
omdetect and omwavelet will produce source lists with
ellipsiodal shapes for objects. This task should be able to produce
ellipses and elliptical annulli for the extraction regions.