next up previous
Next: Basic astronomical image processing Up: Data Reduction Previous: Data reduction procedure

Things to do at the telescope

Biases

Flats

Darks

Checking detector behavior. Image histograms

Determining readout noise and gain.

Consider an ensemble of measurments taken at a light level L. The noise in this ensemble should be $\sigma^2 = LG + rn^2$, where G is the gain and rn is the readout noise in electrons, and $\sigma$ is the noise in electrons.

If you do this at a lot of different light levels, then you can plot $\sigma^2$ vs $L$, and the slope should give you $G$ and the intercept $rn^2$. However, remember that if you compute $\sigma$ from the images, this gives $\sigma$ in DN, so the slope will give you $1/G$. This test is also excellent for checking the basic performancer of a detector. Deviations from nonlinearity can also usually be seen on such a plot.

However, in its most straightforward application the test is very time consuming and hard to analyze: you have to take many exposures at each different light level, and then determine a gain and readout noise for each pixel and look at them all. It is much easier just to use the set of all pixels as your ensemble at each light level. However, you can't do this directly, because each pixel may have a different sensitivity and different fixed pattern noise, so you're not measuring a true ensemble. If there is significant variation of sensitivity than you can't use a whole area at all, because the noise properties will vary across the area. You can avoid these problems by working with the difference between pairs of observations: if the light level is the same in the two images, then you'll be left with an image that only has noise.

Specifically, take a pair of images and form the difference. The expected noise is

\begin{displaymath}\sigma^2 = 2 (L G + rn^2)\end{displaymath}

where $L$ is the light level, $G$ is the gain, and $\sigma$ is the noise in electrons. Since $\sigma(electrons) = G \sigma_{DN}$, we have

\begin{displaymath}\sigma_{DN}^2 = 2 ({L\over G} + {rn^2\over G^2})\end{displaymath}

. You can directly measure $\sigma(DN)$ from your difference image. Make sure to do it over a region which doesn't vary significantly in light level. Now repeat the measurement at a variety of light levels and plot $\sigma^2(DN)$ vs $L$. If you fit a line through this, the slope is $2/G$ and the intercept is $2 rn^2/G^2$. If a straight line doesn't fit the points, then there is some sort of problem, which you should probably track down.

You can abbreviate this test if you just want to get an estimate of the gain and readout noise. First, take a pair of bias frames. These have light level of zero, so the noise from the difference just gives you $\sqrt{2} rn$. (Note that you still need to take a pair in case there is superbias structure). Then take a pair at a high light level; at this level the readout noise is probably negligible, and you can determine the gain from

\begin{displaymath}G = 2 {L\over \sigma_{DN}^2}\end{displaymath}


next up previous
Next: Basic astronomical image processing Up: Data Reduction Previous: Data reduction procedure
Rene Walterbos 2003-04-14