next previous contents

SECTION 5 - PHOTOMETRIC CORRECTION

5.1 OVERVIEW

Because raw images reflect the nonuniformity and nonlinearity of the detector response, a photometric correction is performed on a pixel-by-pixel basis so that pixel DN values are replaced by normalized fluxes (flux numbers, or FN) which are linearly related to the incident photon flux on the vidicon faceplate. For high dispersion images, each pixel value within a circle (as mapped into geometrically correct space) encompassing the useful portion of the image is converted to FN; pixels outside the circle are left with uncorrected (DN) values. In low dispersion, the corrected portion of the image is further restricted to a swath encompassing the spectral order. Detailed data defining these regions and the photometric correction process itself are contained in the following subsections.

5.2 DEFINITION OF ITF

5.2.1 DEFINITION OF ITF

The intensity transfer function (ITF) contains the information needed to map DN to FN on a pixel-by-pixel basis across an IUE image. Given a pixel with a particular DN value, the ITF is used, in the manner described in Section 5.3, to obtain the corresponding linearized FN value.

5.2.2 CHARACTERISTICS OF ITF

The ITF for each camera is defined in geometrically correct space and is generated from a series of geometrically corrected mercury flood-lamp flat-field images at graded exposure levels (11 levels for SWP; 12 levels for LWR and LWP). for the SWP and LWR cameras, each exposure level in the ITF is constructed of an average of several individual images from which reseau marks have been removed by the REMRES program described in Perry and Turnrose (1977). The lowest level is the zero-exposure level ("null", or "pedestal") representing the background left by the vidicon camera flood/erase prep procedure which precedes every normal exposure. The highest level is a several-hundred-second flood-lamp exposure which reaches the vidicon saturation limit of 255 DN over part of the image; less sensitive areas remain unsaturated, although over most of the tube in SWP and LWR, a DN of 220 is reached or exceeded. Holm (1979) and Turnrose (1980) have discussed this aspect of the ITF in the context of low dispersion spectra, where the least sensitive areas of the tube correspond to the longer wavelengths in LWR and the shorter wavelengths in SWP. Table 5-1a, reprinted from Holm (1979), lists the appropriate upper DN limits of the current SWP and LWR ITFs in the vicinity of the low dispersion spectral orders. Table 5-1b, reprinted from Imhoff (1984b), lists the maximum DN levels for an LWP low dispersion spectrum in which the best photometric accuracy of the ITF is preserved. As described in Section 5.3, input image DN levels which exceed the top level of the ITF are now converted to FN by extrapolation.

Table 5-1a:   Approximate Highest DN Values in the SWP and LWR ITFs (for the Low Dispersion Format)
SWP LWR
lambda DN lambda DN
1200197210025O
13002162300245
14002352500240
15002452700225
1600243*2900220
1700239*3100190
1800228*3300180
1900242*  
* Higher levels exist which may include saturated pixels.


Table 5-1b:
  Maximum DN Levels for Best Accuracy

In LWP Low Dispersion Spectra

LWP
lambdaDN
2100205
2300220
2500240
2700245
2900245
3100245
3300245


Intermediate ITF levels are generated from intermediate-length exposures in order to map the camera response through the mid-range of DN values and are spaced sufficiently closely to allow the conversion to FN by linear interpolation between the levels as described in Section 5.3.

5.2.2.1 ITF Exposure Values and Scale Factors

An exposure value (FN value) is associated with each level of the ITF. This exposure value is normalized in an arbitrary (but fixed) way so that the FN associated with the highest level of the ITF is near 20,000 (for SWP the highest ITF flux level is 17632; for LWR, 25220; and for LWP, 19784). The FN associated with level i of the ITF for camera k is

(5-1)
T(k,i) MULT(k,i)
FN(k,i) = -----------------
FACTOR(k)


