next up previous contents
Next: Estimation Up: Basic Concepts Previous: Image sampling

Noise distributions

Besides gross errors which are discussed in Section 2.2.1 the two main sources of noise in a frame come from the detector system N and from photon shot-noise of the image intensity I (see Equation 2.1). It is assumed that the digitalization is done with sufficiently high resolution to resolve the noise. If not, the quantization of output values gives raise to additional noise and errors.

A large number of independent noise sources from different electronics components normally contributes to the system noise of a detector. Using the central limit theorem, the total noise can be approximated by a Gaussian or normal distribution which has the frequency function :

 \begin{displaymath}P_G(x;\mu,\sigma) = \frac{1}{\sigma\sqrt{2\pi}}
\exp \left( -\frac{1}{2} \left[ \frac{x-\mu}{\sigma} \right]^2
\right)
\end{displaymath} (2.3)

where $\mu$ and $\sigma$ are mean and standard deviation, respectively. The photon noise of a source is Poisson distributed with the probability density PP for a given number of photons n :

 \begin{displaymath}P_P(n;\mu) = \frac{\mu^n}{n!} e^{-\mu}
\end{displaymath} (2.4)

where $\mu$ is the mean intensity of the source. It can be approximated with a Gaussian distribution when $\mu$ becomes large. For photon counting devices the number of events is normally so small that Equation 2.4 must be used while Gaussian approximation often can be used for integrating systems (e.g. CCD's).

In the statistical analysis of the probability distribution of data several estimators based on moments are used. The rth moment mr about the mean $\bar{x}$ and its dimensional form $\alpha_r$ are defined as

 \begin{displaymath}m_r = \frac{1}{N} \sum_{i=1}^N (x_i - \bar{x}) \;\;, \;\;
\alpha_r = \frac{m_r}{\sqrt{m_2^r}}.
\end{displaymath} (2.5)

The second moment is the variance while first is always zero. The general shape of a distribution is characterized by the skewness which denotes its asymmetry (i.e. its third moment $\alpha_3$) and the kurtosis showing how peaked it is (i.e. its fourth moment $\alpha_4$). For a normal distribution, these moments are $\alpha_3 = 0$ and $\alpha_4=3$while for a Poisson distribution are $\alpha_3 = 1/\sqrt{\mu}$ and $\alpha_4=3+1/\mu$. Besides these moments other estimators are used to describe a distribution e.g. median and mode. The median of a distribution is defined as the value which has equally many values above and below it while a mode is the local maximum of the probability density function.


next up previous contents
Next: Estimation Up: Basic Concepts Previous: Image sampling
http://www.eso.org/midas/midas-support.html
1999-06-15