Remote Sensing of Environment 113 (2009) 893903

Contents lists available at ScienceDirect

Remote Sensing of Environment

Summary of current radiometric calibration coefcients for Landsat MSS, TM, ETM+, and EO-1 ALI sensors
Gyanesh Chander a,, Brian L. Markham b, Dennis L. Helder c
a b c

SGT, Inc. contractor to the U.S. Geological Survey (USGS) Earth Resources Observation and Science (EROS) Center, Sioux Falls, SD 57198-0001, USA National Aeronautics and Space Administration (NASA) Goddard Space Flight Center (GSFC), Greenbelt, MD 20771, USA South Dakota State University (SDSU), Brookings, SD 57007, USA

a r t i c l e

i n f o

a b s t r a c t
This paper provides a summary of the current equations and rescaling factors for converting calibrated Digital Numbers (DNs) to absolute units of at-sensor spectral radiance, Top-Of-Atmosphere (TOA) reectance, and at-sensor brightness temperature. It tabulates the necessary constants for the Multispectral Scanner (MSS), Thematic Mapper (TM), Enhanced Thematic Mapper Plus (ETM+), and Advanced Land Imager (ALI) sensors. These conversions provide a basis for standardized comparison of data in a single scene or between images acquired on different dates or by different sensors. This paper forms a needed guide for Landsat data users who now have access to the entire Landsat archive at no cost. 2009 Elsevier Inc. All rights reserved.

Article history: Received 10 December 2008 Received in revised form 16 January 2009 Accepted 17 January 2009 Keywords: Landsat MSS, TM, ETM+ EO-1 ALI Radiometric characterization & calibration At-sensor spectral radiance Top-of-atmosphere reectance At-sensor brightness temperature

1. Introduction The Landsat series of satellites provides the longest continuous record of satellite-based observations. As such, Landsat is an invaluable resource for monitoring global change and is a primary source of medium spatial resolution Earth observations used in decision-making (Fuller et al., 1994; Townshend et al.,1995; Goward & Williams,1997; Vogelmann et al., 2001; Woodcock et al., 2001; Cohen & Goward, 2004; Goward et al., 2006; Masek et al., 2008; Wulder et al., 2008). To meet observation requirements at a scale revealing both natural and human-induced landscape changes, Landsat provides the only inventory of the global land surface over time on a seasonal basis (Special issues on Landsat, 1984, 1985, 1997, 2001, 2003, 2004, 2006). The Landsat Program began in early 1972 with the launch of the rst satellite in the series. As technological capabilities increased, so did the amount and quality of image data captured by the various sensors onboard the satellites. Table 1 presents general information about each Landsat satellite. Landsat satellites can be classied into three groups, based on sensor and platform characteristics. The rst group consists of Landsat 1 (L1), Landsat 2 (L2), and Landsat 3 (L3), with the Multispectral Scanner (MSS) sensor and the Return Beam Vidicon (RBV) camera as payloads on a NIMBUS-like platform. The spatial resolution of the

Corresponding author. Tel.: +1 605 594 2554. E-mail address: [email protected] (G. Chander). 1 Work performed under USGS contract 08HQCN0005. 0034-4257/$ see front matter 2009 Elsevier Inc. All rights reserved. doi:10.1016/j.rse.2009.01.007

MSS sensor was approximately 79 m (but often processed to pixel size of 60 m), with four bands ranging from the visible blue to the NearInfrared (NIR) wavelengths. The MSS sensor on L3 included a fth band in the thermal infrared wavelength, with a spectral range from 10.4 to 12.6 m. The L1L3 MSS sensors used a band-naming convention of MSS-4, MSS-5, MSS-6, and MSS-7 for the blue, green, red, and NIR bands, respectively (Markham & Barker, 1983). This designation is obsolete, and to be consistent with the TM and ETM+ sensors, the MSS bands are referred to here as Bands 14, respectively. The second group includes Landsat 4 (L4) and Landsat 5 (L5), which carry the Thematic Mapper (TM) sensor, as well as the MSS, on the Multimission Modular Spacecraft. This second generation of Landsat satellites marked a signicant advance in remote sensing through the addition of a more sophisticated sensor, improved acquisition and transmission of data, and more rapid data processing at a highly automated processing facility. The MSS sensor was included to provide continuity with the earlier Landsat missions, but TM data quickly became the primary source of information used from these satellites because the data offered enhanced spatial, spectral, radiometric, and geometric performance over data from the MSS sensor. The TM sensor has a spatial resolution of 30 m for the six reective bands and 120 m for the thermal band. Because there are no onboard recorders on these sensors, acquisitions are limited to real-time downlink only. The third group consists of Landsat 6 (L6) and Landsat 7 (L7), which include the Enhanced Thematic Mapper (ETM) and the Enhanced Thematic Mapper Plus (ETM+) sensors, respectively. No MSS sensors were included on either satellite. Landsat 6 failed on

894 Table 1 General information about each Landsat satellite. Satellite Landsat Landsat Landsat Landsat Landsat Landsat Landsat EO-1 1 2 3 4 5 6 7 Sensors MSS and MSS and MSS and MSS and MSS and ETM ETM+ ALI RBV RBV RBV TM TM Launch date

G. Chander et al. / Remote Sensing of Environment 113 (2009) 893903

Decommission January 7, 1978 February 25, 1982 March 31, 1983 June 30, 2001 Operational Did not achieve orbit Operational Operational

Altitude km 920 920 920 705 705 705 705

Inclination degrees 99.20 99.20 99.20 98.20 98.20 98.20 98.20

Period min 103.34 103.34 103.34 98.20 98.20 98.20 98.20

Repeat cycle days 18 18 18 16 16 16 16

Crossing time (a.m.) 9:30 9:30 9:30 9:45 9:45 10:00 10:01

July 23, 1972 January 22, 1975 March 5, 1978 July 16, 1982 March 1, 1984 October 5, 1993 April 15, 1999 November 21, 2000

