11.1 Level 0R Product |
Unlike earlier Landsat programs, the Landsat 7 system was not originally designed to produce high level (i.e. Level 1) products for users. The baselined program philosophy was to provide raw data only which would leave the value added domain for commercial companies. A prevailing "wait and see" position by commercial vendors prompted NASA to add a systematic correction capability to ensure product availability. The primary product for users and vendors seeking higher level processing, however, is 0R data - an essentially raw data form that is marginally useful prior to radiometric and geometric correction. This is readily apparent when viewing a simulated 0R image . A Landsat 7 0R product, however, does contain all the ancillary data required to perform these corrections including a calibration parameter file (CPF) generated by the Landsat 7 IAS.
LPS spatially reformats earth imagery and calibration data into Level 0R data. This involves shifting pixels by integer amounts to account for the alternating forward-reverse scanning pattern of the ETM+ sensor, the odd-even detector arrangement within each band, and the detector offsets inherent to the focal plane array engineering design. All LPs 0R corrections are reversible; the pixel shift parameters used are documented in the IAS CPF.
The LPs 0R output is HDF-EOS formatted and archived. Details of the archival format can be found in the Landsat 7 System Wideband DFCB, Vol. 4.
The distance between WRS center points along a path is 161.1 kilometers. A path distance of 90 kilometers before and after a WRS center point defines the standard scene length of 180 km. This length includes 20 scans of overlap with neighboring scenes. The standard WRS scene overlaps neighboring scenes along a path by 5% at the equator and has a width or cross track distance of 185 kilometers.
Landsat 7 browse is framed according to WRS scenes. An ordered scene will cover the same geographic extent observed in the browse with the following caveat. Standard WRS scenes have 375 scans. Partial scenes (less than 375 scans) may exist at the beginning or end of a subinterval due to the fact that imaging events do not always start and end on scene boundaries. Browse and scene metadata for these occurrences accurately reflect their partial scene nature and geographic extent although partials are currently not offered due to complexities associated with level 1 processing.
A user may order a subset of the available bands which will affect the actual file count in a 0R product. In all cases, however, every product includes two PCD files, two MSCD files, three metadata files, the CPF, and the HDF directory. Only the internal calibrator, scan line offset, and earth image file counts are affected by a product possessing less than the full complement of bands.
The product delivered to Landsat 7 data users is packaged in HDF - an open standard selected by NASA for Earth Observing System (EOS) data products. HDF is a self-describing format that allows an application to interpret the structure and contents of a file without outside information. HDF allows Landsat 0R products to be shared across different computer platforms without modification and is supported by a public domain software library consisting of access tools and various utilities.
Product users are directed to the Landsat 7 0R Distribution Product Data Format Control Book, Volume 5 (PDF) for details regarding the HDF design used for the 0R product. Included are references to NCSA-authored documentation. New users should begin with Getting Started with HDF while the HDF User's Guide and HDF Reference Manual are excellent resources for the HDF programmer.
11.2 Level 1R Product |
The 2008 single-product data policy changes at EROS made the Level 1R product option obsolete. The following paragraphs are only relevant only from a historical perspective.
The Level 1R product is a radiometrically corrected 0R product. Radiometric correction is performed using either the CRAM gains in the CPF or gains computed on the fly from the IC data. The choice is available to a user when the product is ordered. The biases used are always calculated from the IC data. Image artifacts such as banding, striping, and scan correlated shift are removed prior to radiometric correction. Radiometric corrections are not reversible. The 1R product geometry is identical to the input Level OR data.
During 1R product rendering image pixels are converted to units of absolute radiance using 32 bit floating point calculations. Pixel values are then multiplied by 100 and converted to 16 bit integers prior to media output. Two digits of decimal precision are thus preserved. One merely divides each pixel value by 100 to convert the 1R image data back to radiance units. The 16 bit 1R product is twice the data volume of an alike 0R product. Note for band 6: A bias was found in the pre-launch calibration by a team of independent investigators post launch. This was corrected for in the LPGS processing system beginning Dec 20, 2000. For data processed before this, the 16 bit image radiances are 0.31 w/m2 ster um too high. See the official announcement for more details.
Two options existed for users when defining the size or spatial extent of a Landsat level 1R product ordered from the LP-DAAC.
A complete scene-sized 1R product consists of 17 data sets derived from the wideband telemetry, an IAS-generated calibration parameter file, a product specific metadata file, a geolocation index generated by EOSDIS Core System (ECS), and an HDF directory. Therefore, if you order a complete (i.e. all bands) scene-based 0R product it will have 21 distinct files. There are two fewer data files than an alike 0R product due to the fact that the multiple PCD and MSCD files are merged into single consensus files. Please reference the 0R file product for individual file descriptions.
A user could order a subset of the available bands which affected the actual file count in a 1R product. In all cases, however, every product included one consensus PCD file, one consensus MSCD files, three metadata files, the CPF, and the HDF directory. Only the internal calibrator, scan line offset, and earth image file counts were affected by a product possessing less than the full complement of bands.
The 1R product was delivered to users only in the HDF format. The HDF 0R and 1R formats are nearly identical. Exceptions include the united PCD and MSCD files and an enhanced product specific metadata file that reflects 1R correction characteristics. Please refer to the Landsat 7 0R Distribution Product Data Format Control Book, Volume 5 (PDF) for details regarding HDF specifics. Additional information unique to the 1R product can be found in the ESDIS Level 1 Product Generation system Output Files DFCB (PDF).
11.3 Level 1G Product |
The 1G product available to users from EROS is a radiometrically and geometrically corrected Level 0R image. The correction algorithms employed model the spacecraft and sensor using data generated by onboard computers during imaging events. Primary inputs are the PCD, which includes the attitude and ephemeris profiles, the definitive ephemeris (if available) and the MSCD. Refined parameters from the CPF, ground control points and a digital elevation model are also used to improve the overall geometric fidelity of the standard level-one terrain-corrected (L1T) product.
The L1T correction process utilizes both ground control points (GCP) and digital elevation models (DEM) to attain absolute geodetic accuracy. The WGS84 ellipsoid is employed as the Earth model for the Universal Transverse Mercator (UTM) coordinate transformation. Associated with the UTM projection is a unique set of projection parameters that flow from the USGS General Cartographic Transformation Package. The end result is a geometrically rectified product free from distortions related to the sensor (e.g. jitter, view angle effects), satellite (e.g. attitude deviations from nominal), and Earth (e.g. rotation, curvature, relief).
Geodetic accuracy of the L1T product depends on the accuracy of the GCPs and the resolution of the DEM used*. The 2005 Global Land Survey is used as the source for GCPs while the primary terrain data is the Shuttle Radar Topographic Mission DEM. Scenes that have a quality scores of 99 and less than 40 percent cloud cover are automatically processed, and any archived scene, regardless of cloud cover, can be ordered through one of two EROS web portals (Product Ordering.)
During L1T processing the 0R image data undergoes two-dimensional resampling according to the following set of parameters:
|
|
which is also expressed as:
Lλ = ((LMAXλ - LMINλ)/(QCALMAX-QCALMIN)) * (QCAL-QCALMIN) + LMINλ
where: | Lλ | = Spectral Radiance at the sensor's aperture in watts/(meter squared * ster * μm) |
Grescale |
= Rescaled gain (the data product "gain" contained in the Level 1 product header or ancillary data record) in watts/(meter squared * ster * μm)/DN |
|
Brescale |
= Rescaled bias (the data product "offset" contained in the Level 1 product header or ancillary data record ) in watts/(meter squared * ster * μm) | |
QCAL |
= the quantized calibrated pixel value in DN | |
LMINλ | = the spectral radiance that is scaled to QCALMIN in watts/(meter squared * ster * μm) | |
LMAXλ | = the spectral radiance that is scaled to QCALMAX in watts/(meter squared * ster * μm) | |
QCALMIN | = the minimum quantized calibrated pixel value (corresponding to LMINλ) in DN |
|
QCALMAX | = the maximum quantized calibrated pixel value (corresponding to LMAXλ) in DN = 255 |
The LMINs and LMAXs are the spectral radiances for each band at digital numbers 0 or 1 and 255 (i.e QCALMIN, QCALMAX), respectively. LPGS used 1 for QCALMIN while NLAPS used 0 for QCALMIN for data products processed before April 5, 2004. NLAPS from that date now uses 1 for the QCALMIN value. Other product differences exist as well. One LMIN/LMAX set exists for each gain state. These values will change slowly over time as the ETM+ detectors lose responsivity. Table 11.2 lists two sets of LMINs and LMAXs. The first set should be used for both LPGS and NLAPS 1G products created before July 1, 2000 and the second set for 1G products created after July 1, 2000. Please note the distinction between acquisition and processing dates. Use of the appropriate LMINs and LMAXs will ensure accurate conversion to radiance units. Note for band 6: A bias was found in the pre-launch calibration by a team of independent investigators post launch. This was corrected for in the LPGS processing system beginning Dec 20, 2000. For data processed before this, the image radiances given by the above transform are 0.31 w/m2 ster um too high. See the official announcement for more details. Note for the Multispectral Scanner (MSS), Thematic Mapper (TM), and Advanced Land Imager (ALI) sensors: the required radiometry constants are tabulated in this PDF file.
Table 11.2 ETM+ Spectral Radiance Range watts/(meter squared * ster * μm) |
||||||||
---|---|---|---|---|---|---|---|---|
Band Number![]() |
Processed Before July 1, 2000 | Proccessed After July 1, 2000 | ||||||
Low Gain | High Gain | Low Gain | High Gain | |||||
LMIN | LMAX | LMIN | LMAX | LMIN | LMAX | LMIN | LMAX | |
1 | -6.2 | 297.5 | -6.2 | 194.3 | -6.2 | 293.7 | -6.2 | 191.6 |
2 | -6.0 | 303.4 | -6.0 | 202.4 | -6.4 | 300.9 | -6.4 | 196.5 |
3 | -4.5 | 235.5 | -4.5 | 158.6 | -5.0 | 234.4 | -5.0 | 152.9 |
4 | -4.5 | 235.0 | -4.5 | 157.5 | -5.1 | 241.1 | -5.1 | 157.4 |
5 | -1.0 | 47.70 | -1.0 | 31.76 | -1.0 | 47.57 | -1.0 | 31.06 |
6 | 0.0 | 17.04 | 3.2 | 12.65 | 0.0 | 17.04 | 3.2 | 12.65 |
7 | -0.35 | 16.60 | -0.35 | 10.932 | -0.35 | 16.54 | -0.35 | 10.80 |
8 | -5.0 | 244.00 | -5.0 | 158.40 | -4.7 | 243.1 | -4.7 | 158.3 |
Where: | ![]() |
![]() |
= Unitless planetary reflectance |
![]() |
= Spectral radiance at the sensor's aperture |
![]() |
= Earth-Sun distance in astronomical units from an Excel file |
![]() |
= Mean solar exoatmospheric irradiances from Table 11.3 |
![]() |
= Solar zenith angle in degrees |
Table 11.3 ETM+ Solar Spectral Irradiances (generated using the Thuillier solar spectrum) |
|
---|---|
Band | watts/(meter squared * μm) |
1 | 1997 |
2 | 1812 |
3 | 1533 |
4 | 1039 |
5 | 230.8 |
7 | 84.90 |
8 | 1362. |
Table 11.4 Earth-Sun Distance in Astronomical Units |
|||||||||
---|---|---|---|---|---|---|---|---|---|
Day of Year |
Distance |
Day of Year |
Distance |
Day of Year |
Distance |
Day of Year |
Distance |
Day of Year |
Distance |
1 |
.98331 |
74 |
.99446 |
152 |
1.01403 |
227 |
1.01281 |
305 |
.99253 |
15 |
.98365 |
91 |
.99926 |
166 |
1.01577 |
242 |
1.00969 |
319 |
.98916 |
32 |
.98536 |
106 |
1.00353 |
182 |
1.01667 |
258 |
1.00566 |
335 |
.98608 |
46 |
.98774 |
121 |
1.00756 |
196 |
1.01646 |
274 |
1.00119 |
349 |
.98426 |
60 |
.99084 |
135 |
1.01087 |
213 |
1.01497 |
288 |
.99718 |
365 |
.98333 |
Where: | ![]() |
T | = Effective at-satellite temperature in Kelvin |
K2 | = Calibration constant 2 from Table 11.5 |
K1 | = Calibration constant 1 from Table 11.5 |
L | = Spectral radiance in watts/(meter squared * ster * µm) |
Table 11.5 ETM+ and TM Thermal Band Calibration Constants | ||
---|---|---|
![]() |
Constant 1- K1 watts/(meter squared * ster * μm) |
Constant 2 - K2 Kelvin |
Landsat 7 | 666.09 | 1282.71 |
Landsat 5 | 607.76 | 1260.56 |
The same two 1R options exist for users when defining the size or spatial extent of a Landsat level 1G product ordered from the LP-DAAC.
The 1G product ordered from the LP-DAAC consists of the corrected image files and descriptive metadata. All other ancillary files delivered with the 0R and 1R products are not included. A user may order a subset of the available bands which affects the actual file count in a 1G product.
The 1G product can be packaged into one of following user-specified output formats:
The number of files comprising an HDF-formatted 1G product will vary according to the number of bands ordered. A product with a full band complement has 11 files - the HDF directory, a metadata file, and a separate file for each band. The HDF directory and metadata files are always present regardless of bands ordered. Please refer to the Landsat 7 0R Distribution Product Data Format Control Book, Volume 5 (PDF) for details regarding band file specifics. The 1R metadata file description can be found in the ESDIS Level 1 Product Generation System Output Files DFCB (PDF).
The HDF format can be specified for any type of 1G product ordered from the LP-DAAC.
A derivative of the Fast Format (Fast-L7) used by EOSAT for Landsat (FAST-B) and Indian Remote Sensing products (Fast-C) was created for Landsat 7. Several differences are worth noting. File names are now included in the administrative record which allows for direct file access. A separate header file now accompanies the panchromatic, thermal and VNIR/SWIR band groups for Landsat 7. For Fast-B and Fast-C all bands were resampled to a common grid cell size thus permitting a single header file. In all likelihood each of the band groups for Landsat 7 will be resampled to a common resolution (i.e. 15, 30, & 60 meters) thus requiring a distinct header file for each.
All critical fields required for product ingest were left unchanged in the Fast L-7A Format. As a consequence Heritage Fast readers residing on user systems can be used for the Landsat 7 Fast formatted product. A full layout of the Fast L-7A Format can be found in the the ESDIS Level 1 Product Generation system Output Files DFCB.
The Fast-L7 format supports all variations of the 1G product.
The TIFF file consists of a number of label (tags) which describe certain properties of the file (such as gray levels, color table, byte format, compression size). After the initial tags comes the image data which may be interrupted by more descriptive tags. GeoTIFF refers to TIFF files which have geographic (or cartographic) data embedded as tags within the TIFF file. The geographic data can then be used to position the image in the correct location and geometry on the screen of a geographic information display.
Baseline TIFF image types can be bilevel, greyscale, palette color, and full color (24 bit). For simplicity's sake the grayscale model was implemented for the Landsat 7 GeoTIFF product. Under this implementation each ordered band is delivered as its own 8 bit greyscale GeoTIFF image. A standard WRS scene possessing the full band complement would thus be comprised of nine separate GeoTIFF images or files. No other files accompany the product. For detailed information regarding the Landsat 7 GeoTIFF implementation please refer to the ESDIS Level 1 Product Generation system Output Files DFCB (PDF). For GeoTIFF details, please download the GeoTIFF Format Specification (PDF) or visit this web site.
At the present time GeoTIFF format cannot be used for the Space Oblique Mercator and Oblique Mercator projections. Products projected into these reference systems must be formatted using HDF or Fast-L7.
An instrument malfunction occurred onboard Landsat 7 on May 31, 2003. The problem was caused by failure of the Scan Line Corrector (SLC), which compensates for the forward motion of the satellite. Subsequent efforts to recover the SLC have not been successful, and the problem is permanent.
The Landsat 7 Enhanced Thematic Mapper Plus (ETM+) is still capable of acquiring useful image data with the SLC turned off, particularly within the central portion of any given scene. Landsat 7 ETM+ will therefore continue to acquire image data in the "SLC-off" mode.
EDC has recently released several Landsat 7 ETM+ SLC-off data products. The first, a gap-present product became available on October 22, 2003. This product release includes all image data acquired by Landsat 7 in SLC-off mode from July 14, 2003 to present, excluding a 2-week interval from 9/3/03 to 9/17/03.
The center of a gap-present SLC-off data product is very similar in quality to previous Landsat 7 data. However, the scene's edges will contain alternating scan lines of missing data (Level 1G) or duplicated data (Level 0Rp or L1R). The precise location of the affected scan lines will vary from scene to scene, and these gaps will not be visible on the browse image preview when ordering SLC-off data. A preliminary report regarding the utility of Landsat 7 SLC-off data is available in PDF form. This report includes input from scientists affiliated with the USGS, NASA, and the Landsat 7 Science Team
The gap-present SLC-off data product is available as a single scene entity in Level 1G terrain corrected (L1T) form. As of November, 2008, the USGS offers all archived Landsat scenes to the public at no charge. Newly acquired Landsat 7 ETM+ SLC-off and Landsat 5 TM images with less than 40 percent cloud cover are automatically processed and made freely available for immediate download. SLC-off data products can be searched and ordered via the Earth Explorer, and Global Visualization L7 Image Browser.
The second product now being offered (as of May 10, 2004) is in 1G form and has the gap areas filled with Landsat 7 data acquired at a similar time of year and prior to the SLC failure. The two scenes are geometrically registered, and a histogram matching technique is applied to the fill pixels which provides the best-expected radiance values for the missing data.
The USGS, in conjunction with NASA, is continuing to research other methods of providing merged data products and will continue to provide information resulting from this work as it becomes available. Various methodoligies have been examined to fill the data gaps with observations acquired during prior or later than the target scene of interest. An exampled of a gap-filled product is illustrated in Figure 11.5.
![]() |
Figure 11.5 Top image: pre-SLC anomaly, middle of image. Middle image: scene after SLC anomaly. Bottom image: scene after SLC anomaly with data gaps filled. |
Source: http://landsathandbook.gsfc.nasa.gov