Method for Calibration of SATI Instruments

Sci report | May 09, 2022

In this work, we suggest a method for calibrating the SATI (Spectral Airglow Temperature Imager) instruments [2] [3] as an alternative to the traditional scheme applied until now involving an integrating photometric sphere. The method shown here uses a known law for the illumination distribution along the instrument field of view. The experimentally obtained flat-field is numerically processed and corrected according to the illumination law with the described scheme. The achieved calibration precision is more than 0.5%. Some of the advantages of the suggested method are the simplicity of its technical performance and portability, which is essential for the cross-calibration of SATI instruments operating at different geographic locations.

Introduction

The cross-calibration and testing of the airglow optical instruments are one of the goals of the developed NATO CLG and INTAS projects. The implementation of this task is also connected to making a low-cost laboratory facility for calibration and testing.

In this paper we present a calibration method for the SATI instruments used in the calibration of SATI-3SZ in Stara Zagora station. The suggested scheme is an alternative to the flat-field calibration by an integrating photometric sphere applied so far [1]. The exposed approach adheres to the specified requirements during technical implementation as well as to mobility.

SATI instrument.
Spectral Airglow Temperature Imager.

Design of the method

The SATI-3SZ calibration includes two stages:

Several laboratory settings have been made for the realization of the first stage:

  1. The SATI calibration follows the optical scheme suggested by Canadian specialists from the Centre for Research in Earth and Space Science (CRESS), York University in Toronto, and is shown in the figure below. A light source illuminates a flat, diffuse reflecting white screen. The scattered light is reflected by the white ring put around the instrument's body towards the concave hemispheric surface and falls within the instrument's field of view.

    SATI calibration with white diffuse hemisphere.
    Calibration scheme using white reflecting hemisphere.

    This scheme has one important disadvantage. It does not ensure equal illuminance along the whole SATI field of view. Therefore, the results of this setting are accepted as insufficiently correct.

  2. SATI-3SZ calibration by a photometric integrating sphere. Due to the specific construction of the sphere we used and the particular illumination scheme (four small halogen lamps arranged in a cross-like way near the sphere inlet), equal luminance of the scattering surface is not yet guaranteed. Precise measurement of the brightness in the radial direction was not performed, which is a disadvantage of the setting. Most integrating spheres are adapted for photometric instrumentation with a central, comparatively small field of view, while the instruments of the SATI or MORTI (Mesopause Oxygen Rotational Temperature Imager) [3] type have an angular field of view with external limits, reaching up to 70°. The last fact raises new strict requirements for the luminance homogeneity of the calibrating spherical surface as well as for the quality of the reflecting coating.

    SATI calibration setup with integrating sphere.
    SATI calibration using photometric integrating sphere.

    The results obtained by the integrating sphere method are compared to the results calculated in our methodology, described in the next section. The relative difference is no more than 2%. Yet, this method requires high-cost and complex laboratory equipment as well as lacks any portability.

A method using a flat scattering surface

Our method uses a flat scattering surface. It has several advantages compared to the previous methods and is very suitable for the cross-calibration of the separate instruments operating around the world. The setting scheme is shown in the figure below.

SATI-3SZ calibration using organic glass diffusor.
Calibration setup using flat opal organic glass screen.

A point-like light source S is located at a distance h (in our case 1500 mm) along the optical axis over the SATI and illuminates a flat scattering screen fixed on top of the instrument. The considerations for the selection of a flat surface are the following:

Due to the specific path of the light rays for the SATI instrument, the separate light beams, forming different falling angles to the interference filter, see the flat surface always under the same angle. The parallel light beams falling on the filter come out of virtual point sources located along the radii of the input field stop, in the front focal plane of the Fresnel lens. These virtual point sources, on their part, are formed by light beams, equal in angular distribution (between 26° and 34°), yet falling on the conic mirror at different heights. In other words, these beams are identical and only translated to one another. Therefore, it is not necessary for the scattering surface to be strictly Lambertian (The Lambertian surface has the property to scatter the light in different angular directions in a specific way, described by the cosine law of the angle to the normal of the surface. As a consequence, the brightness in different directions is measured equal).

SATI instrument optical scheme.
SATI light rays path.

Each homogeneous white surface, sufficiently scattering the passing light, can be used as a diffuse material. In our implementation, we use a 4 mm thick white organic glass flat surface.

Another consideration in support of the flat scattering surface using the light scheme described so far is the well-known distribution of the illuminance E along a given radius in the instrument field of view, set out by Lambert’s cosine law: Е = I cos(θ) / L2. In this way, the illumination EP in point P at a distance r from the optical axis is:

EP = I ⋅ h2 / (r2 + h2)3/2

where I is the intensity of the light source, accepted for convenience to be equal to one. In practice, it should have a sufficiently high intensity value in order to provide the necessary signal-to-noise ratio (SNR) for a typical SATI exposition of 120 s. Because h is 1500 mm and r is almost ten times smaller, the above formula can be simplified using Taylor series, so we yield:

EP ≈ I / h ⋅ (1 - 3/2 ⋅ r2 / h2)

In our case, the maximum difference between the exact and approximate equation is no more than 0.05%.

Practical implementation