launch. The L7 ETM+ sensor has a spatial resolution of 30 m for the six reective bands, 60 m for the thermal band, and includes a panchromatic (pan) band with a 15 m resolution. L7 has a 378 gigabit (Gb) Solid State Recorder (SSR) that can hold 42 min (approximately 100 scenes) of sensor data and 29 h of housekeeping telemetry concurrently (L7 Science Data User's Handbook2). The Advanced Land Imager (ALI) onboard the Earth Observer-1 (EO-1) satellite is a technology demonstration that serves as a prototype for the Landsat Data Continuity Mission (LDCM). The ALI observes the Earth in 10 spectral bands; nine spectral bands have a spatial resolution of 30 m, and a pan band has a spatial resolution of 10 m. The Landsat data archive at the U.S. Geological Survey (USGS) Earth Resources Observation and Science (EROS) Center holds an unequaled 36year record of the Earth's surface and is available at no cost to users via the Internet (Woodcock et al., 2008). Users can access and search the Landsat data archive via the EarthExplorer (EE)3 or Global Visualization Viewer (GloVis)4 web sites. Note that the Landsat scenes collected by locations within the International Ground Station (IGS) network may be available only from the particular station that collected the scene. 2. Purpose Equations and parameters to convert calibrated Digital Numbers (DNs) to physical units, such as at-sensor radiance or Top-Of-Atmosphere (TOA) reectance, have been presented in a sensor-specic manner elsewhere, e.g., MSS (Markham & Barker, 1986, 1987; Helder, 1993), TM (Chander & Markham, 2003; Chander et al., 2007a), ETM+(Handbook2), and ALI (Markham et al., 2004a). This paper, however, tabulates the necessary constants for all of the Landsat sensors in one place dened in a consistent manner and provides a brief overview of the radiometric calibration procedure summarizing the current accuracy of the at-sensor spectral radiances obtained after performing these radiometric conversions on standard data products generated by U.S. ground processing systems. 3. Radiometric calibration procedure The ability to detect and quantify changes in the Earth's environment depends on sensors that can provide calibrated (known accuracy and precision) and consistent measurements of the Earth's surface features through time. The correct interpretation of scientic information from a global, long-term series of remote-sensing products requires the ability to discriminate between product artifacts and changes in the Earth processes being monitored (Roy et al., 2002). Radiometric characterization and calibration is a prerequisite for creating high-quality science data, and consequently, higher-level downstream products.

3.1. MSS sensors Each MSS sensor incorporates an Internal Calibrator (IC) system, consisting of a pair of lamp assemblies (for redundancy) and a rotating shutter wheel. The shutter wheel includes a mirror and a neutral density lter that varies in transmittance with rotation angle. The calibration system output appears as a light pulse at the focal plane that rises rapidly and then decays slowly. This pulse is referred to as the calibration wedge (Markham & Barker, 1987). The radiometric calibration of the MSS sensors is performed in two stages. First, raw data from Bands 13 are decompressed or linearized and rescaled to 7 bits using xed look-up tables. The look-up tables are derived from prelaunch measurements of the compression ampliers. Second, the postlaunch gain and offset for each detector of all four bands are individually calculated by a linear regression of the detector responses to the samples of the in-orbit calibration wedge with the prelaunch radiances for these samples. A reasonable estimate of the overall calibration uncertainty of each MSS sensor at-sensor spectral radiances is 10%, which was the specied accuracy for the sensor (Markham & Barker, 1987). In most cases, the ground processing system must apply an additional step to uncalibrate the MSS data because a number of MSS scenes were archived as radiometrically corrected products. The previously calibrated archived MSS data must be transformed back into raw DNs using the coefcients stored in the data before applying the radiometric calibration procedure. Studies are underway to evaluate the MSS calibration consistency and provide post-calibration adjustments of the MSS sensors so they are consistent over time and consistent between sensors (Helder, 2008). 3.2. TM sensors The TM sensor includes an onboard calibration system called the IC. The IC consists of a black shutter ag, three lamps, a cavity blackbody, and the optical components necessary to get the lamp and blackbody radiance to the focal plane. The lamps are used to calibrate the reective bands, and the blackbody is used to calibrate the thermal band. Historically, the TM radiometric calibration procedure used the detector's response to the IC to determine radiometric gains and offsets on a scene-by-scene basis. Before launch, the effective radiance of each lamp state for each reective band's detector was determined such that each detector's response to the internal lamp was compared to its response to an external calibrated source. The reective band calibration algorithm for in-ight data used a regression of the detector responses against the prelaunch radiances of the eight lamp states. The slope of the regression represented the gain, while the intercept represented the bias. This algorithm assumed that irradiance of the calibration lamps remained constant over time since launch. Any change in response was treated as a change in sensor response, and thus was compensated for during processing. On-orbit data from individual lamps indicated that the lamps were not particularly stable. Because there was no way to validate the lamp radiances once in orbit, the prelaunch measured radiances were the only metrics available for

2 3 4

http://landsathandbook.gsfc.nasa.gov/handbook.html. http://earthexplorer.usgs.gov. http://glovis.usgs.gov.

G. Chander et al. / Remote Sensing of Environment 113 (2009) 893903 Table 2 MSS spectral range, post-calibration dynamic ranges, and mean exoatmospheric solar irradiance (ESUN). MSS sensors (Qcalmin = 0 and Qcalmax = 127) Band Units L1 MSS 1 2 3 4 L2 MSS 1 2 3 4 L3 MSS 1 2 3 4 L4 MSS 1 2 3 4 L5 MSS 1 2 3 4 (NLAPS) 0.4990.597 0.6030.701 0.6940.800 0.8100.989 (NLAPS) 0.4970.598 0.6070.710 0.6970.802 0.8070.990 (NLAPS) 0.4970.593 0.6060.705 0.6930.793 0.8120.979 (NLAPS) 0.4950.605 0.6030.696 0.7010.813 0.8081.023 (NLAPS) 0.4970.607 0.6030.697 0.7040.814 0.8091.036 0.552 0.650 0.759 0.923 3 3 5 3 268 179 148 123 2.086610 1.385830 1.125980 0.944882 3 3 5 3 1824 1570 1249 853.4 0.550 0.650 0.757 0.916 4 4 5 4 238 164 142 116 1.842520 1.259840 1.078740 0.881890 4 4 5 4 1827 1569 1260 866.4 0.545 0.656 0.743 0.896 4 3 3 1 259 179 149 128 2.007870 1.385830 1.149610 1.000000 4 3 3 1 1839 1555 1291 887.9 0.548 0.659 0.750 0.899 8 6 6 3.66667 263 176 152 130.333 2.007870 1.338580 1.149610 0.997373 8 6 6 3.66667 1829 1539 1268 886.6 0.548 0.652 0.747 0.900 0 0 0 0 248 200 176 153 1.952760 1.574800 1.385830 1.204720 0 0 0 0 1823 1559 1276 880.1 Spectral range m Center wavelength LMIN W/(m2 sr m) LMAX Grescale (W/m2 sr m)/DN Brescale W/(m2 sr m) ESUN

895

W/(m2 m)

Note 1: In some cases, the header le may have different rescaling factors than provided here. In these cases, the user should use the header le information that comes with the product. Table 1 (Markham & Barker, 1986, 1987) provides a summary of the band-specic LMIN and LMAX rescaling factors that have been used at different times and by different systems for the ground processing of MSS data.

the regression procedure. Recent studies5 (Thome et al., 1997a, 1997b; Helder et al., 1998; Markham et al., 1998; Teillet et al., 2001, 2004; Chander et al., 2004a) indicate that the regression calibration did not actually represent detector gains for most of the mission. However, the regression procedure was used until 2003 to generate L5 TM data products and is still used to generate L4 TM products. The calibration uncertainties of the L4 TM at-sensor spectral radiances are 10%, which was the specied accuracy for the sensor (GSFC specication, 1981). The L5 TM reective band calibration procedure was updated in 2003 (Chander & Markham, 2003) to remove the dependence on the changing IC lamps. The new calibration gains implemented on May 5, 2003, for the reective bands (15, 7) were based on lifetime radiometric calibration curves derived from the detectors' responses to the IC, cross-calibration with ETM+, and vicarious measurements (Chander et al., 2004a). The gains were further revised on April 2, 2007, based on the detectors' responses to pseudo-invariant desert sites and cross-calibration with ETM+ (Chander et al., 2007a). Although this calibration update applies to all archived and future L5 TM data, the principal improvements in the calibration are for data acquired during the rst eight years of the mission (19841991), where changes in the sensor gain values are as much as 15%. The radiometric scaling coefcients for Bands 1 and 2 for approximately the rst eight years of the mission have also been changed. Along with the revised reective band radiometric calibration on April 2, 2007, an sensor offset correction of 0.092 W/(m2 sr m), or about 0.68 K (at 300 K), was added to all L5 TM thermal band (Band 6) data acquired since April 1999 (Barsi et al., 2007). The L5 TM radiometric calibration uncertainty

of the at-sensor spectral radiances is around 5% and is somewhat worse for early years, when the sensor was changing more rapidly, and better for later years (Helder et al., 2008). The L4 TM reective bands and the thermal band on both the TM sensors continue to be calibrated using the IC. Further updates to improve the thermal band calibration are being investigated, as is the calibration of the L4 TM. 3.3. ETM+ sensor The ETM+ sensor has three onboard calibration devices for the reective bands: a Full Aperture Solar Calibrator (FASC), which is a white painted diffuser panel; a Partial Aperture Solar Calibrator (PASC), which is a set of optics that allows the ETM+ to image the Sun through small holes; and an IC, which consists of two lamps, a blackbody, a shutter, and optics to transfer the energy from the calibration sources to the focal plane. The ETM+ sensor has also been calibrated vicariously using Earth targets such as Railroad Valley (Thome, 2001; Thome et al., 2004) and cross-calibrated with multiple sensors (Teillet et al., 2001, 2006, 2007; Thome et al., 2003; Chander et al., 2004b, 2007b, 2008). The gain trends from the ETM+ sensor are regularly monitored onorbit using the onboard calibrators and vicarious calibration. The calibration uncertainties of ETM+ at-sensor spectral radiances are 5%. ETM+ is the most stable of the Landsat sensors, changing by no more than 0.5% per year in its radiometric calibration (Markham et al., 2004b). The ETM+ radiometric calibration procedure uses prelaunch gain coefcients populated in the Calibration Parameter File (CPF). These CPFs, issued quarterly, have both an effective and version date. The effective date of the CPF must match the acquisition date of the scene. A CPF version is active until a new CPF for that date period supersedes it. Data can be processed with any version of a CPF; the later versions have more rened parameters, as they reect more data-rich post-acquisition analysis. The ETM+ images are acquired in either a low-or high-gain state. The goal of using two gain settings is to maximize the sensors' 8-bit radiometric resolution without saturating the detectors. For all bands,

5 Radiometric performance studies of the TM sensors have also led to a detailed understanding of several image artifacts due to particular sensor characteristics (Helder & Ruggles, 2004). These artifact corrections (such as Scan-Correlated Shift [SCS], Memory Effect [ME], and Coherent Noise [CN]), along with detector-to-detector normalization (Helder et al., 2004), are necessary to maintain the internal consistency of the calibration within a scene.

896

G. Chander et al. / Remote Sensing of Environment 113 (2009) 893903

Table 3 TM spectral range, post-calibration dynamic ranges, and mean exoatmospheric solar irradiance (ESUN).

Note 1: The Q calmin = 0 for data processed using NLAPS. The Q calmin = 1 for data processed using LPGS. Note 2: The LMIN is typically set to a small negative number, so a zero radiance target will be scaled to a small positive DN value, even in the presence of sensor noise (typically 1 DN or less [1 sigma]). This value is usually not changed throughout the mission. Note 3: In mid-2009, the processing of L4 TM data will transition from NLAPS to LPGS. NLAPS used IC-based calibration. The L4 TM data processed by LPGS will be radiometrically calibrated using a new lifetime gain model procedure and revised calibration parameters. Use the header le information that comes with the product and the above rescaling factors will not be applicable. The numbers highlighted in grey are the revised (LMAX =163) post-calibration dynamic ranges for L4 TM Band 1 data acquired between July 16, 1982 (launch), and August 23, 1986. Note 4: The radiometric scaling coefcients for L5 TM Bands 1 and 2 for approximately the rst eight years (19841991) of the mission were changed to optimize the dynamic range and better preserve the sensitivity of the early mission data. The numbers highlighted in grey are the revised (LMAX =169, 333) post-calibration dynamic ranges for L5 TM Band 1 and 2 data acquired between March 1, 1984 (launch), and December 31, 1991 (Chander et al., 2007a).

the low-gain dynamic range is approximately 1.5 times the high-gain dynamic range. Therefore, low-gain mode is used to image surfaces with high brightness (higher dynamic range but low sensitivity), and high-gain mode is used to image surfaces with low brightness (lower dynamic range but high sensitivity).

All of the ETM+ acquisitions after May 31, 2003, have an anomaly caused by the failure of the Scan Line Corrector (SLC), which compensated for the forward motion of the spacecraft so that all the scans were aligned parallel with each other. The images with data loss are referred to as SLC-off images, whereas images collected prior to

Table 4 ETM+ spectral range, post-calibration dynamic ranges, and mean exoatmospheric solar irradiance (ESUN). L7 ETM+ Sensor (Qcalmin = 1 and Qcalmax = 255) Band Units Spectral range m 0.483 0.560 0.662 0.835 1.648 11.335 2.206 0.706 Center wavelength LMIN W/(m2 sr m) 6.2 6.4 5.0 5.1 1.0 0.0 0.35 4.7 293.7 300.9 234.4 241.1 47.57 17.04 16.54 243.1 LMAX Grescale (W/m2 sr m)/DN 1.180709 1.209843 0.942520 0.969291 0.191220 0.067087 0.066496 0.975591 Brescale W/(m2 sr m) 7.38 7.61 5.94 6.07 1.19 0.07 0.42 5.68 ESUN W/(m2 m) 1997 1812 1533 1039 230.8 N/A 84.90 1362

Low gain (LPGS) 1 0.4520.514 2 0.5190.601 3 0.6310.692 4 0.7720.898 5 1.5471.748 6 10.3112.36 7 2.0652.346 PAN 0.5150.896 High Gain (LPGS) 1 2 3 4 5 6 7 PAN

0.4520.514 0.5190.601 0.6310.692 0.7720.898 1.5471.748 10.3112.36 2.0652.346 0.5150.896

0.483 0.560 0.662 0.835 1.648 11.335 2.206 0.706

6.2 6.4 5.0 5.1 1.0 3.2 0.35 4.7

191.6 196.5 152.9 157.4 31.06 12.65 10.80 158.3

0.778740 0.798819 0.621654 0.639764 0.126220 0.037205 0.043898 0.641732

6.98 7.20 5.62 5.74 1.13 3.16 0.39 5.34

1997 1812 1533 1039 230.8 N/A 84.90 1362

G. Chander et al. / Remote Sensing of Environment 113 (2009) 893903 Table 5 ALI spectral range, post-calibration dynamic ranges, and mean exoatmospheric solar irradiance (ESUN). EO-1 ALI Sensor (Qcalmin = 1 and Qcalmax = 32767) Band Units PAN 1P 1 2 3 4 4P 5P 5 7 Spectral range m 0.4800.690 0.4330.453 0.4500.515 0.5250.605 0.6330.690 0.7750.805 0.8450.890 1.2001.300 1.5501.750 2.0802.350 0.585 0.443 0.483 0.565 0.662 0.790 0.868 1.250 1.650 2.215 Center wavelength LMIN W/(m2 sr m) 2.18 3.36 4.36 1.87 1.28 0.84 0.641 1.29 0.597 0.209 784.2 1471 1405 915.5 588.5 359.6 297.5 270.7 91.14 29.61 LMAX Grescale (W/m2 sr m)/DN 0.024 0.045 0.043 0.028 0.018 0.011 0.0091 0.0083 0.0028 0.00091 Brescale W/(m2 sr m) 2.2 3.4 4.4 1.9 1.3 0.85 0.65 1.3 0.6 0.21 ESUN

897

W/(m2 m) 1724 1857 1996 1807 1536 1145 955.8 452.3 235.1 82.38

All EO-1 ALI standard Level 1 products are processed through the EO-1 Product Generation System (EPGS).

the SLC failure are referred to as SLC-on images (i.e., no data gaps exist). The malfunction of the SLC mirror assembly resulted in the loss of approximately 22% of the normal scene area (Storey et al., 2005). The missing data affects most of the image, with scan gaps varying in width from one pixel or less near the center of the image to 14 pixels along the east and west edges of the image, creating a repeating wedge-shaped pattern along the edges. The middle of the scene, approximately 22 km wide on a Level 1 product, contains very little duplication or data loss. Note that the SLC failure has no impact on the radiometric performance with the valid pixels. 3.4. ALI sensor The ALI has two onboard radiometric calibration devices: a lampbased system and a solar-diffuser with variable irradiance controlled by an aperture door. In addition to its onboard calibrators, ALI has the ability to collect lunar and stellar observations for calibration purposes. The ALI radiometric calibration procedure uses a xed set of detector-by-detector gains established shortly after launch and biases measured shortly after each scene acquisition by closing the ALI's shutter. The calibration uncertainties of the ALI at-sensor spectral radiances are 5% (Mendenhall & Lencioni, 2002). The ALI sensor is well-behaved and stable, with changes in the response being less than 2% per year even early in the mission, and averaging, at most, slightly more than 1% per year over the full mission (Markham et al., 2006). 4. Conversion to at-sensor spectral radiance (Qcal-to-L) Calculation of at-sensor spectral radiance is the fundamental step in converting image data from multiple sensors and platforms into a physically meaningful common radiometric scale. Radiometric calibration of the MSS, TM, ETM+, and ALI sensors involves rescaling the raw digital numbers (Q) transmitted from the satellite to calibrated digital numbers (Qcal)6, which have the same radiometric scaling for all scenes processed on the ground for a specic period. During radiometric calibration, pixel values (Q) from raw, unprocessed image data are converted to units of absolute spectral radiance using 32-bit oating-point calculations. The absolute radiance values are then scaled to 7-bit (MSS, Qcalmax = 127), 8-bit (TM and ETM+, Qcalmax = 255), and 16-bit (ALI, Qcalmax = 32767) numbers representing Qcal before output to distribution media. Conversion from Qcal in Level 1 products back to at-sensor spectral radiance (L) requires knowledge of the lower and upper limit of the original rescaling factors. The following equation is used to perform the Qcal-to-L conversion for Level 1 products:

L = or

  LMAX LMIN Qcal Qcalmin + LMIN Qcalmax Qcalmin

L = Grescale Qcal + Brescale where : Grescale = LMAX LMIN Qcalmax Qcalmin   LMAX LMIN Qcalmin Qcalmax Qcalmin 1

Brescale = LMIN where

Spectral radiance at the sensor's aperture [W/(m2 sr m)] L= Qcal= Quantized calibrated pixel value [DN] Qcalmin= Minimum quantized calibrated pixel value corresponding to LMIN [DN] Qcalmax= Maximum quantized calibrated pixel value corresponding to LMAX [DN] LMIN= Spectral at-sensor radiance that is scaled to Qcalmin [W/(m2 sr m)] LMAX= Spectral at-sensor radiance that is scaled to Qcalmax [W/(m2 sr m)] Grescale= Band-specic rescaling gain factor [(W/(m2 sr m))/DN] Brescale= Band-specic rescaling bias factor [W/(m2 sr m)] Historically, the MSS and TM calibration information is presented in spectral radiance units of mW/(cm2 sr m). To maintain consistency with ETM+ spectral radiance, units of W/(m2 sr m) are now used for MSS and TM calibration information. The conversion factor is 1:10 when converting from mW/(cm2 sr m) units to W/(m2 sr m). Tables 2, 3, 4, and 5 summarize the spectral range, post-calibration dynamic ranges7 (LMIN and LMAX scaling parameters and the corresponding rescaling gain [Grescale] and rescaling bias [Brescale] values), and mean exoatmospheric solar irradiance (ESUN) for the MSS, TM, ETM+, and ALI sensors, respectively.

These are the DNs that users receive with Level 1 Landsat products.

7 The post-calibration dynamic ranges summarized in Tables 25 are only applicable to Landsat data processed and distributed by the USGS EROS Center. The IGSs may process the data differently, and these rescaling factors may not be applicable. Special collections, such as the Multi-Resolution Land Characteristics Consortium (MRLC) or Global Land Survey (GLS), may have a different processing history, so the user needs to verify the respective product header information.

898

G. Chander et al. / Remote Sensing of Environment 113 (2009) 893903

Fig. 1. Comparison of the solar reective bands RSR proles of L15 MSS sensors.

Tables 25 give the prelaunch measured (as-built performance) spectral ranges. These numbers are slightly different from the original lter specication. The center wavelengths are the average of the two spectral range numbers. Figs.1 and 2 show the Relative Spectral Response (RSR) proles of the Landsat MSS (Markham & Barker, 1983), TM (Markham & Barker, 1985), ETM+(Handbook2), and ALI (Mendenhall & Parker, 1999) sensors measured during prelaunch characterization. The ETM+ spectral bands were designed to mimic the standard TM spectral bands 17. The ALI bands were designed to mimic the six standard ETM+

solar reective spectral bands 15, and 7; three new bands,1p, 4p, and 5p, were added to more effectively address atmospheric interference effects and specic applications. The ALI band numbering corresponds with the ETM+ spectral bands. Bands not present on the ETM+ sensor are given the p, or prime, designation. MSS spectral bands are signicantly different from TM and ETM+ spectral bands. The post-calibration dynamic ranges are band-specic rescaling factors typically provided in the Level 1 product header le. Over the life of the Landsat sensors, occasional changes have occurred in

G. Chander et al. / Remote Sensing of Environment 113 (2009) 893903

899

Fig. 2. Comparison of the solar reective bands RSR proles of L4 TM, L5 TM, L7 ETM+, and EO-1 ALI sensors.

the post-calibration dynamic range. Future changes are anticipated, especially in the MSS and TM data, because of the possible adjustment of the calibration constants based on comparisons to absolute radiometric measurements made on the ground. In some cases, the header le may have different rescaling factors than provided in the table included here. In these cases, the user should use the header le information that comes with the product. Two processing systems will continue to generate Landsat data products: the Level 1 Product Generation System (LPGS) and the National Land Archive Production System (NLAPS). Starting December 8, 2008, all L7 ETM+ and L5 TM (except Thematic Mapper-Archive

[TM-A]8 products) standard Level 1 products are processed through the LPGS, and all L4 TM and MSS standard Level 1 products are processed through the NLAPS. The Landsat Program is working toward transitioning the processing of all Landsat data to LPGS (Kline, personal
8 A small number of TM scenes were archived as radiometrically corrected products known as TM-A data. The TM-A data are archived on a scene-by-scene basis (instead of intervals). The L4 and L5 TM-A scenes will continue to be processed using NLAPS (with Qcalmin = 0), which attempts to uncalibrate the previously applied calibration and generates the product using updated calibration procedures. Note that approximately 80 L4 TM and approximately 13,300 L5 TM scenes are archived as TM-A data, with acquisition dates ranging between Sept.1982 and Aug. 1990.

900

G. Chander et al. / Remote Sensing of Environment 113 (2009) 893903

Fig. 2 (continued).

communication). In mid-2009, the processing of L4 TM data will transition from NLAPS to LPGS. The scenes processed using LPGS include a header le (.MTL), which lists the LMIN and LMAX values but not the rescaling gain and bias numbers. The scenes processed using NLAPS include a processing history work order report (.WO), which lists the rescaling gain and bias numbers but not the LMIN and LMAX. The sensitivity of the detector changes over time, causing a change in the detector gain applied during radiometric calibration. However, the numbers presented in Tables 25 are the rescaling factors, which are the post-calibration dynamic ranges. The LMIN and LMAX are a representation of how the output Landsat Level 1 data products are scaled in at-sensor radiance units. Generally, there is no need to change the LMIN or LMAX unless something changes drastically on the sensor. Thus, there is no time dependence for any of the rescaling factors in Tables 25. 5. Conversion to TOA reectance (L-to-P) A reduction in scene-to-scene variability can be achieved by converting the at-sensor spectral radiance to exoatmospheric TOA reectance, also known as in-band planetary albedo. When comparing images from different sensors, there are three advantages to using TOA reectance instead of at-sensor spectral radiance. First, it removes the cosine effect of different solar zenith angles due to the time difference between data acquisitions. Second, TOA reectance compensates for different values of the exoatmospheric solar irradiance arising from spectral band differences. Third, the TOA reectance corrects for the variation in the EarthSun distance between different data acquisition dates. These variations can be signicant geographically and temporally. The TOA reectance of the Earth is computed according to the equation: = L d 2 ESUN cos s 2

scene center is typically stored in the Level 1 product header le (.MTL or .WO) or retrieved from the USGS EarthExplorer or GloVis online interfaces under the respective scene metadata (these web sites also contain the acquisition time in hours, minutes, and seconds). The TOA reectance calculation requires the EarthSun distance (d). Table 6 presents d in astronomical units throughout a year generated using the Jet Propulsion Laboratory (JPL) Ephemeris10 (DE405) data. The d numbers are also tabulated in the Nautical Almanac Ofce. The last column of Tables 25 summarizes solar exoatmospheric spectral irradiances (ESUN) for the MSS, TM, ETM+, and ALI sensors using the Thuillier solar spectrum (Thuillier et al., 2003). The Committee on Earth Observation Satellites (CEOS) Working Group on Calibration and Validation (WGCV) recommends11 using this spectrum for applications in optical-based Earth Observation that use an exoatmospheric solar irradiance spectrum. The Thuillier spectrum is believed to be the most accurate and an improvement over the other solar spectrum. Note that the CHKUR solar spectrum in MODTRAN 4.0 (Air Force Laboratory, 1998) was used previously for ETM+(Handbook2) and TM (Chander & Markham, 2003), whereas the Neckel and Lab (Neckel & Labs, 1984) and Iqbal (Iqbal, 1983) solar spectrums were used for MSS and TM solar irradiance values (Markham & Barker,1986). The primary differences occur in Bands 5 and 7. For comparisons to other sensors, users need to verify that the same solar spectrum is used for all sensors. 6. Conversion to at-sensor brightness temperature (L-to-T) The thermal band data (Band 6 on TM and ETM+) can be converted from at-sensor spectral radiance to effective at-sensor brightness temperature. The at-sensor brightness temperature assumes that the Earth's surface is a black body (i.e., spectral emissivity is 1), and includes atmospheric effects (absorption and emissions along path). The at-sensor temperature uses the prelaunch calibration constants given in Table 7. The conversion formula from the at-sensor's spectral radiance to at-sensor brightness temperature is: T= K2   1 ln K L + 1

where Planetary TOA reectance [unitless] = = Mathematical constant equal to ~ 3.14159 [unitless] L= Spectral radiance at the sensor's aperture [W/(m2 sr m)] d= EarthSun distance [astronomical units] ESUN= Mean exoatmospheric solar irradiance [W/(m2 m)] s= Solar zenith angle [degrees9] Note that the cosine of the solar zenith angle is equal to the sine of the solar elevation angle. The solar elevation angle at the Landsat
9 Note that Excel, Matlab, C, and many other software applications use radians, not degrees, to perform calculations. The conversion from degrees to radians is a multiplication factor of pi/180.

where: T= K2= Effective at-sensor brightness temperature [K] Calibration constant 2 [K]

10 11

http://ssd.jpl.nasa.gov/?horizons. CEOS-recommended solar irradiance spectrum, http://wgcv.ceos.org.

G. Chander et al. / Remote Sensing of Environment 113 (2009) 893903 Table 6 EarthSun distance (d) in astronomical units for Day of the Year (DOY). DOY 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60 d 0.98331 0.98330 0.98330 0.98330 0.98330 0.98332 0.98333 0.98335 0.98338 0.98341 0.98345 0.98349 0.98354 0.98359 0.98365 0.98371 0.98378 0.98385 0.98393 0.98401 0.98410 0.98419 0.98428 0.98439 0.98449 0.98460 0.98472 0.98484 0.98496 0.98509 0.98523 0.98536 0.98551 0.98565 0.98580 0.98596 0.98612 0.98628 0.98645 0.98662 0.98680 0.98698 0.98717 0.98735 0.98755 0.98774 0.98794 0.98814 0.98835 0.98856 0.98877 0.98899 0.98921 0.98944 0.98966 0.98989 0.99012 0.99036 0.99060 0.99084 DOY 61 62 63 64 65 66 67 68 69 70 71 72 73 74 75 76 77 78 79 80 81 82 83 84 85 86 87 88 89 90 91 92 93 94 95 96 97 98 99 100 101 102 103 104 105 106 107 108 109 110 111 112 113 114 115 116 117 118 119 120 d 0.99108 0.99133 0.99158 0.99183 0.99208 0.99234 0.99260 0.99286 0.99312 0.99339 0.99365 0.99392 0.99419 0.99446 0.99474 0.99501 0.99529 0.99556 0.99584 0.99612 0.99640 0.99669 0.99697 0.99725 0.99754 0.99782 0.99811 0.99840 0.99868 0.99897 0.99926 0.99954 0.99983 1.00012 1.00041 1.00069 1.00098 1.00127 1.00155 1.00184 1.00212 1.00240 1.00269 1.00297 1.00325 1.00353 1.00381 1.00409 1.00437 1.00464 1.00492 1.00519 1.00546 1.00573 1.00600 1.00626 1.00653 1.00679 1.00705 1.00731 DOY 121 122 123 124 125 126 127 128 129 130 131 132 133 134 135 136 137 138 139 140 141 142 143 144 145 146 147 148 149 150 151 152 153 154 155 156 157 158 159 160 161 162 163 164 165 166 167 168 169 170 171 172 173 174 175 176 177 178 179 180 d 1.00756 1.00781 1.00806 1.00831 1.00856 1.00880 1.00904 1.00928 1.00952 1.00975 1.00998 1.01020 1.01043 1.01065 1.01087 1.01108 1.01129 1.01150 1.01170 1.01191 1.01210 1.01230 1.01249 1.01267 1.01286 1.01304 1.01321 1.01338 1.01355 1.01371 1.01387 1.01403 1.01418 1.01433 1.01447 1.01461 1.01475 1.01488 1.01500 1.01513 1.01524 1.01536 1.01547 1.01557 1.01567 1.01577 1.01586 1.01595 1.01603 1.01610 1.01618 1.01625 1.01631 1.01637 1.01642 1.01647 1.01652 1.01656 1.01659 1.01662 DOY 181 182 183 184 185 186 187 188 189 190 191 192 193 194 195 196 197 198 199 200 201 202 203 204 205 206 207 208 209 210 211 212 213 214 215 216 217 218 219 220 221 222 223 224 225 226 227 228 229 230 231 232 233 234 235 236 237 238 239 240 d 1.01665 1.01667 1.01668 1.01670 1.01670 1.01670 1.01670 1.01669 1.01668 1.01666 1.01664 1.01661 1.01658 1.01655 1.01650 1.01646 1.01641 1.01635 1.01629 1.01623 1.01616 1.01609 1.01601 1.01592 1.01584 1.01575 1.01565 1.01555 1.01544 1.01533 1.01522 1.01510 1.01497 1.01485 1.01471 1.01458 1.01444 1.01429 1.01414 1.01399 1.01383 1.01367 1.01351 1.01334 1.01317 1.01299 1.01281 1.01263 1.01244 1.01225 1.01205 1.01186 1.01165 1.01145 1.01124 1.01103 1.01081 1.01060 1.01037 1.01015 DOY 241 242 243 244 245 246 247 248 249 250 251 252 253 254 255 256 257 258 259 260 261 262 263 264 265 266 267 268 269 270 271 272 273 274 275 276 277 278 279 280 281 282 283 284 285 286 287 288 289 290 291 292 293 294 295 296 297 298 299 300 d 1.00992 1.00969 1.00946 1.00922 1.00898 1.00874 1.00850 1.00825 1.00800 1.00775 1.00750 1.00724 1.00698 1.00672 1.00646 1.00620 1.00593 1.00566 1.00539 1.00512 1.00485 1.00457 1.00430 1.00402 1.00374 1.00346 1.00318 1.00290 1.00262 1.00234 1.00205 1.00177 1.00148 1.00119 1.00091 1.00062 1.00033 1.00005 0.99976 0.99947 0.99918 0.99890 0.99861 0.99832 0.99804 0.99775 0.99747 0.99718 0.99690 0.99662 0.99634 0.99605 0.99577 0.99550 0.99522 0.99494 0.99467 0.99440 0.99412 0.99385 DOY 301 302 303 304 305 306 307 308 309 310 311 312 313 314 315 316 317 318 319 320 321 322 323 324 325 326 327 328 329 330 331 332 333 334 335 336 337 338 339 340 341 342 343 344 345 346 347 348 349 350 351 352 353 354 355 356 357 358 359 360 361 362 363 364 365 366 d

901

0.99359 0.99332 0.99306 0.99279 0.99253 0.99228 0.99202 0.99177 0.99152 0.99127 0.99102 0.99078 0.99054 0.99030 0.99007 0.98983 0.98961 0.98938 0.98916 0.98894 0.98872 0.98851 0.98830 0.98809 0.98789 0.98769 0.98750 0.98731 0.98712 0.98694 0.98676 0.98658 0.98641 0.98624 0.98608 0.98592 0.98577 0.98562 0.98547 0.98533 0.98519 0.98506 0.98493 0.98481 0.98469 0.98457 0.98446 0.98436 0.98426 0.98416 0.98407 0.98399 0.98391 0.98383 0.98376 0.98370 0.98363 0.98358 0.98353 0.98348 0.98344 0.98340 0.98337 0.98335 0.98333 0.98331

K1= L= ln=

Calibration constant 1 [W/(m2 sr m)] Spectral radiance at the sensor's aperture [W/(m2 sr m)] Natural logarithm

The ETM+ Level 1 product has two thermal bands, one acquired using a low gain setting (often referred to as Band 6 L; useful

temperature range of 130350 K) and the other using a high gain setting (often referred to as Band 6H; useful temperature range of 240320 K). The noise equivalent change in temperature (NET) at 280 K for ETM+ high gain is 0.22 and for low gain is 0.28. The TM Level 1 product has only one thermal band (there is no gain setting on the TM sensor), and the thermal band images have a useful temperature

902 Table 7 TM and ETM+ thermal band calibration constants. Constant Units L4 TM L5 TM L7 ETM+ K1 W/(m2 sr m) 671.62 607.76 666.09

G. Chander et al. / Remote Sensing of Environment 113 (2009) 893903 Table A2 Standard Level 1 product specications. K2 Kelvin 1284.30 1260.56 1282.71 Product Type Level 1T (Terrain Corrected) Pixel Size 15/30/60 meters Output format GeoTIFF Resampling Method Cubic Convolution (CC) Map Projection Universal Transverse Mercator (UTM) Polar Stereographic for Antarctica Image Orientation Map (North Up) Distribution File Transfer Protocol (FTP) Download only

range of 200340 K. The NET at 280 K for L5 TM is 0.170.30 (Barsi et al., 2003). 7. Conclusion This paper provides equations and rescaling factors for converting Landsat calibrated DNs to absolute units of at-sensor spectral radiance, TOA reectance, and at-sensor brightness temperature. It tabulates the necessary constants for the MSS, TM, ETM+, and ALI sensors in a coherent manner using the same units and denitions. This paper forms a needed guide for Landsat data users who now have access to the entire Landsat archive at no cost. Studies are ongoing to evaluate the MSS calibration consistency and provide post-calibration adjustments of the MSS sensors so they are consistent over time and consistent between sensors. Further updates to improve the TM and ETM+ thermal band calibration are being investigated, as is the calibration of the L4 TM. Acknowledgments This work was partially supported by the NASA Land-Cover and Land-Use Change (LCLUC) Grant NNH08AI30I. The authors acknowledge the support of David Aaron (SDSU) for digitizing the spectral responses for the MSS sensors and Julia Barsi (SSAI) for generating the EarthSun distance. Special thanks are extended to a number of colleagues for reviewing drafts of this manuscript. The anonymous reviewers' comments were particularly valuable and their efforts are greatly appreciated. Any use of trade, product, or rm names is for descriptive purposes only and does not imply endorsement by the U.S. Government.

References
Air Force Research Laboratory (1998). Modtran Users Manual, Versions 3.7 and 4.0 : Hanscom Air Force Base, MA. Barsi, J. A., Schott, J. R., Palluconi, F. D., Helder, D. L., Hook, S. J., Markham, B. L., Chander, G., & O'Donnell, E. M. (2003). Landsat TM and ETM+ thermal band calibration. Canadian Journal of Remote Sensing, 29(2), 141153. Barsi, J. A., Hook, S. J., Schott, J. R., Raqueno, N. G., & Markham, B. L. (2007). Landsat-5 Thematic Mapper thermal band calibration update. IEEE Transactions on Geoscience and Remote Sensing, 44, 552555. Chander, G., & Markham, B. L. (2003). Revised Landsat-5 TM radiometric calibration procedures, and post-calibration dynamic ranges. IEEE Transactions on Geoscience and Remote Sensing, 41, 26742677. Chander, G., Helder, D. L., Markham, B. L., Dewald, J., Kaita, E., Thome, K. J., Micijevic, E., & Ruggles, T. (2004a). Landsat 5 TM on-orbit absolute radiometric performance. IEEE Transactions on Geoscience and Remote Sensing, 42, 27472760. Chander, G., Meyer, D. J., & Helder, D. L. (2004b). Cross-calibration of the Landsat-7 ETM+ and EO-1 ALI sensors. IEEE Transactions on Geoscience and Remote Sensing, 42(12), 28212831. Chander, G., Markham, B. L., & Barsi, J. A. (2007a). Revised Landsat 5 Thematic Mapper radiometric calibration. IEEE Transactions on Geoscience and Remote Sensing, 44, 490494. Chander, G., Angal, A., Choi, T., Meyer, D. J., Xiong, X., & Teillet, P. M. (2007b). Crosscalibration of the Terra MODIS, Landsat-7 ETM+ and EO-1 ALI sensors using near simultaneous surface observation over Railroad Valley Playa, Nevada test site. In J. J. Butler & J. Xiong (Eds.), Proceedings of SPIE Conference 6677 on Earth Observing Systems XII, SPIE, Vol. 6677. (pp. 66770Y: 112) San Diego, CA. Chander, G., Coan, M. J., & Scaramuzza, P. L. (2008). Evaluation and comparison of the IRS-P6 and the Landsat Sensors. IEEE Transactions on Geoscience and Remote Sensing, 46(1), 209221. Cohen, W. B., & Goward, S. N. (2004). Landsat's role in ecological applications of remote sensing. BioScience, 54(6), 535545. Fuller, R. M., Groom, G. B., & Jones, A. R. (1994). The land cover map of Great Britain: An automated classication of Landsat Thematic Mapper data. Photogrammetric Engineering and Remote Sensing, 60, 553562. Goward, S. N., & Williams, D. L. (1997). Landsat and Earth Systems Science: Development of terrestrial monitoring. Photogrammetric Engineering and Remote Sensing, 63(7), 887900. Goward, S., Irons, J., Franks, S., Arvidson, T., Williams, D., & Faundeen, J. (2006). Historical record of Landsat global coverage: Mission operations, NSLRSDA, and international cooperator stations. Photogrammetric Engineering and Remote Sensing, 72, 11551169. GSFC Specication for the Thematic Mapper System and Associated Test Equipment. (Rev C., January 1981) GSFC 400-8-D-210C, NASA, Goddard Space Flight Center, Greenbelt, MD. Helder, D. L. (1993). MSS radiometric calibration handbook. Report to the U.S. Geological Survey (USGS) Earth Resources Observation and Science (EROS) Center. Helder, D. L., Boncyk, W. C., & Mortt, R. (1998). Absolute calibration of the Landsat Thematic Mapper using the internal calibrator. Proc. IGARSS, Seattle, WA, 1998 (pp. 27162718). Helder, D. L. (2008a). Consistent radiometric calibration of the historical Landsat archive. Proceedings of PECORA Denver, Colorado. Helder, D. L., Markham, B. L., Thome, K. J., Barsi, J. A., Chander, G., & Malla, R. (2008b). Updated radiometric calibration for the Landsat 5 Thematic Mapper reective bands. IEEE Transactions on Geoscience and Remote Sensing, 46(10), 33093325. Helder, D. L., & Ruggles, T. A. (2004a). Landsat Thematic Mapper reectiveband radiometric artifacts. IEEE Transactions on Geoscience and Remote Sensing, 44, 27042716. Helder, D. L., Ruggles, T. A., Dewald, J. D., & Madhavan, S. (2004b). Landsat-5 Thematic Mapper reective-band radiometric stability. IEEE Transactions on Geoscience and Remote Sensing, 44, 27042716. Iqbal, M. (1983). Introduction to solar radiation. New York: Academic Press. Markham, B. L., & Barker, J. L. (1983). Spectral characterization of the Landsat-4 MSS sensors. Photogrammetric Engineering and Remote Sensing, 49(6), 811833. Markham, B. L., & Barker, J. L. (1985). Spectral characterization of the LANDSAT Thematic Mapper sensors'. International Journal of Remote Sensing, 6(5), 697716. Markham, B. L., & Barker, J. L. (1986). Landsat MSS and TM post-calibration dynamic ranges, exoatmospheric reectances and at-satellite temperatures. Earth Observation Satellite Co., Lanham, MD, Landsat Tech. Note 1. Markham, B. L., & Barker, J. L. (1987). Radiometric properties of U.S. processed Landsat MSS Data. Remote Sensing of Environment, 22, 3971. Markham, B. L., Seiferth, J. C., Smid, J., & Barker, J. L. (1998). Lifetime responsivity behavior of the Landsat-5 Thematic Mapper. Proceedings of SPIE Conference 3427 on Earth Observing Systems, SPIE, Vol. 3427. (pp. 420431) San Diego, CA. Markham, B. L., Chander, G., Mortt, R., Hollaren, D., Mendenhall, J. F., & Ong, L. (2004a). Radiometric processing and calibration of EO-1 Advanced Land Imager data. Proceedings of PECORA 16 Global Priorities in Land Remote Sensing South Dakota: Sioux Falls.

Appendix A

Table A1 To maintain consistency, all Landsat scenes are based on the following naming convention. Format Example: LXSPPPRRRYYYYDDDGSIVV L = Landsat X = Sensor S = Satellite PPP = Worldwide Reference System (WRS) Path RRR = WRS Row YYYY = Year DDD = Day of Year GSI=Ground Station Identier Sensor Examples: LM10170391976031AAA01(MSS) LT40170361982320XXX08 (TM) LE70160392004262EDC02 (ETM+)

VV = Version Ground Stations Identiers-Data received at these sites are held at EROS AAA = North American site unknown ASA = Alice Springs, Australia FUI= Fucino, Italy (Historical) GLC = Gilmore Creek, AK, US HOA = Hobart, Australia KIS = Kiruna, Sweden MTI = Matera, Italy EDC=Receiving site unknown PAC = Prince Albert, Canada GNC = Gatineau, Canada LGS = EROS, SD, USA, Landsat 5 data acquired by EROS beginning July 1, 2001 MOR = Moscow, Russia MLK = Malinda, Kenya IKR = Irkutsk, Russia CHM = Chetumal, Mexico XXO = Receiving site unknown XXX = Receiving site unknown

G. Chander et al. / Remote Sensing of Environment 113 (2009) 893903 Markham, B. L., Thome, K., Barsi, J., Kaita, E., Helder, D., Barker, J., & Scaramuzza, P. (2004b). Landsat-7 ETM+ On-Orbit reective-band radiometric stability and absolute calibration. IEEE Transactions on Geoscience and Remote Sensing, 42, 28102820. Markham, B. L., Ong, L., Barsi, J. A., Mendenhall, J. A., Lencioni, D. E., Helder, D. L., Hollaren, D. H., & Mortt, R. M. (2006). Radiometric calibration stability of the EO-1 Advanced Land Imager: 5 Years on-orbit. In R. Meynart, S. P. Neeck, & H. Shimoda (Eds.), Proceedings of SPIE Conference 6361 on Sensors, Systems, and Next-Generation Satellites X, SPIE, Vol. 6361. (pp. 66770U: 112) San Diego, CA. Masek, J. G., Vermote, Huang, C., Wolfe, R., Cohen, W., Hall, F., Kutler, J., & Nelson, P. (2008). North American forest disturbance mapped from a decadal Landsat record. Remote Sensing of Environment, 112, 29142926. Mendenhall, J. A., & Parker, A. C. (1999). Spectral calibration of the EO-1 Advanced Land Imager. Proceedings of SPIE Conference 3750 on Earth Observing Systems IV, SPIE, Vol. 3750. (pp. 109116) Denver, Colorado. Mendenhall, J. A., & Lencioni, D. E. (2002). EO-1 advanced land imager on-orbit radiometric calibration. International Geoscience and Remote Sensing Symposium Toronto, Canada. Nautical Almanac Ofce. The Nautical Almanac for the Year (United States Naval Observatory) (Washington, DC: U.S. Government Printing Ofce). Neckel, H., & Labs, D. (1984). The solar radiation between 3300 and 12500 A. Solar Physics, 90, 205. Roy, D., Borak, J., Devadiga, S., Wolfe, R., Zheng, M., & Descloitres, J. (2002). The MODIS land product quality assessment approach. Remote Sensing of Environment, 83, 6276. Special Issue on Landsat 4 (1984). IEEE Transactions on Geoscience and Remote Sensing, GE-22(3), 160 51(9). Guest Editor: Solomonson, V.V. Special Issue on Landsat Image Data Quality Analysis (LIDQA) (1985). Photogrammetric Engineering and Remote Sensing, 51(9) Guest Editors: Markham, B.L., and Barker, J.L. Special Issue on 25th Anniversary of Landsat (1997). Photogrammetric Engineering and Remote Sensing, 63(7) Guest Editor: Salomonson, V.V. Special Issue on Landsat 7 (2001). Remote Sensing of Environment, 78(12), 1220 Guest Editors: Goward, S.N., and Masek, J.G. Special Issue on Synergistic Utilization of Landsat 7 (2003). Canadian Journal of Remote Sensing, 29(2), 141297 Guest Editor: Teillet, P.M. Special Issue on Landsat Sensor Performance Characterization (2004). IEEE Transactions on Geoscience and Remote Sensing, 42(12), 26872855 Guest Editors: Markham, B.L., Storey, J.C., Crawford, M.M., Goodenough, D.G., Irons, J.R. Special Issue on Landsat Operations: Past, Present and Future (2006). Photogrammetric Engineering and Remote Sensing, 72(10) Guest Editors: Williams, D.L, Goward, S.N., Arvidson, T. Storey, J. C., Scaramuzza, P., & Schmidt, G. (2005). Landsat 7 scan line corrector-off gap lled product development. PECORA 16 Conference Proceedings, Sioux Falls, South Dakota (pp. 2327). Teillet, P. M., Barker, J. L., Markham, B. L., Irish, R. R., Fedosejeves, G., & Storey, J. C. (2001). Radiometric cross-calibration of the Landsat-7 ETM+ and Landsat-5 TM sensors based on tandem data sets. Remote Sensing of Environment, 78(12), 3954.

903

Teillet, P. M., Helder, D. L., Ruggles, T. A., Landry, R., Ahern, F. J., Higgs, N. J., et al. (2004). A denitive calibration record for the Landsat-5 Thematic Mapper anchored to the Landsat-7 radiometric scale. Canadian Journal of Remote Sensing, 30(4), 631643. Teillet, P. M., Markham, B. L., & Irish, R. R. (2006). Landsat cross-calibration based on near simultaneous imaging of common ground targets. Remote Sensing of Environment, 102 (34), 264270. Teillet, P. M., Fedosejevs, G., Thome, K. J., & Barker, J. L. (2007). Impacts of spectral band difference effects on radiometric cross-calibration between satellite sensors in the solar-reective spectral domain. Remote Sensing of Environment, 110, 393409. Thome, K. J. (2001). Absolute radiometric calibration of Landsat 7 ETM+ using the reectance-based method. Remote Sensing of Environment, 78(12), 2738. Thome, K. J., Markham, B. L., Barker, J. L., Slater, P. L., & Biggar, S. (1997a). Radiometric calibration of Landsat. Photogrammetric Engineering and Remote Sensing, 63, 853858. Thome, K. J., Crowther, B. G., & Biggar, S. F. (1997b). Reectance- and irradiance-based calibration of Landsat-5 Thematic Mapper. Canadian Journal of Remote Sensing, 23, 309317. Thome, K. J., Biggar, S. F., & Wisniewski, W. (2003). Cross comparison of EO-1 sensors and other earth resources sensors to Landsat-7 ETM+ using Railroad Valley Playa. IEEE Transactions on Geoscience and Remote Sensing, 41, 11801188. Thome, K. J., Helder, D. L., Aaron, D. A., & Dewald, J. (2004). Landsat 5 TM and Landsat-7 ETM+ absolute radiometric calibration using reectance based method. IEEE Transactions on Geoscience and Remote Sensing, 42(12), 27772785. Thuillier, G., Herse, M., Labs, S., Foujols, T., Peetermans, W., Gillotay, D., Simon, P. C., & Mandel, H. (2003). The solar spectral irradiance from 200 to 2400 nm as measured by SOLSPEC Spectrometer from the ATLAS 123 and EURECA missions.Solar Physics, 214(1), 122 Solar Physics. Townshend, J. R. G., Bell, V., Desch, A. C., Havlicek, C., Justice, C. O., Lawrence, W. E., Skole, D., Chomentowski, W. W., Moore, B., Salas, W., & Tucjer, C. J. (1995). The NASA Landsat Pathnder Humid Tropical Deforestation Project. Proceedings Land Satellite Information in the Next Decade, ASPRS Conference (pp. 7687). Vienna, Virginia. Vogelmann, J. E., Howard, S. M., Yang, L., Larson, C. R., Wylie, B. K., & Van Driel, J. N. (2001). Completion of the 1990's National Land Cover Data Set for the conterminous United States. Photogrammetric Engineering and Remote Sensing, 67, 650662. Woodcock, C. E., Macomber, S. A., Pax-Lenney, M., & Cohen, W. C. (2001). Monitoring large areas for forest change using Landsat: Generalization across space, time and Landsat sensors. Remote Sensing of Environment, 78, 194203. Woodcock, C. E., Allen, A. A., Anderson, M., Belward, A. S., Bindschadler, R., Cohen, W. B., Gao, F., Goward, S. N., Helder, D., Helmer, E., Nemani, R., Oreapoulos, L., Schott, J., Thenkabail, P. S., Vermote, E. F., Vogelmann, J., Wulder, M. A., & Wynne, R. (2008). Free access to Landsat imagery. Science, 320, 1011. Wulder, M. A., White, J. C., Goward, S. N., Masek, J. G., Irons, J. R., Herold, M., Cohen, W. B., Loveland, T. R., & Woodcock, C. E. (2008). Landsat continuity: Issues and opportunities for land cover monitoring. Remote Sensing of Environment, 112, 955969.