Download The Effect of Stacking Multiple Sub

The Effect of Stacking Multiple Sub-exposures on Signal
to Noise Ratio
When making astronomical images with a CCD camera, the question often arises as to
whether it is better to take one long exposure or to stack a number of shorter exposures
(called here sub-exposures). If the sub-exposures are accurately aligned and the
background noise is random and the wanted parts of the image coherently repeatable
across all images then the increase in signal to noise ratio (SNR) is proportional to the
square root of the number of items stacked. However many of the noise sources in a CCD
are more significant in shorter exposures than in longer exposures and thus the resultant
SNR is not as great as might be hoped for. This note examines the amount of loss.
The contribution of the various noise sources is expressed by the well known ‘CCD
equation’, which can be presented in different ways. We will use the following form:
S STAR * t
SNR 
S STAR * t  N SKY * t  N DARK * t  N RO
2
Where the symbols have the following meanings:
SNR
Sstar
t
Nsky
Ndark
Nro
= signal to noise ratio
= signal rate from the wanted part of the image in ADUs/second
= exposure time in seconds
= signal rate from sky background in ADUs/second
= signal rate from thermal noise in ADUS/second
= readout noise in ADUs rms
We assume in the above that we are examining a single pixel in the image. The star noise,
sky noise and dark noise all consist of electrons that arrive at unpredictable intervals of
time, and have what is known as a Poisson distribution. They all have the nature of what is
sometimes called shot noise. The readout noise has a Gaussian distribution.
If we now take a single long exposure of duration TL then the resulting SNR will be:
SNR 
S STAR * TL
S STAR * TL  N SKY * TL  N DARK * TL  N RO
2
------------------- (1)
Similarly the result of one short exposure of duration TS will be:
SNR 
S STAR * TS
S STAR * TS  N SKY * TS  N DARK * TS  N RO
2
Stacking a number M of the short exposures will give
SNR 
S STAR * TS * M
S STAR * TS  N SKY * TS  N DARK * TS  N RO
2
------------------- (2)
To find how many we need to stack, we equate (1) and (2) above and rearrange to give:
2
* TL  N SKY * TL  N DARK * TL  N RO
2
2
M

S
TL S STAR * TS  N SKY * TS  N DARK * TS  N RO
TS
2
STAR


We see that if the readout noise is negligible then M reduces to TL/TS as expected.
Otherwise it may be expressed as follows:

D 
TL 
M 
*
TS 
D 


N RO
TS
2
N RO
TL
2








,
where D  S STAR  N SKY  N DARK
We can see that there is no simple ‘rule of thumb’ to tell us how many sub-exposures to
stack for any required result. Even for a given imaging situation, where D is constant, the
extra number to stack is dependent not only on the proportions of the noise sources but
also on the ratio of TL to TS: there are gains to be had from making the sub exposure
period as long as conveniently possible, usually limited by tracking errors or saturation of
brighter parts of the image.
To help our understanding let us look at a typical image:
This image is a single unguided 10 second subexposure of NGC2903 with a C11
telescope and f6.3 focal reducer. The camera was an Artemis 285. The temperature of the
CCD was approximately -10degC. The limiting visual magnitude was about 4.3 in this
area.
This camera has the following characteristics at this temperature:
N DARK  0.0227 ADU / sec
N RO  18.2 ADUs
From measurements taken from the picture,
S STAR  119.4 ADU / sec for the brightest part of the galaxy core
N SKY  10.93 ADU / sec
 D  119.4  10.93  0.0227  130.35 ADU / sec
To achieve the same result as one 600 second exposure:
18.2  18.2 

130.35 

600 
10
  74.92
M 

18.2  18.2 
10 
130.35 

600 

i.e. we would need to stack 75 sub-exposures rather than the 60 expected by simple
scaling.
Of more interest is to repeat the calculation for objects, such as faint outer arms of
galaxies, that are at or near to the sky noise level:
S STAR  10.93 ADU / sec for such objects in the above image.
And D  10.93  10.93  0.0227  21.88 ADU / sec
18.2  18.2 

 21.88 

600 
10
  147.12

Now M 
18.2  18.2 
10 
 21.88 

600 

i.e. we now need to stack 147 sub-exposures. This is 2.45 times more than the expected
value of 60. For objects that are below the sky noise the increase would be even greater.
Conclusions
Stacking sub-exposures does give a useful increase in SNR, but the amount of that
increase is reduced by the effect of the noise sources. We see that readout noise should
be small compared with sky and dark noise for the best results. Also it pays to keep the
sub-exposures as long as possible, depending on other factors; not only does this
increase the SNR but it also minimises the proportion of time lost downloading the images
from the camera to the computer.
J C Moore
Revised 24th May 2007