IPP Software Navigation Tools IPP Links Communication Pan-STARRS Links

Changes between Version 1 and Version 2 of Stack_Rejection_Discussion


Ignore:
Timestamp:
Jul 31, 2012, 4:38:41 PM (14 years ago)
Author:
watersc1
Comment:

--

Legend:

Unmodified
Added
Removed
Modified
  • Stack_Rejection_Discussion

    v1 v2  
    33== July 31, 2012 ==
    44
    5 Nigel sent around a jpeg of a mean stacked image, ignoring masking and ignoring variances.  This is easy enough to reconstruct, and I did so for this example skycell (skycell.1406.045, stack_id 1033436).  In addition, he overlaid the warp detections with psfQf cuts and stack detections on this image.  I've copied these images to this wiki to keep everything together:
     5Nigel sent around a jpeg of a mean stacked image (http://astro.dur.ac.uk/~nm/mean.jpg), ignoring masking and ignoring variances.  This is easy enough to reconstruct, and I did so for this example skycell (skycell.1406.045, stack_id 1033436).  In addition, he overlaid the warp detections (http://astro.dur.ac.uk/~nm/meandetections.jpg) with psfQf cuts and stack detections (http://astro.dur.ac.uk/~nm/meanstackdetections.jpg) on this image.  I've copied these images to this wiki to keep everything together (I've had to increase the JPEG compression to fit under the wiki filesize limit).
    66
    7 [[Image(mean.jpg,800px)]]
    8 [[Image(meandetections.jpg,800px)]]
    9 [[Image(meanstackdetections.jpg,800px)]]
     7[[Image(mean2.jpg,800px)]]
     8[[Image(meandetections2.jpg,800px)]]
     9[[Image(meanstackdetections2.jpg,800px)]]
    1010
    1111=== Comparison of mean stack to ppStack stack ===
     
    1515This image shows a mosaic of 20/22 inputs around the region of the reflection ghost visible in the mean image.  Only two of these inputs contain this feature, so although it is obvious in the mean stack, it is largely removed in the ppStack result.  We also have masks on this object, which helps rejection on most pixels, as shown by this image which is a mean combined, rejecting all input pixels that have any mask bit set.  This is somewhat more conservative than ppStack, which allows "suspect" pixels to be considered for inclusion if they are consistent with the values from other inputs.
    1616
    17 [[Image(mean_withmasks.png,800px)]]
     17[[Image(mean_withmasks.jpg,800px)]]
    1818
    19 [[Image(ppStack_mean_comp.png,800px)]]
     19[[Image(ppStack_mean_comp.jpg,800px)]]
    2020
    2121Comparing the unconvolved ppStack output (left) to the mean stacked image (right), it is clear that the majority of the non-astronomical sources are eliminated by ppStack.  This includes most of the reflection ghosts, the bright amplifier, the ragged detector edge, and the bad edge cells. 
     
    2929The detections from the meanstackdetections.jpg do seem to more consistently ignore the features visible on the mean image.  I've highlighted two examples where this is not the case:
    3030
    31 [[Image(ppStack_ghost_closeup.png,800px)]]
     31[[Image(ppStack_ghost_closeup.jpg,800px)]]
    3232
    3333ppStack largely rejects the pixels that make up the bright reflection ghost, but a small region shows this printing through somewhat (at the ~few sigma level it seems).  Comparing the shape of this region with the stack detections in this region shows that this bright feature is causing some false detections, but at a substantially reduced rate compared to the warp detections.
    3434
    35 [[Image(ppStack_BEC_closeup.png,800px)]]
     35[[Image(ppStack_BEC_closeup.jpg,800px)]]
    3636
    3737Although this bad edge chip has a mask, some pixels are masked with the "SUSPECT" bit, allowing ppStack to consider them for the stack.  Clearly, it's not being quite strict enough in rejecting these pixels, as some of the bright values pop through (corresponding to a pixel that had only the SUSPECT bit set).  Comparing this to the stack detections shows that these bright pixels do create detections (the pixels are >5 sigma deviations from their neighbors, so this is understandable).  Note however, that these objects should have moment values that are clearly wrong.  A grouping of a few pixels is substantially smaller than a PSF.