where T(k,i) is the effective exposure time (in seconds) for level i and camera k, MULT(k,i) is a scaling factor for level i and camera k, and FACTOR(k) is a level-independent scale factor for camera k. Although in principle MULT(k,i) can vary from level to level, (e.g., to account for flood- lamp variability), in practice MULT(k,i) is used as a function of k only, and lamp drifts during the creation of the ITF are accounted for in adjustments to the effective exposure times T(k,i). Tables
5-2 and 5-3 list the effective exposure times and the scale factors for the three ITFs in use as of November 1981. (These ITFs are designated LWP ITF1, LWR ITF1, and SWP ITF1). With these data and equation 5-1, the exposure value of each ITF level is determined.

Table 5-2:   Effective Exposure Times in ITFs
ITFT (seconds)
Level no.LWPLWRSWP
1 0.00 0.00 0.00
223.0023.0316.84
339.6940.6933.74
460.6280.0868.73
577.90100.7390.91
6102.56118.78105.86
7129.82158.83143.71
8148.50201.49177.45
9179.00244.71215.24
10205.62293.91251.05
11233.79343.33285.00
12329.73420.32 



Table 5-3:   Scale Factors for ITFs
 LWPLWRSWP
MULT17.00 17.00 11.00
FACTOR 0.28333 0.28333 0.1778


As is pointed out in Section 9, the T (in units of 0.01 seconds) and MULT values pertinent to the ITF actually applied to a given image are documented in the processing history portion of the image header labels of the photometrically corrected image and all derived files.

5.2.2.2 Highest Valid ITF DN Levels


Although as mentioned in Section 5.2.2 the highest exposure levels in the ITFs reach the saturation limit of 255 DN over portions of each camera, ITF DN values about 250 present potential hazards due to "hidden" saturation effects. Whitman, Bohlin, and Turnrose (1981) discuss the fact that because of noise in the individual flat-field images of which the ITF is composed and the fact that the production ITFs are composed of averages of several independent (and geometrically corrected) images, even final ITF values several DN below 255 may have been influenced by saturation in some of the raw constituent images.

To alleviate such potential hidden saturation effects, the DN value of 250 has been defined as the maximum valid ITF DN level. This in essence means that DN levels >250 in the ITF are ignored (except in the pathological case (discussed in Section 5.3.2.2.2) in which there are insufficient valid ITF levels for the pixel in question). More specifically, this means that DN > 250 in spectral data images will be converted to FN by extrapolation from lower values, and that ITF extrapolations will be based only on points with DN <= 250; see Section 5.3.2.2.



5.3 CONVERSION OF DN TO FN

Pixels in the raw image are converted from DN to FN by the program PHOTOM only within the portion of the image relevant to the spectral data. In the case of high dispersion, this means that only pixels within a circle slightly smaller than the target boundary are converted; pixels outside the circle are left unconverted. Note that this circle is defined in geometrically correct space, hence the photometrically corrected region will not be bounded by a true circle when viewed in the raw geometry. Table 5-4 lists the center and radius values of these circles for the three operational cameras.



Table 5-4: Center and Radius Values* Defining Region of Photometric Correction
 Center 
CameraLineSampleRadius
LWP400390347
LWR395402350
SWP390390358
*In geometrically correct coordinates

In low dispersion, it is unnecessary to convert the entire region within the circles defined above, so only a swath or band encompassing the spectral orders is converted. This band, oriented parallel to and centered between the large and small-aperture spectra, has a width in the sample direction of 320 pixels. Note that in previous documentation this width has been incorrectly specified as 160 pixels. Although the 160-pixel width would be sufficient to encompass the line-by-line data extracted in low dispersion (see Section 7), the larger area has been retained to allow the conversion of sufficient portions of the tube to facilitate potential studies of the tube background.