A small halogen lamp of 100 W is used as a light source. A round-shaped glass window, 2 mm thick and 10 mm in diameter, is fixed in front of it and is mat-coated on both sides. It plays the role of a point light source. As it is well known, the halogen lamp has a continuous emission spectrum with a shape close to the Planckian curve at Teff ≈ 3000 K. This provides a homogeneous spectrum in the whole operating range of the SATI instrument between 8645 Å and 8680 Å. The difference in the emission intensity between the two wavelengths is no more than 0.25%.

Approximately 250 images at a typical SATI exposition of 120 s are obtained and then averaged. Dark current images are recorded after every eight exposures. The flat-field obtained in this way should be numerically corrected because of the illuminance difference ΔEP1,P2 of the two points P1 and P2 located at a distance of r1 and r2, respectively, along the SATI field of view. ΔE decreases with the increase of h. In the case of SATI-3SZ, h is 1500 mm, so ΔE is less than 1.5%.

The next figure shows that a virtual point source S' located at a distance x from the optical axis in the plane of the input field stop D is formed by the beams coming from the ring-shaped area on the diffusor, locked between the radia r1 and r2.

Point light sources formed at the instrument's entrance.
Virtual light source formed at SATI's entrance stop.

In our case:

r1(x) = 330⋅tan(26°) - 1.110⋅x
r2(x) = 330⋅tan(34°) - 1.203⋅x

Here r1, r2 and x are expressed in mm. On the other side, with known h, the total energetic flux of the source S' can be written as a definite integral taken from r1(x) to r2(x), or:

Ф(х) = ∫r1(x)r2(x) E(r) dr

For convenience, Ф(х) is expressed as a function of radius Rb of a given radial zone of the image in bins (in our case, one bin is composed of 2×2, i.e., four pixels from the CCD sensor) and hence, by the wavelength λ. For this purpose, the following equations are used:

x = 0.312⋅Rb
λ = λ0 ⋅ (1 - Rb2 / (µ2 ⋅ (6252 + Rb2)))1/2

The numerical coefficients 0.312 and 6252 are determined by the specific geometric parameters of SATI-3SZ, such as the focal length of the CCD camera lens, the focal length of the Fresnel lens and bin size. λ0 and μ are the central wavelength of the interference filter and the effective refractive index, respectively. The experimentally obtained flat-field is corrected by the numerical coefficient Ф = Ф(λ) normalized to a unit towards its maximum value. The result is presented in the figure below. It is seen that this profile coincides very precisely with the profile obtained by the photometric sphere, and the relative difference between the two is no more than 2%. Probably, this is due to the inequality of the brightness of the sphere's surface, falling within the SATI field of view.

Flat-field profiles.
Flat-field profiles using different calibration schemes: magenta – Canadian setup, red – photometric sphere, blue – SATI-3SZ.

The second stage of the calibration process is the calculation of the SATI responsivity, expressed in absolute energetic units. It is performed using a Radiometer / Photometer Model 730A, Optronic Laboratories Inc. set in radiometric mode and equipped with a radiometric filter. The receiver measures the energetic brightness in W/cm2•sr. The filter ensures the constancy of the radiometric detector responsivity function in the range from 480 nm to 1000 nm with an accuracy of ±5%. The brightness of the diffusor is determined in 10 points located in a circle with a radius of 185 mm and then is averaged. It is taken into account that the quantity of energy registered from the radiometric detector can be expressed as an integral of the product of the detector's spectral responsivity and the spectral profile of the light source (we accept that the source emits as a black body with an effective temperature of 3000 K). SATI has a frequency range much narrow (8645 Å – 8680 Å) so the energy accepted by the instrument is about 120 times less than the one measured by the radiometer. The average responsivity obtained is 131 R/nm per ADU/s, where ADU is an analog-to-digital unit or simply the count. The last obtained SATI-3SZ profile, expressed in ADU•nm/R•s for a given wavelength λ, is shown in the figure below.

Final averaged SATI-3SZ flad-field calibration profile.
SATI-3SZ flat-field profile.

Conclusion

One possible source of inaccuracies in the proposed method could be the inhomogeneity of the diffuse material, which can be determined experimentally and corrected if necessary. Another source of errors could be the inaccurate centering of the light source relative to the optical axis of the SATI instrument. As final conclusions, the following observations can be made:

Also, in the future, the spectral profile of the light source will be measured with high-enough accuracy. The whole flat-field calibration setting will be realized in a compact mobile version.

Reference

  1. Sargoytchev S., Brown S., Solheim B., Cho Y., Shepherd G., Lopez-Gonzalez M. Spectral airglow temperature imager (SATI) – a ground-based monitor of mesosphere temperature, Applied Optics, 2004.
  2. Shepherd G. Spectral imaging of the atmosphere, International Geophysics Series, Vol. 82, 2002.
  3. Wiens R., Zhang S., Peterson R., Shepherd G. MORTI: a mesopause oxygen rotational temperature imager, Planet. Space Sci., Vol. 39, No. 10, pp. 1363-1375, 1991.

An early version of this report was presented at the Joined NATO-CLG/INTAS Workshop Investigation of Regional Scale Atmospheric Motions and their Influence on the Mesosphere/Thermosphere/Ionosphere Region, September 6-10, 2004, Granada, Spain.