Changes between Version 1 and Version 2 of stack_background_example
- Timestamp:
- Sep 19, 2012, 5:02:29 PM (14 years ago)
Legend:
- Unmodified
- Added
- Removed
- Modified
-
stack_background_example
v1 v2 7 7 I start with an example MD09 stack (MD09.V3.skycell.042.stk.975305.unconv.fits): 8 8 9 [[I MAGE]]9 [[Image(MD09.example.jpg)]] 10 10 11 11 I use psphot to generate a background model image, now displayed with 30x the stretch of the image above: 12 12 13 [[ BACKMDL]]13 [[Image(MD09.example.jpg)]] 14 14 15 15 The 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: 16 16 17 [[I MHIST]]17 [[Image(sky.hist.png)]] 18 18 19 19 The 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. … … 21 21 Now, 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. 22 22 23 [[ LOGIMHIST]]23 [[Image(sky.loghist.png)]]