In August 1982, a further restriction to the area photometrically corrected in the photometric correction in low dispersion is done only within the camera-dependent rectangular area defined by the standard partial read parameters; pixels outside this area are left in DN units. Table 3-1 presented in Section 3.3 lists these parameters, which apply directly to raw-image space. Note that this limitation on the photometric correction is applied to all low dispersion images, whether or not they are partial reads. The only effect on full-frame low dispersion images, however, is a slight truncation of certain corners of the photometrically corrected swath; no extracted spectral data are affected.

5.3.1 SPATIAL MAPPING AND INTERPOLATION

Since the ITFs are defined in geometrically correct space, the conversion of a pixel DN value in a raw image to FN requires first that the corresponding position in the ITF be determined; conceptually, this is done via the geometric mapping function G discussed in Section 4.4. (In practice, for computational simplicity the mapping is actually done in the reverse direction (geom-to-raw), using the G-1 function). Since in general the position in geometrically correct space resulting from this mapping will not coincide with an integral pixel in the ITF, the ITF must be spatially interpolated. To do this, the input pixel DN value is first used to calculate four FN values using the ITF curves of the four pixels surrounding the mapped location; details of such FN assignments are addressed in Section 5.3.2. A flux value representing the final FN corresponding to the input pixel is then computed from the bilinear spatial interpolation of these four neighboring FN values.

Following Lindler (1982a), this process is described mathematically below. Given a pixel with value DN at the location in the raw image (s,l), its position (x,y) in the geometrically corrected ITF is

(5-2)
(x,y) = G(s,l)

where G is defined in Section
4.4. If the function IFIX(x) returns the largest integer <= x, then define

(5-3)
X = IFIX(x)

(5-4)
Y = IFIX(y).

If we now represent the FN-assignment function (described in Section 5.3.2) for a pixel at location (X,Y) with value DN by F(X,Y,DN), then the four neighboring flux values to be used in the spatial interpolation are

(5-5)
F1 = F(X,Y,DN)

F2 = F(X,Y+1,DN)

F3 = F(X+1,Y,DN)

F4 = F(X+1,Y+1,DN).

Finally, the interpolated FN value is calculated as below. If

(5-6)
u = x - X


(5-7)
v = y - Y


then

(5-8)
FN = (1-u)(1-v)F1 + (1-u)vF2 + u(1-v)F3 + uvF 4.

5.3.2 ASSIGNMENT OF FLUX VALUES

In this subsection the techniques used to evaluate the function F(X,Y,DN) defined in Section 5.3.1 are described. F(X,Y,DN) represents the conversion of DN to FN for the specific integer-pixel location (X,Y) in the geometrically correct frame of reference. The evaluation of F(X,Y,DN) proceeds in one of several ways using either interpolation or extrapolation depending on where the DN value falls with respect to the ITF for that pixel.

5.3.2.1 Normal Interpolation Technique

In most instances when the input DN is less than or equal to the maximum valid ITF DN value of 250 (see Section 5.2.2.2), an FN value is assigned by simple linear interpolation between the flux values of the various levels of the ITF. That is, the input DN is compared linearly to the DN corresponding to each level of the ITF, and the associated flux is interpolated accordingly. The majority of the pixels in a well-exposed image will fall below DN = 250 and will thus be converted to FN by this simple interpolation process.

There are three instances, however, in which the input DN may be less than or equal to 250 but extrapolation rather than interpolation will be used. These are when:

All cases in which the input DN exceeds 250 (including saturation) are handled by extrapolation as discussed in Sections 5.3.2.2.2 and 5.3.2.2.3.

5.3.2.2 Extrapolation Techniques

5.3.2.2.1 Negative FN

Extrapolation to negative FN values occurs in all instances in which the input DN falls below the null-level of the ITF. This can result from statistical fluctuations or null drift in which the null level of the particular input image being photometrically corrected falls below the ITF null level as a whole. In this situation FN values are assigned by simple linear extrapolation from the first two points of the ITF (null level and lowest finite-exposure-time level). The current software can accommodate extrapolation down as far as FN = -3488.

5.3.2.2.2 Positive FN

Extrapolation to positive FN values occurs in all instances in which: