Atmospheric correction of high resolution data

The present exercise refers to an introductory application of atmospheric correction performed on remotely sensed data in the area around Providence, Rhode Island, USA. The source data consists of three Landsat 5 TM images (bands 2, 3 and 4 for green [0,6–0,7 μm], red [0,7–0,8 μm] and near-IR [0,8–1,1 μm] respectively), the corresponding metadata text file, and a raster binary layer indicating the locations of the pixels of the Landsat images where lawn grass can be found. The three band layers were combined into one false-color composite image, a basic atmospheric correction was then performed using the metadata attributes and finally the correction method was evaluated by examining the spectral signature of a known type of land cover, i.e. lawn grass.

The atmosphere accounts for distortions in the radiance that reaches the satellite sensor, mainly because of scattering. The two broadly identified scattering mechanisms, Rayleigh scattering (caused by air molecules) and Mie scattering (caused by the aerosols, such as smoke, haze, water vapor and fumes) affect the radiance captured by the satellite sensors and obscure the fine detail in the data stream. The following diagram shows the travel of the electromagnetic radiation from the sun through the atmosphere and back to the satellite sensor.

The effects of the atmosphere in determining various paths for energy to illuminate a (equivalent ground) pixel and to reach the sensor (Richards & Jia, 2006)

The tools available for remote sensing data processing can help reducing the atmospheric distortions and produce “clearer” data for further analysis.

Data and methods

The general steps followed can be summarized as follows:

  1. Data acquisition and exploration
    • Acquisition of the remote sensing data
    • Creation of a three-layer composite from the bands 2, 3 and 4
    • Exploration and close observation of the produced false-color image
  2. Atmospheric correction
    • Exploration of the metadata.txt file
    • Running the ATMOSC module using the values obtained from the metadata file and the DN value haze, after being converted to radiance values, for all three bands
    • Comparison of the results, which is done with the help of histograms for each band, before and after the correction
  3. Evaluation using spectral libraries
    • Extraction of the spectral values for the ‘lawn grass’ locations, as defined by the masking layer, from the atmospherically corrected images, using the EXTRACT module
    • Comparison with the spectral libraries

Source data

The atmospheric correction was performed on three Landsat 5 TM images of bands 2, 3 and 4 of the area around the city of Providence, Rhode Island, USA, acquired on the 9th of September 1987. The resulting false color composite was visually compared to the original false color composite. The correction was assessed using a raster layer with binary value pixels (0 and 1), where ‘1’ indicated an area with lawn grass land cover and ‘0’ an area without lawn grass. The spectral signature of lawn grass was given as well.

Location of the source data on the globe

Atmospheric correction

Many atmospheric correction methods have been proposed for multi-spectral satellite imagery, mainly consisting of image-based methods, methods that use atmospheric modeling and methods that use ground data.

ATMOSC module

IDRISI’s module to perform atmospheric corrections on loaded imagery is ATMOSC, a sophisticated tool, available under “Restoration” menu item. The tool offers four models to perform the corrections:

  • Dark object subtraction model
  • Cos(t) model
  • Full correlation model, and
  • Apparent reflectance model

The model approach used for the atmospheric correction of the present data is the Cos(t), created by Chavez in 1996. The model is based on a technique for approximation, in the events of insufficient data available for proper atmospheric correction. The model incorporates all the elements of the ‘dark object subtraction’ model, which looks for the lowest values in the data (usually at pixels located over deep waters, where reflectance is known to be minimal), regards them as haze and removes them from all the other values. The model requires a number of inputs, all of which can be extracted from the metadata.txt file.

Spectral libraries

The following spectral values have been measured by USGS for ‘lawn grass’ land cover and are normally used to calibrate remote sensing systems. These are the values a “perfect” remote sensing system should capture in ideal conditions and are used as reference values, for the purposes of the present analysis. The bands of particular interest are in bold.

BandSpectral reflectance
14,043227E-02
27,830066E-02
34,706150E-02
46,998996E-01
53,204015E-01
71,464245E-01

The source data were captured approximately three decades ago; it is thus reasonable to expect errors such as the within-band line striping caused by the age of the sensors and a generally lower-quality product. However, such types of corrections fall out of the scope this exercise.

Results

The two images below were created by composing the three band layers and assigning them RGB colors; blue for band 2, green for band 3 and red for band 4. The left composite is the original, unaltered data, where each pixel of each layer carries DN (digital number) values ranging from 0–255 (8-bit), while the right composite carries radiance/reflectance values ranging from 0–1 and is the atmospherically corrected version.

A more precise insight into the data can be offered by the histogram view. The histograms below were adjusted by editing the Graphic View Settings: too high and too low values were excluded and the bar width was set to the most visually proper value. The histograms show the distribution of the values for each image, making it easy to identify the lowest which is used as input to the “haze” cell in ATMOSC.

Extremely large numbers of pixels with very low values are generally regarded as errors. For example, the border of the original images have a DN value of 0 and appear as extremely high bars at the leftmost part of the respective histograms; these values have been excluded.

Histograms of the distribution of pixel values. Top row: atmospherically corrected values in radiance/reflectance values; bottom row: original pixel values in DN (digital numbers)

Lawn grass spectral values

Average spectral values extracted from bands 2, 3 and 4 for pixels without lawn grass and for pixels with lawn grass. The software module used was IDRISI’s EXTRACT, with summary type set to ‘average’. The feature definition image was the lawn grass raster binary image. The results are presented in Table 2, where the anomaly is calculated using the following formula:

n=G1G0SS

where: n the anomaly in %, G1 the average spectral values for pixels with grass, G0 the average spectral values for pixels without grass and S the spectral signature of lawn grass for the specific band.

BandNo grass “0”Grass “1”DifferenceSignatureAnomaly (%)
B20,0005340,0034500,0029160,07830066–96,28%
B30,0188330,0537610,0349280,070615–50,54%
B40,1434220,6316460,4882240,6998996–30,24%

Conclusions

The ATMOSC module managed to produce a clearer image with darker deep waters and clearer differences between vegetation and buildup areas. However, it was unable to remove “bad data” errors, such as the line stripping particularly visible over the Atlantic Ocean.

There is a significant difference between the observed values and the signature values. In fact, the observed values for band 2 are almost half of the ones we should have observed. It is worth considering whether a perfectly atmospherically corrected image would provide us with the actual spectral signature values of a given land cover type. The answer would be no, because there are a number of factors that affect it; in the case of grass, insufficient water availability, for example.