Flat-field correction

From Wikipedia, the free encyclopedia
Jump to: navigation, search
Not to be confused with Petzval field curvature, which refers to focus uniformity.
The brightness variation due to vignetting, as shown here, can be corrected by selectively brightening the perimeter of the image.

Flat-field correction is a technique used to improve quality in digital imaging. The goal is to remove artifacts from 2-D images that are caused by variations in the pixel-to-pixel sensitivity of the detector and/or by distortions in the optical path. It is a standard calibration procedure in everything from pocket digital cameras to giant telescopes.

Flat fielding refers to the process of compensating for different gains and dark currents in a detector. Once a detector has been appropriately flat-fielded, a uniform signal will create a uniform output (hence flat-field). This then means any further signal is due to the phenomenon being detected and not a systematic error.

A flat-field consists of two numbers for each pixel, the pixel's gain and its dark current (or dark frame). The pixel's gain is how the amount of signal given by the detector varies as a function of the amount of light (or equivalent). The gain is almost always a linear variable, as such the gain is given simply as the ratio of the input and output signals. The dark-current is the amount of signal given out by the detector when there is no incident light (hence dark frame). In many detectors this can also be a function of time, for example in astronomical telescopes it is common to take a dark-frame of the same time as the planned light exposure. The gain and dark-frame for optical systems can also be established by using a series of neutral density filters to give input/output signal information and applying a least squares fit to obtain the values for the dark current and gain. 

C = {{(R - D)*m}\over{(F-D)}} = {(R-D)*G}

where:

  • C = corrected image
  • R = raw image
  • F = flat field image
  • D = dark field or dark frame
  • m = average value of (F-D)
  • G = Gain = m\over(F-D)[1]

In this equation, capital letters are 2D matrices, and lowercase letters are scalars. All matrix operations are performed element-by-element.

In order for an astrophotographer to capture a light frame, he or she must place a light source over the imaging instrument's objective lens such that the light source emanates evenly through the users optics. The photographer must then adjust the exposure of their imaging device (CCD or DSLR camera) so that when the histogram of the image is viewed, a peak reaching about 40–70% of the dynamic range (Maximum range of values for a changeable quantity. In this case its pixel value) of the imaging device is seen. The photographer typically takes 15–20 light frames and performs median stacking. Once the desired light frames are acquired, the objective lens is covered so that no light is allowed in, then 15–20 dark frames are taken, each of equal exposure time as a light frame. These are called Dark-Flat frames.

Flat field correction in X-ray imaging[edit]

In X-ray imaging, the acquired projection images generally suffer from fixed-pattern noise, which is one of the limiting factors of image quality. It may stem from beam inhomogeneity, gain variations of the detector response due to inhomogeneities in the photon conversion yield, losses in charge transport, charge trapping, or variations in the performance of the readout. Also, the scintillator screen may accumulate dust and/or scratches on its surface, resulting in systematic patterns in every acquired X-ray projection image. In X-ray Computed Tomography (CT), fixed-pattern noise is known to significantly degrade the achievable spatial resolution and generally leads to ring or band artifacts in the reconstructed images. Fixed pattern noise can be easily removed using flat field correction. In conventional flat field correction, projection images without sample are acquired with and without the X-ray beam turned on, which are referred to as flat fields (F) and dark fields (D). Based on the acquired flat and dark fields, the measured projection images (P) with sample are then normalized to new images (N) according to [2]

 N = {{(P-D)}\over{(F-D)}}

While conventional flat field correction is an elegant and easy procedure that largely reduce fixed-pattern noise, it heavily relies on the stationarity of the X-ray beam, scintillator response and CCD sensitivity. In practice, however, this assumption is only approximately met. Indeed, detector elements are characterized by intensity dependent, nonlinear response functions and the incident beam often shows time dependent non-uniformities, which render conventional FFC inadequate. In synchrotron X-ray tomography, many factors may cause flat field variations: instability of the bending magnets of the synchrotron, temperature variations due to the water cooling in mirrors and the monochromator, or vibrations of the scintillator and other beamline components. The latter is responsible for the biggest variations in the flat fields. To deal with such variations, a dynamic flat field correction procedure can be employed that estimates a flat field for each individual projection. Through principal component analysis of a set of flat fields, which are acquired prior and/or posterior to the actual scan, eigen flat fields can be computed. A linear combination of the most important eigen flat fields can then used to individually normalize each X-ray projection: [2] 

N_j = {{P_j - \bar D}\over{\bar F+\sum_k w_{jk}u_k-\bar D}}

  • N_j = intensity normalized X-ray projection
  • P_j = raw X-ray projection
  • \bar F = mean flat field image (average of flat fields)
  • u_k = k-th eigen flat field
  • w_{jk} = weight of the eigen flat field u_k
  • \bar D = mean dark field (average of dark fields)

See also[edit]

References[edit]

  • V. Van Nieuwenhove, J. De Beenhouwer, F. De Carlo, L. Mancini, F. Marone, and J. Sijbers (2015). "Dynamic intensity normalization using eigen flat fields in X-ray imaging". Optics Express 23 (21): 27975–27989. doi:10.1364/OE.23.027975. 

External links[edit]