IPP Software Navigation Tools IPP Links Communication Pan-STARRS Links

Changes between Version 1 and Version 2 of stack_background_example


Ignore:
Timestamp:
Sep 19, 2012, 5:02:29 PM (14 years ago)
Author:
eugene
Comment:

--

Legend:

Unmodified
Added
Removed
Modified
  • stack_background_example

    v1 v2  
    77I start with an example MD09 stack (MD09.V3.skycell.042.stk.975305.unconv.fits):
    88
    9 [[IMAGE]]
     9[[Image(MD09.example.jpg)]]
    1010
    1111I use psphot to generate a background model image, now displayed with 30x the stretch of the image above:
    1212
    13 [[BACKMDL]]
     13[[Image(MD09.example.jpg)]]
    1414
    1515The mean and median of the background model are 201.6 and 200.6 respectively (according to the `mana` 'stats' command).  The next plot is a histogram of the original image pixel values:
    1616
    17 [[IMHIST]]
     17[[Image(sky.hist.png)]]
    1818
    1919The red line is a Gaussian fit to the full curve (for values from -1000 to 2000) -- the peak is at 214 counts.  The blue line is the fit shifted to have a peak of 200.  The plot below shows the log(Npixel) plot, to show that the small offset in the simple peak fit is due to the asymmetry of the data (the upper end excess is presumably due to the wings of bright stars and the contribution of faint stars and galaxies).  My conclusion: the psphot background model is getting the right answer (200) for the sky in this image (the psphot analysis with attention to a low/high end bias is apparently necessary to get the right answer).  But, why is the stack so far off from 0.0?  This is especially funny since it is the same psphot background model code is used to make the measurement in the first place.
     
    2121Now, in practice, the observed distribution is the result of combining together ~100 exposures.  In order to be this far off, we would need to be off by ~2 counts on average for each of the input images.   
    2222
    23 [[LOGIMHIST]]
     23[[Image(sky.loghist.png)]]