4. Radiometric correction of satellite images: when and why radiometric correction is necessary
Author: Dr Alasdair J.Edwards, University of Newcastle, UK.
Aim of Lesson
To develop your understanding of concepts underlying radiometric correction and how to carry out the radiometric correction of satellite imagery, using two Landsat Thematic Mapper images obtained at different seasons under different atmospheric conditions as examples.
Learning Objectives
- To understand the difference between DN values, radiance, irradiance and reflectance.
- To understand why radiometric correction of imagery is required if (i) you are mapping changes in habitat or other features, or (ii) you are using more than one image in a study.
- To demonstrate how different atmospheric conditions can affect DN values by comparing these in two Landsat TM images acquired during different seasons.
- To understand the basic concepts behind atmospheric correction algorithms.
- To learn how to carry out the process of radiometric correction of two Landsat TM images obtained in different seasons and compare the resultant reflectance values.
How to download the lesson
Why you need radiometric correction
Digital sensors record the intensity of electromagnetic radiation (ER) from each spot viewed on the Earth's surface as a digital number (DN) for each spectral band. The DN values recorded by a sensor are proportional to upwelling ER (spectral radiance), the true units of which are W m-2 sr1 μm-1.
The spectral signature of a habitat (say mangrove) is not transferable if measured in digital numbers. The values are image specific; that is, they are dependent on the viewing geometry of the satellite, sun angle, Earth-Sun distance (which varies seasonally), specific weather conditions, and so on, at the moment the image was taken. It is in general far more useful to convert the DN values to spectral radiances.
However, the spectral radiances obtained from the calibration only account for the spectral radiance measured at the satellite sensor. By the time ER is recorded by a satellite or airborne sensor, it has already passed through the Earth's atmosphere twice (sun to target and target to sensor - see figure 1).
Atmospheric correction uses models (e.g. 5S and 6S models) to attempt to account for the effects of the atmosphere on the signal recorded so that spectral radiances at the satellite can be converted to at surface reflectances.
Lesson overview
This lesson uses two Landsat-5 Thematic Mapper images of the Turks and Caicos Islands, one taken during the northern winter (November) when conditions were relatively clear (see colour composite to right), and one taken during the summer (June) when conditions were much more hazy.
Initially, you compare the raw DN values for 5 pixels (in different habitats) between the two images in band #2 (green) and discover that the average difference between the DN values is almost 30% of their mean value. After radiometric and atmospheric correction of the two images, you find that the average difference between the surface reflectances of the same 5 pixels is less than 1% of their mean value.
The radiometric correction process is broken down into three manageable steps (figure 3), each of which is explained in detail before you are asked to carry it out using a formula document For each step you are given a template formula document with a few lines that you need to complete (see example below in figure 5).
The first step is to use metadata provided with the image to convert the DN values to spectral radiances at the satellite using a formula document.
The second step converts the spectral radiance to top-of-atmosphere reflectances using data on the sun-angle at the time the image was taken (part of image metadata - figure 4), mean solar exoatmospheric irradiance (amount of light from the Sun reaching the satellite), and information on how to calculate the Earth-Sun distance for the date of acquisition.
The third step uses the outputs of a 5S (Simulation of the Satellite Signal in the Solar Spectrum) model to remove the effects of atmospheric scattering and absorption and derive surface reflectances for each pixel in each waveband.