next previous
Up: Deriving object visibilities

7. Dealing with noisy data

  We consider here that the photometers in the fiber interferometer are detector noise limited. Each signal measurement is then affected by an additive, stationary noise, uncorrelated with the data. For each interferogram, a realization of the four noise signals is recorded in the background current sequences. These noise signals have to be taken into account when estimating the photometric signals and the squared coherence factor.

7.1. Estimating the photometric signals

  The measurement tex2html_wrap_inline3186 of the signal produced by photometric detector tex2html_wrap_inline3188 is the sum of tex2html_wrap_inline3190 and the additive noise tex2html_wrap_inline3192:
equation1019

The estimator tex2html_wrap_inline3194 of tex2html_wrap_inline3196 that minimizes the mean quadratic error with the actual signal is obtained by optimal filtering (Press et al. 1988):
equation1025
where tex2html_wrap_inline3198 is the Wiener filter, whose expression is, when the signal and the noise are uncorrelated:
 eqnarray1032

The power spectrum of the background current can be used to estimate tex2html_wrap_inline3200.

The quality of the deconvolution signal tex2html_wrap_inline3202 is critical to achieve a good interferogram correction. What matters is not the average of tex2html_wrap_inline3204 but rather its smallest value: if at the minimum the local signal to noise ratio is too low, then the division is numerically unstable and the interferogram correction process is useless. Therefore it is necessary to adopt a selection scheme: if the minimum value of tex2html_wrap_inline3206 or tex2html_wrap_inline3208 is below a certain rejection threshold, then the interferogram is discarded and does not lead to a fringe visibility measurement. For a photometric signal tex2html_wrap_inline3210, the rejection threshold can be expressed as a multiple of the standard deviation tex2html_wrap_inline3212 of the corresponding background sequence, filtered by tex2html_wrap_inline3214.

There is a trade-off in the choice of the rejection threshold, between rejecting only a few interferograms and having a better correction. For sufficiently good data this choice is not critical, as the correction quality does not significantly depend on the minimum value of the tex2html_wrap_inline3216 as long as this value is greater than 10tex2html_wrap_inline3218. For faint objects or in bad seeing conditions, one may have to use rejection thresholds as low as 3tex2html_wrap_inline3220.

7.2. Estimating the squared coherence factor

We should now seek an estimator of the coherence factor that is not biased by detector noise. The signal M(x) actually measured at the output of the interferometric detector contains an additive noise tex2html_wrap_inline3224:
equation1062
If the data processing described in the preceding sections is applied to M, it leads in the wave number space to
equation1065
and, since the signal and the noise are uncorrelated, their moduli add quadratically:
 equation1073
which yields, in the integral form that is assumed to be valid when there is differential piston:
 equation1082

equation1100

An estimator tex2html_wrap_inline3228 of the noise integral is obtained by applying to the current signal the same treatment that was used to process M. A better estimator is obtained by processing many realizations of the interferometric background current (for example, all the background sequences of a complete batch), and averaging the resulting power spectra. Eq. (58 (click here)) then provides an estimator of S:
equation1104
from which Eq. (45 (click here)) enables us to establish an unbiased estimator of the squared modulus of the coherence factor, averaged over the optical bandpass:
 equation1109
The quantity tex2html_wrap_inline3234 is the final result of the data reduction process on a single interferogram.

7.3. Statistical error and noise sources

 

  figure1118
Figure 7: Result of the data reduction on a batch of 122 interferograms of tex2html_wrap_inline3236Boo. The adopted shape factor was tex2html_wrap_inline3238cm. Because of mediocre seeing conditions, the rejection threshold was set to 3: only 47 interferograms met this level and led to a measurement of the squared coherence factor. For that batch the final estimate of the squared coherence factor is tex2html_wrap_inline3240

A batch of interferograms leads to a collection of n measurements of the fringe visibility. An example is shown in Fig. 7 (click here). The final estimate of the squared coherence factor is the average of the tex2html_wrap_inline3244, and if tex2html_wrap_inline3246 is the standard deviation of the estimator of tex2html_wrap_inline3248, then the standard error is
equation1127
The final error estimate for tex2html_wrap_inline3250 should include the uncertainties on tex2html_wrap_inline3252 and tex2html_wrap_inline3254.

Fluctuations in the measurement of tex2html_wrap_inline3256 can be attributed to three main noise sources:

Another source of ``noise'' results from the incomplete realization of the chromaticism assumption developed in Sect. 3.2 (click here). It will not be discussed here, since in the standard operating conditions of an infrared fiber interferometer (tex2html_wrap_inline3266, bandpass limited by the atmospheric transmission window) its effects are always negligible.

Although we have so far considered only one interferometric output, each scan provides two measurements of tex2html_wrap_inline3268. Coherence factor measurements on each output must be treated separately because each channel has its own instrumental transfer function. But the correlation between the two channels can tell us about the relative importance of the noise sources.

Detector noise is fully uncorrelated between the two interferometric outputs. Conversely, the differential piston perturbations are exactly the same on each output. Deconvolution noise is very strongly correlated since the deconvolution signal is the same for both channels, but not fully correlated because the deconvolution is applied to two separate measurements of the interferometric signal, each one affected by independent noise.

Let tex2html_wrap_inline3270 be the estimator of tex2html_wrap_inline3272 for Channel 1 (Eq. 61 (click here)) and tex2html_wrap_inline3274 for Channel 2. The following holds without making any assumption on the nature of the noise sources:
 equation1148
where tex2html_wrap_inline3276 is the correlated part of the noise for Channel i, and tex2html_wrap_inline3280 is the uncorrelated part. For each interferogram, it is possible to measure the difference tex2html_wrap_inline3282 between the two estimators, i.e.
eqnarray1170
and for a batch of scans the variance of tex2html_wrap_inline3284 is the sum of the variances of the uncorrelated noises:
equation1182
If two photometers of the same type have been employed to measure the interferometric outputs, the statistics of the uncorrelated noise is the same on both channels. It is then possible to determine the variance of the uncorrelated noise:
equation1187
The variance of the correlated noise follows immediately from Eq. (63 (click here)):
equation1193

From what was said above, we can with a good approximation identify tex2html_wrap_inline3286 with the detector noise, and tex2html_wrap_inline3288 with the joint contribution of piston and deconvolution noise. Thus simultaneous measurements of both interferometric outputs make it possible to tell the relative contribution of detector noise to the total noise.

This is useful to work out the optimum fringe speed v, which is the result of a compromise. On the one hand detector noise increases with v, because for a given optical bandpass the fringe signal is spread over a wider frequency bandpass. On the other hand the faster the interferogram is scanned, the more frozen is the seeing for the duration of the scan, and the piston perturbations are smaller. In the tex2html_wrap_inline3294 example, 44% of the noise variance is uncorrelated, which means that no single source dominates the noise.


next previous
Up: Deriving object visibilities

Copyright by the European Southern Observatory (ESO)
web@ed-phys.fr